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.95.0% 163.6 / 3.9 41intergenic (+272/+100)polA/engBDNA polymerase I/putative GTP‑binding protein EngB
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base . (31/7);  total (33/8)
Fisher's exact test for biased strand distribution p-value = 3.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGAGTAACTGAATGGCTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTTCCATCC  >  ACB143/120843‑121106
                                                                                                                                    |                                                                                                                                   
cAGAGTAACTGAATGGCTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCAGACGAAGTCCGGGTgtttt                                                                                                                                 >  4:165504/1‑137 (MQ=255)
          gAATGGCTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGa                                                                                                                      >  1:22659/1‑137 (MQ=255)
                     gAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTCGACCCCGCCCTAGCGGTCCGGACTTGGAACCCGAACTTCCCCCTCACCATACGAAGTCCGGGGTTTTTTTTGACCGGAATTTGCAGGGa                                                                                                           >  5:223512/1‑137 (MQ=255)
                                aTGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTtct                                                                                                >  1:201469/1‑137 (MQ=255)
                                       tATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCAGGGAGCGGGTTTTCTCGCGAAt                                                                                         >  8:290258/1‑137 (MQ=255)
                                       tATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGACCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAAt                                                                                         >  2:295270/1‑137 (MQ=255)
                                             gtgtTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATg                                                                                   >  3:165584/1‑137 (MQ=255)
                                                 tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGa                                                                               >  1:208205/1‑137 (MQ=255)
                                                 tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGa                                                                               >  1:263280/1‑137 (MQ=255)
                                                 tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGa                                                                               >  2:183412/1‑137 (MQ=255)
                                                 tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGa                                                                               >  7:137630/1‑137 (MQ=255)
                                                 tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGa                                                                               >  5:239320/1‑137 (MQ=255)
                                                 tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTGTTTGACCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTAGa                                                                               >  8:300951/1‑137 (MQ=255)
                                                 tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGACCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGa                                                                               >  2:129614/1‑137 (MQ=255)
                                                 tAAGTGTAGGCAGATATCTGGACCCCGCCCTAGCTGTCCGGACTTGGACCAAGAACTTCCCCCTCCCAATACGAAGTCAGTGGTTTTTTTTTCCAGGAATTTTCCGCGAGCTGGTTTTCTCGAGAATGCGATGtggt                                                                               >  5:182445/1‑134 (MQ=255)
                                                           cAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACAccgcc                                                                     >  5:199961/1‑137 (MQ=255)
                                                           cAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGGTTTCTCGCGAATGCGATGTCGAATACAccgac                                                                     >  7:287268/1‑135 (MQ=255)
                                                            aGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACAccgccg                                                                    >  2:215567/1‑137 (MQ=255)
                                                               tatCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCTGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCAt                                                                 <  4:295902/137‑1 (MQ=255)
                                                                  cTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTcg                                                              >  3:125680/1‑137 (MQ=255)
                                                                         ccGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTaa                                                       <  8:171052/137‑1 (MQ=255)
                                                                             ccTAGCGGTCCGGACTTGGACCCCGACCTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAAccc                                                   >  8:180101/1‑137 (MQ=255)
                                                                             ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAAccc                                                   >  5:230124/1‑137 (MQ=255)
                                                                             ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTATTTTTGCCCGGAATTTGCCGGCAGCGGGTTTTCTCGCGAATGCGATCTCGAATACACCGCCCCATTCGCGGGTAAAccc                                                   >  3:254292/1‑137 (MQ=255)
                                                                             ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCAACATACGAAGTCCGGGGTTTTTTTTGACCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATAAACCGCCGCATTAGCGGGTAAAccc                                                   >  7:122602/1‑137 (MQ=255)
                                                                              cTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAAt                                                                                         >  8:11861/1‑98 (MQ=255)
                                                                              cTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAAt                                                                                         <  7:11861/98‑1 (MQ=255)
                                                                               tAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCgg                                                 >  7:13735/1‑137 (MQ=255)
                                                                                     tCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAgg                                           <  7:229381/137‑1 (MQ=255)
                                                                                      ccGGACTTGGACCCCGAACTTCCCCCTCCCCCTACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGgc                                          <  4:89802/137‑1 (MQ=255)
                                                                                          aCTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCaa                                      >  3:8307/1‑137 (MQ=255)
                                                                                          aCTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCaa                                      >  2:291490/1‑137 (MQ=255)
                                                                                          aCTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCAGGAGCGGGTTTTATCTCGAATGCGATGTCGAAGACACAGCAGCATTAGCGGGTAAAACAGGGTCAGGCGCaaa                                     >  5:2953/1‑137 (MQ=255)
                                                                                                  cccGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTaa                              >  3:157469/1‑137 (MQ=255)
                                                                                                  cccGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTaa                              >  5:155688/1‑137 (MQ=255)
                                                                                                       aCTTCCCCCTCCCCATACGAAGTCCGGGGGTTTTTTTTGCCCGGAATTTGCCGGGGCCGGGTTTTCTCGCGATGGCGATGTCGAATACCCCGCCGCATTCCCGGGTAACCCCGGGGCAGGCCCAACCCCATAAACgc                          >  2:255870/1‑136 (MQ=255)
                                                                                                       aCTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGAc                         >  1:217998/1‑137 (MQ=255)
                                                                                                         ttCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCg                       >  3:52934/1‑137 (MQ=255)
                                                                                                             cccTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCAt                                                                 <  8:60246/91‑1 (MQ=255)
                                                                                                             cccTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCAt                                                                 >  7:60246/1‑91 (MQ=255)
                                                                                                                    aaTACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTcc            <  6:230124/136‑1 (MQ=255)
                                                                                                                              tccgggTTTTTTTTTGCCCGTAATTTGCGGGGAGGGGTTTTTCTCGCGAATGCTATGTCGAATACACCGCCGCTTTCGCTGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTtccatc   <  6:62789/134‑1 (MQ=255)
                                                                                                                               ccggggTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTtccatcc  >  4:200988/3‑137 (MQ=255)
                                                                                                                               ccggggTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACACGGGTCAGGAGAAAACAACTAAACGACCGGCTTGTCTTCCAGTtccatcc  >  6:157525/3‑137 (MQ=255)
                                                                                                                                    |                                                                                                                                   
CAGAGTAACTGAATGGCTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTTCCATCC  >  ACB143/120843‑121106

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: