Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,189,451	Δ1 bp	coding (4102/4224 nt)	rpoC →	RNA polymerase subunit beta'

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,189,451	0	G	.	97.8%	459.9 / 1.8	92	coding (4102/4224 nt)	rpoC	RNA polymerase subunit beta'
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base . (50/40);  total (51/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCTCCGACTTCCTGGAAGGCGAACAGGTTGAATACTCTCGCGTCAAGATCGCAAACCGCGAACTGGAAGCGAACGGCAAAGTGGGTGCAACTTACTCCCGCGATCTGCTGGGTATCACCAAAGCGTCTCTGGCAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAACTGACTGAACTAATTGAGTCAAACTCGGCAAGGATTCG  >  NC_000913/4189160‑4189750                                                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                                                            aGCTCCGACTTCCTGGAAGGCGAACAGGTTGAATACTCTCGCGTCAAGATCGCAAACCGCGAACTGGAAGCGAACGGCAAAGTGGGTGCAACTTACTCCCGCGATCTGCTGGGTATCACCAAAGCGTCTATGGCAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGc                                                                                                                                                                                                                                                                                                    >  1:122487/1‑300 (MQ=255)  gCTCCGACTTCCTGGAAGGCGAACAGGTTGAATACTCTCGCGTCAAGATCGCAAACCGCGAACTGGAAGCGAACGGCAAAGTGGGTGCAACTTACTCCCGCGATCTGCTGGGTATCACCAAAGCGTCTATGGCAACCGAGTCCTTCATCTCCGCGGCATCGTTCAAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCAGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCg                                                                                                                                                                                                                                                                                                   >  2:263145/1‑300 (MQ=255)  gCTCCGACTTCCTGGAAGGCGAACAGGTTGAATACTCTCGCGTCAAGATCGCAAACCGCGAACTGGAAGCGAACGGCAAAGTGGGTGCAACTTACTCCCGCGATCTGCTGGGTATCAACAAAGCGTCTATGGCAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCAACAG‑ATCGTATGCg                                                                                                                                                                                                                                                                                                   >  2:521611/1‑300 (MQ=255)        aCTTCCTGGAAGGCGAACAGTTTGAATACTCTCGCGTCAAGATCGCAAACCGCGAACTGGAAGCGAACGGCAAAGTGGGTGCAACTTACTCCCGCGATCTGCTGGGTATCACCAAAGCGTCTCTGGCAACCGAGTCCTTCATCTCCTCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCg                                                                                                                                                                                                                                                                                             <  2:111067/300‑1 (MQ=255)           tCCTGGAAGGCGAACAGGTTGAATACTCTCGCGTCAAGATCGCAAACCGCGAACTGGAAGCGAACGGCAAAGTGGGTGCAACTTACTCCCGCGATCTGCTGGGTATCACCAAAGCGTCTCTGGCAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGtgc                                                                                                                                                                                                                                                                                          <  1:23370/300‑1 (MQ=255)                 aaGGCGAACAGGTTGAATACTCTCGCGTCAAGATCGCAAACCGCGAACTGGAAGCGAACGGCAAAGTGGGTGCAACTTACTCCCGCGATCTGCTGGGTATCACCAAAGCGTCTCTGGCAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGCTCNCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTAGCCGTGCTGCggg                                                                                                                                                                                                                                                                                    >  2:689966/1‑300 (MQ=255)                      gAACAGGTTGAATACTCTCGCGTCAAGATCGCAAACCGCGAACTGGAAGCGAACGGCAAAGTGGGTGCAACTTACTCCCGCGATCTGCTGGGTATCACCAAAGCGTCTCTGGCAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGNCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACC‑NGATCGTATGCGTCGCCGTGCTGCGGGTGAAg                                                                                                                                                                                                                                                                               <  2:666098/300‑1 (MQ=255)                          aGGTTGAATACTCTCGCGTCAAGATCGCAAACCGCGAACTGGAAGCGAACGGCAAAGTGGGTGCAACTTACTCCCGCGATCTGCTGGGTATCACCAAAGCGTCTCTGGCAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTcc                                                                                                                                                                                                                                                                           >  1:596323/1‑300 (MQ=255)                          aGGTTGAATACTCTCGCGTCAAGATCGCAAACCGCGAACTGGAAGCGAACGGCAAAGTGGGTGCAACTTACTCCCGCGATCTGCTGGGTATCACCAAAGCGTCTCTGGCAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCAGTGCTGCGGGTGAAGCTcc                                                                                                                                                                                                                                                                           >  2:769623/1‑300 (MQ=255)                           gggggaggaCTCTCGCGTAAAGATCGCAAACCGCGAACTGGCAGCGAACGGCAAAGTGGGTGCCACCTACTCCCGCGATCTGCTGGGTATCACCAAAGCGTCTCTGGCAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCAATCGCGTACTGACCGCAGNAGACAATGCGGGCAAACGCGACGAACTGAACGGCCTGCACGAGCACGTGATCGTGGGTAGTCTGATCCAGGTAGGTACGCGTTACTCGTATCACCAG‑ACCGTATGTGTCGCCGTGCTGCGGGTGAAGCTCCg                                                                                                                                                                                                                                                                          >  2:221146/9‑300 (MQ=255)                             ttGAATACTCTCGCGTCAAGATCGCAAACCGCGAACTGGAAGCGAACGGCAAAGTGGGTGCAACTTACTCCCGCGATCTGCTGGGTATCACCAAAGCGTCTCTGGCAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGc                                                                                                                                                                                                                                                                        <  1:141/300‑1 (MQ=255)                              tGAATACTCTCGCGTCAAGATCGCAAACCGCGAACTGGAAGCGAACGGCAAAGTGGGTGCAACTTACTCCCGCGATCTGCTGGGTATCACCAAAGCGTCTCTGGCAACCGAGTCCTTCATCTCCGCTGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCt                                                                                                                                                                                                                                                                       <  1:251368/300‑1 (MQ=255)                                  tACTCTCGCGGCAAGAGCGCAAACCGCGAACTGGAAGCGAACGGCAAAGTGGGTGCAACTTACTCCCGCGATCTGCTGGGTATCACCAAAGCGTCTCTGGCAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGNGCTGACCGAAGCAGCCGTTGCGGGCANACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTTGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTacac                                                                                                                                                                                                                                                                   <  2:605213/300‑1 (MQ=255)                                          cgTCAAGATCGCAAACCGCGAACTGGAAGCGAACGGCAAAGTGGGTGCAACTTACTCCCGCGATCTGCTGGGTATCACCAAAGCGTCTCTGGCAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTg                                                                                                                                                                                                                                                           >  1:776838/1‑300 (MQ=255)                                                    gCAAACCGCGAACTGGAAGCGAACGGCAGAGTGGGTGTAACTTACTCCCGCGATCTGCTGGGTATCACCAAAGCGTCTCTGGCAACCGAGTCCTTCAGCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCagaag                                                                                                                                                                                                                                                 <  2:62962/300‑1 (MQ=255)                                                         ccGCGAACTGGAAGCGAACGGCAAAGTGGGTGCAACTTACTCCCGCGATCTGCTGGGTATCACCAAAGCGTCTCTGGCAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCa                                                                                                                                                                                                                                            <  1:562236/300‑1 (MQ=255)                                                                            ggCAAAGTGGGTGCAACTTACTCCCGCGATCTGCTGGGTATCACCAAAGCGTCTCTGGCAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAAc                                                                                                                                                                                                                         <  1:385002/300‑1 (MQ=255)                                                                            ggCAAAGTGGGTGCAACTTACTCCCGCGATCTGCTGGGTATCAACAAAGCGTCTCTGGCAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAAc                                                                                                                                                                                                                         >  1:447699/1‑300 (MQ=255)                                                                                    gggTGCAACTTACTCCCGCAATCTGCTGGGTATCACCAAAGCGTCTATGGCAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAATTGCGCGGCATGAAGGAGAACGTTATCGTGAGTCGTCTGATCCCGGCAGATACCGGTTAAGCGTACCACCAG‑ATCGTATGCGTCGCCATGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGGCTGCAGAAGAAGCATCTGCCAGCATGGCAGAActgctggac                                                                                                                                                                                                                 >  2:166859/1‑297 (MQ=255)                                                                                           aCTTACTCCCGCGATCTGCTGGGTATCACCAAAGCGTTTCTGGCAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTc                                                                                                                                                                                                          <  1:212286/300‑1 (MQ=255)                                                                                           aCTTACTCCCGCGATCTGCTGGGTATCACCAAAGCGTCTCTGGCAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGNGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCACAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTc                                                                                                                                                                                                          <  2:217733/300‑1 (MQ=255)                                                                                                cTCCCGCGATCTGCTGGTTATCACCAGAGCGTCTCTGGCAACCGAGTCCTCCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCATACTGACCGAAGCAGCCGTGGCGGGCAAACGCGACGAACGGCGCGGCCTGAAAGAGAACGTTATCGTGGTTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTAGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGc                                                                                                                                                                                                     <  2:5772/300‑1 (MQ=255)                                                                                                cTCCCGCGATCTGCTGGGTATCACCAAAGCGTCTCTGGCAACCGAGTCCTTCATCTCCGCTGCATCGTTCCAGTAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCTTCGCCGTGCTGCGGTTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGc                                                                                                                                                                                                     <  1:491693/300‑1 (MQ=255)                                                                                                cTCCCGCGATCTGCTGGGTATCACCAAAGCGTCTCTGGCAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGc                                                                                                                                                                                                     <  1:187042/300‑1 (MQ=255)                                                                                                          ctgctgGGTATCACCAAAGCGTCTCTGGCAACCGAGTCCTTCATCTCCGCTGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATa                                                                                                                                                                                           <  1:795795/300‑1 (MQ=255)                                                                                                          ctgctgGGTATCACCAAAGCGTCTCTGGCAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATa                                                                                                                                                                                           <  1:729700/300‑1 (MQ=255)                                                                                                          ctgctgGGTATCAACAAAGCGTCTCTGGCAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTAGTCTGATCCAGGCCGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGAAGGTCTGGGCGGTTCTgcta                                                                                                                                                                                           >  1:245995/1‑297 (MQ=255)                                                                                                            gctgGGTATCACCAAAGCGTCTCTGGCAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAAc                                                                                                                                                                                         <  1:252338/300‑1 (MQ=255)                                                                                                            gctgGGTATCACCAAAGCGTCTATGGCAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAAc                                                                                                                                                                                         >  1:73572/1‑300 (MQ=255)                                                                                                             ctgGGTATCACCAAAGCGTCTCTGGCAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCNGATCCCGGCAGGTACCGGTTACGCGTNCCACCAG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACg                                                                                                                                                                                        <  2:629389/300‑1 (MQ=255)                                                                                                               gggTATCACCAAAGCGTCTCTGGCAACCGAGTCCTTCATCTCCGCGGCATCGTTCCNGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCNTCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAg                                                                                                                                                                                      >  2:848909/1‑300 (MQ=255)                                                                                                               gggTATCACCAAAGCGTCTATGGCAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAATGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGACCGCAGGTCTGGGCGGTTCTGATAACGAg                                                                                                                                                                                      >  2:261736/1‑300 (MQ=255)                                                                                                               gggTATCACCAAAGCGTCTATGGCAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCNTCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGANAACGAg                                                                                                                                                                                      >  2:759919/1‑300 (MQ=255)                                                                                                                     cACCAAAGCGTCTCTGGCAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGTTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCg                                                                                                                                                                                <  1:605533/300‑1 (MQ=255)                                                                                                                      aCCAAAGCGTCTCTGGCAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGt                                                                                                                                                                               >  2:896032/1‑300 (MQ=255)                                                                                                                                             gTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAAc                                                                                                                                                        >  2:772170/1‑300 (MQ=255)                                                                                                                                              tCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAAcc                                                                                                                                                       <  2:580093/300‑1 (MQ=255)                                                                                                                                               ccTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAAccc                                                                                                                                                      >  1:304254/1‑300 (MQ=255)                                                                                                                                                cTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCg                                                                                                                                                     >  1:834465/1‑300 (MQ=255)                                                                                                                                                             gCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGACGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCCGGCAGAGCTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTAGGCGGGttt                                                                                                                                        >  2:880719/1‑300 (MQ=255)                                                                                                                                                                 cGTTCCAGGAGACCCCTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTANCAACAG‑ATCGTATGCGTCGCCGTGCTNCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCGGAAGACGCATCTGCCAGCATGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTGATCGTTAATCCGCACATAACGTAAAAACCCGCTTAGGAGGGttttttt                                                                                                                                    >  2:605487/1‑300 (MQ=255)                                                                                                                                                                 cGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGttttttt                                                                                                                                    >  1:337696/1‑300 (MQ=255)                                                                                                                                                                  gTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTa                                                                                                                                   >  1:590521/1‑300 (MQ=255)                                                                                                                                                                  gTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGAGGGTTTTTTTa                                                                                                                                   >  2:506072/1‑300 (MQ=255)                                                                                                                                                                   ttCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACNAGTAATCGTTAATCCGCAAATAACGTAAAAAACCGCTTCGGCGGGTTTTTTTAt                                                                                                                                  >  2:673032/1‑300 (MQ=255)                                                                                                                                                                   ttCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTAt                                                                                                                                  >  1:555309/1‑300 (MQ=255)                                                                                                                                                                   ttCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAAACCGCTTCGGCGGGTTTTTTTAt                                                                                                                                  >  2:310627/1‑300 (MQ=255)                                                                                                                                                                   ttCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAACCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTAt                                                                                                                                  >  2:150486/1‑300 (MQ=255)                                                                                                                                                                   ttCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGATGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTAt                                                                                                                                  >  1:195623/1‑300 (MQ=255)                                                                                                                                                                    tCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCATAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATg                                                                                                                                 <  2:115307/300‑1 (MQ=255)                                                                                                                                                                             ccACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGG‑TCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGtt                                                                                                                        >  2:13414/1‑300 (MQ=255)                                                                                                                                                                                cTCGCGTTCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAg                                                                                                                     >  1:155479/1‑300 (MQ=255)                                                                                                                                                                                   gcgggCTGACTGAAGCAGCTGTTGCGGGCAAACGCGACGAACGGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCTTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCGGATAACGAGGAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGa                                                                                                                  <  1:425429/296‑1 (MQ=255)                                                                                                                                                                                    cgTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGaa                                                                                                                 >  2:708935/1‑300 (MQ=255)                                                                                                                                                                                        cTGACCTAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAagag                                                                                                             <  1:270040/300‑1 (MQ=255)                                                                                                                                                                                        cTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGTGCTGCGGTTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTGTAGGGAAagag                                                                                                             <  1:884951/300‑1 (MQ=255)                                                                                                                                                                                              gAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTg                                                                                                       >  1:603398/1‑300 (MQ=255)                                                                                                                                                                                               aaGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGt                                                                                                      >  1:528167/1‑300 (MQ=255)                                                                                                                                                                                               aaGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGt                                                                                                      >  2:861206/1‑300 (MQ=255)                                                                                                                                                                                                     ccGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAat                                                                                                <  1:545485/300‑1 (MQ=255)                                                                                                                                                                                                                   cgcgACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCAGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTc                                                                                  >  2:335297/1‑300 (MQ=255)                                                                                                                                                                                                                         gaACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATa                                                                            >  1:244500/1‑300 (MQ=255)                                                                                                                                                                                                                          aaCTGCGCGGCCTGAAATAGAACGTTATCGTGGGTCGTCTGATCCCTGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAGGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATAc                                                                           <  1:496891/300‑1 (MQ=255)                                                                                                                                                                                                                              gcgcGGCCTGAAAAAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGGCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAACAACGTAGAAACACGCTTCGGTGGGTTTTTTTATGGGTGGAGTTTCGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCAt                                                                       >  2:26919/1‑300 (MQ=255)                                                                                                                                                                                                                                  ggCCTGAAAGAGAACGTTATCGTGGTTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATc                                                                   <  1:708397/300‑1 (MQ=255)                                                                                                                                                                                                                                  ggCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATc                                                                   <  1:135344/300‑1 (MQ=255)                                                                                                                                                                                                                                   gCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTCCCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTCCTGAATATTCATAATCa                                                                  >  2:19508/1‑300 (MQ=255)                                                                                                                                                                                                                                   gCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCa                                                                  >  1:343999/1‑300 (MQ=255)                                                                                                                                                                                                                                   gCCTGAAAGAGAACGTTATCGCGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCa                                                                  >  1:537229/1‑300 (MQ=255)                                                                                                                                                                                                                                      tGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGNTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATg                                                               <  1:653327/300‑1 (MQ=255)                                                                                                                                                                                                                                       gAAAGAGAACGTTATCGTGGGTCGTCTGACCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACCGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCNGAACGCAGGTCTGGGCGGTTCTGATANCGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGt                                                              <  2:608314/300‑1 (MQ=255)                                                                                                                                                                                                                                        aaaGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTa                                                             >  2:706835/1‑300 (MQ=255)                                                                                                                                                                                                                                         aaGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCTGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTag                                                            <  2:475107/300‑1 (MQ=255)                                                                                                                                                                                                                                            aggaCGTTATCGCGGGGCGTCTGATCCCTGCAGGTACCAGTTACGCGTACCACCAT‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGTCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTNGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCGGATAACGAGTAATCGTTAAGCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTGTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGGCAGAAGATTTAAGGAATTTCTGAATACTCATAAGCAATGTAgaga                                                         <  2:227712/297‑1 (MQ=255)                                                                                                                                                                                                                                            agaACGTTATCGTGGTTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAgaga                                                         <  2:161770/300‑1 (MQ=255)                                                                                                                                                                                                                                            agaACGTTATCGTGGTTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAgaga                                                         <  1:8072/300‑1 (MQ=255)                                                                                                                                                                                                                                                cgctatnttGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATNGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAACAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGTGGGGAGTTTAGTGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGAc                                                     <  2:748796/292‑1 (MQ=255)                                                                                                                                                                                                                                                 gTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACt                                                    <  1:138090/300‑1 (MQ=255)                                                                                                                                                                                                                                                    aTCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGNGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTNTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAat                                                 >  2:841276/1‑300 (MQ=255)                                                                                                                                                                                                                                                    aTCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGTGCTGTGGGTGAAGCTCAGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAat                                                 >  2:443655/1‑300 (MQ=255)                                                                                                                                                                                                                                                    aTCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAat                                                 >  2:403522/1‑300 (MQ=255)                                                                                                                                                                                                                                                    aTCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGATGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTNGCAGAACTGCTGAACGCAGGTCTGGGNGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTAGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAat                                                 >  2:614284/1‑300 (MQ=255)                                                                                                                                                                                                                                                      cGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAAGACTCATAATCAATGTAGAGATGACTAatat                                               <  2:191148/300‑1 (MQ=255)                                                                                                                                                                                                                                                        tGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATcc                                             <  1:508833/300‑1 (MQ=255)                                                                                                                                                                                                                                                        tGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACGGCTGAACGCAGGTCTGGGCGGTTCTGAGAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATcc                                             <  1:385806/300‑1 (MQ=255)                                                                                                                                                                                                                                                                  gATCCCGGCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAactgac                                   >  1:228711/1‑300 (MQ=255)                                                                                                                                                                                                                                                                        ggCAGGTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCAGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAACTGACTGAACt                             >  2:301063/1‑300 (MQ=255)                                                                                                                                                                                                                                                                             gTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAACTGACTGAACTAATTg                        <  1:159784/300‑1 (MQ=255)                                                                                                                                                                                                                                                                             gTACCGGTTACGCGTACCACCAG‑ATCGTATGCGTCACCGTGCTGCGGTTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAACTGACTGAACTAAttt                        <  2:61215/300‑2 (MQ=255)                                                                                                                                                                                                                                                                                 cGGTTACGCGTACCACCAG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTCGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAAAGCTGACTAATATCCTGAAACTGAATGAACTAATTGAGta                    >  2:237230/1‑299 (MQ=255)                                                                                                                                                                                                                                                                                             ccaccaG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGCATACTCATAATCAATGTAGAGATGACTAATATCCTGAAACTGACTGAACTAATTGAGTCAAATTAGGCAAg        >  2:90847/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                             ccaccaG‑ATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAACTGACTGAACTAATTGAGTCAAACTCGGCAAg        >  1:921434/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                                accagATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAACTGACTGAACTAATTGAGTCAAACTCGGCAAGGa      >  1:464433/5‑300 (MQ=255)                                                                                                                                                                                                                                                                                                    gATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAACTGACTGAACTAATTGAGTCAAACTCGGCAAGGATTCg  >  1:206542/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                                                            AGCTCCGACTTCCTGGAAGGCGAACAGGTTGAATACTCTCGCGTCAAGATCGCAAACCGCGAACTGGAAGCGAACGGCAAAGTGGGTGCAACTTACTCCCGCGATCTGCTGGGTATCACCAAAGCGTCTCTGGCAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAACTGACTGAACTAATTGAGTCAAACTCGGCAAGGATTCG  >  NC_000913/4189160‑4189750

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