Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	640,716	A→G	T322A (ACC→GCC)	ahpF →	alkyl hydroperoxide reductase, AhpF component

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	640,716	0	A	G	100.0%	83.2 / NA	24	T322A (ACC→GCC)	ahpF	alkyl hydroperoxide reductase, AhpF component
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (13/11);  total (13/11)

CGTTGATGAATACGACGTTGATGTGATCGACAGCCAGAGCGCCAGCAAACTGATCCCAGCAGCAGTTGAAGGTGGTCTGCATCAGATTGAAACAGCTTCTGGCGCGGTACTGAAAGCACGCAGCATTATCGTGGCGACCGGTGCAAAATGGCGCAACATGAACGTTCCGGGCGAAGATCAGTATCGCACCAAAGGCGTGACCTACTGCCCGCACTGCGACGGCCCGCTGTTTAAAGGTAAACGCGTAGCGGTTATCGGCGGCGGTAACTCCGG  >  NC_000913/640580‑640852                                                                                                                                         |                                                                                                                                         cGTTGATGAATGCGACGTTGATGTGATCGACAGCCAGAGCGCCAGCAAACTGATCCCAGCAGCAGTTGAAGGTGGTCTGCATCAGATTGAAACAGCTTCTGGCGCGGTACTGAAAGCACGCAGCATTATCGTGGCGGCCg                                                                                                                                       >  2:570443/1‑140 (MQ=255)                 gTTGATGTGATCGACAGCCAGAGCGCCAGCAAACTGATCCCAGCAGCAGTTGAAGGTGGTCTGCATCAGATTGAAACAGCTTCTGGCGCGGTACTGAAAGCACGCAGCATTATCGTGGCGGCCGGTGCAAAATGGCGCaa                                                                                                                       <  2:540342/140‑1 (MQ=255)                       gtgATCGACAGCCAGAGGGCCAGCAAACTGATCCCGGCAGCAGTTGAAGGTGGTCTGCATCAGATTGAAACAGCTTCTGGCGCGGTACTGAAAGCACGCAGCATTATCGTGGCGGCCGGTGCAAAATGGCGCAACATGaa                                                                                                                 <  2:285487/140‑1 (MQ=255)                          aTCGACAGCCAGAGCGCCAGCAAACTGATCCCAGCAGCAGTTGAAGGTGGTCTGCATCAGATTGAAACAGCTTCTGGCGCGGTACTGAAAGCACGCAGCATTATCGTGGCGGCCGGTGCAAAATGGCGCAACATGAACGt                                                                                                              >  1:41356/1‑140 (MQ=255)                               cAGCCAGAGCGCCAGCAAACTGATCCCAGCAGCAGTTGAAGGTGGTCTGCATCAGATTGAAACAGCTTCTGGCGCGGTACTGAAAGCACGCAGCATTATCGTGGCGGCCGGTGCAAAATGGCGCAACATGAACGTTCCgg                                                                                                         >  2:421276/1‑140 (MQ=255)                                aGCCAGAGCGCCAGCAAACTGATCCCAGCAGCAGTTGAAGGTGGTCTGCATCAGATTGAAACAGCTTCTGGCGCGGTACTGAAAGCACGCAGCATTATCGTGGCGGCCGGTGCAAAATGGCGCAACATGAACGTTCCggg                                                                                                        <  2:627485/140‑1 (MQ=255)                                      agCGCCAGCAAACTGATCCCAGCAGCAGTTGAAGGTGGTCTGCATCAGATTGAAACAGCTTCTGGCGCGGTACTGAAAGCACGCAGCATTATCGTGGCGGCCGGTGCAAAATGGCGCAACATGAACGTTCCGGGCg                                                                                                      >  2:595642/1‑136 (MQ=255)                                      agCGCCAGCAAACTGATCCCAGCAGCAGTTGAAGGTGGTCTGCATCAGATTGAAACAGCTTCTGGCGCGGTACTGAAAGCACGCAGCATTATCGTGGCGGCCGGTGCAAAATGGCGCAACATGAACGTTCCGGGCg                                                                                                      <  1:595642/136‑1 (MQ=255)                                             gCAAACTGATCCCAGCAGCAGTTGAAGGTGGTCTGCATCAGATTGAAACAGCTTCTGGCGCGGTACTGAAAGCACGCAGCATTATCGTGGCGGCCGGTGCAAAATGGCGCAACATGAACGTTCCGGGCGAAGATCAGTAt                                                                                           <  1:570443/140‑1 (MQ=255)                                                         cagcagcagTTGAAGGTGGTCTGCATCAGATTGAAACAGCTTCTGGCGCGGTACTGAAAGCACGCAGCATTATCGTGGCGGCCGGTGCAAAATGGCGCAACATGAACGTTCCGGGCGAAGATCAGTATCGCACCAAAGGc                                                                               <  2:95872/140‑1 (MQ=255)                                                                                     gATTGAAATAGCTTCTGACGCGCTACTGAAAGCTCCCAGTATTATCGTTGCGGCCGGTGCAAAATGGCGCAACATGAACGTTCCGGGCGAAGATCAGTATCGCACCAAAGGCGTGACCTACTGCCCGCACTGCGACGGcc                                                   <  2:297784/140‑1 (MQ=255)                                                                                     gATTGAAACAGCTTCTGGCGCGGTACTGAAAGCACGCAGCATTATCGTGGCGGCCGGTGCAAAATGGCGCAACATGAACGTTCCGGGCGAAGATCAGTATCGCACCAAAGGCGTGACCTACTGCCCGCACTGCGACGGcc                                                   >  2:52859/1‑140 (MQ=255)                                                                                     gATTGAAACAGCTTCTGGCGCGGTACTGAAAGCACGCAGCATTATCGTGGCGGCCGGTGCAAAATGGCGCAACATGAACGTTCCGGGCGAAGATCAGTATCGCACCAAAGGCGTGACCTACTGCCCGCACTGCGACGGcc                                                   <  2:58287/140‑1 (MQ=255)                                                                                      aTTGAAACAGCTTCTGGCGCGGTACTGAAAGCACGCAGCATTATCGTGGCGGCCGGTGCAAAATGGCGCAACATGAACGTTCCGGGCGAAGATCAGTATCGCTCCAAAGGCGTGACCCACTGCCCGCACTGCGACGGccc                                                  >  1:297784/1‑140 (MQ=255)                                                                                       ttGAAACAGCTTCTGGCGCGGTACTGAAAGCACGCCGCATTATCGTGGCGGCCGGTGCAAAATGGCGCAACATGAACGTTCCGGGCGAAGATCAGTATCGCACCAAAGGCGTGACCTACTGCCCGCACTGCGACGGCCCg                                                 <  2:228665/140‑1 (MQ=255)                                                                                       ttGAAACAGCTTCTGGCGCGGTACTGAAAGCACGCAGCATTATCGTGGCGGCCGGTGCAAAATGGCGCAACATGAACGTTCCGGGCGAAGATCAGTATCGCACCAAAGGCGAGACCTACTGCCCGCACTGCGACGGCCCg                                                 <  2:41356/140‑1 (MQ=255)                                                                                          aaaCAGCTTCTGGCGCGGTACTGAAAGCACGCAGCATTATCGTGGCGGCCGGTGCAAAATGGCGCAACATGAACGTTCCGGGCGAAGATCAGTATCGCACCAAAGGCGTGACCTACTGCCCGCACTGCGACg                                                      <  2:557797/132‑1 (MQ=255)                                                                                          aaaCAGCTTCTGGCGCGGTACTGAAAGCACGCAGCATTATCGTGGCGGCCGGTGCAAAATGGCGCAACATGAACGTTCCGGGCGAAGATCAGTATCGCACCAAAGGCGTGACCTACTGCCCGCACTGCGACg                                                      >  1:557797/1‑132 (MQ=255)                                                                                                cTTCTGGCGCGGTACTGAAAGCACGCAGCATTATCGTGGCGGCCGGTGCAAAATGGCGCAACATGAACGTTCCGGGCGAAGATCAGTATCGCACCAAAGGCGTGACCTACTGCCCGCACTGCGACGGCCCGCTGTTTaaa                                        >  1:134386/1‑140 (MQ=255)                                                                                                                      aCGCAGCATTATCGTGGCGGCCGGTGCAAAATGGCGCAACATGAACGTTCCGGGCGAAGATCAGTATCGCACCAAAGGCGTGACCTACTGCCCGCACTGCGACGGCCCGCTGTTTAAAGGTAAACGCGTAGCGGTTATcg                  >  1:386826/1‑140 (MQ=255)                                                                                                                      aCGCAGCATTATCGTGGCGGCCGGTGCAAAATGGCGCAACATGAACGTTCCGGGCGAAGATCAGTATCGCACCAAAGGCGTGACCTACTGCCCGCACTGCGACGGCCCGCTGTTTAAAGGTAAACGCGTAGCGGTTATcg                  >  2:650563/1‑140 (MQ=255)                                                                                                                       cGCAGCATTATCGTGGCGGCCGGTGCAAAATGGCGCAACATGAACGTTCCGGGCGAAGATCAGTATCGCACCAAAGGCGTGACCTACTGCCCGCACTGCGACGGCCCGCTGTTTAAAGGTAAACGCGTAGCGGTTATcgg                 >  1:360557/1‑140 (MQ=255)                                                                                                                                    tGGCGGCCGGTGCAAAATGGCGCAACATGAACGTTCCGGGCGAAGATCAGTATCGCACCAAAGGCGTGACCTACTGCCCGCACTGCGACGGCCCGCTGTTGAAAAGTAGACGCGTAGAGGTTATCGGCGGCGGGAACTcc    >  2:217371/1‑140 (MQ=255)                                                                                                                                      gCGGCCGGTGCAAAATGGCGCAACATGAACGTTCCGGGCGAAGATCAGTATCGCACCAAAGGCGTGACCTACTGCCCGCACTGCGACGGCCCGCTGTTTAAAGGTAAACGCGTAGCGGTTATCGGCGGCGGGAACTCCgg  >  1:643967/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                         CGTTGATGAATACGACGTTGATGTGATCGACAGCCAGAGCGCCAGCAAACTGATCCCAGCAGCAGTTGAAGGTGGTCTGCATCAGATTGAAACAGCTTCTGGCGCGGTACTGAAAGCACGCAGCATTATCGTGGCGACCGGTGCAAAATGGCGCAACATGAACGTTCCGGGCGAAGATCAGTATCGCACCAAAGGCGTGACCTACTGCCCGCACTGCGACGGCCCGCTGTTTAAAGGTAAACGCGTAGCGGTTATCGGCGGCGGTAACTCCGG  >  NC_000913/640580‑640852

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