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

AATGGCTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTTCCATCCA  >  ACB122/120854‑121107
                                                                                                                         |                                                                                                                                    
aaTGGCTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGaa                                                                                                                      >  2:92009/1‑137 (MQ=255)
     cTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTgg                                                                                                                 >  1:112965/1‑136 (MQ=255)
     cTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGc                                                                                                                 >  8:65901/1‑137 (MQ=255)
        gTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCgg                                                                                                              >  4:202552/1‑137 (MQ=255)
                 ttCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGtt                                                                                                     >  5:110297/1‑137 (MQ=255)
                     aTGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTtct                                                                                                 >  1:191860/1‑137 (MQ=255)
                            tATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAAt                                                                                          <  7:65901/137‑1 (MQ=255)
                             aTGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTATCGCGAATg                                                                                         >  7:90196/1‑137 (MQ=255)
                                      tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCTGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTTCCGGGAGCGGTTTTTCTAGCGAATGCGATTTCGa                                                                                >  2:52843/1‑137 (MQ=255)
                                      tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAAt                                                                                          >  2:46387/1‑127 (MQ=255)
                                      tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAAt                                                                                          <  7:200791/127‑1 (MQ=255)
                                      tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAAt                                                                                          >  8:200791/1‑127 (MQ=255)
                                      tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAAt                                                                                          <  1:46387/127‑1 (MQ=255)
                                      tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGa                                                                                >  2:17354/1‑137 (MQ=255)
                                      tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGa                                                                                >  8:171950/1‑137 (MQ=255)
                                      tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGa                                                                                >  8:209199/1‑137 (MQ=255)
                                        aGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAAt                                                                              >  8:248589/1‑137 (MQ=255)
                                        aGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAAt                                                                              >  5:184569/1‑137 (MQ=255)
                                        aGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAAt                                                                              >  3:116830/1‑137 (MQ=255)
                                             tGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATTCTATGTCGAATACAcc                                                                         >  3:30709/1‑137 (MQ=255)
                                             tGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACAcc                                                                         >  6:209333/1‑137 (MQ=255)
                                             tGGCAGATATCTGGAACCCGCCCTAGCGGTCCGGACTTGGACCCCGAAATTCCCCCTCCCCATACGAATTCCGGGGTTTTTTTTGCCCGGAATTTGACGGTAGCGGGTTTTCTCGCGAATGCGATGTCGAATACAcc                                                                         >  6:148285/1‑137 (MQ=255)
                                                  gATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACcgccgc                                                                    >  2:229059/1‑137 (MQ=255)
                                                    tatCTGGACCCCGCCCTAGCGGGCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCAt                                                                  <  2:199963/137‑1 (MQ=255)
                                                              ccGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTaa                                                        <  2:191860/137‑1 (MQ=255)
                                                                  ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAAccc                                                    >  1:46811/1‑137 (MQ=255)
                                                                   cTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCg                                                   >  7:34401/1‑137 (MQ=255)
                                                                   cTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCg                                                   >  5:12329/1‑137 (MQ=255)
                                                                   cTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTACGGGGTTTTTTTTGCCAGGAATTTGCCGGGAGCGGGTTTTCTCGCGATTGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCg                                                   >  7:272913/1‑137 (MQ=39)
                                                                    tAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTGTTTTTTGCCCGGAATTTGCCGGGCGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAGCCCgg                                                  >  1:128610/1‑137 (MQ=255)
                                                                           ccGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGgc                                           <  7:248589/137‑1 (MQ=255)
                                                                            cGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGgcg                                          >  1:66250/1‑137 (MQ=255)
                                                                              gACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGTAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCa                                        <  1:92009/137‑1 (MQ=255)
                                                                              gACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGACCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGAAGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCCGgcgcc                                        >  1:150293/1‑136 (MQ=39)
                                                                               aCTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCaa                                       >  5:119670/1‑137 (MQ=255)
                                                                                      ccccGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCGCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCCCCGCATTAGCGGGTAAACCCGGGTCAGGCGCAACCCACTa                                >  2:127093/1‑137 (MQ=255)
                                                                                        ccGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTaaa                              <  8:298878/137‑1 (MQ=255)
                                                                                        ccGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTaaa                              <  8:294712/137‑1 (MQ=255)
                                                                                            aCTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGAc                          >  5:283161/1‑137 (MQ=255)
                                                                                              ttCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCg                        >  1:146721/1‑137 (MQ=255)
                                                                                                cccccTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGc                      <  6:283161/137‑1 (MQ=255)
                                                                                                cccccTACCAATACTAAGTCGTGGTTTTTTTTTGCCGGTAATTTGCGGGGAGCGGGTTTTATGGCTAATGCTATGTCTAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGc                      <  6:253639/137‑1 (MQ=255)
                                                                                                     accccATACGAAGTCCGGGGTTTTTTTTGCCCGTATTTTGCGGGGAGCGGGTTTTCTCGCGATTGCGTTGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTc                 <  1:203853/136‑1 (MQ=255)
                                                                                                     accccATACGAAGTCCGGGGTTTTTTTTGCCCGTAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCTAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTc                 <  3:134310/136‑1 (MQ=255)
                                                                                                           tACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAg           >  5:20690/1‑137 (MQ=255)
                                                                                                               aaGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTtcc       >  4:207925/1‑137 (MQ=255)
                                                                                                                  tCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGTGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTtccatc    <  8:52259/137‑1 (MQ=255)
                                                                                                                     cggggTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGCATGCGATGTCGAATACACCGCAGCATTCGCGGGTAAACCCGGGTCAGGCGCACACCACTAACCGACCGGCTTGTCTTCCAGTtccatcca  >  8:223351/2‑137 (MQ=255)
                                                                                                                         |                                                                                                                                    
AATGGCTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTTCCATCCA  >  ACB122/120854‑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: