Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,410,753 T→C 30.7% F8S (TTT→TCT)  hemY → predicted protoheme IX synthesis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,410,7530TC30.7% 217.3 / 133.1 150F8S (TTT→TCT) hemYpredicted protoheme IX synthesis protein
Reads supporting (aligned to +/- strand):  ref base T (98/6);  new base C (30/16);  total (128/22)
Fisher's exact test for biased strand distribution p-value = 1.27e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 165 total aligned reads displayed.

ACCGCAAGCTGATACTCCGGCAGCCGCGCCGCAAGGAGAATAACCATGCTAAAAGTGTTATTGCTCTTTGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCGATGATTGCCGGCCATCAGGGTTATGTGCTGATCC  >  minE/2410686‑2410821
                                                                   |                                                                    
accGCAAGCTGATACTCCGGCAGCCGCGCCGCAAGGAGAATAACCATGGTAAAAGTGTTATTGCTCTTtgt                                                                   <  1:2107593/71‑1 (MQ=255)
accGCAAGCTGATACTCCGGCAGCCGCGCCGCAAGGAGAATAACCATGCTAAAAGTGTTATTGCTCTTtgt                                                                   <  1:2224039/71‑1 (MQ=255)
accGCAAGCTGATACTCCGGCAGCCGCGCCGCAAGGAGAATAACCATGCTAAAAGTGTTATTGCTCTTtgt                                                                   <  1:592275/71‑1 (MQ=255)
        cTGATACTCCGGCAGCCGCGCCGCAAGGAGAATAACCATGCTAAAAGTGTTATTGCTCTTTGTGTTGCTGa                                                           >  1:1868895/1‑71 (MQ=255)
        cTGATACTCCGGCAGCCGCGCCGCAAGGAGAATAACCATGCTAAAAGTGTTATTGCTCTTTGTGTTGCTGa                                                           >  1:646165/1‑71 (MQ=255)
        cTGATACTCCGGCAGCCGCGCCGCAAGGAGAATAACCATGCTAAAAGTGTTATTGCTCTTTGTGTTGCTGa                                                           >  1:2319751/1‑71 (MQ=255)
                         gcgccgcAAGGAGAATAACCATGCTAAAAGTGTTATTGCTCttt                                                                     >  1:878220/1‑44 (MQ=255)
                         gcgccgcAAGGAGAATAACCATGCTAAAAGTGTTATTGCTCTTTGTGTTGCTGATTGCGGGGATCGTGGtt                                          >  1:2576389/1‑71 (MQ=255)
                         gcgccgcAAGGAGAATAACCATGCTAAAAGTGTTATTGCTCTTTGTGTTGCTGATTGCGGGGATCGTGGtt                                          >  1:2812908/1‑71 (MQ=255)
                         gcgccgcAAGGAGAATAACCATGCTAAAAGTGTTATTGCTCTTTGTGTTGCTGATTGCGGGGATCGTGGtt                                          >  1:1506932/1‑71 (MQ=255)
                         gcgccgcAAGGAGAATAACCATGCTAAAAGTGTTATTGCTCTTTGTGTTGCTGATTGCGGGGATCGTGGtt                                          >  1:150232/1‑71 (MQ=255)
                         gcgccgcAAGGAGAATAACCATGCTAAAAGTGTTATTGCTCTTTGTGTTGCTGATTGCGGGGATCGTGGtt                                          >  1:2864465/1‑71 (MQ=255)
                         gcgccgcAAGGAGAATAACCATGCTAAAAGTGTTATTGCTCTTTGTGTTGCTGATTGCGGGGATCGTGGtt                                          >  1:291352/1‑71 (MQ=255)
                         gcgccgcAAGGAGAATAACCATGCTAAAAGTGTTATTGCTCTTTGTGTTGCTGATTGCGGGGATCGTGGtt                                          >  1:483123/1‑71 (MQ=255)
                         gcgccgcAAGGAGAATAACCATGCTAAAAGTGTTATTGCTCTTTGTGTTGCTGATTGCGGGGATCGTGGtt                                          >  1:2255180/1‑71 (MQ=255)
                         gcgccgcAAGGAGAATAACCATGCTAAAAGTGTTATTGCTCTTTGTGTTGCTGATTGCGGGGATCGTGGtt                                          >  1:48321/1‑71 (MQ=255)
                         gcgccgcAAGGAGAATAACCATGCTAAAAGTGTTATTGCTCTTTGTGTTGCTGATTGCGGGGATCGTGGtt                                          >  1:51380/1‑71 (MQ=255)
                         gcgccgcAAGGAGAATAACCATGCTAAAAGTGTTATTGCTCTTTGTGTTGCTGATTGCGGGGATCGTGGtt                                          >  1:1226965/1‑71 (MQ=255)
                         gcgccgcAAGGAGAATAACCATGCTAAAAGTGTTATTGCTCTTTGTGTTGCTGATTGCGGGGATCGTGGtt                                          >  1:613751/1‑71 (MQ=255)
                         gcgccgcAAGGAGAATAACCATGCTAAAAGTGTTATTGCTCTTTGTGTTGCTGATTGCGGGGATCGTGGtt                                          >  1:1946049/1‑71 (MQ=255)
                         gcgccgcAAGGAGAATAACCATGCTAAAAGTGTTATTGCTCTTTGTGTTGCTGATTGCGGGGATCGTGGtt                                          >  1:69348/1‑71 (MQ=255)
                         gcgccgcAAGGAGAATAACCATGCTAAAAGTGTTATTGCTCTTTGTGTTGCTGATTGCGGGGATCGTGGtt                                          >  1:1971195/1‑71 (MQ=255)
                         gcgccgcAAGGAGAATAACCATGCTAAAAGTGTTATTGCTCTTTGTGTTGCTGATTGCGGGGATCGTGGtt                                          >  1:873863/1‑71 (MQ=255)
                         gcgccgcAAGGAGAATAACCATGCTAAAAGTGTTATTGCTCTTTGTGTTGCTGATTGCGGGGATCGTGGt                                           >  1:1116953/1‑70 (MQ=255)
                         gcgccgcAAGGAGAATAACCATGCTAAAAGTGTTATTGCTCTTTGTGTTGCTGATTGCGGGGATCGTGGt                                           >  1:1029422/1‑70 (MQ=255)
                         gcgccgcAAGGAGAATAACCATGCTAAAAGTGTTATTGCTCTTTGTGTTACTGATTGCGGGGATCGTGGtt                                          >  1:1413152/1‑71 (MQ=255)
                         gcgccgcAAGGAGAATAACCATGCT‑‑AAGTGTTATTGCTCttt                                                                     >  1:1896287/1‑42 (MQ=255)
                                    agaATAACCATGGTAAAAGTGTTATTGCTCTTTGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCgatgat                               >  1:241416/1‑71 (MQ=255)
                                    agaATAACCATGCTAAAAGTGTTATTGCTCTTTGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCgatgat                               >  1:2107759/1‑71 (MQ=255)
                                    agaATAACCATGCTAAAAGTGTTATTGCTCTTTGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCgatgat                               >  1:2103687/1‑71 (MQ=255)
                                    agaATAACCATGCTAAAAGTGTTATTGCTCTTTGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCgatgat                               >  1:2209957/1‑71 (MQ=255)
                                    agaATAACCATGCTAAAAGTGTTATTGCTCTTTGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCgatgat                               >  1:2345615/1‑71 (MQ=255)
                                    agaATAACCATGCTAAAAGTGTTATTGCTCTTTGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCgatgat                               >  1:2608174/1‑71 (MQ=255)
                                    agaATAACCATGCTAAAAGTGTTATTGCTCTTTGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCgatgat                               >  1:2675691/1‑71 (MQ=255)
                                    agaATAACCATGCTAAAAGTGTTATTGCTCTTTGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCgatgat                               >  1:376966/1‑71 (MQ=255)
                                    agaATAACCATGCTAAAAGTGTTATTGCTCTTTGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCgatgat                               >  1:529538/1‑71 (MQ=255)
                                    agaATAACCATGCTAAAAGTGTTATTGCTCTTTGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCgatgat                               >  1:624899/1‑71 (MQ=255)
                                    agaATAACCATGCTAAAAGTGTTATTGCTCTTTGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCgatgat                               >  1:745090/1‑71 (MQ=255)
                                    agaATAACCATGCTAAAAGTGTTATTGCTCTTTGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCgatgat                               >  1:969797/1‑71 (MQ=255)
                                    agaATAACCATGCTAAAAGTGTTATTGCTCTTTGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCgatgat                               >  1:1360276/1‑71 (MQ=255)
                                    agaATAACCATGCTAAAAGTGTTATTGCTCTTTGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCgatgat                               >  1:121291/1‑71 (MQ=255)
                                    agaATAACCATGCTAAAAGTGTTATTGCTCTTTGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCgatgat                               >  1:1125690/1‑71 (MQ=255)
                                    agaATAACCATGCTAAAAGTGTTATTGCTCTTTGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCgatgat                               >  1:1795499/1‑71 (MQ=255)
                                    agaATAACCATGCTAAAAGTGTTATTGCTCTTTGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCgatgat                               >  1:1264899/1‑71 (MQ=255)
                                    agaATAACCATGCTAAAAGTGTTATTGCTCTTTGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCgatgat                               >  1:1318266/1‑71 (MQ=255)
                                    agaATAACCATGCTAAAAGTGTTATTGCTCTTTGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCgatgat                               >  1:1811599/1‑71 (MQ=255)
                                    agaATAACCATGCTAAAAGTGTTATTGCTCTTTGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCgatgat                               >  1:1850729/1‑71 (MQ=255)
                                    agaATAACCATGCTAAAAGTGTTATTGCTCTTTGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCgatgat                               >  1:1863315/1‑71 (MQ=255)
                                    agaATAACCATGCTAAAAGTGTTATTGCTCTTTGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCgatgat                               >  1:1563122/1‑71 (MQ=255)
                                    agaATAACCATGCTAAAAGTGTTATTGCTCTTTGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCgatgat                               >  1:1073980/1‑71 (MQ=255)
                                    agaATAACCATGCTAAAAGTGTTATTGCTCTTTGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCgatgat                               >  1:1921308/1‑71 (MQ=255)
                                    agaATAACCATGCTAAAAGTGTTATTGCTCTTTGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCgatgat                               >  1:1420004/1‑71 (MQ=255)
                                    agaATAACCATGCTAAAAGTGTTATTGCTCTTTGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCgatgat                               >  1:1434664/1‑71 (MQ=255)
                                                     aGTGTTATTGCTCTTTGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCGATGATTGCCGGCCATCAGGGtt              <  1:2177482/71‑1 (MQ=255)
                                                     aGTGTTATTGCTCTCTGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCGATGATTGCCGGCCATCAGGGtt              <  1:840633/71‑1 (MQ=255)
                                                     aGTGTTATTGCTCTCTGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCGATGATTGCCGGCCATCAGGGtt              <  1:494847/71‑1 (MQ=255)
                                                     aGTGTTATTGCTCTCTGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCGATGATTGCCGGCCATCAGGGtt              <  1:473125/71‑1 (MQ=255)
                                                     aGTGTTATTGCTCTCTGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCGATGATTGCCGGCCATCAGGGtt              <  1:1428693/71‑1 (MQ=255)
                                                     aGTGTTATTGCTCTCTGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCGATGATTGCCGGCCATCAGGGtt              <  1:1485684/71‑1 (MQ=255)
                                                     aGTGTTATTGCTCTCTGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCGATGATTGCCGGCCATCAGGGtt              <  1:1146354/71‑1 (MQ=255)
                                                     aGTGTTATTGCTCTCTGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCGATGATTGCCGGCCATCAGGGtt              <  1:2295916/71‑1 (MQ=255)
                                                     aGTGTTATTGCTCTCTGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCGATGATTGCCGGCCATCAGGGtt              <  1:2199399/71‑1 (MQ=255)
                                                     aGTGTTATTGCTCTCTGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCGATGATTGCCGGCCATCAGGGtt              <  1:844744/71‑1 (MQ=255)
                                                     aGTGTTATTGCTCTCTGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCGATGATTGCCGGCCATCAGGGtt              <  1:2063600/71‑1 (MQ=255)
                                                     aGTGTTATTGCTCTCTGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCGATGATTGCCGGCCATCAGGGtt              <  1:2013349/71‑1 (MQ=255)
                                                        gtTATTGCTCTTTGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCGATGATTGCCGGCCATCAGGGTTAtg           >  1:722067/1‑71 (MQ=255)
                                                        gtTATTGCTCTTTGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCGATGATTGCCGGCCATCAGGGTTAtg           >  1:2841660/1‑71 (MQ=255)
                                                        gtTATTGCTCTTTGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCGATGATTGCCGGCCATCAGGGTTAtg           >  1:884712/1‑71 (MQ=255)
                                                        gtTATTGCTCTTTGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCGATGATTGCCGGCCATCAGGGTTAtg           >  1:891461/1‑71 (MQ=255)
                                                        gtTATTGCTCTTTGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCGATGATTGCCGGCCATCAGGGTTAtg           >  1:931623/1‑71 (MQ=255)
                                                        gtTATTGCTCTTTGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCGATGATTGCCGGCCATCAGGGTTAtg           >  1:201839/1‑71 (MQ=255)
                                                        gtTATTGCTCTTTGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCGATGATTGCCGGCCATCAGGGTTAtg           >  1:2677925/1‑71 (MQ=255)
                                                        gtTATTGCTCTTTGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCGATGATTGCCGGCCATCAGGGTTAtg           >  1:2149719/1‑71 (MQ=255)
                                                        gtTATTGCTCTTTGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCGATGATTGCCGGCCATCAGGGTTAtg           >  1:2331697/1‑71 (MQ=255)
                                                        gtTATTGCTCTTTGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCGATGATTGCCGGCCATCAGGGTTAtg           >  1:1673910/1‑71 (MQ=255)
                                                        gtTATTGCTCTTTGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCGATGATTGCCGGCCATCAGGGTTAtg           >  1:2430607/1‑71 (MQ=255)
                                                           aTTGCTCTCTGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCGATGATTGCCGGCCAt                     >  1:1190986/1‑58 (MQ=255)
                                                           aTTGCTCTCTGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCGATGATTGCCGGCCATCAGGGTTATGTGc        >  1:103965/1‑71 (MQ=255)
                                                           aTTGCTCTCTGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCGATGATTGCCGGCCATCAGGGTTATGTGc        >  1:1998140/1‑71 (MQ=255)
                                                           aTTGCTCTCTGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCGATGATTGCCGGCCATCAGGGTTATGTGc        >  1:935534/1‑71 (MQ=255)
                                                           aTTGCTCTCTGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCGATGATTGCCGGCCATCAGGGTTATGTGc        >  1:2065984/1‑71 (MQ=255)
                                                           aTTGCTCTCTGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCGATGATTGCCGGCCATCAGGGTTATGTGc        >  1:1981573/1‑71 (MQ=255)
                                                           aTTGCTCTCTGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCGATGATTGCCGGCCATCAGGGTTATGTGc        >  1:1914558/1‑71 (MQ=255)
                                                           aTTGCTCTCTGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCGATGATTGCCGGCCATCAGGGTTATGTGc        >  1:1803782/1‑71 (MQ=255)
                                                           aTTGCTCTCTGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCGATGATTGCCGGCCATCAGGGTTATGTGc        >  1:724070/1‑71 (MQ=255)
                                                           aTTGCTCTCTGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCGATGATTGCCGGCCATCAGGGTTATGTGc        >  1:1265381/1‑71 (MQ=255)
                                                           aTTGCTCTCTGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCGATGATTGCCGGCCATCAGGGTTATGTGc        >  1:163409/1‑71 (MQ=255)
                                                           aTTGCTCTCTGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCGATGATTGCCGGCCATCAGGGTTATGTGc        >  1:511686/1‑71 (MQ=255)
                                                           aTTGCTCTCTGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCGATGATTGCCGGCCATCAGGGTTATGTGc        >  1:2721572/1‑71 (MQ=255)
                                                           aTTGCTCTCTGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCGATGATTGCCGGCCATCAGGGTTATGTGc        >  1:384958/1‑71 (MQ=255)
                                                           aTTGCTCTCTGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCGATGATTGCCGGCCATCAGGGTTATGTGc        >  1:2788210/1‑71 (MQ=255)
                                                           aTTGCTCTCTGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCGATGATTGCCGGCCATCAGGGTTATGTGc        >  1:1561522/1‑71 (MQ=255)
                                                           aTTGCTCTCTGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCGATGATTGACGGCCATCa                   >  1:447518/1‑60 (MQ=255)
                                                                  tttGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCGATGATTGCCGGCCATCAGGGTTATGTGCTGAt  >  1:1078970/1‑68 (MQ=255)
                                                                  tttGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCGATGATTGCCGGCCATCAGGGTTATGTGCTGATcc  >  1:2537735/1‑70 (MQ=255)
                                                                  tttGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCGATGATTGCCGGCCATCAGGGTTATGTGCTGATc   >  1:2495022/1‑69 (MQ=255)
                                                                  tttGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCGATGATTGCCGGCCATCAGGGTTATGTGCTGATc   >  1:1095813/1‑69 (MQ=255)
                                                                  tttGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCGATGATTGCCGGCCATCAGGGTTATGTGCTGATc   >  1:1585069/1‑69 (MQ=255)
                                                                  tttGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCGATGATTGCCGGCCATCAGGGTTATGTGCTGATc   >  1:2637111/1‑69 (MQ=255)
                                                                  tttGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCGATGATTGCCGGCCATCAGGGTTATGTGCTGATc   >  1:399402/1‑69 (MQ=255)
                                                                   |                                                                    
ACCGCAAGCTGATACTCCGGCAGCCGCGCCGCAAGGAGAATAACCATGCTAAAAGTGTTATTGCTCTTTGTGTTGCTGATTGCGGGGATCGTGGTTGGCCCGATGATTGCCGGCCATCAGGGTTATGTGCTGATCC  >  minE/2410686‑2410821

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: