Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,620,256	C→T	G218G (GGC→GGT)	deoB →	phosphopentomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,620,256	0	C	T	94.3%	183.1 / ‑3.8	75	G218G (GGC→GGT)	deoB	phosphopentomutase
Reads supporting (aligned to +/- strand):  ref base C (3/1);  new base T (30/36);  total (38/37)
Fisher's exact test for biased strand distribution p-value = 3.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAAGAAACTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCTACGCCAACTGCGGTATCACCAAAAAAGTGAAAGCGACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTC  >  NC_000913/4620010‑4620499                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      tAATCTACCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGCGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAACCTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCCACACAATCGGTCGGGTGTA‑CGCCCATCCGTTTATCGGTGAc                                                                                                                                                                                                                                                    >  2:587606/1‑250 (MQ=255)  aaTCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGGCGTGT‑TATAGCTCGTCCGTTTATCGGTGACa                                                                                                                                                                                                                                                   >  1:78078/1‑250 (MQ=255)   atcggCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATGCGGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACaaa                                                                                                                                                                                                                                                 <  1:607517/247‑1 (MQ=255)                      aaccggCCACTCTGCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTa                                                                                                                                                                                                                              <  2:316872/246‑1 (MQ=255)                        aCTGCCACTCTTCCGGTACGGTCATTCTGGAACAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTGCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCTGTAACTTCCAGCGTAcc                                                                                                                                                                                                                            <  2:526310/251‑1 (MQ=255)                           gCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCgg                                                                                                                                                                                                                          <  1:319516/250‑1 (MQ=255)                            cccgacgctcttccgaactaTTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTa                                                                                                                                                                                                                        <  2:29512/232‑1 (MQ=255)                             cACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTa                                                                                                                                                                                                                        >  1:301911/1‑250 (MQ=255)                                       gTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGGCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGAc                                                                                                                                                                                                             <  1:26666/251‑1 (MQ=255)                                         aCGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCt                                                                                                                                                                                                           <  1:301163/251‑1 (MQ=255)                                           ggTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAACTCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACAACCGTCGTGC‑TATCGCTCGTCCGGTTATCGGTGACCAAGCCGGTAACATCCCGCGTACCGGTCACCGTCACGAACTgg                                                                                                                                                                                                         >  2:251500/1‑251 (MQ=255)                                              cATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTg                                                                                                                                                                                                      <  1:130413/251‑1 (MQ=255)                                               aTTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTa                                                                                                                                                                                                                        >  1:29512/1‑232 (MQ=255)                                                  cTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTgg                                                                                                                                                                                                 >  1:94773/1‑250 (MQ=255)                                                      aTCAGCTGGGCGAAGAGCACATGAAAAACGGCAAGCCGATTTTCTATACCTCCGCTGACTACGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAAATCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCCGTCGGTT‑TATCGCTCGTCCCTTATTCGGTGACAAAACCGGTCACTTCCCACGTACCGGTAACCGtaccgaccgggaagtggcgacg                                                                                                                                                                                              >  2:675156/1‑230 (MQ=255)                                                             gggCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAACCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCGG‑TTGAGCCGCCCGCac                                                                                                                                                                                      >  1:465348/1‑251 (MQ=255)                                                                  aaGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATAGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCaccgac‑c                                                                                                                                                                                >  1:478383/1‑251 (MQ=255)                                                                      gCACATGAAAACCGGCAAGACGATTTTCTATACCTCCGCTGACTCCGTGATCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGGGCGAAATCGCCCGAGAAGAGCTGACAAACGGCGGCTACAAACTCGGTCGGGG‑TAACGCACGAACGTTTAACGGGGACCAAGCAGGTAAATTCCAGCGTACCGGTAAACGTGACGACCTGGCTG‑GGGGGCCTCCAGCaccgac‑cgggc                                                                                                                                                                            >  2:9619/1‑248 (MQ=255)                                                                       cacaTGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGtgct                                                                                                                                                                           >  1:539849/1‑251 (MQ=255)                                                                         caTGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAAACGTCACGACCTGGCTG‑TTGGGCCGCCAGCCACGAC‑CGtgctgg                                                                                                                                                                         >  1:133902/1‑250 (MQ=255)                                                                          aTGAAACCCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCa                                                                                                                                                                        <  1:57819/251‑1 (MQ=255)                                                                             aaaaCCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGGGACAAAGACGGTAACTTCCAACGTACCGGGAACCGTCACGACCGGGCTG‑GTGGGCCCCCAGGacccgccctggtgcaga                                                                                                                                                                       >  2:558897/1‑236 (MQ=255)                                                                                aCCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTTACAAAGCCGGTAACTTCCAGCGTACCGGGAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACt                                                                                                                                                                  >  1:537280/1‑251 (MQ=255)                                                                                aCCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACt                                                                                                                                                                  <  1:651832/251‑1 (MQ=255)                                                                                aCCGGAAAGCCTATTTTCGATCCCTCTGCTGATTCCGTGTTCCAGATTGCCTTCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTTCGAAATCGCCCTTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCGCGTCCGTTTACCGGTGACAAAGCCGGTAACTTCCAGCGTTCCGTTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACt                                                                                                                                                                  <  2:106666/251‑1 (MQ=255)                                                                                 ccGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGGCACGCCCTGGCTG‑TTGGGCCGCCAGCACCGAA‑CGGGCTGCAGAAACTg                                                                                                                                                                 >  1:557871/1‑251 (MQ=255)                                                                                    gCAAGCCGATTTTCTATTCCTCCGCGGACGCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGAAAAAGCCGTTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGtt                                                                                                                                                              <  1:47343/251‑1 (MQ=255)                                                                                          cGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTTTCGGGGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCCGCCCCCAC‑CGCGCTGCAGAAACTGGTtgatta                                                                                                                                                         >  2:235200/1‑248 (MQ=255)                                                                                             ttttCTATACCTCCGCTGACTCCGTGTTCCCGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGaaaa                                                                                                                                                      >  1:520621/1‑250 (MQ=255)                                                                                               ttCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAAATGGTTGATGAAAAac                                                                                                                                                    >  1:90776/1‑250 (MQ=255)                                                                                                         ccGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTGCAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAggtg                                                                                                                                         <  1:235200/251‑1 (MQ=255)                                                                                                          cGCTGACTCCGTGTTCCAGTTTGCTGGACACGAAGGAACCGTTGGCCTGGATAAACACTACGAAATGTGCGAAATCGCCCGTGAAGAGCGGCCCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTGGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAggtgg                                                                                                                                        <  2:172989/251‑1 (MQ=255)                                                                                                            cgtactccgtttgtcAGATTGCCTGCCCAGAAGGAACCTTTGGCCTGGATAAACTCTACGTACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGtt                                                                                                                                      <  2:564585/237‑1 (MQ=255)                                                                                                             tGACTCCGTGTTCCAGATTGCCTGCCATGAAGAAACTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGA‑TATCGCTCGTCCGTTTATCGGCGACAAAGCCGGACACTTCCAGCGTACCGGGAAACGGCACGACCTGGCTGGTTG‑GCCGCCAGCCCCGACACGGCCTAC‑GAAACCTGTTGATGAAAAACACGGGCAggagggtt                                                                                                                                     >  2:332042/1‑245 (MQ=255)                                                                                                                cTCCGTGTTCCAGATTGCCTGCCATGAAGAAACTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAac                                                                                                                                                    <  1:235364/233‑1 (MQ=255)                                                                                                                cTCCGTGTTCCAGATTGCCTGCCATGAAGAAACTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGAGT‑TATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGAACTGGCTG‑TTGGGCCGCCAGCACCGAC‑CGGGCTGCAGAAACTGGGTGGTGTAAAacagaacggaagagaggcgg                                                                                                                                  >  2:235364/1‑234 (MQ=255)                                                                                                                  ccGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGTTTCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTc                                                                                                                                <  1:587606/251‑1 (MQ=255)                                                                                                                       ttCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGGGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGGGCCGCCAGCAACGAC‑CGTGCTGCAGAAACCGGTTGATGAAAAAAACGGTAAGGGGGGTTCTGTCGGTa                                                                                                                            >  1:371230/1‑250 (MQ=255)                                                                                                                         ccAGATGGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTTTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTaaa                                                                                                                          <  2:94773/250‑1 (MQ=255)                                                                                                                           aGATTGCCTGCCACGAAGGAACCTTGGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACGTCCAGCGTACCGTTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAAtt                                                                                                                       <  1:319198/251‑1 (MQ=255)                                                                                                                              ttGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTtgatga                                                                                                                                                         >  1:360305/1‑214 (MQ=255)                                                                                                                                gcctgccACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGcaga                                                                                                                  <  1:626593/251‑4 (MQ=255)                                                                                                                                              aacctttggccTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGTCGGCTACAACATCGGTCAAGT‑TATCGCTCGACCGTTTATCGGAGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGATGCATAAACTGGTTGATGAAAAAAAAGGACAGGGGGTTTCTGTCGGTAAAATTGAAGACATATACGCCAACt                                                                                                    >  1:397843/11‑251 (MQ=255)                                                                                                                                                  tttggccTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGATACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCAAGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGGGCTGCCGAAACTGGTTGATGAAAAACACGGTCAGGGGGTTTCTGTCGGTAAAATTGCAGACCTCTACGCACACTGCgg                                                                                                >  1:208218/7‑251 (MQ=255)                                                                                                                                                    tggccTGGATAAACCCTACGAACTGGGCGAAATCGCCCGTGTAGAGCTGACCAACGGCGGCTACAACAGCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGt                                                                                               <  2:345062/246‑1 (MQ=255)                                                                                                                                                      gccTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGACGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCACCACCGCC‑CGTGCTGCAGAAACTGGTTGATGGAAAACACGGTCAGGTGGTTTCTGTAGGTAAAATTTCAGACATCTAAGACAACTCCGGtgtc                                                                                            >  1:105480/3‑248 (MQ=255)                                                                                                                                                       ccTGGATTAACTCTCCGAACTGTGCGAAATCGCGCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGATAAAGCCGGTAATGTCCAGCGTACCGGTAACCGTCACGACCTGGTTG‑TTGAGCCGCCATCACCGAC‑CGTGCTTCAGAAACTGGTTGATGAAAAACACTGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCa                                                                                           <  1:617375/250‑1 (MQ=255)                                                                                                                                                       ccTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCa                                                                                           <  1:446750/250‑1 (MQ=255)                                                                                                                                                           gtaaaaccctccGAACTGTTCGAAACTGCCAGTGAAGTGCTGCCCAACGGCGGCTACAACTTCGTCCGTTT‑GATCGCTCGTCCGTTTATCGGTTCCAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACAACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCaa                                                                                       <  2:172842/240‑1 (MQ=255)                                                                                                                                                           gATAAACTCTACGAACTGGGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCaa                                                                                       <  1:501122/251‑1 (MQ=255)                                                                                                                                                              aaaCTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCaaga                                                                                     >  1:513447/1‑248 (MQ=255)                                                                                                                                                                            gtgCGAAATCGCCCGTGAAGAGCTGTCCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTg                                                                      <  2:363562/251‑1 (MQ=255)                                                                                                                                                                            gtgCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACt                                                                       <  1:231169/250‑1 (MQ=255)                                                                                                                                                                                     cGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGGGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGAc                                                              <  1:168044/250‑1 (MQ=255)                                                                                                                                                                                      gCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGAcgc                                                            >  1:156971/1‑251 (MQ=255)                                                                                                                                                                                           tGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTaa                                                                                                                           >  1:235136/1‑183 (MQ=255)                                                                                                                                                                                           ggAAGAGCTGACCAACGGCGGCTACAACATCTGTCGTGT‑TATCGCTCGTCCGTTTACCGGTGACAAAGCCGGTACCTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTaa                                                                                                                           <  2:235136/182‑1 (MQ=255)                                                                                                                                                                                             aaGAGCTGACCAACGGCGGCTACAACATCGGTCGGGT‑TATCGCTCGTCTGTGTATCGGTGAAAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGGTGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGttt                                                     <  2:388344/251‑1 (MQ=255)                                                                                                                                                                                                   tGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGcc                                               >  1:386476/1‑251 (MQ=255)                                                                                                                                                                                                        aaCGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCAt                                          <  1:294656/251‑1 (MQ=255)                                                                                                                                                                                                             cggcTACAACATCGGCCGTGT‑TATCGCTCGCCCGTTTATCGGCGACAAAGCCGGTAAATTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCCCCGAC‑AGCGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGCGGTTTCTGTCGGTAAAATTGCAGACCTCTACGCCAACTGCGGTACCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGAAGCCACCAGCaa                                       >  2:293571/1‑249 (MQ=255)                                                                                                                                                                                                               gcTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAag                                     >  1:315824/1‑249 (MQ=255)                                                                                                                                                                                                                cTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAag                                     >  1:34968/1‑248 (MQ=255)                                                                                                                                                                                                                    aacatCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTAACGGTAACCGTCACGAACTGGCTG‑TTGAGCCGCCAGCACCGCC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGATTCTGTCGGCAAACTTGCAGACATCTACGCCACCTGCGGCATCACCACGAAAGTGACAGCAACTGGCCTGGACGCGCTGTTTGACACCACCATCAAAGAGCTGa                               >  1:286307/4‑250 (MQ=255)                                                                                                                                                                                                                         cgggcgttt‑tATCGCGCGTCCGTTTAGCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGACACGTCCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAaagaag                         <  2:46657/243‑1 (MQ=255)                                                                                                                                                                                                                                 t‑tATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGAt                 >  2:167342/1‑251 (MQ=255)                                                                                                                                                                                                                                    atctcTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCg                                                                                                 <  2:343004/165‑1 (MQ=255)                                                                                                                                                                                                                                    atctcTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCg                                                                                                 >  1:343004/5‑169 (MQ=255)                                                                                                                                                                                                                                    agcGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCAGCTAGCACCGAC‑CGTGCTTCAGAAACTGGTTGATGAAAAACACGGTCAGGTTGTTTCTTTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAATAAAGTGAAAGCAACTGGCCTGGACGCGTTTTTTGACGCCTCCATCAAAGAGATGAAAGAAGCGGGTGATaa               <  1:605667/249‑1 (MQ=255)                                                                                                                                                                                                                                      cGCTCGTCCGTTTATCGGTGACAAAGCCGTTAACTTCCAGCGTACCTGTAACCTTCACGACCTGGCTT‑GTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCAAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGTGATGAAAGAAGCGGGTGATAac              <  1:400470/250‑1 (MQ=255)                                                                                                                                                                                                                                      cGCTCGTCCGTTTACCGGTATCAAAGCCGGTAGCTTCCAGCGTGCCGGAAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGGAACTGGTTGATGAAAAACACGGTCAGGTGTTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTTGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAaca             <  1:285304/251‑1 (MQ=255)                                                                                                                                                                                                                                          cGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCAt         >  1:562635/1‑251 (MQ=255)                                                                                                                                                                                                                                            tCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCg       >  1:591204/1‑251 (MQ=255)                                                                                                                                                                                                                                                tttATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCATACATCTACGCCAACTGCGGCATCACCAAGAAAGAGAAGGCAACTGGCCTGGACGCGCTGTTTGACGCCCCCATCAAAGAGATGAAAGAATCggtgaataaccccaacgcctt   >  2:256456/1‑232 (MQ=255)                                                                                                                                                                                                                                                  tATCGGTGAAAAAGCCGGAAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTc  <  2:156971/250‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      TAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAAGAAACTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCTACGCCAACTGCGGTATCACCAAAAAAGTGAAAGCGACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTC  >  NC_000913/4620010‑4620499

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