Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP006881	3,725,082	A→C	100%	A29A (GCT→GCG)	U712_19205 ←	Uncharacterized protein ywcI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006881	3,725,082	0	A	C	100.0%	164.5 / ‑6.6	45	A29A (GCT→GCG)	U712_19205	Uncharacterized protein ywcI
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (21/23);  total (22/23)
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CTTTTATAAAGCGCGGCGGCTTTCGCTTTTTGTTTTTCTTCCTCATTCACAGCCATCCGCCACTCATAAAACGAAACAAGTCCCGCGGCCTGCACCGCCGCTTTTCTTTGGTATGTACGCTGAGCCGAATAGACAGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTT  >  CP006881/3724948‑3725207                                                                                                                                       |                                                                                                                              cttTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACCGCCGCCGTTGTTTTGTATTTTAGCAGAGCCGAATAGACCGccg                                                                                                                           >  5:173606/1‑137 (MQ=255)      ataAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATACACCGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGc                                                                                                                       <  2:125917/138‑1 (MQ=255)       taAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCg                                                                                                                     <  7:218522/139‑1 (MQ=255)       taAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCg                                                                                                                     <  3:190632/139‑1 (MQ=255)       aaaaGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCAACCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCg                                                                                                                     <  1:138551/138‑1 (MQ=255)               ggcggcTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGttttt                                                                                                             <  8:247543/139‑1 (MQ=255)                    cTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTc                                                                                                        <  5:175822/139‑1 (MQ=255)                              ctgcttttCTTCCTCATTCACCGCCAGCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCTCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCAc                                                                                               <  8:304392/134‑1 (MQ=255)                                           tCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTGTTTGGTATTTACGCGGAGCCGCATAGACCGCGGCAGCCGCCGGTTTTTGGTCCGGCGTCACCCAGGTGATCaaa                                                                                  >  3:62099/1‑138 (MQ=255)                                                    gCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATcc                                                                        <  3:17932/139‑1 (MQ=255)                                                     ccATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCca                                                                       <  5:283664/139‑1 (MQ=255)                                                     ccATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCca                                                                       <  1:121353/139‑1 (MQ=255)                                                     ccATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGGCACCCAGTTCATCAAAAAGACCATCca                                                                       <  7:28091/139‑1 (MQ=255)                                                       aTCCGCCATTCATAACACGAATCGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGGTTTTCTGTCCGGCGGCACCCAGTTCATCAAAAAGACCATCcaca                                                                     <  8:351555/139‑1 (MQ=255)                                                       aTCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCcaca                                                                     <  4:224426/139‑1 (MQ=255)                                                               ttCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCGTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGGTCATCACAAAGGCCATCCACACCCgga                                                               >  1:231886/2‑135 (MQ=255)                                                               ttCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCGTTCTTTGGTATGTAAGCGGCGCCGAATAGACCGCGGCCGCCCCAGTTTTTCTGTCCGGCGTCCCCCACTTCATCAAAAAGACCCACCACACCCCTaa                                                              >  5:132908/2‑138 (MQ=255)                                                                   tAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATAAAAAATACCATCCACACCCGTAAAGATa                                                         <  3:174108/139‑1 (MQ=255)                                                                   tAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGATGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATa                                                         <  4:126397/139‑1 (MQ=255)                                                                   tAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCATTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATa                                                         <  7:162821/139‑1 (MQ=255)                                                                        cGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGc                                                    <  6:173606/139‑1 (MQ=255)                                                                           aaCGAGTCCTGCGGCCTGCACAGCCGCCGTTCTTTGGTATGTACCCTGAGCCGAATAGACACCGGCAGCCCCAGTTTTTATGTCAGTCGTCCCCCAGGTAATCAAAAATACCAACCACACCCGTAAAGATACAAGCaaa                                                 >  7:300107/1‑139 (MQ=255)                                                                                          cTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGaaaaa                                  <  6:137721/139‑1 (MQ=255)                                                                                          cTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGaaaaa                                  <  3:299626/139‑1 (MQ=255)                                                                                           tGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGCATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGGAAAGATAAAAGCCAACCTTTGATTGaaaaa                                  >  4:66078/1‑138 (MQ=255)                                                                                                  ccgcCTTTGTTTGGTATTTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTGCGGCGTCACCCAGTTCATGAAAAAGACCCTCCACCCCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTc                          >  2:108810/1‑139 (MQ=255)                                                                                                  ccgcCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTGCGGCGTCACCCCGTTCATCAAAAAGCCCATCCACACCCGGAAAGATACAATCAAAAGTTTCATAGAAAACCATCCCTc                          >  4:249724/1‑139 (MQ=255)                                                                                                  ccgcCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTc                          >  1:34390/1‑139 (MQ=255)                                                                                                                    tACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTggcagg         >  6:141687/1‑138 (MQ=255)                                                                                                                    tACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTggcagg         >  7:64812/1‑138 (MQ=255)                                                                                                                      cGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTggcaggca       <  3:350149/138‑1 (MQ=255)                                                                                                                      cGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTggcaggca       <  8:64812/138‑1 (MQ=255)                                                                                                                         tGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAGAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCGGGCAGCCt   >  5:326215/1‑139 (MQ=255)                                                                                                                         tGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGCCCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCt   >  8:256164/1‑139 (MQ=255)                                                                                                                         tGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCca                                                                       >  7:279690/1‑71 (MQ=255)                                                                                                                         tGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCca                                                                       <  8:279690/71‑1 (MQ=255)                                                                                                                         tGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCt   >  8:20780/1‑139 (MQ=255)                                                                                                                         tGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCt   >  1:315556/1‑139 (MQ=255)                                                                                                                         tGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCt   >  1:56439/1‑139 (MQ=255)                                                                                                                         tGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGGCTCCCAGTGCATCAAAAAGACCATCCCCCCCCGTAAAGATACATGCAAACGTATCATAGAGAAACATCCCTCCGCTTCTTTTTtgccgggcgcgg   >  4:45067/1‑133 (MQ=255)                                                                                                                          gAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGa                                                                              <  7:323465/63‑1 (MQ=255)                                                                                                                          gAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGa                                                                              >  8:323465/1‑63 (MQ=255)                                                                                                                          gAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAACCGTTTCATAGAAAAACTTCCCTCCGCGACTTTTTTGCCGGCAGCCtt  >  6:52254/1‑139 (MQ=255)                                                                                                                          gAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCGTCTTTTTGGCAGGCAGccgt  >  8:162106/1‑137 (MQ=255)                                                                                                                          gAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCAACCCCACCCCTAAAGATACAATCAAACGTTTCATAGAAAAACATCCCACCGCGGATTTTTTGCAGGCCGCCtt  >  7:132778/1‑139 (MQ=255)                                                                                                                                       |                                                                                                                              CTTTTATAAAGCGCGGCGGCTTTCGCTTTTTGTTTTTCTTCCTCATTCACAGCCATCCGCCACTCATAAAACGAAACAAGTCCCGCGGCCTGCACCGCCGCTTTTCTTTGGTATGTACGCTGAGCCGAATAGACAGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGCCTT  >  CP006881/3724948‑3725207

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: ‑