Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP006881	3,725,049	T→C	K40K (AAA→AAG)	U712_19205 ←	Uncharacterized protein ywcI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006881	3,725,049	0	T	C	100.0%	102.4 / NA	29	K40K (AAA→AAG)	U712_19205	Uncharacterized protein ywcI
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (4/25);  total (4/25)

GGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTTTGTTTTTCTTCCTCATTCACAGCCATCCGCCACTCATAAAACGAAACAAGTCCCGCGGCCTGCACCGCCGCTTTTCTTTGGTATGTACGCTGAGCCGAATAGACAGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACAT  >  CP006881/3724918‑3725178                                                                                                                                    |                                                                                                                                  ggCGCTGCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCGTTCGCTTTCTGCTTTTCTTCCTCATTCACCGCCATCCGCCTTTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCt                                                                                                                              <  4:204767/137‑1 (MQ=255)  gCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTg                                                                                                                           <  5:120628/139‑1 (MQ=255)      ttCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTAt                                                                                                                       <  3:290352/139‑1 (MQ=255)            gAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACg                                                                                                                  <  8:217636/138‑1 (MQ=255)            gAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACg                                                                                                                  <  7:145799/138‑1 (MQ=255)            caaGGCTTGCAAACATCTTCTTTTATAACGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCAATAAAAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACg                                                                                                                  <  8:131106/137‑1 (MQ=255)             aaGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCGTCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCt                                                                                                                <  8:6704/139‑1 (MQ=255)               ggCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGa                                                                                                              <  8:156703/139‑1 (MQ=255)                gCTTGCAAACTTCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAg                                                                                                             <  8:88631/139‑1 (MQ=255)                gCTTGCAAACCTCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTGCTTCCTCATTCACAGCCCTCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAg                                                                                                             <  2:107973/139‑1 (MQ=255)                gCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAg                                                                                                             <  1:87363/139‑1 (MQ=255)                gCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAg                                                                                                             <  4:262049/139‑1 (MQ=255)                 cttcaaaaCATATTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGc                                                                                                            <  1:181358/134‑1 (MQ=255)                                                  cTTTCGCTTTCTGCTTTTCTTCCTCATTCACCGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTc                                                                           <  1:174103/139‑1 (MQ=255)                                                  cTTTCGCCTTCAGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTc                                                                           <  1:292766/139‑1 (MQ=255)                                                                 tttCTACCACATTCACCGCCACCCCCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGt                                                             <  3:20143/138‑1 (MQ=255)                                                                   tcttcCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTtc                                                           <  5:3623/138‑1 (MQ=255)                                                                                   ccATCCGCCATTCATAAAACGAATCGAGTCCTGCGGCCTGCACAGCCGCCTTTCATTGGTATGTACGCGGAGCCGAATAGACCGCGGCAGCCGCAGGTTTTCTGTCCGGCGTCACGCAGTTCATCAAAAAGACCATCca                                          <  3:141654/139‑1 (MQ=255)                                                                                   ccATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCca                                          <  4:71878/139‑1 (MQ=255)                                                                                   ccATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCca                                          <  3:45594/139‑1 (MQ=255)                                                                                   ccATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCca                                          <  1:246315/139‑1 (MQ=255)                                                                                     aTCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCa                                                                   >  7:22698/1‑112 (MQ=255)                                                                                     aTCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCa                                                                   <  8:22698/112‑1 (MQ=255)                                                                                            tttCATAAACCTAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTTGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTaa                                 <  1:362444/137‑1 (MQ=255)                                                                                                                     ggCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGaa        <  4:43982/139‑1 (MQ=255)                                                                                                                         tGCACAGCCGCCGTTCTTTGGTATGTACGCTGAGCCGAATAGCCCGCGGCACCCGCCGTTTTTGGGTCCGGCGTCGCCCAGTTGCTCAAAAAGACCCTCCCCACCCGTAAAGATACAAGCCAACGTTTGATAGTAAAAc    >  1:37583/1‑139 (MQ=255)                                                                                                                           cacaGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACAt  >  6:361274/1‑139 (MQ=255)                                                                                                                                ccgcCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCca                                          <  3:162424/94‑1 (MQ=255)                                                                                                                                ccgcCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTATGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCca                                          >  4:162424/1‑94 (MQ=255)                                                                                                                                    |                                                                                                                                  GGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTTTGTTTTTCTTCCTCATTCACAGCCATCCGCCACTCATAAAACGAAACAAGTCCCGCGGCCTGCACCGCCGCTTTTCTTTGGTATGTACGCTGAGCCGAATAGACAGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACAT  >  CP006881/3724918‑3725178

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