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.2% 177.2 / 3.9 43intergenic (+272/+100)polA/engBDNA polymerase I/putative GTP‑binding protein EngB
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base . (29/11);  total (29/14)
Fisher's exact test for biased strand distribution p-value = 9.06e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACTGAATGGCTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTTCCATCCA  >  ACB122/120849‑121107
                                                                                                                              |                                                                                                                                    
aaCTGAATGGCTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGttttttttgcc                                                                                                                            >  7:178892/1‑134 (MQ=255)
             gTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCgg                                                                                                              >  4:90637/1‑137 (MQ=255)
              tGAAGCCCTGCGATGGTCCTATGTTGTGTTAAGTGTTGTCAGATATCTGGACCCCGCCCTTGCGGTCCGGACTTGGACCCCGAACTTCACCCTGCCCTTACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCggg                                                                                                             <  6:192035/137‑1 (MQ=255)
                    ccTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCggg                                                                                                       <  4:276339/137‑1 (MQ=255)
                           tGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTtcta                                                                                                >  7:230599/1‑136 (MQ=255)
                                cTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGaa                                                                                           >  4:126606/1‑137 (MQ=255)
                                 tATGTTGTGTTAAGTGTTGGCATATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAAt                                                                                          <  8:237875/137‑1 (MQ=255)
                                 tATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAAt                                                                                          <  6:128864/137‑1 (MQ=255)
                                  aTGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCTGACTTTGACCCCGAACTTCCCCCACCCCATACGAATTCAGGTGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGTTTTTCTCGCGAATg                                                                                         >  6:258334/1‑137 (MQ=255)
                                     ttgtGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGTAGCGGGTTTTCTCGCGAATGCGa                                                                                      >  2:72772/1‑137 (MQ=255)
                                      ggtgtTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGAt                                                                                     <  3:26310/136‑1 (MQ=255)
                                        tgtTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGt                                                                                   >  2:273924/1‑137 (MQ=255)
                                           tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAAtt                                                                                                                    <  7:42057/101‑1 (MQ=255)
                                           tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAAtt                                                                                                                    >  8:42057/1‑101 (MQ=255)
                                           tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGa                                                                                >  4:118206/1‑137 (MQ=255)
                                           tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGa                                                                                >  4:326161/1‑137 (MQ=255)
                                           tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGa                                                                                >  5:325742/1‑137 (MQ=255)
                                           tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGACCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGa                                                                                >  5:117666/1‑137 (MQ=255)
                                           tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑‑TTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGaa                                                                               >  8:156164/1‑137 (MQ=255)
                                           tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCAGGGGTTTTTTTTGACCGGAATTAGCCGGGAGCGGGTTTTCTCGCGAATACCATGTCGa                                                                                >  6:313627/1‑137 (MQ=255)
                                             aGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAAt                                                                              >  4:305134/1‑137 (MQ=255)
                                             aGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAAt                                                                              >  5:138558/1‑137 (MQ=255)
                                                gtTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATaca                                                                           >  5:137369/1‑137 (MQ=255)
                                                     cAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACAccgcc                                                                      >  7:178801/1‑137 (MQ=255)
                                                          atCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACTAATTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGCATACACCGCCGCAtt                                                                 >  6:161939/1‑137 (MQ=255)
                                                            cTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTcg                                                               >  2:123851/1‑137 (MQ=255)
                                                                       ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAAccc                                                    >  6:89398/1‑137 (MQ=255)
                                                                       ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAAccc                                                    >  6:347827/1‑137 (MQ=255)
                                                                       ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAAccc                                                    >  3:355966/1‑137 (MQ=255)
                                                                       ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTAGCGGGTATAccc                                                    >  1:104879/1‑137 (MQ=255)
                                                                       ccTAGCGGTCCGGACTTGGACCCCGAACTACCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACAGCCGCATTCGCGGGTAAAccc                                                    >  3:15330/1‑137 (MQ=255)
                                                                        cTAGCGGTCCGGACTTGGACCCCGAACATCCCCCTCCCCAAACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCg                                                   <  7:108846/137‑1 (MQ=255)
                                                                              gTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAg                                             >  5:120302/1‑137 (MQ=255)
                                                                                ccGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGgc                                           <  6:325742/137‑1 (MQ=255)
                                                                                 cGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGgcg                                          >  3:273053/1‑137 (MQ=255)
                                                                                        gcaCCCCGAAGTTCCCCCTCCCATTACGATGTCGGGGGTTTTTTTTGCCCGTATTTTGCCGGGAGCGGGTTTTCTCGCTATTGTGTTTTGTATTACACCGCCGCTTTCGCGGTTAAACCCGGGTCAGGCGCAAAcaa                                   <  8:66961/135‑3 (MQ=255)
                                                                                                 aCTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGAc                          >  8:221923/1‑137 (MQ=255)
                                                                                                 aCTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGAc                          >  1:77318/1‑137 (MQ=255)
                                                                                                     cccccTCCCCATCCGAAGTCCGGGTTTTTTTTTGCCCGTATTTTGCCGGGTGGGGGTTTTCTCGCTATTGCGATGTCGAATACACCGCCGAATTCGCGGGTAAACCGGGGTCAGGCGCAAACCACTAAACGACCGGc                      <  5:89398/137‑1 (MQ=255)
                                                                                                     cccccTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGc                      <  5:347827/137‑1 (MQ=255)
                                                                                                     cccccACCCCATACGAAGTCCGGGGTTTTTTTTGCCCGTAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGc                      <  4:238871/137‑1 (MQ=255)
                                                                                                        ccTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTg                   >  2:206820/1‑137 (MQ=255)
                                                                                                          tCCCAATACGAAGTCCGGGGTTTTTTTTGCCCGTAATTTGCCGTGAGCGGTTTTTCTCGCGAATGCGATGTCGAATCCACCGCCTCATTCGCGGGTAAACCCGGGGCAGGCGCAAACCACTAAACGACCGGCTTGTc                 <  4:24204/137‑1 (MQ=255)
                                                                                                          cccccATACGATGTCGGGGTTTTTTTTTGCCCGTATTTTGCGGGGAGCGGGTTTTCTCGCGAATGCGATGTCTAAAACACCGCCGCATGCGCGGGTAAACCGGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTc                 <  1:123851/136‑1 (MQ=255)
                                                                                                                          cggggTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTtccatcca  >  4:223569/2‑137 (MQ=255)
                                                                                                                          cggggTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACACGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTtccatcca  >  1:338940/2‑137 (MQ=255)
                                                                                                                              |                                                                                                                                    
AACTGAATGGCTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTTCCATCCA  >  ACB122/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: