Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,311,303 A→G E72G (GAA→GGA)  tonB → Ton complex subunit TonB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,311,3030AG84.6% 29.5 / ‑5.8 13E72G (GAA→GGA) tonBTon complex subunit TonB
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base G (10/1);  total (12/1)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGAAAAAGGTACAGGAGCAGCCAAAACGCGATGTCA  >  NC_000913/1311180‑1311431
                                                                                                                           |                                                                                                                                
cGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCTCAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAAcc                                                                                                                         >  2:314942/1‑133 (MQ=255)
 ggTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGGGGTAGAACCAGGACCGGAACCt                                                                                                                        >  1:231322/1‑133 (MQ=255)
                     tACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAAGCTGGGCGGATCCCAGAACCGcc                                                                                                    >  1:142984/1‑133 (MQ=255)
                              ctgcgcAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGGACCGGAACCTGAGCCGGTCCGCGGACCGCCCAAAGAAGc                                                                                           >  1:202774/1‑133 (MQ=255)
                                 cgcAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGGACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCAcc                                                                                        >  2:358953/1‑133 (MQ=255)
                                 cgcAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGCAGAGCCAGGACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCAcc                                                                                        >  2:46877/1‑133 (MQ=255)
                                                tCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGGACCGGAACCTGAGCCTATCCCCGAACCGACAAAAGAAGCACTGGTGCTCATTGaaaa                                                                         >  2:180853/1‑133 (MQ=255)
                                                                               gCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGGACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAg                                          >  2:356594/1‑133 (MQ=255)
                                                                                 cacaAGCCGTCCAGCAGCCACCGGTGCCAGTGTTAGAGCCACGACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGcc                                        <  2:142984/133‑1 (MQ=255)
                                                                                      gCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGGACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGGAAAGCCGAAGCCGAAATCTAAGCCAAAACCTAAGCCGGTGa                                   >  2:171408/1‑133 (MQ=255)
                                                                                         gTTCAGCCGTAACCGGAGCCGGTGGTTGCGCCTGGACCGGATCCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGaaaa                                <  2:231243/133‑1 (MQ=255)
                                                                                          ttCAGCCGCCCCCGGCGCCGGTGGTAGAGCCAAGACCGGAAGCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGaaaaa                               <  2:202774/133‑1 (MQ=255)
                                                                                          ttCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGGACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAACCCGAAGCCAAAACCTAAGCCAAAACCGAAGCCGGTGaaaaa                               >  2:182736/1‑133 (MQ=255)
                                                                                                             ggtggtAGAGCCAGGACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGAAAAAGGTACAGGAGCAGCCaaaa            >  2:438796/1‑133 (MQ=255)
                                                                                                                       ccAGGACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGAAAAAGGTACAGGAGCAGCCAAAACGCGATGTCa  >  1:37111/1‑133 (MQ=255)
                                                                                                                           |                                                                                                                                
CGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGAAAAAGGTACAGGAGCAGCCAAAACGCGATGTCA  >  NC_000913/1311180‑1311431

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: