Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,375,528	0	C	A	23.2%	121.0 / 11.2	56	D220E (GAC→GAA)	dcuD	putative transporter DcuD
Reads supporting (aligned to +/- strand):  ref base C (29/14);  new base A (1/12);  total (30/26)
Fisher's exact test for biased strand distribution p-value = 2.22e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TGGGTATGTTGCTGATGGTCACCTTATTTCCGACGCTGGTGAGTCTGGGAGTAAGTCGTCTCTCTGCGGTGGCAGTTATCGCAACCACGATGTCCATTGAGTGGGGGATTCTGGAAACGAACTCCATTTTTGCTGCCCAGGTAGCGGGAATGAAAATTGCCACATACTTCTTCCACTACCAGCTTCCGGTCGCCTCTTGCGTCATTATCTCGGTGGCGATCTCCCACTTTTTCGTGCAACGCGCTTTTGACAAAAAAGATAAAAATATCAATCACGAACAGG‑CAGAGCAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTTTGTTGAGTTCTTCCGCAAGCATAACTTGCGCGAAACAATGGACGATG  >  NC_000913/3375278‑3375772                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                      tGGGTATGTTGCTGATGGTCACCTTATTTCCGACGCTGGTGAGTCTGGGAGTAAGTCGTCTCTCTGCGGTGGCAGTTATCGCAACCACGATGTCCATTGAGTGGGGGATTCTGGAAACGAACTCCATTTTTGCTGCCCAGGTAGCGGGAATGAAAATTGCCACATACTTCTTCCACTACCAGCTTCCGGTCGCCTCTTGCGTCATTATCTCGGTGGCGATCTCCCACTTTTTCGTGCAACGCGCTTTTGAc                                                                                                                                                                                                                                                       >  1:168568/1‑251 (MQ=255)  gggTATGTTGCTGATGGTCGCCTTATTTCCGACGCTGGGGAGTCTGGGAGTAAGTCGTCTCTCTGCTGTGGCAGTTATCGCAACCACGATGTCCATTGAGTGGGGGATTCTGGAAACGAACTCCATTTTTGCTGCCCAGGTAGCGGGAATGAAAATTGCCACATACTTCTTCCACTACCAGCTTCCGGTCTCCTCTTGCGTCATTATCTCGGTGGCGATCTCCCACTTTTTCGTGCAACGCGCTTTTGACa                                                                                                                                                                                                                                                      <  2:300610/251‑1 (MQ=255)   ggTATGTTGCTGATGGTCACCTTATTTCCGACGCTGGTGAGTCTGGGAGTAAGTCGTCTCTCTGCGGTGGCAGTTATCGCAACCACGATGTCCATTGAGTGGGGGATTCTGGAAACGAACTCCATTTTTGCTGCCCAGGTAGCGGGAATGAAAATTGCCACATACTTCT‑CCACTACCAGCTTCCGGTCGCCTCTTGCGTCATTATCTCGGTGGCGATCTCCCACTTTTTCGTGCAACGCGCTTTTGACaaa                                                                                                                                                                                                                                                    <  1:214885/251‑1 (MQ=255)    gTATGTTGCTGATGGTCACCTTATTTCCGACGCTGGTGAGTCTGGGAGTAAGTCGTCTCTCTGCGGGGGCAGTTATCGCAACCACGATGTCCATTGAGTGGGGGATTCTGGAAACGAACTCCATTTTTGCTGCCCAGGTAGCGGGAATGAAAATTGCCACATACTTCTTCCACTACCAGCTTCCGGTCGCCTCTAGCGTCATTATCTCGGTGGCGATCTCCCACTTTTTCGTGCAACGCGCTTTTGACaa                                                                                                                                                                                                                                                     >  1:129366/1‑250 (MQ=255)      aTGTTGCTGATGGTCACCTTATTTCCGACGCTGGTGAGTCTGGGAGTAAGTCGTCTCTCTGCGGTGGCAGTTATCGCAACCACGATGTCCATTGAGTGGGGGATTCTGGAAACGAACTCCATTTATGCTGCCCAGGTAGCGGGAATGAAAATTGCCACATACTTCTTCCACTACCAGCTTCCGGTCGCCTCTTGCGTCATTATCTCGATGGCGATCTCCCACTTTTTCGTGCAACGCGCTTTTGTCaaaa                                                                                                                                                                                                                                                   >  2:330742/1‑250 (MQ=255)                        ttATTTCCGACGCTGGTGAGTCTGGGAGTAAGACGTCTCTCTGCGGTGGCAGTTATCGCAACCACGATGTCCATTGAGTGGGGGATTCTGGAAACGAACTCCATTTTTGCTGCCCAGGTAGCGGGAATGAAAATTGCCACATACTTCTTCCACTACCAGCTTCCGGTCGCCTCTTGCGTCATTATCTCGGTGGCGATCTCCCACTTTTTCGTGCAACGCGCTTTTGACAAAAAAGATAAAAATatcaatc                                                                                                                                                                                                                                 >  1:358232/1‑250 (MQ=255)                          aTTTGCGACGCTGGTGAGTCTGGGAGTAAGCCGTCTCTCTGCGGTGGCAGTTATCGCAACCACGATGTCCATTGAGTGGGGGATTCTGGAAACGAACTCCATTTTTGCTGCCCAGGTAGCGGGAATGAAAACTGCCTCATACTTCTGCCACTACCAGCATCCGGTCGCCTCTTGCGTCATTATCTCGGTGGCGATCTCCCTCTTTTTCGTGCAACGCGCTTTTGACAAAAAAGATAAAAAGatcaatcgcg                                                                                                                                                                                                                              <  2:192916/251‑4 (MQ=255)                          aTTTCCGACGCTGGTGAGTCTGGGAGTAAGTCGTCTCTCTGCGGTGGCAGTTATCGCAACCACGATGTCCATTGAGTGGGGGATTCTGGAAACGAACTCCATTTTTGCTGCCCAGGTAGCGGGAATGAAAATTGCCACATACTTCTTCCACTACCAGCTTCCGGTCGCCTCTTGCGTCATTATCTCGGTGGCGATCTCCCACTTTTTCGTGCAACGCGCTTTTGACAAAAAAGATAAAAATATCAATCAc                                                                                                                                                                                                                               <  1:386107/250‑1 (MQ=255)                             tCCGACGCTGGTGAGTCTGGGAGTAAGTCGTCTCTCTGCGGTGGCAGTTATCGCAACCACGATGTCCATTGAGTGGGGGATTCTGGAAACGAACTCCATTTTTGCTGCCCAGGTAGCGGGAATGAAAATTGCCACATACTTCTTCCACTACCAGCTTCCGGTCGCCTCTTGCGTCATTATCTCGGTGGCGCTCTCCCACTTTTTCGTGCAACGCGCTTTTGACAAAAAAGATAAAAATATCAATCACGaa                                                                                                                                                                                                                            >  2:259480/1‑250 (MQ=255)                                 aCGCTGGTGAGTCTGGGAGTAAGTCGTCTCTCTGCGGTGGCAGTTATCGCAACCACGATGTCCATTGAGTGGGGGATTCTGGAAACGAACTCCATTTTTGCTGCCCAGGTAGCGGGAATGAAAATTGCCACATACTTCTTCCACTACCAGCTTCCGGTCGCCTCTTGCGTCATTATCTCGGTGGCGATCTCCCACTTTTTCGTGCAACGCGCTTTTGACAAAAAAGATAAAAATATCAATCACGAACAGG‑c                                                                                                                                                                                                                      >  1:499303/1‑251 (MQ=255)                                  cGCTGGTGAGTCTGGGAGTAAGTCGTCTCTCTGCGGTTGCAGTGATCGTAACCATGATGTCCATTGAGTGGGGGAGTCTGGAAACGAACTCCATTTTTGCTGCCCAGGTAGCGGGAATGAAAATTGCCACATGCTTCTTCCACTACCAGCTTCCGGTCGCCTCTTGCGTCATTATCTCGGTGGCGATCTCCCACTTTTTCGTGCAACGCGCTTTTGACAAAAAAGATAAAAAGATCAATCACGAACAGG‑Ca                                                                                                                                                                                                                     <  2:213908/251‑1 (MQ=255)                                        ggAGTCTGGGAGTAAGTCGTCTCTCTGCGGTGGCAGTTATCGCAACCACGATGTCCATTGAGTGGGGGATTCTGGAAACGAACTCCATTTTTGCTGCCCAGGTAGCGGGAATGAAAATTGCCACATACTTCTTCCACTACCAGCTTCCGGTCGCCTCTTGCGTCATTATCTCGGTGGCGATCTCCCACTTTTTCGTGCAACGCGCTTTTGACAAAAAAGATAAAAATATCAATCACGAACAGG‑CAGAGCaa                                                                                                                                                                                                               <  1:35074/250‑1 (MQ=255)                                                 gAGTAAGTCGTCTCTCTGCGGTGGCAGTTATCGCAACCACGATGTCCATTGAGTGGGGGATTCTGGAAACGAACTCCATTTTTGCTGCCCAGGTAGCGGGAATGAAAATTGCCACATACTTCTTCCACTACCAGCTTCCGGTCGCCTCTTGCGTCATTATCTCGGTGGCGATCTCCCACTTTTTCGTGCAACGCGCTTTTGACAAAAAAGATAAAAATATCAATCACGAACAGG‑CAGAGCAAAAAGCtctc                                                                                                                                                                                                      >  2:378740/1‑251 (MQ=255)                                                              ctctGCGGTGGCAGTTATCGCAACCACGATGTCCATTGAGTGGGGGATTCTGGAAACGAACTCCATTTTTGCTGCCCAGGTAGCGGGAATGAAAATTGCCACATACTTCTTCCACTACCAGCTTCCGGTCGCCTCTTGCGTCATTATCTCGGTGGCGATCTCCCACTTTTTCGTGCAACGCGCTTTTGACAAAAAAGATAAAAATATCAATCACGAACAGG‑CAGAGCAAAAAGCTCTCGATAATGTcc                                                                                                                                                                                            <  1:105894/248‑1 (MQ=255)                                                                ctGCGGTGGCAGTTATCGCAACCACGATGTCCATTGAGTGGGGGATTCTGGAAACGAACTCCATTTTTGCTGCCCAGGTAGCGGGAATGAAAATTGCCACATACTTCTTCCACTACCAGCTTCCGGTCGCCTCTTGCGTCATTATCTCGGTGGCGATCTCCCACTTTTTCGTGCAACGCGCTTTTGACAAAAAAGATAAAAATATCAATCACGAACAGG‑CAGCGCAAAAAGCTCTCGATAATGTCccgccg                                                                                                                                                                                       >  2:418232/1‑251 (MQ=255)                                                                  gCGGTGGCAGTTATCGCAACCACGATGTCCATTGAGTGGGGGATTCTGGAAACGAACTCCATTTTTGCTGCCCAGGTAGCGGGAATGAAAATTGCCACATACTTCTTCCACTACCAGCTTCCGGTCGCCTCTTGCGTCATTATCTCGGTGGCGATCTCCCACTTTTTCGTGCAACGCGCTTTTGACAAAAAAGATAAAAATATCAATCACGAACAGG‑CAGAGCAAAAAGCTCTCGATAATGTCccgcc                                                                                                                                                                                        <  1:160467/248‑1 (MQ=255)                                                                     gTGGCAGTTATCGCAACCACGATGTCCATTGAGTGGGGGATTCTGGAAACGAACTCCATTTTTGCTGCCAAGGTAGCGGGAATGAAAATTGCCACATACTTCTTCCACTACCAGCTTCCGGTCGCCTCTTGCGTCATTATCTCGGTGGCGATCTCCCACTTTTTCGTGCAACGCGCTTTTGACAAAAAAGATAAAAATATCAATCACGAACAGG‑CAGAGCAAAAAGCTCTCGCTAATGTCCCGCCGctct                                                                                                                                                                                   >  1:272194/1‑250 (MQ=255)                                                                       ggCAGTTATCGCAACCACGATGTCCATTGAGTGGGGGATTCTGGAAACGAACTCCATTTTTGCTGCCCAGGTAGCGGGAATGAAAATTGCCACATACTTCTTCCACTACCAGCTTCCGGTCGCCTCTTGCGTCATTATCTCGGTGGCGATCTCCCACTTTTTCGTGCAACGCGCTTTTGACAAAAAAGATAAAAATATCAATCACGAACAGG‑CAGAGCAAAAAGCTCTCGATAATGTCCCGCCGCTCTAtt                                                                                                                                                                                >  2:361522/1‑251 (MQ=255)                                                                             tATCGCAACCACGATGTCCATTGAGTGGGGGATTCTGGAAACGAACTCCATTTTTGCTGCCCAGGTAGCGGGAATGAAAATTGCCACATACTTCTTCCACTACCAGCTTCCGGTCGCCTCTTGCGTCATTATCTCGGTGGCGATCTCCCACTTTTTCGTGCAACGCGCTTTTGACAAAAAAGATAAAAATCTCAATCACGAACAGG‑CAGCGCAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGcc                                                                                                                                                                           >  2:271535/1‑250 (MQ=255)                                                                              aTCGCAACCACGATGTCCATTGAGTGGGGGATTCTGGAAACGAACTCCATTTTTGCTGCCCAGGTAGCGGGAATGAAAATTGCCACATACTTCTTCCACTACCAGCTTCCGGTCGCCTCTTGCGTCATTATCTCGGTGGCGATCTCCCACTTTTTCGTGCAACGCGCTTTTGACAAAAAAGATAAAAATATCAATCACGAACAGG‑CAGAGCAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCAt                                                                                                                                                                         >  2:234374/1‑251 (MQ=255)                                                                                        cGATGTCCCTTGAGTGGGGGATTCTGGAAACGAACTCCATTTTTGCTGCCCAGGTAGCGGGAATGGAAATTTCCACATACTTCTTCCCCTACCAGCTTCCGGTCGCCTCTTGCGTCATTATCTCGGTGGCGATCTCCTACTTTTTCTTGCAACGCGCTTTTGAAAAAAAAGATAAAAATATCAATCACGAACAGG‑CAGAGCAAAAAGCACTCGAAAATGGCCCGCCACTCAATCACGCCATTTTACCtgtg                                                                                                                                                               >  2:443456/1‑251 (MQ=255)                                                                                                    aGTGGGGGATTCTGGAAACGAACTCCATTTTTGCTGCCCAGGTAGCGGGAATGAAAATTGCCACATACTTCTTCCACTACCAGCTTCCGGTCGCCTCTTGCGTCATTATCTCGGTGGCGATCTCCCACTTTTTCGTGCAACGCGCTTTTGACAAAAAAGATAAAAATATCAATCACGAACAGG‑CAGAGCAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATc                                                                                                                                                   >  1:144614/1‑251 (MQ=255)                                                                                                              tctgaaaacaatccccATTTTTGCTGCCCAGGTAGGGGGAATGAAATTTGCCACATACTTCTTCCGCTCCCAGCTTCCGGTCGCCTCTTGCGTCATTATCGCGGTGGCGATCTCCCACTTTTTCGTGCAACGCGCTTTTGACAAAAAAGATAAAAATATCAACCTCGAACAGG‑CAGAGCAAAACGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATgctc                                                                                                                                          <  1:292536/236‑1 (MQ=255)                                                                                                                 ggAAACGAACTCCATTTTTGCTGCCCAGGTAGCGGGAATGAAAATTGCCACATACTTCTTCCACTACCAGCTTCCGGTCGCCTCTTGCGTCATTATCTCGGTGGCGATCTCCCACTTTTTCGTGCAACGCGCTTTTGACAAAAAAGATAAAAATATCAATCACGAACAGG‑CAGAGCAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATgctcggct                                                                                                                                      >  1:343632/1‑251 (MQ=255)                                                                                                                    aaCGAACTCCATTTTTGCTGCCCAGGTAGCGGGAATGAAAATTGCCACATACTTCTTCCACTACCAGCTTCCGGTCGCCTCTTGCGTCATTATCTCGGTGGCGATCTCCCACTTTTTCGTGCAACGCGCTTTTGACAAAAAAGATAAAAATATCAATCACGAACAGG‑CAGAGCAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGc                                                                                                                                   >  2:510644/1‑251 (MQ=255)                                                                                                                       gAACTCCATTTTTGCTGCCCAGGTAGCGGGAATGCAAATTGCAACATACTTCTTCCACTACCAGCTTCCGGTCGCCTCTTGCGTCATTATCTCGGTGGCGATCTCCCACTTTTTCGTGCAACGCGCTTTTGACAACAAAGATAAAAATATCAATCACGAACAGG‑CAGAGCAACAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGACGTTAATCCTGATGCTCGGCTCGCTGt                                                                                                                                >  2:286922/1‑251 (MQ=255)                                                                                                                       gAACTCCATTTTTGCTGCCCAGGTAGCGGGAATGAAAATTGCCACATACTTCTTCCACTACCAGCTTCCGGTCGCCTCTTGCGTCATTATCTCGGTGGCGATCACCCACTTTTTCGTGCAACGCGCTTTTGACACAAATGATAAAAATATCAATCACGAACAGG‑CAGAGCAAAAAGCTCTCGATCATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTACTCCTGATGCTCGGCTCGCTGt                                                                                                                                >  2:54648/1‑251 (MQ=255)                                                                                                                          cTCCATTTTTGCTGCCCAGGTAGCGGGTATGAAAATTGCCACATACTTCTTCCACTACCAGCTTCCGGGCGCCTCTTGCGTCATTATCTCGGTGCCGATCTCACACTTTTTCGTGCAACGCGCTTTTGACAAAAAAGATAAAGATATCAATCACGAACAGC‑CAGAGCAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCAATTTTACCTATGATGCCGTTAATACTGAAGCTCGGCTCGCTTTTcc                                                                                                                             >  2:470917/1‑251 (MQ=255)                                                                                                                           tCCATTTTTGCTGCCCAGGTAGCGGGAATGCAAATTGCCAAATACTTCTTCCACTACCAGCTTCCGGTCGCCTCTTGCGTCATTATCTCGGTGGCGATCTCCCACTTTTTCGTGCAACACGCTTTTGACAAAAAAGATAAAAATATCAATCACGAACAGG‑CAGAGCAAAAAGCTCTCGATAATGACACGCCGCTCTAGTACGCCATTTTACCTGTGATTCCGTTAATCCTGATGCTCGTCTCGCTGtttct                                                                                                                            >  2:377725/1‑248 (MQ=255)                                                                                                                              aTTTTTGCTGCCCGGGTAGGGGGAATGAAAATGGCCACATCCTTCTTCCACTCCCAGCTTCCGGTCGCCTCTTGCGTCATGATCGCGTGGGCGCTCTCCCATTTTTTCGGGCACCGCGCTTTTGAAAAAAAAGATAAAAATATCAACCCCGAAAAGG‑AAGAGAAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTc                                                                                                                           <  1:128101/249‑1 (MQ=255)                                                                                                                                                ggggaaaTGAAAAGTGCCCCCTTCTTCTTCCCCTCCCTGCTTCCGGTCCCCTCTTGCGTCATTATCTCGGTGGGGACCCCCCCCTTTTCCGTGAAACGCGTTTTTGACAAAAAAGAAAAAAATATCAATCCCGAACAGG‑CAGAGCAAAAAGCTCTCAATAATGCCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCa                                                                                                       <  1:2470/246‑1 (MQ=255)                                                                                                                                                   ggAATGAAAATTGCCACATACTTCTTCCACTACCAGCTTCCGGTCGCCTCTTGCGTCATTATCTCGGTGGCGATCTCCCACTTTTTCGTGCAACGCGCTTTTGACAAAAAAGATAAACATATCAATCACGAACAGG‑CAGAGCAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTc                                                                                                    >  2:446031/1‑251 (MQ=255)                                                                                                                                                             ttGCCACATACTTCTTCCACTACCAGCTTCCGGTCGCCTCTTGCGTCATTATCTCGGTGGCGATCTCCCACTTTTTCGTGCAACGCGCTTTTGACACAAAAGATAAAAATATCAATCACGAACAGG‑CAGCGCAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCACCGTCGGGCTGATGCCGTCAGCATTGCAt                                                                                          >  2:63832/1‑251 (MQ=255)                                                                                                                                                                      aCTTCTCCCACTACCACCTTCCGGTCGCCTCTTGCGTCATGACCTCGGGGTCTATCTCCCTCTTTTTGGTGCACCGCCCTTTTGAAAAAAAAGATAAAAATATCAATCACGAACAGG‑CAGGGCAAAAAGCCCTCTATAGTGTTCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCtggtggtg                                                                                 <  1:448699/251‑1 (MQ=255)                                                                                                                                                                         tctgccccTACCAGCTTCCGGTCGCCTCTTGCGTCATTATCTCGGTGGCGATCTCCCTCTTTTTCGTGCAACGCGCTTTTGACAAAAAAGATAAAAATATCAACCACGAACAGG‑GAGAGCAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCtggtggtgg                                                                                <  1:474644/242‑1 (MQ=255)                                                                                                                                                                         tctccctcTACCAGCTTCCGGTCGCCTCTTGCGTCATTACCTCGGTGGCGTCCTCCCCCTTTTTCGTGCAACGCGCTTTTGAAAAAAAAGATAAAAATATCAATCCCGAACAGG‑CAGAGAAAAAAGTTCCCGATAATGCCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCtggtggtg                                                                                 <  1:301157/241‑1 (MQ=255)                                                                                                                                                                         tctccccctcccAGCTTCCGGTCGCCTCTTGCGTCATTATCTCGGTGGCGATCTCCCACTTTTTCGTGCAACGCGCTTTTGAAAAAAAAAAAAAAAATATCAATCACGAACAGG‑CAGAGCAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTggtggtggt                                                                               <  1:234374/240‑1 (MQ=255)                                                                                                                                                                                 tACCAGCTTCCGGTCGCCTCTTGCGTCATTATCTCGGTGGCGATCTCCCACTTTTTCGTGCAACGCGCTTTTGACAAAAAAGATAAAAATATCAATCACGAACAGG‑CAGAGCAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACt                                                                      >  2:300216/1‑251 (MQ=255)                                                                                                                                                                                     aGCTTCCGGTCGCCTCTTGCGTCTTTATCCCGGTGGCGATCTCCCACTTTTCCGTGCACCGCGCTTTTAACAAAAAAGATAAAAATATCAATCACGAACAGG‑CAGAGAAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGt                                                                  <  1:268218/251‑1 (MQ=255)                                                                                                                                                                                         tCCGGTCGCCTCTTGCGTCATTATCTCGGTGGCGATCTCCCACTTTTTCGTGCACCGCGTTTTTGAAAAAAAAGATAAAAATATCAATCACGAACAGG‑CAGAGCAAAAATCTCTCTATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGttt                                                                <  1:257450/249‑1 (MQ=255)                                                                                                                                                                                           gggtggcctttgccgtaaTTATCTGGGTGGCGATCTCCCCTTTTTTGGTGAAACGCGCTTTGGAAAAAAAATAAAAAAATATCAACCACGACCGGG‑AAGACCAAAAAGCTCTGTATAATTTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTGCCTCGCCCACGTCGGGTTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGAcg                                                            <  2:11515/234‑2 (MQ=255)                                                                                                                                                                                            ggtcgcctgtagcgtagtttcctcGGTGCCCATCTCCCACTTTTTGGGGCACCGCGTTTTTGAAAAAAAAGAAAAAAATATCAA‑CATGGA‑AGGACAGAGCAAAACGCTCGCGATATTGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACtgt                                                          <  2:101602/230‑1 (MQ=255)                                                                                                                                                                                            ggTCGCCTCTTGCGTCATTATCTCGGTGGCGATCTCCCACTTTTTCGTGCAACGCGCTTTTGACAAAAAAGATAAAAATATCAATCACGAACAGG‑CAGAGCAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACtg                                                           >  1:321096/1‑251 (MQ=255)                                                                                                                                                                                                  ctctTGCGTCATTATCTCGGTGGCGATCTCCCACTTTTTCGTGCAACGCGCTTTTGACAAAAAAGATAAAAATATCAATCACGAACAGG‑CAGAGCAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTAAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACg                                                      >  1:476811/1‑250 (MQ=255)                                                                                                                                                                                                       gCGTCATTATCTCGGTGGCGATCTCCCACTTTTTCGTGCAACGCGCTTTTGACAAAAAAGATAAAAATATCAATCACGAACAGG‑CAGAGCAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGCTGCAGTCAGAACTGCATCTGTTGGTGGTGCTGTTACTGAGTTTGACTGTGACGAtgttt                                                >  2:181165/1‑251 (MQ=255)                                                                                                                                                                                                       gCGTCATTATCTCGGTGGCGATCTCCCACTTTTTCGTGCAACGCGCTTTTGACAAAAAAGAAAAAAATATCAATCACGAACAGG‑CAGAGAAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGAtgtt                                                 <  1:319457/250‑1 (MQ=255)                                                                                                                                                                                                        cGTCATTATCTCGGTGGCGATCTCCCTCTTTTTCGTGCACCGCGCTTTTGACAAAAAAGATAAAAATATCAATCACGAACAGG‑CAGAGAAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGAtgtttg                                               <  1:418232/251‑1 (MQ=255)                                                                                                                                                                                                        cGTCATTATCCCGGTGGCGATCTCCCACTTTTCCGTGCACCGCGCTTTTGAAAAAAAAGAAAAAAATATCAATCACGAACAGG‑CAGAGCAAAAAGCTCTCGAAAATGCCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGAtgtttg                                               <  1:361239/251‑1 (MQ=255)                                                                                                                                                                                                             ttATCTCGGTGGCGATCTCCCACTTTTTCGTGCAACGCGCTTTTGACAAAAAAGATAAAAATATCAATCACGAACAGG‑CAGAGCAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTttgttg                                            >  1:198210/1‑249 (MQ=255)                                                                                                                                                                                                                 cgcGGTGGCACTCTCCCACTTTTTCGGCCACCGCTTTGTTGAAAAAAAAGAAAAAAATATCAAACAAGAACAGG‑CAAATAAAAAACCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTTTGTTGAGtt                                        <  2:168568/247‑1 (MQ=255)                                                                                                                                                                                                                   cgggggCGATCCCCCACTTTTTCGTGCAACGCGTTTTTGAAAAAAAAGATAAAAATATCAATCACGAACAGG‑CAGAGCAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTTTGTTGAGttc                                       <  1:330742/244‑1 (MQ=255)                                                                                                                                                                                                                     ggggcaatcacccACTTTTTGGCGAACCGCGCTTGTAAAAAAAAAGAAAAAAAAAGCAACCACGAAAGGG‑AAGAGCAAAACGTTCCAGATAAGGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTTTGTTGAGTTCTTCCg                                  <  2:144614/241‑1 (MQ=255)                                                                                                                                                                                                                       ggCGATCTCCCATTTTTTCGTGCAACGCGCTTTTGAAAAAAAAGATAAAAATATCAATCACGAACGGG‑CAGAGCAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTTTGTTGAGTTCTTCCGCa                                <  1:210126/251‑1 (MQ=255)                                                                                                                                                                                                                           aTCTCCCACTTTTTCGTGCAACGCGCTTTTGACAAAAAAGATAAAAATATCAATCACGAACAGG‑CAGAGCAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTTTGTTGAGTTCTTCCGCAAGc                             >  1:22877/1‑250 (MQ=255)                                                                                                                                                                                                                                          gTGCAACGCGCTTTTGACAAAAAAGATAAAAATATCAATCACGAACAGG‑CAGAGCAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTTTGTTGAGTTCTTCCGCAAGCATAACTTGCGCGAAAc             >  1:295008/1‑251 (MQ=255)                                                                                                                                                                                                                                             cAACGCGCTTTTGACAAAAAAGATAAAAATATCAATCACGAACAGG‑CAGAGCAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTTTGTTGAGTTCTTCCGCAAGCATAACTTGCGCGAAACAAt          >  1:65220/1‑251 (MQ=255)                                                                                                                                                                                                                                                     ttttGACAAAAAAGATAAAAATATCAATCACGAACAGG‑CAGAGCAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTTTGTTGAGTTCTTCCGCAAGCATAACTTGCGCGAAACAATGGACGAtg  >  1:496837/1‑251 (MQ=255)                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                      TGGGTATGTTGCTGATGGTCACCTTATTTCCGACGCTGGTGAGTCTGGGAGTAAGTCGTCTCTCTGCGGTGGCAGTTATCGCAACCACGATGTCCATTGAGTGGGGGATTCTGGAAACGAACTCCATTTTTGCTGCCCAGGTAGCGGGAATGAAAATTGCCACATACTTCTTCCACTACCAGCTTCCGGTCGCCTCTTGCGTCATTATCTCGGTGGCGATCTCCCACTTTTTCGTGCAACGCGCTTTTGACAAAAAAGATAAAAATATCAATCACGAACAGG‑CAGAGCAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTTTGTTGAGTTCTTCCGCAAGCATAACTTGCGCGAAACAATGGACGATG  >  NC_000913/3375278‑3375772

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