Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,733,178	T→C	P1364P (CCT→CCC)	lhr →	putative ATP‑dependent helicase Lhr

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,733,178	0	T	C	100.0%	81.7 / NA	24	P1364P (CCT→CCC)	lhr	putative ATP‑dependent helicase Lhr
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (9/15);  total (9/15)

GCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGTCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGACTCCGGGCGAATTATGCGAGGTCGTTTTGTA  >  NC_000913/1733109‑1733253                                                                      |                                                                            gctggcgGAAAATATGCTCGACCGCTACGGCATCATCAGTCGTCAGGCGGTGATAGCCGAAAATATCCCCGGCGGGTTTcc                                                                  <  2:298946/81‑1 (MQ=255) gctggcgGAAAATATGCTCGACCGCTACGGCATCATCAGTCGTCAGGCGGTGATAGCCGAAAATATCCCCGGCGGGTTTc                                                                   >  1:298978/1‑80 (MQ=255)             tatGCTCGACCGCTACGGCATCATCAGTCGTCAGGCGGTGATAGCCGAAAATATCCCCGGCGGGTTTCCATCGATGCAAAc                                                      <  1:298947/81‑1 (MQ=255)                 cTCGACCGCTACGGCATCATCAGTCGTCAGGCGGTGATAGCCGAAAATATCCCCGGCGGGTTTCCATCGATGCAAACGCtt                                                  >  1:298975/1‑81 (MQ=255)                 cTCGACCGCTACGGCATCATCAGTCGTCAGGCGGTGATAGCCGAAAATATCCCCGGCGGGTTTCCATCGATGCAAACGCt                                                   <  2:298949/80‑1 (MQ=255)                 cTCGACCGCTACGGCATCATCAGTCGTCAGGCGGTGATAGCCGAAAATATCCCCGGCGGGTTTCCATCGATGCAAACGCt                                                   <  2:298948/80‑1 (MQ=255)                   cGACCGCTACGGCATCATCAGTCGTCAGGCGGTGATAGCCGAAAATATCCCCGGCGGGTTTCCATCGATGCAAACGCTTTg                                                >  2:298976/1‑81 (MQ=255)                    gACCGCTACGGCATCATCAGTCGTCAGGCGGTGATAGCCGAAAATATCCCCGGCGGGTTTCCATCGATGCAAACGCttt                                                 <  2:298950/79‑1 (MQ=255)                            cGGCATCATCAGTCGTCAGGCGGTGATAGCCGAAAATATCCCCGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGt                                         <  2:298951/79‑1 (MQ=255)                              gCATCATCAGTCGTCAGGCGGTGATAGCCGAAAATATCCCCGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATgg                                     <  1:298952/81‑1 (MQ=255)                               catcatCAGTCGTCAGGCGGTGATAGCCGAAAATATCCCCGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGa                                    <  1:298953/81‑1 (MQ=255)                                 tcatcaGTCGTCAGGCGGTGATAGCCGAAAATATCCCCGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAg                                  <  1:298954/81‑1 (MQ=255)                                  catcaGTCGTCAGGCGGTGATAGCCGAAAATATCCCCGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGa                                 >  2:298989/1‑81 (MQ=255)                                   atcaGTCGTCAGGCGGTGATAGCCGAAAATATCCCCGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGAc                                <  1:298955/81‑1 (MQ=255)                                      aGTCGTCAGGCGGTGATAGCCGAAAATATCCCCGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGACTcc                             <  1:298956/81‑1 (MQ=255)                                          gtcAGGCGGTGATAGCCGAAAATATCCCCGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGACTCCGGGc                         >  2:298984/1‑81 (MQ=255)                                              ggCGGTGATAGCCGAAAATATCCCCGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGACTCCGGGCGAAt                     >  2:298982/1‑81 (MQ=255)                                               gCGGTGATAACCGAAGATATCCCCGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGACTCCGGGCGAAtt                    <  1:298957/81‑1 (MQ=255)                                                        gCCGAAAATATCCCCGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGACTCCGGGCGAATTATGCGAGGt           <  2:298958/81‑1 (MQ=255)                                                        gCCGAAAATATCCCCGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGACTCCGGGCGAATTATGCGAGGt           >  2:298979/1‑81 (MQ=255)                                                         ccGAAAATATCCCCGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGACTCCGGGCGAATTATGCGAGGTc          >  1:298985/1‑81 (MQ=255)                                                          cGAAAATATCCCCGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGACTCCGGGCGAATTATGCGAGGTCg         >  1:298980/1‑81 (MQ=255)                                                          cGAAAATATCCCCGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGACTCCGGGCGAATTATGCGAGGTCg         <  2:298959/81‑1 (MQ=255)                                                                  tCCCCGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGACTCCGGGCGAATTATGCGAGGTCGTTTTGTa  <  2:298960/80‑1 (MQ=255)                                                                      |                                                                            GCTGGCGGAAAATATGCTCGACCGCTACGGCATCATCAGTCGTCAGGCGGTGATAGCCGAAAATATCCCTGGCGGGTTTCCATCGATGCAAACGCTTTGTCGAAGTATGGAAGACTCCGGGCGAATTATGCGAGGTCGTTTTGTA  >  NC_000913/1733109‑1733253

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