Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,151,829	C→T	K52K (AAG→AAA)	exbB ←	Ton complex subunit ExbB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,151,829	0	C	T	97.0%	103.3 / ‑4.4	33	K52K (AAG→AAA)	exbB	Ton complex subunit ExbB
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (16/16);  total (16/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

GTTCCAGCTCATTCTGCGCTTCATTGAGCAAATGCAGGCTTAAGCTTTTGCTACCAAAATCAGCGGCGATATCGTTGGCCTGGTTTAAGGAACGCGCTTCAGCCAGCAGTTGCTGCTCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCACGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTC  >  NC_000913/3151696‑3151964                                                                                                                                      |                                                                                                                                        gTTCCAGCTCATTCTGCGCTTCATTGAGCAAATGCAGGCTTAAGCTTTTGCTACCAAAATCAGCGGCGATATCGTTGGCCTGGTTTAAGGAACGCGCTTCAGCCAGCAGTTGCTGCTCGCGCTTAAGGCGACGTTTCTGa                                                                                                                                   >  1:197855/1‑140 (MQ=255)            ttCTGCGCTTCATTGAGCAAATGCAGGCTTAAGCTTTTGCTACCAAAATCAGCGGCGATATCGTTGGCCTGGTTTAAGGAACGCGCTTCAGCCAGCAGTTGCTGCTCGCGCTTAAGGCGACGTTTCTGATTGAAGAActc                                                                                                                        <  1:155164/140‑1 (MQ=255)               tGCGCTTCATTGAGCAAATGCAGGCTTAAGCTTTTGCTACCAAAATCAGCGGCGATATCGTTGGCCTGGTTTAAGGAACGCGCTTCAGCCAGCAGTTGCTGCTCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTAc                                                                                                                     <  2:28387/140‑1 (MQ=255)                      cATTGAGCAAATGCAGGCTTAAGCTTTTGCTACCAAAATCAGCGGCGATATCGTTGGCCTGGTTTAAGGAACGCGCTTCAGCCAGCAGTTGCTGCTCGCGCTTAAGGCGACGTTTCTGATTGAAGAACTCTACGCTCTTa                                                                                                              >  2:57069/1‑140 (MQ=255)                      cATTGAGCAAATGCAGGCTTAAGCTTTTGCTACCAAAATCAGCGGCGATATCGTTGGCCTGGTTTAAGGAACGCGCTTCAGCCAGCAGTTGCTGCTCGCGCTTAAGGCGACGTTTCTGATTGAAGAACTCTACGCTCTTa                                                                                                              <  2:136480/140‑1 (MQ=255)                       aTTGAGCAAATGCAGGCTTAAGCTTTTGCTACCAAAATCAGCGGCGATATCGTTGGCCTGGTTTAAGGAACGCGCTTCAGCCAGCAGTTGCTGCTCGCGCTTAAGGCGACGTTTCTGATTGAAGAACTCTACGCTCTTAc                                                                                                             >  1:236975/1‑140 (MQ=255)                           aGCAAATGCAGGCTTAAGCTTTTGCTACCAAAATCAGCGGCGATATCGTTGGCCTGGTTTAAGGAACGCGCTTCAGCCAGCAGTTGCTGCTCGCGCTTAAGGCGACGTTTCTGATTGAAGAACTCTACGCTCTTACTgaa                                                                                                         >  1:52909/1‑140 (MQ=255)                           aGCAAATGCAGGCTTAAGCTTTTGCTACCAAAATCAGCGGCGATATCGTTGGCCTGGTTTAAGGAACGCGCTTCAGCCAGCAGTTGCTGCTCGCGCTTAAGGCGACGTTTCTGACTGAAGAACTCTACGCTCTTACTgaa                                                                                                         <  2:209623/140‑1 (MQ=255)                                     ggCTTAAGCTTTTGCTACCAAAATCAGCGGCGATATCGTTGGCCTGGTTTAAGGAACGCGCTTCAGCCAGCAGTTGCTGCTCGCGCTTAAGGCGACGTTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGcc                                                                                               <  2:236975/140‑1 (MQ=255)                                       cTTAAGCTTTTGCTACCAAAATCAGCGGCGATATCGTTGGCCTGGTTTAAGGAACGCGCTTCAGCCAGCAGTTGCTGCTCGCGCTTAAGGCGACGTTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCa                                                                                             <  1:57069/140‑1 (MQ=255)                                         tAAGCTTTTGCTACCAAAATCAGCGGCGATATCGTTGGCCTGGTTTAAGGAACGCGCTTCAGCCAGCAGTTGCTGCTCGCGCTTAAGGCGACGTTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAgg                                                                                           <  2:197855/140‑1 (MQ=255)                                                ttGCTACCAAAATCAGCGGCGATATCGTTGGCCTGGTTTAAGGAACGCGCTTCAGCCAGCAGTTGCTGCTCGCGCTTAAGGCGACGTTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTAc                                                                                    >  2:265350/1‑140 (MQ=255)                                                ttGCTACCAAAATCAGCGGCGATATCGTTGGCCTGGTTTAAGGAACGCGCTTCAGCCAGCAGTTGCTGCTCGCGCTTAAGGCGACGTTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTAc                                                                                    >  2:267349/1‑140 (MQ=255)                                                 tGCTACCAAAATCAGCGGCGATATCGTTGGCCTGGTTTAAGGAACGCGCTTCAGCCAGCAGTTGCTGCTCGCGCTTAAGGCGACGTTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACg                                                                                   <  1:110989/140‑1 (MQ=255)                                                                    gATATCGTTGGCCTGGTTTAAGGAACGCGCTTCAGCCAGCAGTTGCTGCTCGCGCTTAAGGCGACGTTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCaa                                                                <  2:52909/140‑1 (MQ=255)                                                                            tGGCCTGGTTTAAGGAACGCGCTTCAGCCAGCAGTTGCTGCTCGCGCTTAAGGCGACGTTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCACg                                                        <  2:18447/140‑1 (MQ=255)                                                                            tGGCCTGGTTTAAGGAACGCGCTTCAGCCAGCAGTTGCTGCTCGCGCTTAAGGCGACGTTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCACg                                                        <  2:35688/140‑1 (MQ=255)                                                                                   gTTTAAGGAACGCGCTTCAGCCAGCAGTTGCTGCTCGCGCTTAAGGCGACGTTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCACGCACTTaa                                                 >  2:271865/1‑140 (MQ=255)                                                                                             cgcgCTTCAGCCAGCAGTTGCTGCTCGCGCTTAAGGCGACGTTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCACGCACTTAACGACAATATc                                       >  2:322889/1‑140 (MQ=255)                                                                                             cgcgCTTCAGCCAGCAGTTGCTGCTCGCGCTTAAGGCGACGTTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCACGCACTTAACGACAATATc                                       >  2:36681/1‑140 (MQ=255)                                                                                             cgcgCTTCAGCCAGCAGTTGCTGCTCGCGCTTAAGGCGACGTTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCACGCACTTAACGACAATATc                                       >  2:91376/1‑140 (MQ=255)                                                                                             cgcgCTTCAGCCAGCAGTTGCTGCTCGCGCTTAAGGCGACGTTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCACGCACTTAACGACAATATc                                       >  2:237984/1‑140 (MQ=255)                                                                                             cgcgCTTCAGCCAGCAGTTGCTGCTCGCGCTTAAGGCGACGTTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCACGCACTTAACGACAATATc                                       >  1:388223/1‑140 (MQ=255)                                                                                                gcTTCAGCCAGCAGTTGCTGCTCGCGCTTAAGGCGACGTTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCACGCACTTAACGACAATATCGGc                                    <  1:267349/140‑1 (MQ=255)                                                                                                   tCAGCCAGCAGTTGCTGCTCGCGCTTAAGGCGACGTTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCACGCACTTAACGACAATATCGGCGTg                                 <  1:237984/140‑1 (MQ=255)                                                                                                      gccagcAGTTGCTGCTCGCGCTTAAGGCGACGTTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCACGCACTTAACGACAATATCGGCGTGCTg                              <  1:91376/140‑1 (MQ=255)                                                                                                              ttGCTGCTCGCGCTTAAGGCGACGTTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCACGCACTTAACGACAATATCGGCGTGCTGatacatac                      >  1:385016/1‑140 (MQ=255)                                                                                                              ttGCTGCTCGCGCTTAAGGCGACGTTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCACGCACTTAACGACAATATCGGCGTGCTGat                            <  2:128093/134‑1 (MQ=255)                                                                                                              ttGCTGCTCGCGCTTAAGGCGACGTTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCACGCACTTAACGACAATATCGGCGTGCTGat                            >  1:128093/1‑134 (MQ=255)                                                                                                                gctgctCGCGCTTAAGGCGACGTTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCACGCACTTAACGACAATATCGGCGTGCTGATACATAccc                    >  2:375731/1‑140 (MQ=255)                                                                                                                     tCGCGCTTAAGGCGACGTTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCACGCACTTAACGACAATATCGGCGTGCTGataca                         <  2:317412/130‑1 (MQ=255)                                                                                                                     tCGCGCTTAAGGCGACGTTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCACGCACTTAACGACAATATCGGCGTGCTGataca                         >  1:317412/1‑130 (MQ=255)                                                                                                                                  gACGTTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCACGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTc  <  2:385016/140‑1 (MQ=255)                                                                                                                                      |                                                                                                                                        GTTCCAGCTCATTCTGCGCTTCATTGAGCAAATGCAGGCTTAAGCTTTTGCTACCAAAATCAGCGGCGATATCGTTGGCCTGGTTTAAGGAACGCGCTTCAGCCAGCAGTTGCTGCTCGCGCTTAAGGCGACGCTTCTGATTGAAGAACTCTACGCTCTTACTGAAGAAGATTGCCCAGGTGACTACGGAGGCCAAAATAAGCCCAATCATCACGCACTTAACGACAATATCGGCGTGCTGATACATACCCCAGACGGAAAGGTCCGTC  >  NC_000913/3151696‑3151964

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