Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,896,475	0	A	C	34.5%	82.7 / 11.6	55	G27G (GGT→GGG)	ygcR	putative flavoprotein
Reads supporting (aligned to +/- strand):  ref base A (25/11);  new base C (0/19);  total (25/30)
Fisher's exact test for biased strand distribution p-value = 2.47e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TTCGGCAATATGGCGAGCGACAAATTCCGGCGCAAAGCGCAGATCTGCCGCCGTCTCCAGCAATACAGCTTCCTCAAACCCAAGAGCCATGAGATAGCGCAGCCAGTGCAATGCCCGTTCATCCCCCATGCTCAACGCGGTTAAAGACATCGGCGTACCGTTTTTTCTCTGCGCCAACAGCAGCGCGGCAGCGGCCTGTT‑CATCAGCACCGAGTAAACTTCGCAGTAGCGAAATATCCGGTCCGCTTTT‑ACCCTGAGCCGCCGCCTGCCACTCTTTTTCCGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCATTCCCCACTGCGCCAGCACCCACGGCAGCAGGCCAGTACTT  >  NC_000913/2896226‑2896716                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                  ttCGGCAATATGGCGAGCGACAAATTCCGGCGCAAAGCGCAGATCTGCCGCCGTCTCCAGCAATACAGCTTCCTCAAACCCAAGAGCCATGAGATAGCGCAGCCAGTGCAATGCCCGTTCATCCTCCATGCTCAACGCGGTTAAAGACATCGGCGTACCGTTTTTTCTCTGCGCCAACAGCAGCGCGGCAGCGGCCTGTT‑CATCAGCACCGAGTAAACTTCGCAGTAGCGAAATATCCGGTCCGCTTTT‑Ac                                                                                                                                                                                                                                                  >  2:393492/1‑251 (MQ=255)   cGGCAATATGGCGAGCGACAAATTCCGGCGCAAAGCGCAGATCTGCCGCCGTCTCCAGCAATACAGCTTCCTCAAACCCAAGAGCCATGAGATAGCGCAGCCAGTGCAATGCCCGTTCATCCCCCATGCTCAACGCGGTTAAAGACATCGGCGTACCGTTTTTTCTCTGCGCCAACAGCAGCGCGGCAGCGGCCTGTT‑CATCAGCACCGAGTAAACTTCGCAGTAGCGAAATATCCGGTCCGCTTTT‑Acc                                                                                                                                                                                                                                                 <  1:18373/250‑1 (MQ=255)                     aaaaTGCCGGCGCAAGGCGCAGATCTGCCGCCGTCTCCAGCAATACAGCTTCCTCAAACCCAAGAGCCATGAGATAGCGCAGCCAGTGCAATGCCCGTTCATCCCCCATGCTCAACGCGGTTAAAGACATCGGCGTACCGTTTTTTCTCTGCGCCAACAGCAGCGCGGCAGCGGCCTGTT‑CATCAGCACCGAGTAAACTTCGCAGTAGCGAAATATCCGGTCCGCTTTT‑CCCCTGAGCCGCCGCCTGCCAc                                                                                                                                                                                                                              <  2:453806/250‑1 (MQ=255)                      aaaTTCCGGCGCAAAGCGCAGATCTGCCGCCGTCTCCAGCAATACAGCTTCCTCAACCCCAAGAGCCATGAGATAGCGCAGCCAGTGCAATGCCCGTTCATCCCCCATGCTCAACGCGGTTAAAGACATCGGCGTACCGTTTTTTCTCTGCGCCAACAGCAGCGCGGCAGCGGCCTGTT‑CATCAGCACCGAGTAAACTTCGCAGTAGCGAAATATCCGGTCCGCTTTT‑ACCCTGAGCCGCCGCCTGCCAct                                                                                                                                                                                                                             <  1:514498/251‑1 (MQ=255)                            cGGCGCAAAGCGCAGATCTGCCGCCGTCTCCAGCAATACAGCTTCCTCAAACCCAAGAGCCATGAGATAGCGCAGCCAGTGCAATGCCCGTTCATCCCCCATGCTCAACGCGGTTAAAGACATCGGCGTACCGTTTTTTCTCTGCGCCAACAGCAGCGCGGCAGCGGCCTGTT‑CATCAGCACCGAGTAAACTTCGCAGTAGCGAAATATCCGGTCCGCTTTT‑CCCCTGAGCCGCCGCCTGCCACTCttttt                                                                                                                                                                                                                       <  1:96520/251‑1 (MQ=255)                             ggCGCAAAGCGCAGACCTGCCGCCGTCCCCAGCAATACAGTTTCCTCAACCCCAAGAGCCATGAGATAGCGCAGCCAGTGCAATGCCCGTTCATCCCCCATGCTCAACGCGGTTAAAGACATCGGCGTACCGTTTTTTCTCTGCGCCAACAGCAGCGCGGCAGCGGCCTGTT‑CATCAGCACCGAGTAAACTTCGCAGTGGCGAAATATCCGGTCCGCTTTT‑ACCCTGAGCCGCCGCCTGCCACTCGTTTTc                                                                                                                                                                                                                      <  2:561273/251‑1 (MQ=255)                                 cAAAGCCCAGATCCGCCGCCTCCTCCAGCAATACAGCTTCCTCAAACCCAAGAGCCATGGGATAGCGCAGCCAGTGCAATTCCCGTTCACCCCCCATGCTCAACGCGGTTAAACACATCGGCGTACCTTTTTTTCTCTGCGCCAACAGCAGCGCGGCAGCGGCCTGTT‑CATCAGCCCCGATTAAACTTCGCAGTCGCGAAATATCCGCTCCGCTTTT‑CCCCTGAGCCGCCGCCTGCCACTCTTTTTCCGcc                                                                                                                                                                                                                  <  2:429450/251‑1 (MQ=255)                                  aaaGCGCAGATCTGCCGCCGTCTCCAGCAATACAGCTTCCTCAAACCCAAGAGCCATGAGATAGCGCAGCCAGTGCAATGCCCGTTCATCCCCCATGCTCAACGCGGTTAAAGACATCGGCGTACCGTTTTTTCTCTGCGCCAACAGCAGCGCGGCAGCGGCCTGTT‑CATCAGCACCGAGTAAACTTCGCAGTAGCGAAATATCCGGTCCGCTTTT‑ACCCTGAGCCGCCGCCTGCCACTCTTTTTCCGcc                                                                                                                                                                                                                  >  1:277592/1‑250 (MQ=255)                                           atctccgcccatctccagatgtACAGCTTCCTCAACCACCAGAGCCACGAGTGAGCGCGCCCAGTGCTATGCCCGTTCACCCCCCAGGCTCACGGCGGTTAAAGACGCCGGCGTCCCGTTTTTTCTCTGCGCAAACAGCAGAGCGGCAGCGGCCTTTT‑CATCAGCCCCGAGTAAACTTCGCAGTAGCGAAATATCCGCTCATCTTTT‑ACCCGGAGCACCCGCCTGCCACTCTTTTTCCGCCAGCATTCCg                                                                                                                                                                                                         <  2:344195/229‑1 (MQ=255)                                              tGCCGCCGTCTCCAGCAATACAGCTTCCTCAAACCCAAGAGCCATGAGATAGCGCAGCCAGTGCAATGCCCGTTCATCCCCCATGCTCAACGCGGTTAAAGACATCGGCGTACCGTTTTTTCTCTGCGCCAACAGCAGCGCGGCAGCGGCCTGTT‑CATCAGCACCGAGTAAACTTCGCAGTAGCGAAATATCCGGTCCGCTTTT‑ACCCTGAGCCGCCGCCTGCCACTCTTTTTCCGCCAGCATTCCGGCAt                                                                                                                                                                                                     >  1:60771/1‑251 (MQ=255)                                                   ccgTCTCAAGAAATACAGCTCCCTCAACCCCAAGAGCCATGAGATAGCGCGCCCAGTGCAATGCCCGTTCACCCCCCATGCTCAACGCGGTTAAAGACACCGGCGTACCGTTTTTTCTCTGCGCCAACAGCAGCGCGGCAGCGGCCTGTT‑CATCAGCACCGAGTAAACTTCGCAGTAGCGAAATATCCGGTCCTCTTTT‑CCCCTGAGCCGCCGCCTGCCACTCTTTTTCCGCCAGCATTCCGGCATCCGGt                                                                                                                                                                                                <  2:528705/251‑1 (MQ=255)                                                     gtcccaaGCAATACGGCTCCCACAACCCCAAGAGCCAGGAGAGAGCGCAGCCAGGGCAAAGCCCGTTCACCCCCCAGGCTCACCGGGGTTAAAGAAATCGTCGTACCGTGTTTTCTCTGCGCCAACAGCAGCGGGCCAGGGGCCGGTT‑CATCAACACCGAGTAAACTCCGCTTTCCCAAAATACCCGGACCCTTTTT‑CCCCTGCGCCGCCGCCTGCCACTCTTTTTCCGCCAGCATTCCGGCATCCGGTTc                                                                                                                                                                                              <  2:351436/245‑1 (MQ=255)                                                                 aCAGCTTCCTCAAACCCAAGAGCCATGAGATAGCGCAGCCAGTGCAATGCCCGTTCATCCCCCATGCTCAACGCGGTTAAAGACATCGGCGTACCGTTTTTTCTCTGCGCCAACAGCAGCGCGGCAGCGGCCTGTT‑CATCAGCACCGAGTAAACTTCGCAGTAGCGAAATATCCGGTCCGCTTTT‑ACCCTGAGCCGCCGCCTGCCACTCTTTTTCCGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACg                                                                                                                                                                                   <  1:569401/250‑1 (MQ=255)                                                                      ttCCTCAAACCCAAGAGCCATGAGATAGCGCAGCCAGTGCAATGCCCGTTCATCCCCCATGCTCAACGCGGTTAAAGACATCGGCGTACCGTTTTTTCTCTGCGCCAACAGCAGCGCGGCAGCGGCCTGTT‑CATCAGCACCGAGTAAACTTCGCAGTAGCGAAATATCCGGTCCGCTTTT‑ACCCTGAGCCGCCGCCTGCCACTCTTTTTCCGCCAGCATTCCGGCATCCGGTTCGGCCTTAAACGCTAAc                                                                                                                                                                              >  1:117114/1‑250 (MQ=255)                                                                            aaaCCCAAGAGCCATGAGATAGCGCAGCCAGTGCAATGCCCGTTCATCCCCCATGCTCAACGCGGTTAAAGACATCGGCGTACCGTTTTTTCTCTGCGCCAACAGCAGCGCGGCAGCGGCCTGTT‑CATCAGCACCGAGTAAACTTCGCAGTAGCGAAATATCCGGTCCGCTTTT‑ACCCTGAGCCGCCGCCTGCCACTCTTTTTCCGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGt                                                                                                                                                                       >  2:343550/1‑251 (MQ=255)                                                                            aaaCCCAAGAGCCATGAGATAGCGCAGCCAGTGCAATGCCCGTTCATCCCCCATGCTCAACGCGGTTAAAGACATCGGCGTACCGTTTTTTCTCTGCACCAACAGCAGCGCGGCAGCGGCCTGTT‑CATCAGCACCGAGTAAACTTCGCAGTAGCGAAATATCCGGTCCGCTTTT‑ACCCTGAGCCGCCGCCTGCCACTCTTTTTCCGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGt                                                                                                                                                                       >  2:134693/1‑251 (MQ=255)                                                                                ccAAGAGCCATGAGATAGCGCAGCCAGTGCAATGCCCGTTCATCCCCCATGCTCAACGCGGTTAAAGACATCGGCGTACCGTTTTTTCTCTGCGCCAACAGCAGCGCGGCAGCGGCCTGTT‑CATCAGCACCGAGTAAACTTCGCAGTAGCGAAATATCCGGTCCGCTTTT‑ACCCTGAGCCGCCGCCTGCCACTCTTTTTCCGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCAt                                                                                                                                                                   <  1:163796/251‑1 (MQ=255)                                                                                          tGAGATAGCGCAGCCAGTGCAATGCCCGTTCATCCCCCATGCTCAACGCGGTTAAAGACATCGGCGTACCGTTTTTTCTCTGCGCCAACAGCAGCGCGGCAGCGGCCTGTT‑CATCAGCACCGAGTAAACTTCGCAGTAGCGAAATATCCGGTCCGCTTTT‑ACCCTGAGCCGCCGCCTGCCACTCTTTTTCCGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCAcctcc                                                                                                                                                         >  1:213254/1‑251 (MQ=255)                                                                                          tGAGATAGCGCAGCCAGTGCAATGCCCGTTCATCCCCCATGCTCAACGCGGTTAAAGAAATCGGCGTACCGTTTTTTCTCTGCGCCAACAGCAGCGCGGCAGCGGCCTGTT‑CATCAGCACCGAGTAAACTTCGCAGTAGCGAAATATCCGGTCCGCTTTT‑ACCCTGAGCCGCCGCCTGCCACTCTTTTTCCGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCAcctcc                                                                                                                                                         >  2:349977/1‑251 (MQ=255)                                                                                           gagaTAGCGCAGCCAGTGCAATGCCCGTTCATCCCCCATGCTCAACGCGGTTAAAGACATCGGCGTACCGTTTTTTCTCTGCGCCAACAGCAGCGCGGCAGCGGCCTGTT‑CATCAGCACCGAGTAAACTTCGCAGTAGCGAAATATCCGGTCCGCTTTT‑ACCCTGAGCCGCCGCCTGCCACTCTTTTTCCGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCAcctcct                                                                                                                                                        >  2:116049/1‑251 (MQ=255)                                                                                             gaTAGGGCGGCCAGTACAATGCCCGTTCACCCCCCATGCTCAACGCGGTTAAAAACAGCGGCGTACCGTTTTTTCTCTGCGCCAACAGCAGCGCGGCAGCGGCCTGTC‑CACCAGCACCGAGTAAACTTCGCAGTAGCGAAATATCCGGTCCGCTTTT‑CCCCGGAGCCGCCGCCTGCCACTCTTTTTCCGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTAc                                                                                                                                                      <  2:18630/251‑1 (MQ=255)                                                                                                aGCGCCGCCAGTGCAATGCCCGTTCACCCCCCATGCTCACCGCGGTTAAAGACATCGGCGTACCGTTTTTTCTCTGCGCCAACAGCAGCGCGGCAGCGGCCTGGT‑CATCAGCCCCGATAAAACTTCGCAGTAGCGAAATATCCGGTCCGCTTTT‑CCCCTGCGCCGCCGCCTGCCACTCTTTTTCCGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACAt                                                                                                                                                    <  2:35951/250‑1 (MQ=255)                                                                                                  cgcAGCCAGTGCAATGCCCGTTCATCCCCCATGCTCAACGCGGTTAAAGACATCGGCGTACCGTTTTTTCTCTGCGCCAACAGCAGCGCGGCAGCGGCCTGTT‑CATCAGCACCGAGTAAACTTCGCAGTAGCGAAATATCCGGTCCGCTTTT‑ACCCTGAGCCGCCGCCTGCCACTCTTTTTCCGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTc                                                                                                                                                 >  2:501386/1‑251 (MQ=255)                                                                                                    cAGCCAGTGCAATGCCCGTTCATCCCCCATGCTCAACGCGGTTAAAGACATCGGCGTACCGTTTTTTCTCTGCGCCAACAGCAGCGCGGCAGCGGCCTGTT‑CATCAGCACCGAGTAAACTTCGCAGTAGCGAAATATCCGGTCCGCTTTT‑ACCCTGGGCCGCCGCCTGCCACTCTTTTTCCGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTccc                                                                                                                                               <  1:246055/251‑1 (MQ=255)                                                                                                            gCAATGCCCGTTCATCCCCCATGCTCAACGCGGTTAAAGACATCGCCGTACCTTTTTTTCTCTGCGCCAACAGCAGCGCGGCAGCGGCCTGTT‑CATCAGCACCGAGTAAACTTCGCAGTAGCGAAATATCCGGTCCGCTTTT‑ACCCTGAGCCGCCGCCTGCCACTCTTTTTCCGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGcca                                                                                                                                        <  1:351220/250‑1 (MQ=255)                                                                                                                tGCCCGTTCATCCCCCGTGCGCAACGGGGTTAAAGCCATCGGCGTCCCGTTTTTTCGCTGCCCAAACAGCAGCGCGGCAGCGGCCTGTT‑CACAAGCACCGAGTAAACTCCGCAGTAGCGAAATATCCGGTCCGCTTTT‑CCCCTGAGCCCCGGCCTGCCACTCTTGTTCCGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAg                                                                                                                                    <  1:571930/250‑1 (MQ=255)                                                                                                                          ccccccATGCCCAACGCGTTTAAAGACACCGGCGTACCGTTTTTTCTCTGCCCCAACAGCAGCGCGGCAGCGGCCTGTT‑CATCAGCACCGAGAAAACTTCGCAGTAGCGAAATACCCGGTCCGCTTTT‑ACCCGGAGCCGCCGCCTGCCACTCTTTTTCCGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCtt                                                                                                                         <  1:95364/250‑1 (MQ=255)                                                                                                                                aTGCTCAACGCGGTTAAAGACATCGGCGTACCGTTTTTTCTCTGCGCCAACAGCAGCGCGGCAGCGGCCTGTT‑CATCAGCACCGAGTAAACTTCGCAGTAGCGAAATATCCGGTCCGCTTTT‑ACCCTGAGCCGCCGCCTGCCACTCTTTTTCCGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGttt                                                                                                                   >  1:97510/1‑251 (MQ=255)                                                                                                                                 ggCTCACCGCGGTTAAAGACATCGGCGTACCGTTTTTTCTCTGCGCCAACAGCAGCGCGGCAGCGGCCTGTT‑CATCACCACCGAGTAAACTTCGCAGTAGCGAAATATCCGGTCCGCTTTT‑CCCCTGAGCCGCCGCCTGCCACTCTTTTTCCGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTc                                                                                                                  <  1:100888/250‑1 (MQ=255)                                                                                                                                         gcgGTTAAAGACATCGGCGTACCGTTTTTTCTCTGCGCCAACAGCAGCGCGGCAGCGGCCTGTT‑CATCAGCACCGAGTAAACTTCGCAGTAGCGAAATATCCGGTCCGCTTTT‑ACCCTGAGCCGCCGCCTGCCACTCTTTTTCCGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGccc                                                                                                          >  2:158908/1‑251 (MQ=255)                                                                                                                                            gTTAAAGACATCGGCGTACCGTTTTTTCTCTGCGCCAACAGCAGCGCGGCAGCGGCCTGTT‑CATCAGCACCGAGTAAACTTCGCAGTAGCGAAATATCCGGTCCGCTTTT‑ACCCTGAGCCGCCGCCTGCCACTCTTTTTCCGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCAt                                                                                                        >  1:579516/1‑250 (MQ=255)                                                                                                                                                          ggacctttttttccgcggcgcCAACAGCGGCGCGGCAGCGGC‑TGTTTCACCAGCCCCGAGTAAACTTCGCAGTAGCAAAATACCCGGCCCCCTTTT‑CCCCGGAGCCGCCGCCTGCCACTCTTTTTCCGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTTATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCacaaca                                                                                         <  2:297032/234‑1 (MQ=255)                                                                                                                                                                   ttttCTCTGCGCCAACAGCAGCGCGGCAGCGGCCTGTT‑CATCAGCACCGAGTAAACTTCGCAGTAGCGAAATATCCGGTCCGCTTTT‑ACCCTGAGCCGCCGCCTGCCACTCTTTTTCCGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAaaccaacc                                                                                >  1:110754/1‑251 (MQ=255)                                                                                                                                                                   ttttCTCTGCCCAAACAGCAGCGCGGCAGCGGCCTGTT‑CATCAGCCCCGAGTAACCTCCGCAGTAGCGAAATATCCGGTCCGCTTTT‑CCCCTGAGCCGCCGCCTGCCACTCTTTTTCCGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAaaccaac                                                                                 <  1:345660/250‑1 (MQ=255)                                                                                                                                                                         cggcgcCAACAGTAGCGCGGCAGGGGCTTGTT‑CATCAGCACCGAGTAATCTTCGCAGTAGCAAAATTCCCGGCCCGCTTTT‑CCCCTGGCCCGCCGCCTGCCACTCTTTTTCCGCCAGCATGCCGGCATCCGGTTCGGCTTTAAACGCGAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCaacaac                                                                          <  2:249410/249‑1 (MQ=255)                                                                                                                                                                               cAACAGCAGCGCGGCAGCGGCCTGTT‑CATCAGCACCGAGTAAACTTCGCAGTAGCGAAATATCCGGTCCGCTTTT‑ACCCTGAGCCGCCGCCTGCCACTCTTTTTCCGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTc                                                                    >  1:76819/1‑251 (MQ=255)                                                                                                                                                                                 aCAGCAGCGCGGCAGCGGCCTGTT‑CATCAGCACCGAGTAAACTTCGCAGTAGCGAAATATCCGGTCCGCTTTT‑ACCCTGAGCCGCCGACTGCCACTCTTTTTCCGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGc                                                                  >  1:26797/1‑251 (MQ=255)                                                                                                                                                                                     cagcGCGGCAGCGGCCTGTT‑CATCAGCACCGAGTAAACTTCGCAGTAGCGAAATATCCGGTCCGCTTTT‑ACCCTGAGCCGCCGCCTGCCACTCTTTTTCCGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAgg                                                              >  1:396580/1‑251 (MQ=255)                                                                                                                                                                                         gcgGCAGCGGCCTGTT‑CATCAGCACCGAGTAAACTTCGCAGTAGCGAAATATCCGGTCCGCTTTT‑ACCCTGAGCCGCCGCCTGCCACTCTTTTTCCGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCCGGAGc                                                           >  2:41234/1‑250 (MQ=255)                                                                                                                                                                                           ggCAGCGGCCTGTT‑CATCAGCCCCGGGTAAACTTCGCAGAAGCGAAATATCCGGTCCGCTTTT‑CCCCTGGGCCGCCGCCTGCCACTCTTTTTCCGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAAt                                                        <  1:175510/251‑1 (MQ=255)                                                                                                                                                                                            gCAGCGGCCTGTT‑CATCAGCACCGAGTAAACTTCGCAGTAGCGAAATATCCGGTCCGCTTTT‑ACCCTGAGCCGCCGCCTGCCACTCTTTTTCCGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATa                                                       <  1:415893/251‑1 (MQ=255)                                                                                                                                                                                                 ggCCTGTT‑CATCAGCACCGAGTAAACTTCGCAGTAGCGAAATATCCGGTCCGCTTTT‑ACCCTGAGCCGCCGCCTGCCACTCTTTTTCCGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTg                                                  >  1:55393/1‑251 (MQ=255)                                                                                                                                                                                                  gCCTGTT‑CATCAGCACCGCGTAAACTTCGCAGTAGCGAAATATCCGGTCCGCTTTT‑ACCCTGAGCCGCCGCCTGCCACTCTTTTTCCGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTgc                                                 >  2:94083/1‑250 (MQ=255)                                                                                                                                                                                                          caccaGCCCCGATTAACCTTCGCAGTGGCGAAATTCCCGGCCCGCTTTT‑CCCCTGCGCCCCCGCCTGCCCCTCTTTTTCCGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCCCCTCCTCCATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCAt                                          <  2:97510/248‑1 (MQ=255)                                                                                                                                                                                                           cacaGCACCGATAAAACTTCGCAGGAGCAAAATATCCG‑TCCCCTTTTTCCCCTGACCCGCCCCCTGCCACTCTTTTTCCGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCAtt                                         <  2:355944/249‑1 (MQ=255)                                                                                                                                                                                                            tcaGCCCCGAGAAACCTTCGCAGTAGCGAAATATCCGGTCCGTTTTT‑CCCCTGAGCCCCCGCCTGCCACTCTTTTTCCGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCAt                                          <  1:25432/249‑1 (MQ=255)                                                                                                                                                                                                                  ccaaGTAACCTCCGCAGTAGGAAAATACCCGGTCCGCTTTT‑CCCCGGACCCGCCGCCTGCCACTCTTTTTCCGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCATTCCCCACt                                  <  2:110754/248‑1 (MQ=255)                                                                                                                                                                                                                     aGTAAACTTCGCAGTAGCGAAATATCCGTTCCGCTTTT‑CCCCTGAGCCGCCGCCTGCCACTCTTTTTCCGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCATTCCCCACTgcg                               <  1:307768/251‑1 (MQ=255)                                                                                                                                                                                                                          aCTTCGCAGTAGCGAAATATCCGGTCCGCTTTT‑ACCCTGAGCCGCCGCCTGCCACTCTTTTTCCGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCATTCCCCACTGCGCCAGc                          >  2:133826/1‑251 (MQ=255)                                                                                                                                                                                                                             ccGCAGTAGCAAAATATCCGGTCCGCTTTT‑ACCCGGAGCCGCCGCCTGCCACTCTTTTTCCGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCATTCCCCACTGCGCCAGCAcc                       <  1:167834/250‑1 (MQ=255)                                                                                                                                                                                                                                  gTAGCAAAATACCCGGTCCGTTTTT‑CCCCTGAGCCGCCGCCTGCCACTCTTTTTCCGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCATTCCCCACTGCGCCAGCACCCACgg                  <  1:452562/251‑1 (MQ=255)                                                                                                                                                                                                                                   tAGCGAAATATCCGGTCCGCTTTT‑ACCCTGAGCCGCCGCCTGCCACTCTTTTTCCGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCATTCCCCACTGCGCCAGCACCCACGgc                 >  1:142551/1‑251 (MQ=255)                                                                                                                                                                                                                                              ccGGTCCGCTTTT‑ACCCTGAGCCGCCGCCTGCCACTCTTTTTCCGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCATTCCCCACTGCGCCAGCACCCACGGCAGCAGGCCAGt      >  1:187559/1‑251 (MQ=255)                                                                                                                                                                                                                                                ggTCCGCTTTT‑ACCCTGAGCCGCCGCCTGCCACTCTTTTTCCGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCATTCCCCACTGCGCCAGCACCCACGGCAGCAGGCCAGTAc    >  1:575242/1‑251 (MQ=255)                                                                                                                                                                                                                                                  tCCGCTTTT‑ACCCTGAGCCGCCGCCTGCCACTCTTTTTCCGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCATTCCCCACTGCGCCAGCACCCACGGCAGCAGGCCAGTACtt  >  1:319814/1‑251 (MQ=255)                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                  TTCGGCAATATGGCGAGCGACAAATTCCGGCGCAAAGCGCAGATCTGCCGCCGTCTCCAGCAATACAGCTTCCTCAAACCCAAGAGCCATGAGATAGCGCAGCCAGTGCAATGCCCGTTCATCCCCCATGCTCAACGCGGTTAAAGACATCGGCGTACCGTTTTTTCTCTGCGCCAACAGCAGCGCGGCAGCGGCCTGTT‑CATCAGCACCGAGTAAACTTCGCAGTAGCGAAATATCCGGTCCGCTTTT‑ACCCTGAGCCGCCGCCTGCCACTCTTTTTCCGCCAGCATTCCGGCATCCGGTTCGGCTTTAAACGCTAACAGAATGTTCATTCGCACCTCCTACATTTCCCGCCGCCACCAGCGGGAGGGCTTTAGTTTCTGGTGCCCATAGCCAGGTCACAACAAAACCAACCAACAACACTGTCGCCAGGAGCAATAAGGTGACTTGCATTCCCCACTGCGCCAGCACCCACGGCAGCAGGCCAGTACTT  >  NC_000913/2896226‑2896716

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