Predicted mutation
evidence seq id position mutation annotation gene description
RA ACB122 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
*ACB122120,9750T.95.5% 179.5 / 2.8 45intergenic (+272/+100)polA/engBDNA polymerase I/putative GTP‑binding protein EngB
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base . (28/14);  total (28/17)
Fisher's exact test for biased strand distribution p-value = 1.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTGAATGGCTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTTCCATCCAGTC  >  ACB122/120850‑121110
                                                                                                                             |                                                                                                                                       
aCTGAATGGCTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGttttttttgcc                                                                                                                               <  8:300990/136‑4 (MQ=255)
aCTGAATGGCTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGttttttttgc                                                                                                                                <  8:1717/135‑3 (MQ=255)
    aaTGGCTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGttttttttggacgg                                                                                                                            >  3:199496/1‑129 (MQ=255)
               aaGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAACTCCGGGG‑GGTTTTTTGCCCGGCATTTGCCgg                                                                                                                 >  7:231666/1‑134 (MQ=255)
               aaGCCCTTCGATGCTCCTATGTTGAGTGAAGTGTTGGCAGATATCTGGACCCCGCCCTATCGGTCCGGACTTGGACCCCGAACTTCCCACTCCCCATACGAAGTCCGGGG‑GTTTTTTTGCCCGGAATTTGCCgg                                                                                                                 >  7:58789/1‑134 (MQ=255)
                       cGATGCTCCTATGTTGGGTTAAGTGTTGGCAGATATCGGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGt                                                                                                         <  7:275762/134‑1 (MQ=255)
                        gaagCTCCTATCTTCTGTTAAGTCTATCCAGAACTCTGGACCCCGCCCTAGCGGTCCGGACGTGGAGCCCGCACTTCCGCCTCCCCCTACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGt                                                                                                         <  2:155413/130‑1 (MQ=255)
                         aTGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGtttt                                                                                                      >  1:49610/1‑135 (MQ=255)
                           gCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTtct                                                                                                    <  1:223357/135‑1 (MQ=255)
                               cTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGATCGGGTTTTCTcgct                                                                                                >  7:245709/1‑134 (MQ=255)
                               cTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTcgagt                                                                                               >  5:330974/1‑133 (MQ=255)
                               cTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGa                                                                                               <  4:212493/136‑1 (MQ=255)
                                  tgttgtGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAAt                                                                                             <  8:81395/135‑1 (MQ=255)
                                    ttgtGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGc                                                                                           <  7:201951/135‑1 (MQ=255)
                                    ttgtGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACCTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGc                                                                                           <  3:170686/135‑1 (MQ=255)
                                     tgtgTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGc                                                                                           <  1:189300/134‑1 (MQ=255)
                                         ttAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATAAGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTc                                                                                     >  5:145854/1‑136 (MQ=255)
                                          tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAAt                                                                                             >  4:314396/1‑127 (MQ=255)
                                          tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAAt                                                                                             <  3:314396/127‑1 (MQ=255)
                                          tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGt                                                                                      >  2:97258/1‑134 (MQ=255)
                                          tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCACGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTc                                                                                     >  4:269660/1‑135 (MQ=255)
                                            aGTGTTGGCAGATATCTGGCCCCCGCCCTAGCGGTCCGGACTTGGACCCCGCACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGa                                                                                   >  1:183094/1‑135 (MQ=255)
                                            aGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGaa                                                                                  >  8:203894/1‑136 (MQ=255)
                                            aGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGaa                                                                                  >  3:291366/1‑136 (MQ=255)
                                            aGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGa                                                                                   >  1:298023/1‑135 (MQ=255)
                                            aGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCACCCCAAACGAAGTACGGGG‑TTTTTTTTGCCCGGAATTAGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCg                                                                                    >  2:101819/1‑134 (MQ=255)
                                                    cAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACAccgc                                                                          >  5:158556/1‑136 (MQ=255)
                                                    cAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TGTTTTTTGCCCGGAATTTGCCGGGAGCGGGTGTTCTCGCGATTGCGCTGTCGAATACCccgc                                                                          >  7:126810/1‑136 (MQ=255)
                                                      gATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAAt                                                                                             <  6:311700/115‑1 (MQ=255)
                                                      gATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAAt                                                                                             >  5:311700/1‑115 (MQ=255)
                                                           cTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGATGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATCTCGAATACACCGCCGCATTc                                                                   >  7:221552/1‑136 (MQ=255)
                                                             ggACCCCCCCCTCGCGGGCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACTAAGTCC‑GGGTTTTTTTTTGCCCGGAATTTGCCGGTAGCGGGTTTTCTCGCGAATGCGATGTCTAAAACACCGCCGCATTcgc                                                                 <  4:212791/136‑1 (MQ=255)
                                                                      ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGCGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGGAAAcc                                                        >  5:229695/1‑136 (MQ=255)
                                                                      ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTaagc                                                         >  1:231617/1‑133 (MQ=255)
                                                                      ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAAcc                                                        >  7:185674/1‑136 (MQ=255)
                                                                      ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCAATCGCGGGTAAAcc                                                        >  5:301418/1‑136 (MQ=255)
                                                                      ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑‑TTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAAccc                                                       >  5:141749/1‑136 (MQ=255)
                                                                          gCGGGCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCgg                                                     <  3:304719/135‑1 (MQ=255)
                                                                                 ggACTTGGACCCCGAACTTCGCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGgc                                              <  2:23491/135‑1 (MQ=255)
                                                                                 ggACTTGGACCCCGAACTTCCCCCTCCCCATCCGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGgc                                              <  1:24446/135‑1 (MQ=255)
                                                                                   aCTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCa                                           <  6:158556/136‑1 (MQ=255)
                                                                                               aaCTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCTATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAACCCACTAAAc                                >  8:116435/1‑135 (MQ=255)
                                                                                                aCTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGa                              >  6:232575/1‑136 (MQ=255)
                                                                                                aCTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGa                              >  4:292168/1‑136 (MQ=255)
                                                                                                 cTTCCCCCTCCCCATACGAAGT‑CGGGGGTTTTTTTTGCCGGTATTTTGCCGGGAGGGGGTTTTCTCGCTAATGCGATGTCTAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACAACTAAACGa                              <  2:231617/135‑1 (MQ=255)
                                                                                                       ccTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCACACCACTAAACGACCGGCtt                       >  3:93587/1‑136 (MQ=255)
                                                                                                            ccATACTACGT‑CGGGGTTTTTTTTTCCCCGTATTTTGCGGGGAGCGGTTTTGCCCGCGATTGGGATGTCTAATACACCGCCGCATTCGGGTTTAACCCCGGGTCAGGCGCAAACCACAAAACGACCGGCTTGTCtt                  <  8:125024/136‑1 (MQ=255)
                                                                                                            ccATACGAAGTCCGGGG‑TTTTTTTTGCACGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCCTTCGCGGGTAACCCCGGGTCAGGCGCAAACCACTAAACGCCCGGCTTGTCtg                  >  6:269587/1‑135 (MQ=255)
                                                                                                                        ccggggTTTTTTTTGCCGGTAATTTGCCGGGAGCGGGTTTTCTCGCAAATGCTATGTCTAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTtccatc       <  7:203894/134‑1 (MQ=255)
                                                                                                                        ccggggTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTtccatc       >  8:259171/3‑136 (MQ=255)
                                                                                                                             tttttttttGCCGGGATTTTGCGGGGAGTGGGTTTTCTGGCGAATGCTATGTCGATTCCCCGCCCGAATTGGCGGTTAAACGGGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTTCCATCCAGTc  <  8:317904/136‑1 (MQ=255)
                                                                                                                             |                                                                                                                                       
ACTGAATGGCTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTTCCATCCAGTC  >  ACB122/120850‑121110

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: ‑