Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 996,989 A→T 10.0% Y1191F (TAT→TTT)  mukB → cell division protein MukB

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606996,9890A→T10.0% 2.9 ‑4.5 ‑0.0403572mukBcell division protein MukB
Reads supporting (aligned to +/- strand):  new base (4/0):  ref base (10/26):  total (14/26)
Fisher's exact test strand distribution p-value = 1.10e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.51e-01

GGAGCGTAAAATTCAGTTCTTCGTGGCGGTTTATCAGCATCTGCGTGAACGTATTCGTCAGGATATTA  >  REL606/996957‑997024
                                |                                   
ggAGCGTAAAATTCAGTTCTTCGTGGCGGTTTATCa                                  <  1:2348382/36‑1 (MQ=255)
  aGCGTAAAATTCAGTTCTTCGTGGCGGTTTTTTAGc                                >  1:1354969/1‑36 (MQ=37)
  aGCGTAAAATTCAGTTCTTCGTGGCGGTTTATCAGc                                <  1:5897263/36‑1 (MQ=255)
  aGCGTAAAATTCAGTTCTTCGTGGCGGTTTATCAGc                                <  1:1930887/36‑1 (MQ=255)
    cGTAAAATTCAGTTCTTCGTGGCGGTTTATCAGCAt                              <  1:5040643/36‑1 (MQ=255)
    cGTAAAATTCAGTTCTTCGTGGCGGTTTATCAGCAt                              <  1:6457368/36‑1 (MQ=255)
     gTAAAATTCAGTTCTTCGTGGCGGTTTATCAGCATc                             <  1:2953354/36‑1 (MQ=255)
      tAAAATTCAGTTCTTCGTGGCGGTTTATCAGCATCt                            <  1:4453239/36‑1 (MQ=255)
       aaaaTTCAGTTCTTCGTGGCGGTTTATCAGCATCTg                           <  1:4105752/36‑1 (MQ=255)
           ttCAGTTCTTCGTGGCGGTTTTTTAGCCTTTTCGTg                       >  1:2107048/1‑36 (MQ=16)
           ttCAGTTCTTCGTGGCGGTTTATCAGCATCTGCgag                       <  1:5855980/36‑3 (MQ=255)
           ttCAGTTCTTCGTGGCGGTTTATCAGCATCTGCGTg                       <  1:486018/36‑1 (MQ=255)
            tCAGTTCTTCGTGGCGGTTTATCAGCATCTGCGTGa                      <  1:6651987/36‑1 (MQ=255)
              aGTTCTTCGTGGCGGTTTATCAGCATCTGCGTGAAc                    <  1:7295931/36‑1 (MQ=255)
              aGTTCTTCGTGGCGGTTTATCAGCATCTGCGTGAAc                    <  1:382250/36‑1 (MQ=255)
               gTTCTTCGTGGCGGTTTATCAGCATCTGCGTGAACg                   <  1:1187474/36‑1 (MQ=255)
                ttcttcGTGGCGGTTTATCAGCATCTGCGTGAACGt                  <  1:7484111/36‑1 (MQ=255)
                ttcttcGTGGCGGGTTTTTATCATTTGCGTGGACGt                  >  1:5493091/1‑36 (MQ=9)
                 tcttcGTGGCGGTTTTTTATCATCTTCGTGAACGTa                 >  1:5327630/1‑36 (MQ=19)
                    tcGTGGCGGTTTATCAGCATCTGCGTGAACGTATTc              <  1:3893041/36‑1 (MQ=255)
                     cGTGGCGGTTTATCAGCATCTGCTTGAACGTATTCg             <  1:1807914/36‑1 (MQ=255)
                     agTGGCGTTTTATCAGCATCTGCTTGAACGTATTCg             <  1:3364419/35‑1 (MQ=25)
                      ggggCGGTTTAGCAGCATCTGCGGGAACGTATTCGt            <  1:3403747/34‑1 (MQ=37)
                       tGGCGGTTTATCAGCATCTGCGTGAACGTATTCGTc           <  1:3338931/36‑1 (MQ=255)
                        ggCGGTTTATCAGCATCTGCGTTAACGTATTCGTCa          >  1:2794877/1‑36 (MQ=255)
                        ggCGGTTTATCAGCATCTGCGTGAACGTATTCGTCa          <  1:1092624/36‑1 (MQ=255)
                        ggCAGTTTATCAGCATCTGCGTGAACGTATTCGTCa          <  1:859395/36‑1 (MQ=255)
                          cGGTTTATCAGCATCTGCGTGAACGTATTCGTCAgg        <  1:2355908/36‑1 (MQ=255)
                          cGGTTTATCAGCATCTGCGTGAACGTATTCGTAAgg        >  1:3684118/1‑36 (MQ=255)
                           ggTTTATCAGCATCTGCGTGAACGTATTCGTCAGGa       >  1:7588501/1‑36 (MQ=255)
                            gTTTATCAGCATCTGCGTGAACGTATTCGTCAGGat      <  1:5359645/36‑1 (MQ=255)
                            gTTTATCAGCATCTGCGTGAACGTATTCGTCAGGat      <  1:6130110/36‑1 (MQ=255)
                             tttATCAGCATCTGCGTGAACGTATTCGTCAGGata     >  1:2636443/1‑36 (MQ=255)
                             tttATCAGCATCTGCGTGAACGTATTCGTCAGGata     >  1:3279918/1‑36 (MQ=255)
                             tttATCAGCATCTGCGTGAACGTATTCGTCAGGata     >  1:4510313/1‑36 (MQ=255)
                             tttATCAGCATCTGCGCGAACGTATTCGTCAGGata     >  1:1003693/1‑36 (MQ=255)
                              ttATCAGCATCTGCGTGAACGTATTCGTCAGGatat    >  1:3434910/1‑36 (MQ=255)
                               tATCAGCATCTGCGTGAACGTATTCGTCAGGAtatt   >  1:6076571/1‑36 (MQ=255)
                               tATCAGCATCTGCGTGAACGTATTCGTCAGGAtatt   >  1:7007043/1‑36 (MQ=255)
                               tATCAGCATCTGCGTGAACGTATTCGTCAGGAtatt   <  1:7495454/36‑1 (MQ=255)
                                aTCAGCCTCTGCGTGAACGTATTCGTCAGGAtatta  <  1:1551699/36‑1 (MQ=255)
                                aTCAGCATCTGCGTGAACGTATTCGTCAGGAtatta  <  1:4054229/36‑1 (MQ=255)
                                aTCAGCATCTGCGTGAACGTATTCGTCAGGAtatta  >  1:158237/1‑36 (MQ=255)
                                aTCAGCATCTGCGTGAACGTATTCGTCAGGAtatta  <  1:1002630/36‑1 (MQ=255)
                                |                                   
GGAGCGTAAAATTCAGTTCTTCGTGGCGGTTTATCAGCATCTGCGTGAACGTATTCGTCAGGATATTA  >  REL606/996957‑997024

Base quality scores:  ATCG  < 3 ≤  ATCG  < 11 ≤  ATCG  < 22 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG