Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	460,772	T→G	S629A (TCG→GCG)	lon →	Lon protease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	460,772	0	T	G	96.7%	77.2 / NA	30	S629A (TCG→GCG)	lon	Lon protease
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base G (16/13);  minor base A (1/0);  total (17/13)
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

GCGTTTCGACTATGGTCGCGCGGATAACGAAAACCGTGTCGGTCAGGTAACCGGTCTGGCGTGGACGGAAGTGGGCGGTGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGA  >  NC_000913/460633‑460904                                                                                                                                            |                                                                                                                                     gCGTTTCGACTATGGTCGCGCGGATAACGAAAACCGTGTCGGTAAGGTAACCGGTCTGGCGTGGACGGAAGTGGGCGGTGACTTGCTGACCATTGAAACAGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACAGGTACGc                                                                                                                                   >  2:507895/1‑143 (MQ=255)                tCGCGCGGATAACGAAAACCGTGTCGGTCAGGTAACCGGTCTGGCGTGGACGGAAGTGGGCGGTGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTACAGGCAAGCTGACCTATACCGGTGGGCTCGGCGAAGtgctgc                                                                                                                    >  2:134699/1‑143 (MQ=255)                     aggATAACGAAAACCGTGTCGGTCAGGTAACCGGTCTGGCGTGGACGGAAGTGGGCGGTGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTGCGCTCGGCGAAGTGATGCAGGAg                                                                                                               >  2:609997/2‑143 (MQ=255)                         tAACGAAAACCGTGTCGGTCAGGTAACCGGTCTGGCGTGGACGGAAGTGGGCGGTGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTGCGCTCGGCGAAGTGATGCAGGAGTCCa                                                                                                           <  1:395818/143‑1 (MQ=255)                            cGAAAACCGTGTCGGTCAGGTAACCGGTCTGGCGTGGACGGAAGTGGGCGGTGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTGCGCTCGGCGAAGTGATGCAg                                                                                                                  >  2:627276/1‑133 (MQ=255)                            cGAAAACCGTGTCGGTCAGGTAACCGGTCTGGCGTGGACGGAAGTGGGCGGTGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTGCGCTCGGCGAAGTGATGCAg                                                                                                                  <  1:627276/133‑1 (MQ=255)                             gAAAACCGTGTCGGTCAGGTAACCGGTCTGGCGTGGACGGAAGTGGGCGGTGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACCGACCTATACCGGTGCGCTCGGCGAAGTGAGGCAGGAGTCCATTCa                                                                                                       >  2:786818/1‑143 (MQ=255)                                aaCCGTGTCGGTCAGGTAACCGGTCTGGCGTGGACGGAAGTGGGCGGTGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTGCGCTCGGGGAAGTGATGCAGGAGTCCATTGAggc                                                                                                    >  2:230073/1‑143 (MQ=255)                                   cgggtCGGTCAGGTGCGCGGTCTGGCGTGGGCGGAAGTGGGCGGTGACTTGCTGACCATTGAATCCGCAGGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTGCGCTCGGCGAAGTGATGCAGGAGTCCATTCAggcggc                                                                                                 <  2:167915/140‑1 (MQ=255)                                               gTAACCGGTCTGGCGTGGACGGAAGTGGGCGGTGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTGCGCTCGGCGAAGTGCTGCAGGAGTCCATTCAGGCGGCGTTAACggtggt                                                                                     >  2:810290/1‑143 (MQ=255)                                                              tGGACGGAAGTGGGCGGTGACTTGCTGACCATTGAAACCGCACTGGTTCCGGGTAAAGGCAAACTGACCTATACCGGTGCGCTCGGCGAAGTGATGCAGGAGTCCAGCCGGGCGGCGCGAACGGTGGTTCGTGCGCGGGCGGa                                                                      >  2:748067/1‑143 (MQ=255)                                                                gACGGAAGTGGGCGGTGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTGCGCTCGGCGAAGTGATGCAGGAGGCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGaaa                                                                    >  1:169736/1‑143 (MQ=255)                                                                   tggaGTGGGCGGTGACTTGCTGAGCGTGTGGACCGGGGGGGTCGCGGTTAAATGCAAACTGTCCTATACCGGTGCGCCCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGAGTTAACGGTGGTTCGTGTGCGTGCGGAAAAAc                                                                 <  2:150167/140‑1 (MQ=255)                                                                     aaGTGGGCGGTGACTTGCTTACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTGCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTg                                                               <  1:134699/143‑1 (MQ=255)                                                                              gTGACTTGCTGACCATTGAACCCGCCTGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTGCGCCCGGCGACGTGATGCACGAGTCCATTCAGGCGGCGTTCACGGTGGTCCGTGCGCGCGCGCGAAAACCGCGGATCGAc                                                      >  2:60609/1‑143 (MQ=255)                                                                               tGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTGCGCTCGGCGAAGTGATGCAGGAGTCAATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAGAGACTGGGGGTCAAcc                                                     >  1:817923/1‑143 (MQ=255)                                                                                   ttGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTGCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAGAACTGGGAAGGCACCCTGa                                                 >  2:249797/1‑143 (MQ=255)                                                                                           cATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTGCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACg                                         >  1:599357/1‑143 (MQ=255)                                                                                           cATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTGCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTCCGTGCGCGTGCGGCACAACTGGGGCTGAACCCTGATTTGTACg                                         >  1:358220/1‑143 (MQ=255)                                                                                                aaaCCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCCGAGCGCTCGGCGGCGTGATGCAGGACACCATTCAGGCCGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGCGTCAACCCGGATTTTTAAGaacca                                    >  2:458905/1‑140 (MQ=255)                                                                                                   ccGCAGGGGTTGCGGGTAAAGGCAAACCGACCGACACCGGTGCGGGCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGATAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGt                                 <  2:336079/143‑1 (MQ=255)                                                                                                                gggTAAATGCAACCTGACCTATACCGGTGCGCGCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACAtccacgtcc                    <  1:846467/143‑1 (MQ=255)                                                                                                                 ggTAAAGGCAAACTGACCTCTCCCGGTGCGCTCGGCCAAGTGCTGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACAtccacgtcca                   <  1:796759/143‑1 (MQ=255)                                                                                                                      aGGCAAACTGACCTATACCGGTGCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTAc              <  2:357798/143‑1 (MQ=255)                                                                                                                        gCAAACTGACCTATCCCGGTGCGCTCGGCGAAGTGATGCAGGAGAACGTTCCCTCCCCGTGAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCg            <  1:208113/143‑1 (MQ=255)                                                                                                                          aaaCTGATCTATACCGGTGCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGA‑TTTTACGAAAAACGTGACATCCACGTCCACGTACCGGaa         <  1:291690/143‑1 (MQ=255)                                                                                                                            aCTGACCTATACCGGTGCGCTCGGCGAAGTGATGCAGGAGGCCATTCAGGCGGCGCTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGAGCAACCCTGATTTTTACGAAGAACGTGACATCCACGAGCGCGTACCGGagg        >  1:588897/1‑141 (MQ=255)                                                                                                                             ccgACCTATACCGGTGCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAgg       <  1:786818/141‑1 (MQ=255)                                                                                                                               gACCTATACCGGTGCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTg     <  1:507895/143‑1 (MQ=255)                                                                                                                                  cTATACCGGTGCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAGCGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGCCAAACGTGACATCCACGTCGACGTACCGCGAGGTGCGa  >  1:813187/1‑143 (MQ=255)                                                                                                                                            |                                                                                                                                     GCGTTTCGACTATGGTCGCGCGGATAACGAAAACCGTGTCGGTCAGGTAACCGGTCTGGCGTGGACGGAAGTGGGCGGTGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGA  >  NC_000913/460633‑460904

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