Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,940,095	C→T	M33I (ATG→ATA)	lpxM ←	Lipid A biosynthesis myristoyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,940,095	0	C	T	100.0%	83.6 / NA	26	M33I (ATG→ATA)	lpxM	Lipid A biosynthesis myristoyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (12/14);  total (12/14)

TTCTCCGGCCCGCGTATTGCCAACTCAGCCATCATTGCCATCGCTTGCGGCGCGGTGGCAAACATCTCATCAACAATCGCTTCGCGTTCAGCTTCACTACGTTCTGGAAAGCAGAGCGACAGATTGATTAACGCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCCATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGG  >  NC_000913/1939865‑1940333                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                               ttCTCCGGCCCGCGTATTGCCAACTCAGCCATCATTGCCATCGCTTGCGGCGCGGTGGCAAACATCTCATCAACAATCGCTTCGCGTTCAGCTTCACTACGTTCTGGAAAGCAGAGCGACAGATTGATTAACGCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGcc                                                                                                                                                                                                                             <  1:146946/250‑1 (MQ=255)      cGGCCCGCGTATTGCCAACTCAGCCATCATTGCCATCGCTTGCGGCGCGGTGGCAAACATCTCATCAACAATCGCTTCGCGTTCAGCTTCACTACGTTCTGGAAAGCAGAGCGACAGATTGATTAACGCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACTCCCAGCCATGCTc                                                                                                                                                                                                                       >  1:281861/1‑251 (MQ=255)          ccGCGTATTGCCAACTCAGCCATCATTGCCATCGCTTGCGGCGCGGTGGCAAACATCTCATCAACAATCGCTTCGCGTTCAGCTTCACTACGTTCTGGAAAGCAGAGCGACAGATTGATTAACGCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTcccc                                                                                                                                                                                                                    >  1:266885/1‑250 (MQ=255)              gTATTGCCAACTCAGCCATCATTGCCATCGCTTGCGGCGCGGTGGCAAACATCTCATCAACAATCGCTTCGCGTTCAGCTTCACTACGTTCTGGAAAGCAGAGCGACAGATTGATTAACGCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGt                                                                                                                                                                                                                 >  1:253322/1‑249 (MQ=255)                aTTGCCAACTCAGCCATCATTGCCATCGCTTGCGGCGCGGTGGCAAACATCTCATCAACAATCGCTTCGCGTTCAGCTTCACTACGTTCTGGAAAGCAGAGCGACAGATTGATTAACGCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAgcg                                                                                                                                                                                                             <  2:253322/251‑1 (MQ=255)                     cAACTCAGCCATCATTGCCATCGCTTGCGGCGCGGTGGCAAACATCTCATCAACAATCGCTTCGCGTTCAGCTTCACTACGTTCTGGAAAGCAGAGCGACAGATTGATTAACGCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGt                                                                                                                                                                                                        <  1:186759/251‑1 (MQ=255)                             ccATCATTGCCATCGCTTGCGGCGCGGTGGCAAACATCTCATCAACAATCGCTTCGCGTTCAGCTTCACTACGTTCTGGAAAGCAGAGCGACAGATTGATTAACGCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAg                                                                                                                                                                                                <  1:42917/251‑1 (MQ=255)                                  atTGCCATCGCTTGCGGCGCGGTGGCAAACATCTCATCAACAATCGCTTCGCGTTCAGCTGCACTACGTTCTGGAAAGCAGAGCGACAGATTGATTAACGCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGAttt                                                                                                                                                                                           <  2:214712/251‑1 (MQ=255)                                        agcGCTTGCTGCGCGTTGTCAAACATCTCACCAACAATCGCTTCGCGTTCAGCTTCACTACGTTCTGGAAAGCAGAGCGACAGATTGATTAACGCACGACGGCCTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCaaa                                                                                                                                                                                     <  2:77949/249‑1 (MQ=255)                                        aTCGCTTGCGGCGCGGTGGCAAACATCTCATCAACAATCGCTTCGCGTTCAGCTTCACTACGTTATGGAAAGCAGAGCGACAGATTGATTAACGCACGACGGCGTGAGCTTTTTCCCAGTCGCCCGTCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCaac                                                                                                                                                                                     >  2:160233/1‑250 (MQ=255)                                          cGCTTGCGGCGCGGTGGCAAACATCTCATCAACAATCGCTTCGCGTTCAGCTTCACTACGTTCTGGAAAGCAGAGCGACAGATTGATTAACGCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCAAAc                                                                                                                                                                                    <  1:27825/250‑1 (MQ=255)                                            cTTGCGGCGCGGTGGCAAACATCTCATCAACAATCGCTTCGCGTTCAGCTTCACTACGTTCTGGAAAGCAGAGCGACAGATTGATTAACGCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCAAACTc                                                                                                                                                                                  >  2:39258/1‑250 (MQ=255)                                              tGCGGCGCGGTGGCAAACAGCTCATCAACAATCGCTTCGCGTTTAGCTTCACGACGTTCTGGAAATCAGATCGTCTGATTGATTAACGCACGACGGCGGGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCTGTGCCAGAATGGGATCACGGAACTTTGTCGGCGTTAAAGCGGTGCCCGCTGTCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCAAACTCAg                                                                                                                                                                                <  1:188545/250‑1 (MQ=255)                                                 gacgcgGCGGCAAACATCTCATCAACAATCGCTTCGCGTTCAGCTTCACTACGTTCTGGAAAGCAGAGCGACAGATTGATTAACGCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCAAACTCCGGAAt                                                                                                                                                                            >  2:47637/3‑251 (MQ=255)                                                                            aTCGCTTCGCGTTCAGCTTCACTACGTTCTGGAAAGCAGAGCGACAGATTGATTAACGCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCGTTTc                                                                                                                                                 >  1:39899/1‑251 (MQ=255)                                                                                              tCACTACGTTCTGGAAAGCAGAGCGACAGATTGATTAACGCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGttt                                                                                                                                 <  1:67670/249‑1 (MQ=255)                                                                                                                       aCAGATTGATTAACGCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTg                                                                                                       >  2:140860/1‑250 (MQ=255)                                                                                                                                          acgGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACtt                                                                                   <  2:197785/251‑1 (MQ=255)                                                                                                                                                          cccAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTGCGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATa                                                                    <  1:72854/250‑1 (MQ=255)                                                                                                                                                                 gtcCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCAAACTCAGGAATGTATTCGCTATTATTTTTTTTc                                                                                                                                                      <  2:230496/161‑1 (MQ=255)                                                                                                                                                                 gtcCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCAAACTCAGGAATGTATTCGCTATTATTTTTTTTc                                                                                                                                                      >  1:230496/1‑161 (MQ=255)                                                                                                                                                                                ccAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGt                                             >  1:168007/1‑251 (MQ=255)                                                                                                                                                                                          cAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCg                                   >  2:270183/1‑251 (MQ=255)                                                                                                                                                                                           aGAATGGGATCAAGGAACTTTGGCGGCGTTAAAGCGATACCCGCTATCGATGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCgc                                  >  2:291594/1‑250 (MQ=255)                                                                                                                                                                                                              ttGGCGGCGTTAAAGCGTTACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGTTGGCGAAAGGATTTATCAAACTCAGGAATGTATTCGCTATTATTTGTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCa                <  1:73787/250‑1 (MQ=255)                                                                                                                                                                                                                           aGCGATACCCGCTATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCgg  <  1:190975/251‑1 (MQ=255)                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                               TTCTCCGGCCCGCGTATTGCCAACTCAGCCATCATTGCCATCGCTTGCGGCGCGGTGGCAAACATCTCATCAACAATCGCTTCGCGTTCAGCTTCACTACGTTCTGGAAAGCAGAGCGACAGATTGATTAACGCACGACGGCGTGAGCTTTTTCCCAGTCGTCCGGCAAAACGTCCCAGCCGTGCCAGAATGGGATCACGGAACTTTGGCGGCGTTAAAGCGATACCCGCCATCGCTGCTACGCCCAGCCATGCTCCCCAGTAGCGCGGGTGGCGAAAGGATTTATCAAACTCAGGAATGTATTCGCTATTATTTTTTTTCGTTTCCATGCTTTTCCAGTTTCGGATAAGGCAAAAATCAATCTGGTGATAGTGTAGCGGCGCAACTTGCCCCGCACCAAATAAAAAAAGCCGGTACTGACTGCGTACCGGCTGCGAATGGATGTTAATTAATCAAACCGTAGCTGCGG  >  NC_000913/1939865‑1940333

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑