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	100.0%	114.9 / NA	41	S629A (TCG→GCG)	lon	Lon protease
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (17/24);  total (17/24)

GCGTTTCGACTATGGTCGCGCGGATAACGAAAACCGTGTCGGTCAGGTAACCGGTCTGGCGTGGACGGAAGTGGGCGGTGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGA  >  NC_000913/460633‑460910                                                                                                                                            |                                                                                                                                           gCGTTTCGACTATGGTCGCGCGGATAACGAAAACCGTGTCGGTCAGGTAACCGGTCTGGCGTGGACGGAAGTGGGCGGTGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGGAGAGGCAAACTGACCTGTACCGGTgcgc                                                                                                                                         >  1:596296/1‑143 (MQ=255)         actagcgTCTCGCTGATAACGAACACCGGGGCGGTCAGGTAGCCGGTCTGGCGTGGTGGGATGTGGGCGCTGACTCGCGAACCATTGAAACCGCTTGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTGCGCTCGGCGaa                                                                                                                                 <  2:603746/137‑1 (MQ=255)               gTCGCGCGGATAACGAAAACCGTGTCGGTCAGGTAACCGGTCTGGCGTGGACGGAAGTGGGCGGTGACTTGCTGACCATTGAAACCGCATGGGTTCCGGGTAGAGGCAAAGTCACCTATACCGGGGCGCGCGGCGACgtggtg                                                                                                                           >  2:330213/1‑143 (MQ=255)               gTCGCGCGGATAACGAAAACCGTGTCGGTCAGGTAACCGGTCTGGCGTGGACGGAAGTGGGCGGTGACTCGCTGACCATTGAAACCGCATGTGTTCCGGGTACAGGCAAACTGACCTATACCGGTGCGCTCGGCGAAGTGATg                                                                                                                           >  2:615689/1‑143 (MQ=255)                tCGCGCGGATAACGAAAACCGTGTCGGTCAGGTACCCGGTCTGGCGTGGACGGAAGTGGGAGGTGACTTGCTGACCATTGAAACCGCGAGTGTTCTGGGTAAAGGCAGACTGACCTATACCGGTGCGCTCGGCGAAGTGATGc                                                                                                                          >  1:129497/1‑143 (MQ=255)                    gcgGATACCGGAAACCTTGCCGGTCAGGTAACCGGTCTGGCGTGGGCGGAAGTGGGCGGTGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAAATGACCTATACCGGTGCGCTCGGCGAAGTGATGCAGGa                                                                                                                      <  2:265345/143‑1 (MQ=255)                      ggATGACGAAAGCCGTGTCGGTCAGGTAACCGGTCTGGCGTGGACGGAAGTGGGCGGTGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTGCGCTCGGCGAAGTGATGCAGGAGt                                                                                                                    <  1:164146/143‑1 (MQ=255)                                 gcccggtCGGTCAGGTAACCGGTCTGGCGTGGACGGAAGTGGGCGGTGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTGCGCTCGGCGAAGTGATGCAGGAGTCCATGCAggcg                                                                                                         <  2:323318/138‑1 (MQ=255)                                            cAGGTAACCGGGCCGGCGGGGACGGAAGTGGGCGGTGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTGCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACggt                                                                                              <  1:597341/143‑1 (MQ=255)                                                 acccGGTCGGGCGAGCACGGAACGGGGCGGTGACTTGCTGACCCTTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGCTGCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTc                                                                                         <  1:91868/141‑1 (MQ=255)                                                     ggTCTGGCGTTGACTGGTGCGGGCGGCGACGTGCCGACCATGGAAACCCCATGGGTTCCGGGTAAAAGCCAGCGGCCCTATACGGGTGCGCTGGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTcgtgc                                                                                     <  2:401422/143‑1 (MQ=255)                                                              tGGACGGAAGCGGGCGGTTGCTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATCCCGGTGCGCTCGGCGAAGTGATGCCGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGa                                                                            <  1:100278/143‑1 (MQ=255)                                                                  cGGGAGTGGGTGGTGCCTTGCTTACCATGGGAGCCGGATGGGGTGCGTGTAAAGGCAAACTGCCCTATACCGGTGCGCTCGGCGAAGTGACGCAGGAGGCCATTCAGGCGGCGTTAACGGTGGGTCGTGCGCGGGCGGaaaaa                                                                        <  2:163049/143‑1 (MQ=255)                                                                   ggAAGTGGGCGGTGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTCTACCGGTGCGCTCGGCCAAGGGACGCCGGAGGCCGTTCAGGCGGCGTGAACGGTGCTTCGTGCGCGTGCggccatcc                                                                       >  1:150578/1‑137 (MQ=255)                                                                      ggTGGGCGGTGACTTGCGGACCAGTGGGAGCGCATGCGTTCCGGGTAAAGGCAAACTGACCTATACCGGTGCGCTTGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTgg                                                                    <  2:235386/142‑1 (MQ=255)                                                                      aGTGGGCGGTGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTGCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTgg                                                                    <  1:765472/143‑1 (MQ=255)                                                                      aGTGGGCGGTCACGTGCTGACTACTGAGATGGCAGGGGTGCCGGGTAAATGCAAAGTGACCTATACGGGTGCGCGCGGCGAGGTGATGCAGGAGTCCAGTCAGGTGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTgg                                                                    <  2:786519/143‑1 (MQ=255)                                                                       gTGGGCGGTGACTTTCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTGCGCTCGGCTAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTggg                                                                   <  1:511748/143‑1 (MQ=255)                                                                        tgctcGGTGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTGCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTgggg                                                                  <  1:137909/139‑1 (MQ=255)                                                                         gggCGGTGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTGCGCTCGGCGAAGTGATGCAGGAATCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGACAAACTagggc                                                                 >  2:38949/1‑142 (MQ=255)                                                                             ggTGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTGCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGAACaa                                                             >  2:71738/1‑143 (MQ=255)                                                                                      cTGACCATTGAAACCGCATGTCTTCCGGGTAAAGGCAAACTGACCTATACCGGTGCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCCGCGTTAACGGTGGTTCGTGCGCGTGCGGAAACGCTGGGGATCAACCCTgcgtt                                                    >  2:285876/1‑139 (MQ=255)                                                                                       tGACCAGTGAAACCGCATGGGTTCCGGGTAAAGGCAAACTGACCTATACCGGTGCGCTCGGCGAAGTGATGCAGGAGGCCATTCAGGCGGCGTTAACGGTGGTGCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGAtttt                                                   <  1:150713/143‑1 (MQ=255)                                                                                         aCCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTGCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTa                                                 <  1:381537/143‑1 (MQ=255)                                                                                         aCCATTGAAACCGCATGTGTTCCGGGTAAAGGAAAACTGACCTATACCGGTGCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTa                                                 <  1:308509/143‑1 (MQ=255)                                                                                         aCCATCGGAACCGCATGTGGTCCGGGTAAAGGCAAACTGACCTATACCGGCGCGCTCGGCGAAGGGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTa                                                 <  1:631672/143‑1 (MQ=255)                                                                                           cATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTGCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAGAATGGGGATCAACCCTGAtgtcgccg                                               >  2:656397/1‑136 (MQ=255)                                                                                            aTTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTGCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGTTAAACTGGGGATCAACCCTGATTTTTACGa                                              >  1:146172/1‑143 (MQ=255)                                                                                            aTTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTGCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGa                                              >  1:72278/1‑143 (MQ=255)                                                                                                aacccGCATGTGTTCCGGGTAAAGGCAAACTGACCTATCCCGGTGCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGaaaaa                                          <  1:31911/140‑1 (MQ=255)                                                                                                 caCCCCATCTGCTCCGGGTAAAGGCACACTGACCTATACCGGTGCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCTTCCGCGTGCGGAAAAACTGGGTATCAACCCTGATTTTTACGAAAAAc                                         <  1:332948/142‑1 (MQ=255)                                                                                                          tgtTTCGGGTAAAGTTTATCTGTCATATACCGGTGCCAACGGCGAAGGGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACAtcc                                <  1:328628/143‑1 (MQ=255)                                                                                                            ttCCGGGTAAAGGCAAACTGACCTATACCGGTGCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACAtccac                              >  2:362061/1‑143 (MQ=255)                                                                                                               cGGGTAAAGGCAAACTGACCTATACCGGTGCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACAtccacgtc                           >  2:424718/1‑143 (MQ=255)                                                                                                               cGGGTAAAGGCAAACTGACCTATACCGGTGCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGGCAACAGTGGTTCGTGCGAGCCCCGAAAGACGGGGGATCAACCCTGATTTTTACGAAAAGCGGGACAtccacggc                           >  2:432513/1‑141 (MQ=255)                                                                                                                   tAAAGGCAAACTGACCTATACCGGTGCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGTGTTTAGGGTGGTCCGAGCAGGTGAGAACATAACGCAGATCAACCCTGATTTCTACGCACAACCCGACAgccacgtccacg                       >  2:565825/1‑143 (MQ=255)                                                                                                                      aGGCGAACTGTCCTATCCCGGTGCGCCTGGCGAAGTGATGCAGGAAATAATGCTTGCGACATTCAAAGTTTTGCTGGCTAGTGCGGAAAACCTGGGGATCTCCACTGATTTTTACGAAAAACGTGACATCCACGTCCACGTAc                    <  2:565766/143‑1 (MQ=255)                                                                                                                        gCAAACTGACCTATACCGGTGCGCTCGGCGAAGTGATGCAGGAGCCCATTCAGGCGGCGTCAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGAGCAACCCTGCTTTTTCCGCAACACGTGACATCCACGCCCACGTACCg                  >  2:113633/1‑143 (MQ=255)                                                                                                                             cTGACCTATACCGGTGCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAgg             <  1:71738/143‑1 (MQ=255)                                                                                                                              tGACCTAGTCCGGTGCGCTCCGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGt            <  2:445608/143‑1 (MQ=255)                                                                                                                                        cGGTGCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGGTTTTTACGAAAAACGTGACATCCACGTCCGCGTACCGGAAGGTGCGACGCCgg  >  1:504859/1‑142 (MQ=255)                                                                                                                                        cGGTGCGCTCGGCGAAGTGATGCAGGAGTCCATCCAGGCGGCGTTAACGGTGGTCCGGGCGCGTGCGGAACAACTGGGGATCAACCCTGATTGCTCCGAACAACGTGACCTCCCCGTCCACGTGCCGGCAGGTGCGACGCCGa  >  2:399215/1‑143 (MQ=255)                                                                                                                                            |                                                                                                                                           GCGTTTCGACTATGGTCGCGCGGATAACGAAAACCGTGTCGGTCAGGTAACCGGTCTGGCGTGGACGGAAGTGGGCGGTGACTTGCTGACCATTGAAACCGCATGTGTTCCGGGTAAAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCACGTCCACGTACCGGAAGGTGCGACGCCGA  >  NC_000913/460633‑460910

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