Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,446,813	G→A	6.4%	P221S (CCG→TCG)	aroC ←	chorismate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,446,813	0	G	A	6.4%	193.4 / 3.0	62	P221S (CCG→TCG)	aroC	chorismate synthase
Reads supporting (aligned to +/- strand):  ref base G (28/30);  new base A (2/2);  total (30/32)
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.98e-01

CTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGCAAAAAAACGGATTTTGCTCGACCTGCGACCAGTCTTTGATATCCAGCGGAATGTCGCCCATCTGGGTCAGGCAGCCACGGATTTCAATACCAAATTTCTCGGCGAGATATTTTTTGGCAATAGCTCCTGCCGCCACGCGCATGGCGGTTTCGCGGGCGGAAGAACGTCCAC  >  NC_000913/2446522‑2447106                                                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                                                      cTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCa                                                                                                                                                                                                                                                                                                >  1:531579/1‑299 (MQ=255)         ctcaccTCAATTTTGTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTCTGTTTCAGCGCCATATAGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCgg                                                                                                                                                                                                                                                                                                       <  2:317885/283‑1 (MQ=255)                cAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAAc                                                                                                                                                                                                                                                                                  >  1:726698/1‑298 (MQ=255)                 aaCTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCCTCCAGGCGGTCAAAGACCGGCTCGCCaa                                                                                                                                                                                                                                                                                               >  2:278547/1‑284 (MQ=255)                                  ggTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATAGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAACAACGGt                                                                                                                                                                                                                                                                <  1:708170/298‑1 (MQ=255)                                        tGGTACGACCCGGCACGGTAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAACAACGGTGACTTTa                                                                                                                                                                                                                                                         >  1:353179/1‑299 (MQ=255)                                                cccGGCACGGTAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAACAACGGTGACTTTAGCGCCGAt                                                                                                                                                                                                                                                 <  1:277795/299‑1 (MQ=255)                                                              tgtGGGGGTCGGTCTCAGCGCCAGATGGGGAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGACACCGTCTTTGGTGATGTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAACAACGGTGACTTTAGCGCCGa                                                                                                                                                                                                                                                  <  2:849867/282‑1 (MQ=255)                                                                           ttCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGcgcac                                                                                                                                                                                                                       <  1:371312/298‑1 (MQ=255)                                                                            tCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAACAAAGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGcgcacgc                                                                                                                                                                                                                     >  1:982725/1‑299 (MQ=255)                                                                                    atatGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTc                                                                                                                                                                                                             <  1:201632/299‑1 (MQ=255)                                                                                        gggCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGACACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAACGCCGTCGCCACTTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCCGGAACGCCCCCGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAacgccc                                                                                                                                                                                                                       >  2:458498/1‑283 (MQ=255)                                                                                          gCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTa                                                                                                                                                                                                        >  1:518949/1‑298 (MQ=255)                                                                                            aaTGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAAc                                                                                                                                                                                                      <  1:314424/298‑1 (MQ=255)                                                                                                     ctgcCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGAttt                                                                                                                                                                                            >  1:230169/1‑299 (MQ=255)                                                                                                            ctgctgATACCGCCGAGAATGCCGCCCTCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCggg                                                                                                                                                                                     <  1:321593/299‑1 (MQ=255)                                                                                                                  aTACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCgggggatcg                                                                                                                                                                               >  1:689381/1‑294 (MQ=255)                                                                                                                   tACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCAGCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCGGGGTCCggg                                                                                                                                                                              <  1:527426/299‑1 (MQ=255)                                                                                                                    aCCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAAATCGTCTAACGAGTCGATTTTGTCGGGGTCCGGGc                                                                                                                                                                             >  1:981882/1‑299 (MQ=255)                                                                                                                       gccgAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGCTTTTGTCGGGGTCCGGGCaaa                                                                                                                                                                          >  1:33114/1‑299 (MQ=255)                                                                                                                        ccgAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCg                                                                                                                                                                                       >  2:270173/1‑285 (MQ=255)                                                                                                                         cgggactGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCgg                                                                                                                                                                                      <  2:594239/279‑1 (MQ=255)                                                                                                                          gagaATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGCTTTTGTCGGGGTCCGGGCaaaaaa                                                                                                                                                                       >  1:516954/1‑299 (MQ=255)                                                                                                                                       gcATGGTTGCGCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGCAAAAAAACGGATTTTGCTc                                                                                                                                                          <  1:718510/299‑1 (MQ=255)                                                                                                                                               gCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGCAAAAAAAAGGAATTTGGTCGGcccgtc                                                                                                                                                  >  1:595947/1‑295 (MQ=255)                                                                                                                                                ctctGGAACCCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGCAAAAAAACGGAtt                                                                                                                                                                <  2:451994/284‑1 (MQ=255)                                                                                                                                                ctctGGAAACCGTCTGTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGCAAAAAAACGGAttt                                                                                                                                                               <  2:138783/285‑1 (MQ=255)                                                                                                                                                 tctGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGCAAAAAAACGGAttt                                                                                                                                                               >  2:166611/1‑284 (MQ=255)                                                                                                                                                 tctGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGACTCGCCAAGTCCGGCAGGAACGCCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGCAAAAAAACGGATTTTGCTCGACCTgaga                                                                                                                                                 >  1:49178/1‑298 (MQ=255)                                                                                                                                                         cctgcTTTGGTGATGTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTGGCCAAGTTCAACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAGCGCCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGCAAAAAAACGGATTTTGCTCGACCTGCGACCAGTCttt                                                                                                                                        <  1:989830/295‑1 (MQ=255)                                                                                                                                                           ggcTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGACTCGCCAAGTCCGGCAGGAACGCCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGCAAAAAAACGGATTTTGCTCGACCt                                                                                                                                                     <  2:40735/282‑1 (MQ=255)                                                                                                                                                                      ttccATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACGCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGCAAAAAAACGGATTTTGCTCGACCTGCGACCAGTCtt                                                                                                                                         <  2:415809/282‑1 (MQ=255)                                                                                                                                                                           tCGCGGTTCTGGCTGCCGCGCAGCGCCACCACTTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGCAAAAAAACGGATTTTGCTCGACCTGCGACCAGTCTTTGat                                                                                                                                     <  2:346355/284‑1 (MQ=255)                                                                                                                                                                              cgGTTCTGGCCGCCGCGCCGCGCCACCACGGCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGCAAAAAAGCGGATTTTGCTCGACCTGCGACCAGTCTTTGATATcc                                                                                                                                 <  2:450196/285‑1 (MQ=255)                                                                                                                                                                                      gCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGCAAAGAAACGGAGTTTGCTCGACCTGCGACCAGTCTTTGATATCCAGCGGAAt                                                                                                                         >  2:841142/1‑285 (MQ=255)                                                                                                                                                                                      gCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGACTCGCCAAGTCCGGCAGGAACGCCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGCAAAAAAACGGATTTTGCTCGACCTGCGACCAGTCTTTGATATCCAGCGGAAt                                                                                                                         >  2:512947/1‑285 (MQ=255)                                                                                                                                                                                        tGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCCTCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAACAACGGTGCCTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGACTTTGTCGGGGTCCGGGCAAAAAAACGGATTTTGCTCCACCTGCGCCCAGTCTTTGATATCCAGCGGACTGt                                                                                                                       >  2:651352/1‑285 (MQ=255)                                                                                                                                                                                        tGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGCAAAAAAACGGATTTTGCTCGACCTGCGACCAGTCTTTGATATCCAGCGGAATGt                                                                                                                       <  2:26567/285‑1 (MQ=255)                                                                                                                                                                                            gcgcAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCCACTCGTCTAACGCGTCGATTTTGTAGGGGTCCGGGCAAAAAAACGGATTTTGCTCGACCTGCGACCAGTCTTTGATATCCAGCGGCATGTCGc                                                                                                                    >  2:748577/1‑284 (MQ=255)                                                                                                                                                                                              gcAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATTAACTCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGCAACAAAACGGATTTTGCTCGACCTGATACCAGTCTTTGATATCCAGCGGAATGTCGCCCATCTGGGTCAGGCa                                                                                                    >  1:218261/1‑298 (MQ=255)                                                                                                                                                                                                      accacGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGCAAAAAAACGGATTTTGCTCGACCTGCGACCAGTCTTTGATATCCAGCGGAATGTCGCCCATCTGGGTc                                                                                                         <  2:841268/285‑1 (MQ=255)                                                                                                                                                                                                                      agCCAATTTCCACGCCTTTCACCGCGTGGATGCTCATCAGCGCATGGGCGAGGTCAGCATCCAGGCGTTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGCAAAAAAACGGATTTTGCTCGACCTGCGACCAGTCTTTGATATCCAGCGGAATGTCGCCCATCTGGGTCAGGCAGCCACGGAttt                                                                                         <  2:625878/284‑1 (MQ=255)                                                                                                                                                                                                                         cAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAACAACGGTGACTTTAGTGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGCAAAAAAACGGATTTTGCTCGACCTGCGACCAGTCTTTGATATCCAGCGGAATGTCGCCCATCTGGGTCAGGCAGCCACGGATTTCaa                                                                                      <  2:532128/285‑1 (MQ=255)                                                                                                                                                                                                                         cAATTTCCACGACGTTCACCGCGGTGATGCTCAGCAGCGCATGGGCGAGGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGCAAAAAAACGGATTTTGCTCGACCTGCGACCAGTCTTTGATATCCAGCGGAATGTCGCCCATCTGGGTCAGGCAGCCACGGATTTCaa                                                                                      <  2:461863/285‑1 (MQ=255)                                                                                                                                                                                                                                  cGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGCAAAAAAACGGATTTTGCTCGACCTGCGACCAGTCTTTGATATCCAGCGGAATGTCGCCCATCTGGGTCAGGCAGCCACGGATTTCAATACCAAAtt                                                                             >  2:276233/1‑285 (MQ=255)                                                                                                                                                                                                                                     cTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGCAAAAAAACGGATTTTGCTCGACCTGCGACCAGTCTTTGATATCCAGCGGAATGTCGCCCATCTGGGTCAGGCAGCCACGGATTTCAATACCAAATTTCCCGGCGAGATAtttt                                                            >  1:25742/1‑299 (MQ=255)                                                                                                                                                                                                                                       ttCACCTCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGGCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGCAAAAAAACGGATTTTGCTCGACCTGCGACCAGTCTTTGATATCCAGCGGAATGTCGCCCATCTGGGTCAGGCAGCCACGGATTTCAATACCAAATTTCTCGGCGAGATAtttttt                                                          <  1:981731/299‑1 (MQ=255)                                                                                                                                                                                                                                          aCCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGCAAAAAAACGGATTTTGCTCGACCTGCGACCAGTCTTTGATATCCAGCGGAATGTCGCCCCTCTGGGTCAGGCAGCCACGGATTTCAATACCAAATTTCCCGGCg                                                                     >  2:257429/1‑285 (MQ=255)                                                                                                                                                                                                                                            cgcgTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGCAAAAAAACGGATTTTGCTCGACCTGCGACCAGTCTTTGATATCCAGCGGAATGTCGCCCATCTGGGTCAGGCAGCCACGGATTTCAATACCAAATTTCTCGGCGAGATATTTTTTGGCaa                                                     <  1:276038/299‑1 (MQ=255)                                                                                                                                                                                                                                             gcgTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCCCTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCACTCTTTTTTCAGCGCACGCATCAACTCGTCTACCGCGTCGTTTTTGTCGGGGTCCGGGCAAAACAACGGATTTTGCTCGACCTTCGACCAGTCCTTGATATACAGCGGAATGTCGCCCATCTGGGACAGGCAGCCACGGAAATCTAGACCAAATTTCAAGGCgaga                                                                  >  2:391399/1‑285 (MQ=255)                                                                                                                                                                                                                                                  gATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGCAAAAAAACGGATTTTGCTCGACCTGCGACCAGTCTTTGATATCCAGCGGAATGTCGACCATCTGGGGCAGGCAGCCACGGATTTCAATACCAAATTTCTCGGCGAGATATTTTTTGGCAATAGCt                                                >  1:692303/1‑298 (MQ=255)                                                                                                                                                                                                                                                    tGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGCAAAAAAACGGATTTTGCTCGACCTGCGACCAGTCTTTGATATCCAGCGGAATGTCGCCCATCTGGGTCAGGCAGCCACGGATTTCAATACCAAATTTCTCGGCGAGATATTTTTTGGCAATAGCTcc                                              <  1:318993/298‑1 (MQ=255)                                                                                                                                                                                                                                                     gCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGCAAAAAAACGGATTTTGCTCGACCTGCGACCAGTCTTTGATATCCAGCGGAATGTCGCCCATCTGGGTCAGGCAGCCACGGATTTCAATACCAAATTTCTCGGCGAGATAtttttt                                                          <  2:305846/285‑1 (MQ=255)                                                                                                                                                                                                                                                     gCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGCAAAAAAACGGATTTTGCTCGACCTGCGACCAGTCTTTGATATCCAGCGGAATGTCGCCCATCTGGGTCAGGCAGCCACGGATTTCAATACCAAATTTCTCGGCGAGATAttttt                                                           >  2:677587/1‑284 (MQ=255)                                                                                                                                                                                                                                                       tcatcaGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGCAAAAAAACGGATTTTGCTCGACCTGCGACCAGTCTTTGATATCCAGCGGAATGTCGCCCATCTGGGTCAGGCAGCCACGGATTTCAATACCAAATTTCTCGGCGAGATATTTTTTgg                                                        <  2:198081/285‑1 (MQ=255)                                                                                                                                                                                                                                                         atcaGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCATGCATCAACTCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGCAAAAAAACGGATTTTGCTCGACCTGCGACCAGTCTTTGATATCCAGCGGAATGTCGCCCATCTGGGTCAGGCAGCCACGGATTTCAATACCAAATTTCTCGGCGAGATATTTTTTGGCAATAGCTCCTgccg                                         >  1:42830/1‑298 (MQ=255)                                                                                                                                                                                                                                                               gcATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGACTCGCCAAGTCCGGCAGGAACGCCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGCAAAAAAACGGATTTTGCTCGACCTGCGACCAGTCTTTGATATCCAGCGGAATGTCGCCCATCTGGGTCAGGCAGCCACGGATTTCAATACCAAATTTCTCGGCGAGATATTTTTTGGCAATAGCTCCTGCCGCCAcgc                                   <  1:593509/298‑1 (MQ=255)                                                                                                                                                                                                                                                                      cGATGTCCGCATCCAGGCGGGCAAAGACCTGCTCGCCAGGTCCTGCATGAACTCCACGGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGAGCGCATCACCTCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGCAAAAAAACGGATTTTGCTCGACCTGCGACCAGTCTTTGATATCCAGCGGAATGTCGCCCATCTGGGTCAGGCAGCCACGGATTTCAATACCAAATTTCTCGGCGAGATATTTTTTGGCAATAGCTCCTgccg                                         <  2:20854/285‑1 (MQ=255)                                                                                                                                                                                                                                                                                    aGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAACAACGGTGACTTTAGCGCCGATGGCGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGCTTTTGTCGGGGTCCGGGCAAAAAAACGGATTTTGCTCGACCTGCGACCAGTCTTTGATATCCTGCGGAATGTCGCCCATATGGGTCAGGCAGCCACGGATTTCACTACCACATTTCTCGGCGAGATATTTTTTGGCAATAGCTCCTGCCGCCACGCGCATGGCg                           >  2:173905/1‑285 (MQ=255)                                                                                                                                                                                                                                                                                           cAAAGACAGGCTCGGCAAGTCCGGCAGGAACGCCACTGGCAACAACGGTGACTTTAGCGCCGAGGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGCAAAAAAACGGATTTTGCTCGACCTGCGACCAGTCTTTGATATCCAGCGGAATGTCGCCCATCTGGGTCAGGCAGCCACGGATTTCAATACCAAATTTCTCGGCGAGATATTTTTTGGCAATAGCTCCTGCCGCCACGCGCATGGCGGTTTCGCGGGCGGAAGAACGt      <  1:758245/299‑1 (MQ=255)                                                                                                                                                                                                                                                                                             aaGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCTATTTTGTCGGGGTCCGGGCAAAAAAACGGATTTTGCTCGACCTGCGACCAGTCTTTGATATCCAGCGGAATGTCGCCCATCTGGGTCAGGCAGCCACGGATTTCAATACCAAATTTCTCGGCGAGATATTTTTTGGCAATAGCTCCTGCCGCCACGCGCATGGCGGTTTCGCGGGCGGAAGAACGTcc    <  1:297604/299‑1 (MQ=255)                                                                                                                                                                                                                                                                                              aGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGCAAAAAAACGGATTTTGCTCGACCTGCGACCAGTCTTTGATATCCAGCGGCATGTCGCCCATCTGGGTCAGGCAGACACGGATTTCAATACCAAATTTCTCGGCGAGATATTTTTTGGCAATAGCTCCTGCCGCCACGCGCATGGCGGTTTCgcgg                  >  2:773055/1‑284 (MQ=255)                                                                                                                                                                                                                                                                                                 ccGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGCAAAAAAACGGATTTTGCTCGACCTGCGACCAGTCTTTGATATCCAGCGGAATGTCGCCCATCTGGGTCAGGCAGCCACGGATTTCAATACCAAATTTCTCGGCGAGATATTTTTTGGCAATAGCTCCTGCCGCCACGCGCATGGCGGTTTCGCGGGCGGAAGAACGTCCAc  <  1:312300/297‑1 (MQ=255)                                                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                                                      CTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGCAAAAAAACGGATTTTGCTCGACCTGCGACCAGTCTTTGATATCCAGCGGAATGTCGCCCATCTGGGTCAGGCAGCCACGGATTTCAATACCAAATTTCTCGGCGAGATATTTTTTGGCAATAGCTCCTGCCGCCACGCGCATGGCGGTTTCGCGGGCGGAAGAACGTCCAC  >  NC_000913/2446522‑2447106

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑