Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,415,866	0	A	T	21.9%	151.6 / 34.7	76	I374F (ATC→TTC)	pta	phosphate acetyltransferase
Reads supporting (aligned to +/- strand):  ref base A (34/23);  new base T (11/5);  total (47/29)
Fisher's exact test for biased strand distribution p-value = 5.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTGGATCGAATCTCTGACTGCCACTTCTGAGCGCAGCCGTCGTCTGTCTCCGCCTGCGTTCCGTTATCAGCTGACTGAACTTGCGCGCAAAGCGGGCAAACGTATCGTACTGCCGGAAGGTGACGAACCGCGTACCGTTAAAGCAGCCGCTATCTGTGCTGAACGTGGTATCGCAACTTGCGTACTGCTGGGTAATCCGGCAGAGATCA‑ACCGTGTTGCAGCGTCTCAGGGTGTAGAACTGGGTGC  >  NC_000913/2415624‑2416107                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                   cAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTGGATCGAAtct                                                                                                                                                                                                                                            <  1:381597/251‑1 (MQ=255)    ccATGAGCGGCGTAGAAATCGGTGCCCGGGTGCTGGCTGGCGGTTACGAAATGGGCGAGCGCATTTCTAATCTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGGGAACACCAACACCTGGCAGTCCTCTCTGAGGCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCGCGAACGTATGGAGAAAGTTCAGGAATGCGTTGCTAACTACATCAACGGTTGCTGGATCGAATCTCTg                                                                                                                                                                                                                                         <  2:215731/251‑1 (MQ=255)         aaCGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGCTACGGAATCGCCGCGCGCATTTCTTAAAGGTGCGAACGTGCTTTCGCTACCAGCCTGCCGGGATTTATGGAGAACAACAACACCCGGCCGGCCTTTCTTAGCCTGCAGAGCTTTAAACTGGAAAATCCGGGTGACGATCCCGGACGTATCGCGAAAGTTCAGGAATACGTTGCTAAATCCATCAACGCTGACTGGGTCGAATTTCTGGCTGc                                                                                                                                                                                                                                    >  1:295792/1‑251 (MQ=255)                tAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGGCTTTCTTAAATGTGCGAAACTGCTTTTGCTACCGGCCTGCCGGGATTTATTGTGAACAACAACACCTGGCCGACCTCTCTGAGGCTGCAGAGCCTCAACCTGGGAGGTCCCGTTGACGATCACGAACGTCACGCGAAAGTTCTGGAAAACGTTGCTAACTACATCAACGCTGACTGGATCGAATCTCTGACTGCCACctct                                                                                                                                                                                                                             >  1:414369/1‑251 (MQ=255)                 aGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTGCCTGTTGACGAGCACGAGCGTATCGAGAAAGTTCAGGAATACGTTGGTAACTACATCAACGCTGACTGGATCGAATCTCTGACTGCCACTTCt                                                                                                                                                                                                                             <  1:184917/250‑1 (MQ=255)                 aGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAAGGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGGGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTGGATCGAATCTCTGACTGCCACTTCt                                                                                                                                                                                                                             <  2:307684/250‑1 (MQ=255)                             ccTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTGGATCGAATCTCTGACTGCCACTTCTGAGCGCAGCcgtc                                                                                                                                                                                                                <  1:134020/251‑1 (MQ=255)                                      gACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGGATTTATGGGGGACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCCGGAATACGTTGCTAACTACATCAACGCTGACTGGATCGAATCTCTGACTTCCACTTCTGAGCGCAGCCGTCGTCTGtctc                                                                                                                                                                                                       >  1:87816/1‑251 (MQ=255)                                        cTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGCAATACGTTGCTAACTACATCAACGCTGACTGGTTCGAATCTCTGACTGCCACTTCTGAGCGCAGCCGTCGTCTGTCTCCg                                                                                                                                                                                                     >  1:389466/1‑251 (MQ=255)                                              gTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACAGCCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTGGATCGAATCTCTGACTGCCACTTCTGAGCGCAGCCGTCGTCTGTCTCCGCCTGCg                                                                                                                                                                                               >  2:369320/1‑251 (MQ=255)                                                tccgagaTGGGCGCCCGCATTCTTAAACGGTGAGAACGTGCTTGCGCTACCGGCCGGCCGGTATTTATGGGGAACACCAACCCCGGGCGGACCTCTCGGAGCCTGCAGATCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTGGATCGAATCTCTGACTGCCACTTCTGAGCGCAGCCGTCGTCTGTCTCCGCCTGCGt                                                                                                                                                                                              <  2:312188/244‑1 (MQ=255)                                                    aaaTGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCACCCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTGGTTCGAATCTCTGACTGCCACTTCTGAGCGCAGCCGTCGTCTGTCTCCGCCTGCGTTCCGt                                                                                                                                                                                         >  2:292880/1‑251 (MQ=255)                                                    aaaTGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTAGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTGGATCGAATCTCTGACTGCCACTTCTGAGCGCAGCCGTCGTCTGTCTCCGCCTGCGTTCCGt                                                                                                                                                                                         <  1:269823/251‑1 (MQ=255)                                                     aaTGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTGGATCGAATCTCTGACTGCCACTTCTGAGCGCAGCCGTCGTCTGTCTCCGCCTGCGTTCCGtt                                                                                                                                                                                        >  1:130749/1‑251 (MQ=255)                                                                gcATTTCTAAACTGTGCCAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGTCAGACCTCTCTGAGCCTGCAGAGCTTCAACCGGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTGGATCGAATCTCTGACTGCCACTTCTGAGCGCAGCCGGCGTCTGTCTCCGCCTGCGTTCCGTTATCAGctgact                                                                                                                                                                             <  1:158361/251‑1 (MQ=255)                                                                        aaaCTGTGCGACCGTGCTTTCGCTCCCGCCCTGCCGTTATTTATGGTGACCACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGACCGTATCGAGAAAGTTCAGGAATCCTTTGCTAACTACATCACCGCTGCCTGGATCGATTCTCTGACTGCCACTTCTGAGCGCAGCCGTCGTCTGTCTCCGCCTGCGTTCCGTTATCAGCTGACTGAACt                                                                                                                                                                        <  2:105010/248‑1 (MQ=255)                                                                        aaaCTGTGAGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGCGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCCGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGCATCGAGAAAGTTCAGGAATACGTTGCTAACTACATTCACGCTGACTGGATCGAATATCTGCATGCCTCTTCTGAGCGCAGCCGTCGTCTGTCTCCCCCTGCGTTCCGCTATCAGCTGACTGAACttcc                                                                                                                                                                     >  2:231253/1‑249 (MQ=255)                                                                         caCTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGCGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTGGATCGAATATCTGACTGCCAATTCTGAGCGCAGCCGTCGTCTGTCTCCGCCTGCGTTCCGTTGTCAGCTGACTGAACTTggg                                                                                                                                                                    >  2:251643/2‑249 (MQ=255)                                                                         aaCTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTGGTTCGAATCTCTGACTGCCACTTCTGAGCGCAGCCGTCGTCTGTCTCCGCCTGCGTTCCGTTATCAGCTGACTGAACTTgcg                                                                                                                                                                    >  2:105947/1‑251 (MQ=255)                                                                            gtggCGAACGTGCTTTCCCTCCCGGCCTGCCGGTATTTATGGTGAACCCCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCGGGAATCCGTTGCTAACTACATCAACGCTGACTGGTTCGAATCTCTGACTGCCACTTCTTAGCGCAGCCGTCGTCTGTCTCCGCCTGCGTTCCGTTATCAGCTGACTGAACTTgcgc                                                                                                                                                                   <  2:262380/246‑1 (MQ=255)                                                                               gCGAACGTGCTTTCGCTCCCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTGGATCGAATCTCTGACTGCCACTTCTGAGCGCAGCCGTCGTCTGTCTCCGCCTGCGTTCCGTTATCAGCTGACTGAACTTGCGCGCaaa                                                                                                                                                              <  2:87816/251‑1 (MQ=255)                                                                                cGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGCGAACACCAACACCTGGCAAACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTGGTTCGAATCTCTGACTGCCACTTCTGAGCGCAGCCGTCGTCTGTCTCCGACTGCGTCCCGTTATCACCTGACTGAACTTgcgcgcca                                                                                                                                                               >  2:393333/1‑247 (MQ=255)                                                                                 gAACGTGCTTTCGGTACCGGCATGCCGGAATATATGGCGCCAACCAAAACCTGACACCCCTCTACGAGCCTGCCGAGCTTCAACATGGACGTTCCGGTTGACGATCACGAACGTAGCGAGTCAGCTCCGGAATACGTTGCTAACTACATCCACGCTGACTGGATCGAATCTCTGACTGCCACTTCTGCGCGCAGCCGTCCTCTGTCTCCGCCTGCGTTCCGTCATCAGCGGACTGCACTTGCGCCCaacgc                                                                                                                                                            >  2:31735/1‑248 (MQ=255)                                                                                 gAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTGGATCGAATCTCTGACTGCCACTTCTGAGCGCAGCCGTCGTCTGTCTCCGCCTGCGTTCCGTTATCAGCTGACTGAACTTGCGCGCAAAGc                                                                                                                                                            >  1:271311/1‑251 (MQ=255)                                                                                  aaCGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTGGTTCGAATCTCTGACTGCCACTTCTGAGCGCAGCCGTCGTCTGTCTCCGCCTGCGTTCCGTTATCAGCTGACTGAACTTGCGCGCAAAGCg                                                                                                                                                           >  2:181266/1‑251 (MQ=255)                                                                                     gTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTGGATCGAATCTCTGACTGCCACTTCTGAGCGCAGCCGTCGTCTGTCTCCGCCTGCGTTCCGTTATCAGCTGACTGAACTTGCGCGCAAAGCGGGc                                                                                                                                                        >  1:195841/1‑251 (MQ=255)                                                                                                  cGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTGGATCGAATCTCTGACTGCCACTTCTGAGCGCAGCCGTCGTCTGTCTCCGCCTGCGTTCCGTTATCAGCTGACTGAACTTGCGCGCAAAGCGGGCAAACGTATCGTAc                                                                                                                                           >  1:19487/1‑251 (MQ=255)                                                                                                      cTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTACCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATAAGTTGCTAACTACATCAACGATGACTGGATCGAATCTCTGACTGCCACTTCTGAGCGCAGCCGTCGTCTGTCTCCGCCTGCGTTCCGTTATCAGCTGACTGAACTTGCGCGCAAAGCGGGCAAACGAATCGTCCTGcc                                                                                                                                       >  2:283227/1‑251 (MQ=255)                                                                                                        gCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTGGATCGAATCTCTGACTGCCACTTCTGAGAGCAGCCGTCGTCTGTCTCCGCCTGCGTTCCGTTATCAGCTGACTGAACTTGCGCGCAAAGCGGGCAAACGTATCGTACTGCCgg                                                                                                                                     >  1:311871/1‑251 (MQ=255)                                                                                                        gCCGGTATTTATGGTGAACACCAACACCGGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTGGATCGAATCTCTGACTGCCACTTCTGGGCGCAGCCGTCGTCTGTCTCCGCCTGCGTTCCGTTATCAGCTTACTGAACTTGCGCGCAAAGCGGGCAAACGTATCGTACTGCCgg                                                                                                                                     <  1:275778/251‑1 (MQ=255)                                                                                                           ggTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGCAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGCCTGGATCGAATCTCTGACTGCCACTTCTGAGCGCAGCCGTCGTCTGTCTCCGCCTGCGTTCCGTTATCAGCTGACTGAACTTGCGCGCACAGCGGGCAAACGTATCGTACTGCCGGa                                                                                                                                    >  2:265312/1‑249 (MQ=255)                                                                                                              aTTTATGGTGAACCCCAACACCGGGCAGCCCTCTCTGGGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTGGAGCGAATCTCTGACTGCCACTTCTGAGCGCAGCCGTCGTCTGTCTCCGCCTGCGTTCCGTTATCAGCTGACTGAACTTGCGCGCAAAGCGGGCAAACGTATCGTACTGCCGGAAGGt                                                                                                                                <  2:414369/250‑1 (MQ=255)                                                                                                               tttATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTGGATCGAATCTCTGACTGCCACTTCTGAGCGCAGCCGTCGTCTGTCTCCGCCTGCGTTCCGTTATCAGCTGACTGAACTTGCGCGCAAAGCGGGCAAACGTATCGTACTGCCGGAAGGTg                                                                                                                               >  1:100631/1‑250 (MQ=255)                                                                                                                 tATGGTGCACACCAACCACTGGCAGACCGCTCTGAGCCTGCAGGGCTTCAACCTGGAAGTTCCGGTTGCCGATCCCGAACGTATCGGGAAATTTCAGGAATACGTTGCTAGCTACATCAACGCTGACTGTACCGAGTCTCTGACTGCCGCTTCTGAGCGCGGCGGTCGGCTGTCTCCGCCTGCGTTCCGTTATCAGCTGACGGAAGTTTCGCGCAATGCGGGCTAACGTATCGTACTGcgtgagggtgac                                                                                                                             <  2:194949/250‑12 (MQ=255)                                                                                                                  aTGGTGAACACCAACACCTGGCAGAAATCTCTGAGCCTGCAGAGCTTCACCCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAACGTTCAGGAATACGTTGCTAACTCCATCCACGCTGACTGGATCGAATCTCTGACTGCCACTTCTGAGCGCAGCCGTCGTCTGTCTCCGCCTGCGTTCCGTTATCAGCTGACTGAACTTGCGCGCAATGCGGGCAAACGTATCGTACTGCCGGAAGGAGACgc                                                                                                                           >  2:131996/1‑250 (MQ=255)                                                                                                                          caccaacaccTGGCAGACATCTCTGATCCTGCAGAGCTTCAACCGGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAGAGGTCAGGAATACGTTGCTAACTACCTCAACGCTGACTGGATCGAATCTCTGACTGACATTTCTGAGCGCAGCCGTCGTCTGTCGCCGCCTGCGATCCGTTATCAGCTGACTGAACTTGCGCGCCACTAGGAAAAGCGTCTCGTACTACCGGACGGTGACGAAcctggt                                                                                                                    >  2:512/1‑246 (MQ=255)                                                                                                                                ccccTGGCAGACCTCTCGAAGCCGGCAGGGCTTCAACCTGGAAGTTCCGGTGGCCGTTCACGAACGTTTCGAGAAAGTTAAGGAATACGTAGCTAACTACTTCAACGCTGACTGGGTCGAATCTCTGACTGCCACTTCCGAGCGCAGCCGTCGTCTGTCTCCGCCTGCGTTCCGTTATCATCTGTCTGAACTTGCGCGCAAAGCGGGCAAACGTATCGTACTGCCGTAAGGTGACGAACCGCGTACCGTTa                                                                                                             <  2:272267/249‑1 (MQ=255)                                                                                                                                     ggctgaactgtctGGGCCTGCAGTGCTTCAACCGGGAAGTTCCGGTTGACGATCACGAGCGTATCGAGAAAGTTCAGAAATACGTGGCTAACTACATCATCGCTGACTGTATCGAATCTCTGACTGCCACTTCTGAGCGCAGCCGGCGTCTGTCTCCGCCTGCGTTCCGTTATCAGCTGACTGAACGTGCGCGCAAAGCGGAGAAACGTATCGTACTTCCGGAAGGTGACTAACCGCGTACCGTTAAagca                                                                                                        <  2:422541/241‑1 (MQ=255)                                                                                                                                        aGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTGGTTCGAATCTCTGACTGCCACTTCTGAGCGCAGCCGTCGTCTGTCTCCGCCTGCGTTCCGTTATCAGCTGACTGAACTTGCGCGCAAAGCGGGCAAACGTATCGTACTGCCGGAAGGTGACGAACCGCGTACCGTTAAAGCAGcc                                                                                                     >  1:15497/1‑251 (MQ=255)                                                                                                                                        aGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTGGATCGAATCTCTGACTGCCACTTCTGAGCGCAGCCGTCGTCTGTCTCCGCCTGCGTTCCGTTATCAGCTGACTGAACTTGCGCGCAAAGCGGGCAAACGTATCGTACTGCCGGAAGGTGACGAACCGCGTACCGTTAAAGCAGcc                                                                                                     <  2:100631/251‑1 (MQ=255)                                                                                                                                        aGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTGGATCGAATCTCTGACTGCCACTTCTGAGCGCAGCCGTCGTCTGTCTCCGCCTGCGTTCCGTTATCAGCTGACTGAACTTGCGCGCAAAGCGGGCAAACGTATCGTACTGCCGGAAGGTGACGAACCGCGTACCGTTAAAGCAGcc                                                                                                     >  1:5353/1‑251 (MQ=255)                                                                                                                                                        gCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTGGATCGAATCTCTGACTGCCACTTCTGAGCGCAGCCGTCGTCTGTCTCCGCCTGCGTTCCGTTATCAGCTGACTGAACTTGCGCGCAAAGCGGGCAAACGTATCGTACTGCCGGAAGGTGACGAACCGCGTACCGTTAAAGCAGCCGCTATCTGTGCTGAAc                                                                                     <  1:99135/251‑1 (MQ=255)                                                                                                                                                        gCAGAGCTTCAACCTGGAAGTTCCGGTGGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTGGTTCGAATCTCTGACTGCCACTTCTGAGCGCAGCCGTCGTCTGTCTCCGCCTGCGTTCCGTGATCAGCTGACTGAACTTGCGCGCAAAGCGTGCAAACGTATCGTACTGCGGGAAGGTGACGAACCGCGTACCGTTAAAGCAGCCGCTATCTGTGCTGAAc                                                                                     <  1:177488/251‑1 (MQ=255)                                                                                                                                                        cgagggctccaacgcggacgttgcGGTTGACGACCACGAACGAATCGAGAAAGTTAGGGAAGACGTTGCTAACTAAAACAGCGCTGACTTGTCCGAATCTCTGCGTTCCACTTCGGAGCGCAGCCGTCGTCTGTCTCCGCCTGCGTTCCGTTATCAGCTGACTGAACTTGCGCGCAAAGCGGGCAAACGTATCGTACTGCCGGAAGATGACGAACCCCATACCGTTAAAGTATCCGCTTTCCGTGCTGaa                                                                                      <  2:430109/227‑1 (MQ=255)                                                                                                                                                                tcaccctgtaacttacGGTTGACGATCACGAACGATTCAGGAAAGTGCAGGAATACGTTGCTAACTACATCACCGTGGACGGGAGCGAACCTCTGACTGCCACTCCTGAGCGCAGCCGTCGTCTGTCTCCGCCTGCGTTCCGTGATCAGCTGACGGAACTTGCGCGCAAAGCGGGCAAACGTATCGTACTGCCGGAAGGTGACGAACCGCGTACCGTTAAAGCAGCCGCTATCTGTGCTGAACGTGGTAt                                                                              <  2:119878/235‑1 (MQ=255)                                                                                                                                                                tCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTGGATCGAATCTCTGACTGCCACTTCTGAGCGCAGCCGTCGTCTGTCTCCGCCTGCGTTCCGTTATCAGCTGACTGAACTTGCGCGCAAAGCGGGCAAACGTATCGTACTGCCGGAAGGTGACGAACCGCGTACCGTTAAAGCAGCCGCTATCTGTGCTGAACGTGGTa                                                                               <  1:184235/249‑1 (MQ=255)                                                                                                                                                                  aaCCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTGGATCGAATCTCTGACTGCCACTTCTGAGCGCAGCCGTCGTCTGTCTCCGCCTGCGTTCCGTTATCAGCTGACTGAACTTGCGCGCAAAGCGGGCAAACGTATCGTACTGCCGGAAGGTGACGAACCGCGTACCGTTAAAGCAGCCGCTATCTGTGCTGAACGTGGTATCGc                                                                           >  2:95616/1‑251 (MQ=255)                                                                                                                                                                   aCCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGGTCAGGAATACGTTGCTAACTACATCAACGCTGACTGGTTCGAATCTCTGACTGCCACTTCTGAGCGCAGCCGTCGTCTGTCTCCGCCTGCGTTCCGTTATCAGCTGACTGAACTTGCGCGCAAAGCGCGCAAACGTATCGTACTGCCGGAAGGTGACGAACTGCGTACAGTTAAAGCAGCCGCTATCTGTGCTGAACGTgtgatcgc                                                                           >  1:345866/1‑243 (MQ=255)                                                                                                                                                                     cTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAACGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTGGATCGAATCTCTGACTGCCACTTCTGCGCGCCGCCGTCGTCTGTCTCCGCCTGCGTTCCGTTATCAGCTGACTGAACTTGCGCGCACAGCGGGCAAACGTATCGTACTGCCGGAAGGTGACGAACCGCGTACCGTTAAAGCAGCCGCTATCTGTGCTGAACGTGGCATCGCAAc                                                                        >  2:277149/1‑251 (MQ=255)                                                                                                                                                                         aaGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATCCGTTGCTCACTACATCAACGCTGACTGGATCGACTCTCTGACTGCCACTTCTGAGCGCAGCCGTCGTCTGTCTCCGCCTGCGTTCCGTTATCAGCTGACTGAACTTGCGCGCAAAGCGGGCAAACGTATCGTACTGCCGGAAGGTGACGAACCGCGTACCGTTAAAGCAGCCGCTATCTGTGCTGCACGTGGTATCGCAACTTGc                                                                    >  1:289935/1‑251 (MQ=255)                                                                                                                                                                         aaGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCGGACTGGATCGAATCTCTGACTGCCACTTCTGAGCGCAGCCGTCGTCTGTCTCCGCCTGCGTTCCGTTATCAGCTGACTGAACTTGCGCGCAAAGCGGGCAAACGTATCGTACTGCCGGAAGGTGACGAACCGCGTACCGTTAAAGCAGCCGCTATCTGTGCTGAACGTGGTATCGCAACTTGc                                                                    <  1:265312/251‑1 (MQ=255)                                                                                                                                                                          aGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTGGTTCGAATCTCTGACTGCCACTTCTGAGCGCAGCCGTCGTCTGTCTCCGCCTGCGTTCCGTTATCAGCTGACTGAACTTGCGCGCAAAGCGGGCAAACGTATCGTACTGCCGGAAGGTGACGAACCGCGTACCGTTAAAGCAGCCGCTATCTGTGCTGAACGTGGTATCGCAACTTGCg                                                                   <  1:292655/251‑1 (MQ=255)                                                                                                                                                                                  ttGACGCTCACGAACGTATCGAGAAAGTTAAGGAATACGATGCTAACTACATCCACGCTGACTGGCTCGCATCTCTGAATGCCCCTTCTGCGCGCCGCCGTCGTCTGTCTCCGCCTGCGTTCCGTTATCAGCTGACTGAACTTGCGCGCAAAGCGGGCAAACGTATCGTACTGCCGGAAGGTGCCGAACCGCGTCCCGTTAAAGCAGCCGCTATCTGTGCTGAACGTGGTATCGCAACTTGCGTActgctg                                                           >  1:62041/1‑251 (MQ=255)                                                                                                                                                                                   tGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTGGATCGAATCTCTGACTGCCACTTCTGAGCGCAGCCGTCGTCTGTCTCCGCCTGCGTTCCGTTATCAGCTGACTGAACTTGCGCGCAAAGCGGGCAAACGTATCGTACTGCCGGAAGGTGACGAACCGCGTACCGTTAAAGCAGCCGCTATCTGTGCTGAACGTGGTATCGCAACTTGCGTACTGCTgg                                                          >  1:315901/1‑251 (MQ=255)                                                                                                                                                                                    gACGATCACGAACGTCTCGAGAACGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTGGATCGAATCTCTGACTGCCACTTCTGAGCGCAGCCGTCGTCTGTCTCCGCCTGCGTTCCGTTATCAGCTGACTGAACTTGCGCGCAAAGCGGGCAAACGTATCGTACTGCCGGAAGGTGACGAACCGCGTACCGTTAAAGCAGCCGCTATCTGTGCTGAACGTGGTATCGCAACTTGCGTACTGCTggg                                                         >  2:273176/1‑251 (MQ=255)                                                                                                                                                                                           aCGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTGGATCGAATCTCTGACTGCCACTTCTGAGCGCAGCCGTCGTCTGTCTCCGCCTGCGTTCCGTTATCAGCTGACTGAACTTGCGCGCAAAGCGGGCAAACGTATCGTACTGCCGGAAGTTGACGAACCGCGTACCGTTAAAGCAGCCGCTATCTGTGCTGAACGTGGTATCGCAACTTGCGTACTGCTGGGTAATCCg                                                  >  1:231032/1‑251 (MQ=255)                                                                                                                                                                                           aCGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTGGATCGAATCTCTGACTGCCACTTCTGAGCGAAGCCGTCGTCTGTCTCCGCCTGCGTTCCGTTATCAGCTGACTGAACTTGCGCGCAAAGCGGGCAAACGTATCGTACTGCCGGAAGGTGACGAACCGCGTACCGTTAAAGCAGCCGCTATCTGTGCTGAACGTGGTATCGAAACTTGCGTACTGCTGGGTAAtacg                                                  >  1:136581/1‑248 (MQ=255)                                                                                                                                                                                            cgcaCGTATCGAGAAAGTTCAGGACTACGTTGCTAACTACATCAACGCTGACTGGATCGAATCTCTGACTGCCACTTCTGAGCGCAGCCGTCGTCTGTCTCCGCCTGCGTTCCGTTATCAGCTGACTGAACTTGCGCGCAAAGCGGGCAAACGTATCGTACTGCCGGAAGGTGACGACCCGCGTACCGTTCAAGCAGCCGCTATCTGTGCTGAACGTGGTATCGCAACTTGCGTACTGCTGGGTAATCCgg                                                 >  2:343474/4‑251 (MQ=255)                                                                                                                                                                                             gAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTGGATCGAATCTCTGACTGCCACTTCTGAGCGCAGCCGTCGTCTGTCTCCGCCTGCGTTCCGTTATCAGCTGACTGAACTTGCGCGCAAAGCGGGCAAACGTATCGTACTGCCGGAAGGTGACGAACCGCGTACCGTTAAAGCAGCCGCTATCTGTGCTGAACGTGGTATCGCAACTTGCGTACTGCTGGGTAATCCGGc                                                >  1:345994/1‑251 (MQ=255)                                                                                                                                                                                                   aTCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTGGTTCGAATCTCTGACTGCCACTTCTGAGCGCAGCCGTCGTCTGTCTCCGCCTGCGTTCCGTTATCAGCTGACTGAACTTGCGCGCAAAGCGGGCAAACGTATCGTACTGCCGGAAGGTGACGAACCGCGTACCGTTAAAGCAGCCGCTATCTGGGCTGAACGTGGTATCGCAACTTGCGTACTGCTGGGTAATCCGGCag                                              >  2:431247/1‑247 (MQ=255)                                                                                                                                                                                                   aTCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTGGATCGAATCTCTGACTGCCACTTCTGAGCGCAGCCGTCGTCTGTCTCCGCCTGCGTTCCGTTATCAGCTGACTGAACTTGCGCGCAAAGCGGGCAAACGTATCGTACTGCCGGAAGGTGACGAACCGCGTACCGTTAAAGCAGCCGCTATCTGTGCTGAACGTGGTATCGCAACTTGCGTACTGCTGGGTAATCCGGCAGAGAt                                          <  1:357660/251‑1 (MQ=255)                                                                                                                                                                                                       agaAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTGGATCGAATCTCTGACTGCCACTTCTGAGCGCAGCCGTCGTCTGTCTCCGCCTGCGTTCCGTTATCAGCTGACTGAACTTGCGCGCAAAGCGGGCAAACGTATCGTACTGCCGGAAGGTGACGAACCGCGTACCGTTAAAGCAGCCGCTATCTGTGCTGAACGTGGTATCGCAACTTGCGTACTGCTGGGTAATCCGGCAGAGATCA‑Ac                                     <  1:95616/251‑1 (MQ=255)                                                                                                                                                                                                          aaGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTGGATCGAATCTCTGACTGCCACTTCTGAGCGCAGCCGTCGTCTGTCTCCGCCTGCGTTCCGTTATCAGCTGACTGAACTTGCGCGCAAAGCGGGCAAACGTATCGTACTGCCGGAAGGTGACGAACCGCGTACCGTTAAAGCAGCCGCTATCTGTGCTGAACGTGGTATCGCAACTTGCGTACTGCTGGGTAATCCGGCAGAGATCA‑ACCgt                                  >  2:293584/1‑251 (MQ=255)                                                                                                                                                                                                            gTTCAGGAATACGTTGCTAACTACATCAACGCTGACTGGATCGAATCTCTGACTGCCACTTCTGAGCGCAGCCGTCGTCTGTCTCCGCCTGCGTTCCGTTATCAGCTGACTGAACTTGCGCGCAAAGCGGGCAAACGTATCGTACTGCCGGAAGGTGACGAACCGCGTACCGTTAAAGCAGCCGCTATCTGTGCTGAACGTGGTATCGCAACTTGCGTACTGCTGGGTAATCCGGCAGAGATCA‑ACCgtgt                                >  2:108871/1‑251 (MQ=255)                                                                                                                                                                                                             ttCATGAATACGTTTCTAACTACATCAACGCTAACTGGTGCGAATCTCTGACTGCCACTGCTTATCGCAGCCGTCGTCTGTCTCCGCCTGCGTTCCGTTATCAGCTGACTGAACTGGCACGCAAAGCGGGCAAACGTATCGTACTGCCGGAAGGTGACGa                                                                                                                           <  1:419021/160‑1 (MQ=255)                                                                                                                                                                                                             ttCAGGACTACGTTGCTAACAACATCAACGCTGTCTGGTTGGAATATCTGACTGCCACTTCTGAGCGCAGCCGTCGTCTGTCTCCGCCTGCGTTCAGTTCTCAGCTGACTGAACTTGCGCGCAAAGCGGGCAAACGTATCGTAATGCCGGAAGGTGACGa                                                                                                                           >  2:419021/1‑160 (MQ=255)                                                                                                                                                                                                                         ttGCTAACTACATCAACGCTGACTGGATCGAATCTCTGACTGCCACTTCTGAGCGCAGCCGTCGTCTGTCTCCGCCTGCGTTCCGTTATCAGCTGACTGAACTTGCGCGCAAAGCGGGCAAACGTATCGTACTGCCGGAAGGTGACGAACCGCGTAcc                                                                                                                 <  1:154446/158‑1 (MQ=255)                                                                                                                                                                                                                         ttGCTAACTACATCAACGCTGACTGGATCGAATCTCTGACTGCCACTTCTGAGCGCAGCCGTCGTCTGTCTCCGCCTGCGTTCCGTTATCAGCTGACTGAACTTGCGCGCAAAGCGGGCAAACGTATCGTACTGCCGGAAGGTGACGAACCGCGTACCGTTAAAGCAGCCGCTATCTGTGCTGAACGTGGTATCGCAACTTGCGTACTGCTGGGTAATCCGGCAGAGATCA‑ACCGTGTTGCAGCGtctc                     >  2:44415/1‑249 (MQ=255)                                                                                                                                                                                                                         ttGCTAAATACCTCAAAGCTGACTAGATCGAATCTCTGACTGCCACTTATGAGCGCAGCCGTCGTCTGTCTCCGCCTGCGTTCCGTTATCAGCTGACTGAACTTGCGCGCACAGCGGGCAAACGTATCGTACTGCCGGAAGGTGACGAACCGCGTAcc                                                                                                                 >  2:154446/1‑158 (MQ=255)                                                                                                                                                                                                                           gCTAACTACATCAACGCTGACTGGATCGAATCTCTGACTGCCACTTCTGAGCGCAGCCGTCGTCTGTCTCCGCCTGCGTTCCGTTATCAGCTGACTGAACTTGCGCGCAAAGCGGGCAAACGTATCGTACTGCCGGAAGGTGACGAACCGCGTACCGTTAAAGCAGCCGCTATCTGTGCTGAACGTGGTATCGCAACTTGCGTACTGCTGGGTAATCCGGCAGAGATCA‑ACCGTGTTGCAGCGTCTCCggg                 >  2:109043/1‑251 (MQ=255)                                                                                                                                                                                                                              aaCTACATCAACGCTGACTGGATCGAATCTCTGACTGCCACTTCTGAGCGCAGCCGTCGTCTGTCTCCGCCTGCGTTCCGTTATCAGCTGACTGAACTTGCGCGCAAAGCGGGCAAACGTATCGTACTGCCGGAAGGTGACGAACCGCGTACCGTTAAAGCAGCCGCTATCTGTGCTGAACGTGGTATCGCAACTTGCGTACTGCTGGGTAATCCGGCAGAGATCA‑ACCGTGTTGCAGCGTCTCAGGgtgt              <  1:9338/251‑1 (MQ=255)                                                                                                                                                                                                                              aaCTACATCAACGCTGACTGGATCGAATCTCTGACTGCCACTTCTGAACGCAGCCGTCGTCTGTCTCCGCCTGCGTTCCGTTATCAGCTGACTGAACTTGCGCGCAACGCGAGCAAACGTATCGTACTGCCGGAAGGTGACGAACCGCGTACCGTTAAAGCAGCCGCTATCTGTGCTGAACGTGGTATCGCAACTTGCGTACTGCTGGGTAATCCGGCAGAGCTCA‑ACCGTGTTGCAGCGTCTCAGGgtg               >  1:280768/1‑250 (MQ=255)                                                                                                                                                                                                                                 tACATCAACGCTGACTGGATCGAATCTCTGACTGCCACTTCTGAGCGCAGCCGTCGTCTGTCTCCGCCTGCGTTCCGTTATCAGCTGACTGAACTTGCGCGCAAAGCGGGCAAACGTATCGTACTGCCGGAAGGTGACGAACCGCGTACCGTTAAAGCAGCCGCTATCTGTGCTGAACGTGGTATCGCAACTTGCGTACTGCTGGGTAATCCGGCAGAGATCA‑ACCGTGTTGCAGCGTCTCAGGgtgg              >  2:307544/1‑247 (MQ=255)                                                                                                                                                                                                                                        aCGCTGACTGGTTCGAATCTCTGACTGCCACTTCTGAGCGCAGCCGTCGTCTGTCTCCGCCTGCGTTCCGTTATAAGCTGACTGAACTTGCGCGCAAAGCGGGCAAACGTATCGTACTGCCGGAAGGTGACGAACCGCGTACCGTTAAAGCAGCCGCTATCTGTGCTGAACGTGGGATCGCAACTTGCGCACTGCTGGGTAATCCGGCAGAGATCA‑ACCGTGTTGCAGCGTCTCAGGGTGTAGAACTgggg    >  2:407556/1‑250 (MQ=255)                                                                                                                                                                                                                                          gcggACTGGATCGAATCTCTGACTGCCAATTGTGCGCGCAGCCGTCGTCTGTCTCCGCCTTCGTTCCGTTATCAGCTGACTGAACTTGCGCGCAATGCGGGCAAACGTATCGTACTGCCGGAAGTTGACGAACCGCGTACCGTTAAAGCAGCCGCTATCTGTGCTGAACGTGGTATCGCACCTTGCGGACTGCTGGGTAATTCGGTAGACA‑CATACCGTGGTGCAGCGTCCCTGTGTGTAGAAATGGGTgg  >  2:2612/4‑250 (MQ=255)                                                                                                                                                                                                                                            tGACTGGTTCGAATCTCTGACTGCCACTTCTGAGCGCAGCCGTCGTCTGTCTCCGCCTGCGTTCCGTTATCAGCTGACTGAACTTGCGCGCAAAGCGGGCAAACGTATCGTACTGCCGGAAGGTGACGAACCGCGTACCGTTAAAGCAGCCGCTATCTGTGCTGAACGTGGTATCGCAACTTGCGTACTGCTGGGTAATCCGGCAGAGATCA‑ACCGTGTTGCAGCGTCTCAGGGTGTAGAACTGGGTGc  >  1:218926/1‑249 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                   CAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTGGATCGAATCTCTGACTGCCACTTCTGAGCGCAGCCGTCGTCTGTCTCCGCCTGCGTTCCGTTATCAGCTGACTGAACTTGCGCGCAAAGCGGGCAAACGTATCGTACTGCCGGAAGGTGACGAACCGCGTACCGTTAAAGCAGCCGCTATCTGTGCTGAACGTGGTATCGCAACTTGCGTACTGCTGGGTAATCCGGCAGAGATCA‑ACCGTGTTGCAGCGTCTCAGGGTGTAGAACTGGGTGC  >  NC_000913/2415624‑2416107

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