Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	13,539	A→G	L459L (CTA→CTG)	dnaK →	chaperone protein DnaK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	13,539	0	A	G	97.4%	252.0 / ‑5.9	77	L459L (CTA→CTG)	dnaK	chaperone protein DnaK
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (32/43);  minor base T (0/2);  total (32/45)
Fisher's exact test for biased strand distribution p-value = 5.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

TGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGT‑TACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTAGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACC  >  NC_000913/13299‑13783                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     tggtgTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGT‑TACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGGGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCCGTTCAACCTGGATGGTATc                                                                                                                                                                                                                                             >  1:56442/1‑250 (MQ=255)       tCTGACTGGTGACGTGAAAGACGTACTGCTTCTGGACGT‑TACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCg                                                                                                                                                                                                                                       <  1:514502/250‑1 (MQ=255)         tgaagtgtgtcgtgaaaGACGTACTGCGTCTGGA‑GTGTACCACGCTGTCTCTGGGTATCTAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCTCGCCCATTGGTTCTCTACCGCTGAAGTCAACCAGTCTGTGGTAACCATCCATGTGCTGCATGGTGATCTGAAATTTGCGTCTGATAATAAATCTCGGGTTCATTTCACCCTTGATGGTATCATCCCGGc                                                                                                                                                                                                                                     <  2:256202/236‑1 (MQ=255)              ggTGACGTGACAGACGTACTGCTGCTGGACGT‑TCCCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAACACACCCCTATCCCGACCAAGCACAGCAAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCCTGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAACTCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACcgcg                                                                                                                                                                                                                               >  2:584609/1‑251 (MQ=255)              ggTGACGTGAAAGACGTACTGCTGCTGGACGT‑TACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACcgc                                                                                                                                                                                                                                <  2:162301/250‑1 (MQ=255)                  aCGTGAAAGACGTACTGCTGCTGGACGT‑TACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAACACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTACCCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACCAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGc                                                                                                                                                                                                                           >  2:27652/1‑251 (MQ=255)                   cgtgaaaGACGTACTGCTGCTGGACGT‑TACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCa                                                                                                                                                                                                                          <  1:220883/247‑1 (MQ=255)                    gtgaaaGACGTACTGCTGCTGGACGT‑TACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCAt                                                                                                                                                                                                                         <  2:210367/248‑1 (MQ=255)                      gaaaGACGTACTGCTGCTGGACGT‑TACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAACCACCACTATCCAGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACACCCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATCACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGc                                                                                                                                                                                                                       >  2:300983/2‑251 (MQ=255)                       aaaGACGTACTGCTGCTGGACGT‑TACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGcc                                                                                                                                                                                                                      >  2:487514/1‑251 (MQ=255)                         aGACGTACTGCTGCTGGACGT‑TACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGc                                                                                                                                                                                                                    <  2:495190/251‑1 (MQ=255)                                         gACGT‑TACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCtt                                                                                                                                                                                                    <  1:203574/251‑1 (MQ=255)                                                 ccccGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATc                                                                                                                                                                                             <  1:156753/251‑1 (MQ=255)                                                            ctGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCTCTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGAc                                                                                                                                                                                    <  1:577129/249‑1 (MQ=255)                                                              gggTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACgg                                                                                                                                                                                  <  2:191034/249‑1 (MQ=255)                                                              gggTATCGAAACCATGAGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACgg                                                                                                                                                                                  >  1:191034/1‑249 (MQ=255)                                                                    cGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGc                                                                                                                                                                          >  1:437083/1‑251 (MQ=255)                                                                     gAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGc                                                                                                                                                                          >  1:17743/1‑250 (MQ=255)                                                                        cccATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTGCGATATCGATGCTGACGGTATCCTGCACGt                                                                                                                                                                      <  2:119372/250‑1 (MQ=255)                                                                               gCGGTGTGATGACGACGCGGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGTCACCGCGTGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCcgcg                                                                                                                                                               <  2:79662/251‑1 (MQ=255)                                                                                 ggTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGaa                                                                                                                                                             <  1:503415/251‑1 (MQ=255)                                                                                      gatgaCGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGAt                                                                                                                                                         >  2:172958/1‑250 (MQ=255)                                                                                         gacgacGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGGCAGGTGTTCTCGACCGCTGAAGTCAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGCTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATaaa                                                                                                                                                      <  1:309762/250‑1 (MQ=255)                                                                                          acgacgCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATaaaaa                                                                                                                                                    <  1:559805/251‑1 (MQ=255)                                                                                            gacgCTGATCGCGAAAAACACCACTATCCCGACCACGCACAGCCCGGTGTTCTCTACAGCTGAAGACAACCAGTCTGCGGAAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAAc                                                                                                                                                   >  2:39696/1‑250 (MQ=255)                                                                                               ccTGATCGCGAAAAACACCACTAGCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCTGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCg                                                                                                                                               <  2:420784/250‑1 (MQ=255)                                                                                                  gatagagAAACACACCACTATCCCGACCAAGCACCGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGCGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGCTGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGCAAGATAAAAACAGCGGt                                                                                                                                             >  2:429172/7‑250 (MQ=255)                                                                                                  gATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGt                                                                                                                                             >  1:73204/1‑250 (MQ=255)                                                                                                     cgcgAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAag                                                                                                                                         <  1:300983/251‑1 (MQ=255)                                                                                                         aaaatCACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCCGTCTGCGGTAACCATCCATGTGCGGCAGTGTGAACGTAAACGTGCGGCTGATAACAAATCGCTGGGTCAGTGCAACCTGGATGGTATCAACTCGGCACCGCGCGGCATGCCGGTGATCGAAGTTATCTGCGATATCGATGGTGACGGTATCCTGCACGTGTTCGCGAAAGATAGAATCAGCGGTAAAGAGCa                                                                                                                                     <  2:67117/251‑1 (MQ=255)                                                                                                          aaaaCACCACTATCCCGTCCAAGCACAGCCGAGTGTTCTCTACCGGTTAAGACAACCAGTCTGCGGGAACCATCCATGTGCTGCAGGGTGAACGTAATCGTGCGGCTGATAACAAATCTCTGGGTCATTTCAACCTGGATGGTATCAGCCCGGCTTCGCGCGTCATGCCGCAGATCGAAGTTACCTTCGATATCGATTCTGACGGGTTCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCag                                                                                                                                    <  1:50524/251‑1 (MQ=255)                                                                                                           aaaCACCACTACCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGAGAAAAACAGCGGTAAAGAGCaga                                                                                                                                   <  2:123096/251‑1 (MQ=255)                                                                                                            aaCACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCGGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCagaa                                                                                                                                  <  2:209473/251‑1 (MQ=255)                                                                                                             acacCACTATCCCGACCAAGCACAGCCATGTTTTCTCTACCGCTGAAGACATCCAGTCTGCGGTAACCAGCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAAGCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCTTCACCGCTCGGCATGCCGCAGATCGAAGTTTCCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGATCagaag                                                                                                                                 <  1:578205/251‑1 (MQ=255)                                                                                                               accacTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCaga                                                                                                                                   >  1:437291/1‑247 (MQ=255)                                                                                                                     aTCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCAcc                                                                                                                           <  1:15055/249‑1 (MQ=255)                                                                                                                      ccccGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAAGCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGTAGAAGATCTCCATc                                                                                                                        <  2:477395/250‑1 (MQ=255)                                                                                                                       cccGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATc                                                                                                                        <  2:469150/250‑1 (MQ=255)                                                                                                                         cGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCGGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGGTCAACCTGGATGGTATCAACCCGGCACCGCGGGGCATGCCGCAGAGCGAAGTTACCTTCGATATCGATGCTGACTGTATCCTGCACGTTTCCGCGAAATATAAAAACAGCGGTAAAGAGTGGGAGATCACCATCacg                                                                                                                     <  2:281174/251‑3 (MQ=255)                                                                                                                           aCCAAGAACCGCCAGGTGCTCTCTACCGCTGAAGACAACCAGCCTGAGGTAAACATCCATGAGCTGCAGGGTGACCGTAAACGTGCGGCTGCTAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCATCCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATACTGCACGTTTCCGCGCAAGATAAAAACAGCGGTAACGAGCAGAAGATCACCATCAAGGc                                                                                                                   >  1:399495/1‑251 (MQ=255)                                                                                                                             cAAGCCCAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATTTGCTGCAGGGGGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGAGGGTATCAACCCGGCGCCGCGCGGCATGCCGCAGATCGAAGGTACCTTCGATAGCGATGCTGACGGTATCCGGCACGTTTCCGCTAAAGATAAAAACAGCTGTAAAGAGCAGAAGATCACCATCAAGGctt                                                                                                                 <  2:578417/251‑1 (MQ=255)                                                                                                                               aGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTTTGCGGTAACCATCCAGGTGCTGCAGGGTGAACGTAAACGTGCGGCGGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGcttct                                                                                                               <  2:128911/251‑1 (MQ=255)                                                                                                                               aGCACAGCCAGGTGATCTCTACCGCTGAAGCCAACCAGTCTGCGGTCACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTATGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGCCGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGcttct                                                                                                               >  1:596906/1‑251 (MQ=255)                                                                                                                                 cactaccAGGTGTTCTCTGCCGCTGAAGACAACGAGACGGTGGTAACCATGCATGTGCTTGAGGGTGAACGGAAACGTGGGGCTGATAACAAAGCTCCGGGTCAGTTCAACCCGGATGGTATCAACCCGGCTACGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGAAGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGAAGAAGAACACCATCAAGGCttcttc                                                                                                             <  2:175415/246‑1 (MQ=255)                                                                                                                                   caGCCAGGTGTTCTCTACAGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTg                                                                                                           <  1:470433/251‑1 (MQ=255)                                                                                                                                    aGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGTTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTgg                                                                                                          <  1:145847/251‑1 (MQ=255)                                                                                                                                      ccAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTc                                                                                                        >  1:150394/1‑251 (MQ=255)                                                                                                                                       catgtgatctcTACCTCCGAAGACAACCAGTCTGCTGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCTGATGATAACAAATCTCTGGGTCAGTTCAACCTTGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCGTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAACTGTAAAGAGCAGAAGATCGTCAGCAAGGCTTCTTCTGGTCt                                                                                                       <  2:62058/244‑1 (MQ=255)                                                                                                                                        aGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTg                                                                                                      >  2:422257/1‑251 (MQ=255)                                                                                                                                           tgtTCTCTACCGCTGAAGACACCCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCCACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAAc                                                                                                   >  2:265310/1‑251 (MQ=255)                                                                                                                                                    ccGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGa                                                                                           >  1:135453/1‑250 (MQ=255)                                                                                                                                                      gCTGAAGACAACCAGTCTGCGGTAACCCTCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAAt                                                                                        >  2:219544/1‑251 (MQ=255)                                                                                                                                                       cTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTg                                                                                                           <  2:258496/231‑1 (MQ=255)                                                                                                                                                       cTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTg                                                                                                           >  1:258496/1‑231 (MQ=255)                                                                                                                                                           aGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATcc                                                                                      >  1:512292/1‑248 (MQ=255)                                                                                                                                                             aCAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGaaa                                                                                 <  1:414852/251‑1 (MQ=255)                                                                                                                                                                aCCAGTCTGCGGTAACCATCCATGTGCTGCAGGTTGAACGTAAACGTGCGGCTGAGAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAATATCACCATCAAGGCTTCGTCTGGTCTGAACGAAGATGAAATCCAGAAAATg                                                                              <  1:408606/251‑1 (MQ=255)                                                                                                                                                                   aGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTa                                                                           <  1:27652/251‑1 (MQ=255)                                                                                                                                                                    gTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTAc                                                                          <  1:474120/251‑1 (MQ=255)                                                                                                                                                                        gCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTGCAACCTGGATGGTATCAACCCGGCGCCGCGCGTCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTAcgc                                                                        <  1:584609/249‑1 (MQ=255)                                                                                                                                                                             aaCCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAg                                                                 <  1:151366/251‑1 (MQ=255)                                                                                                                                                                              aCCATCCATGTACTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGACAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGcgccgc                                                                   >  2:529099/1‑248 (MQ=255)                                                                                                                                                                                                     aCGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGt                                          >  1:421565/1‑250 (MQ=255)                                                                                                                                                                                                                        aaCAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTAc                       >  1:540824/1‑250 (MQ=255)                                                                                                                                                                                                                         aCAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTAc                       >  1:162699/1‑249 (MQ=255)                                                                                                                                                                                                                           aaaTCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGAc                   >  2:222997/1‑251 (MQ=255)                                                                                                                                                                                                                            aaTCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACt                  <  1:603376/251‑1 (MQ=255)                                                                                                                                                                                                                             aTCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCATGCACGTTTCCGCGAAAGATAAAAACAGCGGAAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGAAGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTc                 >  2:541349/1‑251 (MQ=255)                                                                                                                                                                                                                              tctcTGGGTCCGTTCAACCTGGATGGTGTCAACCCGGCACCGCGCGGCAGGCCGCAGATCGAAGTTACCTTCGATATCGCTGCTGACGGTATCCTGCACGTTTCCGCGGAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGGAATCCAGAAACTGGTACGCGACGGAGAAGCTAACGCCGaa                                                   >  2:404803/1‑216 (MQ=255)                                                                                                                                                                                                                              tctcTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGaa                                                   <  1:404803/216‑1 (MQ=255)                                                                                                                                                                                                                                tctGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTCAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGCAGATGAAATCCAGAAAATGGTAAGCGACGCAGAAGCTAACGCCGAAGCTGCACGTAAGTTTGAAGAGCTGGTACAGACTCGc               >  2:344219/1‑250 (MQ=255)                                                                                                                                                                                                                                 cggggTCAGTTCAACCTGGATGGTATCAACCCGGCGCCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCGGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAACTGTTACAGACTCGCaa             <  2:421565/249‑1 (MQ=255)                                                                                                                                                                                                                                       cAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAATAGTTGGTACATACGCGCAACCAgg        <  2:117266/250‑1 (MQ=255)                                                                                                                                                                                                                                        aGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGc      <  1:425813/251‑1 (MQ=255)                                                                                                                                                                                                                                           tCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGAc   >  1:40296/1‑251 (MQ=255)                                                                                                                                                                                                                                            cAACCTGGATGGTATCACCCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGAcc  >  2:298508/1‑251 (MQ=255)                                                                                                                                                                                                                                            cAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGAcc  <  1:21415/251‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     TGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGT‑TACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTAGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACC  >  NC_000913/13299‑13783

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