Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,940,308	A→G	L102P (CTG→CCG)	hsrA ←	putative multidrug or homocysteine efflux system

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,940,308	0	A	G	100.0%	115.5 / NA	32	L102P (CTG→CCG)	hsrA	putative multidrug or homocysteine efflux system
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (16/16);  total (16/16)

CACCAGACCCGGCATGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGTATAACTG  >  NC_000913/3940175‑3940454                                                                                                                                      |                                                                                                                                                   caccaGACCCGGCATGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCGGCTGTGGTAGCGAAtt                                                                                                                                    >  2:132747/1‑150 (MQ=255)                 gCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCGGCTGTGGTAGCGAATTAGAAAGTGCGCAGGcc                                                                                                                    >  1:241588/1‑150 (MQ=255)                   gACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCGGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCaga                                                                                                                 <  2:431972/151‑1 (MQ=255)                      aaaaTTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCGGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGaa                                                                                                              <  1:247096/151‑1 (MQ=255)                         aTTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCGGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAAccc                                                                                                           <  1:220346/151‑1 (MQ=255)                         aTTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCGGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAAccc                                                                                                           <  1:624810/151‑1 (MQ=255)                             aaTACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCGGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAAtg                                                                                                       >  2:309104/1‑151 (MQ=255)                                aCTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCGGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAAtgtg                                                                                                     >  2:552140/1‑150 (MQ=255)                                      agaagTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCGGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGAc                                                                                              <  2:408011/151‑1 (MQ=255)                                                tACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCGGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCaa                                                                                     <  2:491878/150‑1 (MQ=255)                                                 aCGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCGGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAgg                                                                                   <  1:613170/151‑1 (MQ=255)                                                   gAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCGGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAggg                                                                                  >  2:314416/1‑150 (MQ=255)                                                    aGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCGGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGt                                                                                 >  2:327686/1‑150 (MQ=255)                                                    aGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCGGCTGTGGTAGCGAAGTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTa                                                                                <  1:412416/151‑1 (MQ=255)                                                      gATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCGGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTaa                                                                               >  2:626346/1‑150 (MQ=255)                                                       aTAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCGGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTaaa                                                                              >  1:62001/1‑150 (MQ=255)                                                           gcgcgcAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCGGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGc                                                                         >  1:171298/1‑151 (MQ=255)                                                              cgcAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCGGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGAc                                                                      >  1:34636/1‑151 (MQ=255)                                                               gcAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCGGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGAcg                                                                     >  2:442292/1‑151 (MQ=255)                                                                        ccAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCGGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGaa                                                            <  1:275970/151‑1 (MQ=255)                                                                           gCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCGGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCg                                                         <  2:174149/151‑1 (MQ=255)                                                                               aGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCGGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCg                                                     <  2:39103/151‑1 (MQ=255)                                                                                  aaCAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCGGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCg                                                     >  1:357124/1‑148 (MQ=255)                                                                                         atcatcatTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCGGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATc                                           <  2:212944/151‑1 (MQ=255)                                                                                           catcatTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAGCCGGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATcca                                         <  1:5398/151‑2 (MQ=255)                                                                                                  gcgcCGCCTATCCCCTGAATAACCCGGAAGACAACCGGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCg                                  <  2:53011/151‑1 (MQ=255)                                                                                                       gccTATCCCCTGAATAACCCGGAAGACAACCGGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAAt                              >  2:316488/1‑150 (MQ=255)                                                                                                              cccTGAATAACCCGGAAGACAACCGGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCg                      >  2:5158/1‑151 (MQ=255)                                                                                                                      aaCCCGGAAGACAACCGGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTc              >  1:548715/1‑151 (MQ=255)                                                                                                                         ccGGAAGACAACCGGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGc           <  1:608963/151‑1 (MQ=255)                                                                                                                          cGGAAGACAACCGGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCg          <  1:132743/151‑1 (MQ=255)                                                                                                                                  aaCCGGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGTATAACtg  >  2:640314/1‑151 (MQ=255)                                                                                                                                      |                                                                                                                                                   CACCAGACCCGGCATGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGTATAACTG  >  NC_000913/3940175‑3940454

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