Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP006881	3,724,808	G→A	100%	L7L (CTA→TTA)	U712_19200 ←	SacPA operon antiterminator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006881	3,724,808	0	G	A	100.0%	144.4 / NA	43	L7L (CTA→TTA)	U712_19200	SacPA operon antiterminator
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (16/27);  total (16/27)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AGATCTTTTCCACTTTATTCATAGGGATGAGGTCATTTTTCTTTTTTTGAAAAGCGATTCCTGGTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAGTACTTTATAGATTTTCAAGAGGCTATCTCTCCCGCCAATTAGCTTTTGACTGGGATTTTACCTCCTTTTGACCGCCATGTAAAGATGCGGGCCAGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGC  >  CP006881/3724688‑3724934                                                                                                                         |                                                                                                                               aGATCTTTTCCACTTTATTCATAGGGATGAGATCATTTTTCTTTTTTTGAAAAGCGATTCCCGGTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATACTTTATAGATTTTCaa                                                                                                              <  8:113563/139‑1 (MQ=255) aGATCTTTTCAACTTTATTCATAGGGATGAGATCATTTTTCTTTTTTTGAAAAGCGATTCCCGGTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATACTTTATAGATTTTCaa                                                                                                              <  4:149702/139‑1 (MQ=255)        ttCCACTTTATTCATAGGGATGAGATCATTTTTCTTTTTTTGAAAAGCGATTCCCGGTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATACTTTATAGATTTTCAAGAGGTTa                                                                                                       >  8:27985/1‑139 (MQ=255)        ttCCACTTTATTCATAGGGATGAGATCATTTTTCTTTTTTTGAAAAGCGATTCCCGGTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATACTTTATAGATTTTCAAGAGGTTa                                                                                                       >  1:277642/1‑139 (MQ=255)        ttCCACTTTATTCATAGGGATGAGATCATTTTTCTTTTTTTGAAAAGCGATTCCCGGTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATACTTTATAGATTTTCAAGAGGTTa                                                                                                       >  1:53274/1‑139 (MQ=255)        ttCCACTTTATTCATAGGGATGAGATCATTTTTCTTTTTTTGAAAAGCGATTCCCGGTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGACGCATTGTTGTTTAATACTTTATAGATTTTAATGAGGTTa                                                                                                       >  3:232778/1‑139 (MQ=255)                  ttCATAGGGATGAGATCATTTTTCTTTTTTTGAAAAGCGATTCCCGGTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATACTTTATAGATTTTCAAGAGGTTATCTCTCCCGc                                                                                             >  4:93142/1‑139 (MQ=255)                      tAGGGATGAGATCATTTTTCTTTTTTTGAAAAGCGATTCCCGGTACCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATACTTTATAGATTTTCAAGAGGTTATCTCTCCCGCCAAt                                                                                         <  8:119468/139‑1 (MQ=255)                       aGGGATGAGATCATTTTTCTTTTTTTGAAAAGCGATTCCCGGTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATACTTTATAGATTTTCAAGAGGTTATCTCTCCCGCCAAtt                                                                                        >  7:19568/1‑139 (MQ=255)                           aTGAGATCATTTTTCTTTTTTTGAAAAGCGATTCCCGGTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATACTTTATAGATTTTCAAGAGGTTATCTCTCCCGCCAATTTGCt                                                                                    <  1:253650/139‑1 (MQ=255)                           aTGAGATCATTTTTCTTTTTTTGAAAAGCGATTCCCGGTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATACTTTATAGATTTTCAAGAGGTTATCTCTCCCGCCAATTTGCt                                                                                    <  1:249933/139‑1 (MQ=255)                            tGAGATCATTTTTCTTTTTTTGAAAAGCGATTCCCGGTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATACTTTATAGATTTTCAAGAGGTTAtct                                                                                                    >  3:168434/1‑122 (MQ=255)                            tGAGATCATTTTTCTTTTTTTGAAAAGCGATTCCCGGTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATACTTTATAGATTTTCAAGAGGTTAtct                                                                                                    <  4:168434/122‑1 (MQ=255)                            tGAGATCATTTTTCTTTTTTTGAAAAGCGATTCCCGGTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATACTTTATAGATTTTCAAGAGGTTATCTCTCCCGCCAATTTGCtt                                                                                   >  2:244309/1‑139 (MQ=255)                            tGAGATCATTTTTCTTTTTTTGAAAAGCGATTCCCGGTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATACTTTATAGATTTTCAAGAGGTTATCTCTCCCGCCAATTTGCtt                                                                                   >  6:199999/1‑139 (MQ=255)                            tGAGATCATTTTTCTTTTTTTGAAAAGCGATTCCCGGTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATACTTTATAGATTTTCAAGAGGTTATCTCTCCCGCCAATTTGCtt                                                                                   >  5:151360/1‑139 (MQ=255)                                     ttttCTTTTTTTGAAAAGCGATTCCCGGTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATACTTTATAGATTTTCAAGAGGTTATCTCTCCCGCCAATTTGCTTTTGACt                                                                             <  3:12038/136‑1 (MQ=255)                                      tttCTTTTTTTGAAAAGCGATTCCCGGTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATACTTTATAGATTTTCAAGAGGTTATCTCTCCCGCCAATTTGCTTTTGACTggg                                                                          <  5:145936/138‑1 (MQ=255)                                             ttttGAAAAGCGATTCCCGGTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATACTTTATAGATTTTCAAGAGGTTATCTCTCCCGCCAATTTGCTTTTGACTGGGATTTTAc                                                                   <  4:76974/138‑1 (MQ=255)                                                       cgcttCCCGGTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATACTTTATAGATTTTCAAGAGGTTATCTCTCCCGCCAATTTGCTTTTGACTGGGATTTTACCTCCTTTTGAc                                                        <  4:209084/136‑1 (MQ=255)                                                         aTTCCCGGTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATACTTTATAGATTTTCAAGAGGTTATCTCTCCCGCCAATTTGCTTTTGACTGGGATTTTACCTCCTTTTGACCg                                                      <  5:46846/139‑1 (MQ=255)                                                         aTTCCCGGTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATACTTTATAGATTTTCAAGAGGTTATCTCTCCCGCCAATTTGCTTTTGACTGGGATTTTACCTCCTTTTGACCg                                                      <  5:199999/139‑1 (MQ=255)                                                         aTTCCCGGTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATACTTTATAGATTTTCAAGAGGTTATCTCTCCCGCCAATTTGCTTTTGACTGGGATTTTACCTCCTTTTGACCg                                                      <  1:244309/139‑1 (MQ=255)                                                         aTTCCCGGTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATACTTTATAGATTTTCAAGAGGTTATCTCTCCCGCCAATTTGCTTTTGACTGGGATTTTACCTCCTTTTGACCg                                                      <  1:140282/139‑1 (MQ=255)                                                         aTTCCCGGTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATACTTTATAGATTTTCAAGAGGTTATCTCTCCCGCCAATTTGCTTTTGACTGGGATTTTACCTCCTTTTGACCg                                                      <  2:53274/139‑1 (MQ=255)                                                         aTTCCCGGGCCCATCACGATTTTATCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATACTTTATAGATTTTCAAGAGGTTATCTCTCCCGCCAATTTGCTTTTGACTGGGATTTTACCTCCTTTTGACCg                                                      <  6:282652/139‑1 (MQ=255)                                                            cccggTCCCATGACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATACTTTATAGATTTTCAAGAGGTTATCTCTCCCGCCAATTTGCTTTTGACTGGGATTTTACCTCCTTTTGACCGcc                                                    <  2:277642/135‑1 (MQ=255)                                                            cccggTCCCATCACGTTTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATACTTTATAGATTTTCAAGAGGTTATCTCTCCCGCCAATTTGCTTTTGACTGGGATTTTACCTCCTTTTGACCGcc                                                    <  3:242217/135‑1 (MQ=255)                                                            cccggTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATACTTTATAGATTTTCAAGAGGTTATCTCTCCCGCCAATTTGCTTTTGACTGGGATTTTACCTCCTTTTGACCGcc                                                    <  5:271479/135‑1 (MQ=255)                                                            cccggTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATACTTTATAGATTTTCAAGAGGTTATCTCTCCCGCCAATTTGCTTTTGACTGGGATTTTACCTCCTTTTGACCGcc                                                    <  7:95370/135‑1 (MQ=255)                                                            cccggTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATACTTTATAGATTTTCAAGAGGTTATCTCTCCCGCCAATTTGCTTTTGACTGGGATTTTACCTCCTTTTGACCGcc                                                    <  7:27985/135‑1 (MQ=255)                                                            cccggTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATACTTTATAGATTTTCAAGAGGTTATCTCTCCCGCCAATTTGCTTTTGACTGGGATTTTACCTCCTTTTGACCGCCa                                                   >  8:95271/4‑139 (MQ=255)                                                                 tCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATACTTTATAGATTTTCAAGAGGTTATCTCTCCCGCCAATTTGCTTTTGACTGGGATTTTCCGTCCTTTTGACCGCCAGGTa                                               >  7:116669/1‑138 (MQ=255)                                                                           aTTTTTTCCTGATAATCCTATTTTATTAAATCCGCATTGTTGTTTAATACTTTATAGATTTTCAAGAGGTTATCTATACCGCCAATTTGCTTTTGACTGTGATTTTACCTCCTTTTGACCGCCATGTAAAGATGCGGGc                                    <  7:201718/139‑1 (MQ=255)                                                                               tttCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATACTTTATAGATTTTCAAGAGGTTATCTCTCCCGCCAATTTGCTTTTGACTGGGATTTTACCTCCTTTTGACCGCCATGTAAAGATGCGGGcccg                                 <  3:93142/138‑3 (MQ=255)                                                                                          acagctTTTATTAAAGCCGCATTGTTGTTTAATACTTTATAGATTTTCAAGAGGTTATCTCTCCCGCCAATTTGCTTTTGACTGGGATTTTACCTCCTTTTGACCGCCATGTAAAGATGCGGGCCCGCTCTCTCCTGCt                     >  3:236843/5‑139 (MQ=255)                                                                                                 ttattaAAGCCGCATTGTTGTTTAATACTTTATAGATTTTCAAGAGGTTATCTCTCCCGCCAATTTGCTTTTGACTGGGATTTTACCTCCTTTTGACCGCCATGTAAAGATGCGGGCCCGCTCTCTCCTGCTTTGGCg               <  4:232778/138‑1 (MQ=255)                                                                                                  tattaAAGCCGCATTGTTGTTTAATACTTTATAGATTTTCAAGAGGTTATCTCTCCCGCCAATTTGCTTTTGACTGGGATTTTACCTCCTTTTGACCGCCATGTAAAGATGCGGGCCCGCTCTCTCCTGCTTTGGCGa              <  5:216827/138‑1 (MQ=255)                                                                                                  tattaAAGCCGCATTGTTGTTTAATACTTTATAGATTTTCAAGAGGTTATCTCTCCCGCCAATTTGCTTTTGACTGGGATTTTACCTCCTTTTGACCGCCATGTAAAGATGCGGGCCCGCTCTCTCCTGCTTTGGCGa              <  2:103484/138‑1 (MQ=255)                                                                                                  aattaAAGCCGAATTGTTGTTTAATACTTTATAGATTTTCAAGAGGTAATCTCTCCCGCCAATTTGCTTTTGACTGGGATTTTACCTCCTTTTGACCGCCATGTAAAGATGCGGGCCCGCTCTATCCTGCTTTGGCGa              <  8:116669/137‑1 (MQ=255)                                                                                                         gCCGCATTGTTGTTTAATACTTTATAGATTTTCAAGAGGTTATCTCTCCCGCCAATTTGCTTTTGACTGGGATTTTACCTCCTTTTGACCGCCATGTAAAGATGCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTg       >  6:191012/1‑138 (MQ=255)                                                                                                             cATTGTTGTTTAATACTTTATAGATTTTCAAGAGGTTATCTCTCCCGCCAATTTGCTTTTGACTGGGATTTTACCTCCTTTTGACCGCCATGTAAAGATGCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGc  >  6:276360/1‑139 (MQ=255)                                                                                                             cATTGTTGTTTAATACTTTATAGATTTTCAAGAGGTTATCTCTCCCGCCAATTTGCTTTTGACTGGGATTTTACCTCCTTTTGACCGCCATGTAAAGATGCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGc  >  4:212744/1‑139 (MQ=255)                                                                                                                         |                                                                                                                               AGATCTTTTCCACTTTATTCATAGGGATGAGGTCATTTTTCTTTTTTTGAAAAGCGATTCCTGGTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAGTACTTTATAGATTTTCAAGAGGCTATCTCTCCCGCCAATTAGCTTTTGACTGGGATTTTACCTCCTTTTGACCGCCATGTAAAGATGCGGGCCAGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGC  >  CP006881/3724688‑3724934

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