Predicted mutation
evidence seq id position mutation annotation gene description
RA ACB143 120,983 (T)9→8 intergenic (+280/+92) polA → / ← engB DNA polymerase I/putative GTP‑binding protein EngB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB143120,9750T.90.4% 196.2 / 13.2 52intergenic (+272/+100)polA/engBDNA polymerase I/putative GTP‑binding protein EngB
Reads supporting (aligned to +/- strand):  ref base T (1/4);  new base . (30/17);  total (31/21)
Fisher's exact test for biased strand distribution p-value = 1.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACTGAATGGCTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTTCCATCCA  >  ACB143/120849‑121107
                                                                                                                              |                                                                                                                                    
aaCTGAATGGCTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGttttttttgac                                                                                                                            >  1:121586/1‑134 (MQ=255)
     aaTGGCTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTT‑CCCGGaa                                                                                                                      >  7:174975/1‑137 (MQ=255)
     aaTGGCTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGaa                                                                                                                      >  1:140544/1‑137 (MQ=255)
          cTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAAttttc                                                                                                                 >  1:184242/1‑135 (MQ=255)
          cTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGc                                                                                                                 >  1:283543/1‑137 (MQ=255)
            ggTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACACCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCg                                                                                                               >  7:341167/1‑137 (MQ=255)
                          aTGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTtct                                                                                                 <  8:342236/137‑1 (MQ=255)
                                cTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGaa                                                                                           >  6:174700/1‑137 (MQ=255)
                                 tATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAAt                                                                                          <  1:328637/137‑1 (MQ=255)
                                 tATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAAt                                                                                          <  4:259327/137‑1 (MQ=255)
                                 aaTGTAGTGTTAAGTGTTCGCAGATATCTCGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACATCCCCCTCCCCATACGAAGTCCGGTG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAAt                                                                                          <  3:45890/136‑1 (MQ=255)
                                   tgttgtGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCAGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGc                                                                                        >  7:313077/1‑137 (MQ=255)
                                     ttgtGTTAAGTTTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCCTACGAAGTCAGGGG‑TTTTTTTTGCCCGGAATTTTCCGGGAGCGGGTTTTCTCGCGAATGCGa                                                                                      >  8:40458/1‑137 (MQ=255)
                                          ttAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTAGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCg                                                                                 >  5:312830/1‑137 (MQ=255)
                                           tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGa                                                                                >  7:44862/1‑137 (MQ=255)
                                           tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGa                                                                                >  3:247148/1‑137 (MQ=255)
                                           tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGa                                                                                >  1:300338/1‑137 (MQ=255)
                                           tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATCTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGa                                                                                >  8:184622/1‑137 (MQ=255)
                                           tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTGTTTTTGCCCGGAATTTGCCGGGAGCGCGTTTTCTCGCGAATGCGAGGTCGa                                                                                >  2:120966/1‑137 (MQ=255)
                                                     cAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACAccgcc                                                                      >  3:8316/1‑137 (MQ=255)
                                                         cagcTGGACCCCGCCCTAGCGGGCCGGACTTGGACCCCGAACTTCCCCCTCCCCCTACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCAt                                                                  <  7:140533/134‑1 (MQ=255)
                                                         aatCTGGACCCCCCCCTCGCGGTCCGGACTTGCACCCCCACCTGCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCACGGAATTTGCCGGGAGCGGGTTTTCTCGCGATTGCGATGTCTAATACACCGCCGCAt                                                                  <  4:175770/136‑1 (MQ=255)
                                                            cTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTcg                                                               >  7:243933/1‑137 (MQ=255)
                                                                   ccGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTaa                                                        <  2:300338/137‑1 (MQ=255)
                                                                   ccGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTaa                                                        <  6:129879/137‑1 (MQ=255)
                                                                    cGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTaaa                                                       >  3:333639/1‑137 (MQ=255)
                                                                       cttAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGGAAAccc                                                    >  8:321108/3‑137 (MQ=255)
                                                                       ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCAt                                                                  <  8:246647/123‑1 (MQ=255)
                                                                       ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCAt                                                                  >  7:246647/1‑123 (MQ=255)
                                                                       ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAAccc                                                    >  2:9172/1‑137 (MQ=255)
                                                                       ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAAccc                                                    >  1:343408/1‑137 (MQ=255)
                                                                       ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAAccc                                                    >  6:127986/1‑137 (MQ=255)
                                                                       ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAAccc                                                    >  5:337621/1‑137 (MQ=255)
                                                                       ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAAccc                                                    >  2:11570/1‑137 (MQ=255)
                                                                       ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAAccc                                                    >  4:129096/1‑137 (MQ=255)
                                                                         tAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCgg                                                  >  2:277180/1‑137 (MQ=255)
                                                                                  ggACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGgcgc                                         >  8:228106/1‑137 (MQ=255)
                                                                                    aCTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCaa                                       >  6:136757/1‑137 (MQ=255)
                                                                                        gcaCCCCGAACTTCCCCCTCCCAAAACGAAGTCC‑GGGTTTTTTTTTGCCGGTATTTTCCCGGGAGCGTTTTTTCTCGCTAATGCGTTGTCTAATACACCGCCGCATTCGCGGTTAAACCCGGGTCAGGCGAAAACCa                                   <  4:198004/135‑1 (MQ=255)
                                                                                                 ccTTCCCCCTCCCCATACGAAGTCCGGGG‑ATATTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAAAGCGATGTGGAATACACCGCCCCATTCGCGGGTATACCCGGGTCACGCGCAAACCACAAAACGAc                          >  8:90125/2‑137 (MQ=255)
                                                                                                    tCCCCCTCCCCATACGATGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCgg                       <  7:291414/137‑1 (MQ=255)
                                                                                                    tCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCgg                       <  4:247148/137‑1 (MQ=255)
                                                                                                     cccccTCCCCATAGTAAGTCCGGGG‑TTTTTTTTGCCCGTAATTTGCCGGGTGCGGGTTTTCTCGCTAATGCTATGTCTAATACACCGCCGCATTCGCGGGTAAACCGGGGTCAGGCGCAAACCACTAAACGACCGGc                      <  3:194029/137‑1 (MQ=255)
                                                                                                     cccccTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGTAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCTATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGc                      <  2:141380/137‑1 (MQ=255)
                                                                                                     cccccACCCCATACGAAGT‑CGGGGTTTTTTTTTGCCCGTAATTTGCCGGGTGGGGGTTTTCTCGCTATTGGTATGTCGAAAACACCGCCGCATTCGCGGGTAAACCGGGGTCAGGCGAAAACCACTAAACGACCGGc                      <  4:8316/137‑1 (MQ=255)
                                                                                                          tCCCCATTCGATGT‑CGGGGTTTTTTTTTGCCGGTAATTTGCGGGGTGGGGGTTTTCTCGGTATTGCTTTGTGGAATACACGGCCGGATTCGGGGTTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTc                 <  2:184242/137‑1 (MQ=255)
                                                                                                          tCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTc                 <  7:228106/137‑1 (MQ=255)
                                                                                                          cccccATCCGAAGTCCGGGG‑TTTTTTTTGCCGGTAATTTGCCGGGAGCGGGTTTTCTCGCTAATGCGATGTCTAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTc                 <  2:283543/136‑1 (MQ=255)
                                                                                                          cccccATCCGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAAGGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTc                 <  2:56476/136‑1 (MQ=255)
                                                                                                          cccccATAGGATGTCC‑GGGTTTTTTTTTGCCGGTATTTTGCCGGGAGCGGGTTTTCTCGCTAATGCTATGTCTAAACCACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTc                 <  2:303210/136‑1 (MQ=255)
                                                                                                          cccccATACGCGGTCCGGGG‑TTTTTTTTGCCCGTAATTTGCGGGGAGCGGGTTTTCTCGCTAATGGGATGTCGAATACCCCGCCGCATTCGCGGGTAAACCGGGGGCAGGCGCAAACCACTAAACGACCGGCTTGTc                 <  8:83077/136‑1 (MQ=255)
                                                                                                              cATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTcc             >  8:354477/1‑137 (MQ=255)
                                                                                                                  cGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGGATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAACCCACTAACCGACCGGCTTGTCTTCCAGtt         >  3:252183/1‑137 (MQ=39)
                                                                                                                          cggggTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTtccatcca  >  1:157282/2‑137 (MQ=255)
                                                                                                                          cggggTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTATCGAGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCATGAGCAAACACATAAACGACCGGCTTGTCTTCAAGTtcaatcca  >  7:123108/2‑137 (MQ=255)
                                                                                                                              |                                                                                                                                    
AACTGAATGGCTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTTCCATCCA  >  ACB143/120849‑121107

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