Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,830,800	C→T	A115V (GCG→GTG)	yicH →	AsmA family protein YicH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,830,800	0	C	T	97.5%	258.7 / NA	82	A115V (GCG→GTG)	yicH	AsmA family protein YicH
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (47/32);  minor base G (0/2);  total (48/34)
Fisher's exact test for biased strand distribution p-value = 1.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACATATCAGCGCATGGGTTTCCGAGAATAGCGACTATCATCTGGCCTTCGGGGCGATGGATCACCGTTTTTCCGCGCCATCTCATATCGTGCTGGAGAACGTCACGTTTGGTCGTGATGGCCAGCCCGCGACCCTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGCGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAG‑TGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTG  >  NC_000913/3830554‑3831044                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      aaCAAATCAGCGCATGGGTTTCCGAGATTAGCGACTATCATCTGGCCTTCGGGGCGATGGATCACCGTTTTTCCGCGCCATCTCATCTCGTGCTGGAGAACGTCGCGTTTGGGCGGGAGGGCGAGCCCGCGACCCTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTTGGCAACTGACCGAACCACGCCATGTCGAGACCATCCGGCTGGAAAACAAGACGCTGAATCTCACCGACCAGACCGtgccg                                                                                                                                                                                                                                                   <  2:52036/251‑1 (MQ=255)    atatCAGCGCATGGGTTTCCGAGAATAGCGACTATCATCTGGCCTTCGGGGCGATGGATCACCGTTTTTCCGCGCCATCTCATATCGTGCTGGAGAACGTCACGTTTGGTCGTGATGGCCAGCCCGCGACCCTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGCAAACGGGACGCTGAATCTCACCGACCCGACCGTGCCGCtc                                                                                                                                                                                                                                                >  1:393916/1‑250 (MQ=255)     tatCAGCGCATGGGTTTCCGAGAATAGCGACTGTAATCTGGCATTCGGGGCGATGGATCACCGTTTTTCCGCGCCATCTCATATCGTGCTGGAGAACGCCACGTTTGGTCGTGATGGCCAGCCCGCGACCCTGGTGGCCACAAGTGTCGCAATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCCTGTCGATACCATCCTGCTGGAAAACGGGACGATGACTCTCACCGACCAGACCGTGCCGCTAc                                                                                                                                                                                                                                               >  1:364400/1‑251 (MQ=255)            gcATGGGTTTCCGAGAATAGCGCCTATAATCTGGCCTTAGGGACGATGGATAACCGTTTTTCCGAGACATCTAATCTCGTGCTGGAGACCGTCACGTTTGGTCGTGATGGCAAGCCCGCGACCCTGGTGGCCAACAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGCCCGAACCCCGCCATGTCGAAACCATCCTGCTGGAAAACGGGACGCTGACTCCCACCGACCAGCACGTGCCGCTACCGTTCa                                                                                                                                                                                                                                         >  2:49059/1‑250 (MQ=255)                       cGAGAATAGCGACTATCATCTGGCCTTCGGGGCGATGGATCACCGTTTTTCCGCGCCATCTCATATCGTGCTGGAGAACGTCACGTTTGGTCGTGATGGCCAGCCCGCGACCCTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGTGCCGCTACCGTTCAAAGCCGATCGTc                                                                                                                                                                                                                             <  1:179596/251‑1 (MQ=255)                       cGAGAAAAGCGACTATCATCGGGTCTGCGGGGCGATGGATCACCGTTTTTCCGCGCCATCTCATATCGTTCTGGAGAACGTCACGTTTGGTCGTGATGGCCAGCCCGCGACCCTGTTGGCCAAAAGTGTCGACATTGCGCTAAGCAGGCGGCAACGGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCCTGCCGCTACCGTTCAAACCCGCTCGTc                                                                                                                                                                                                                             <  2:237041/251‑1 (MQ=255)                           aaTAGCGACTATCATCTGGCCTTCGGGGCGATGGATCACCGTTTTTCCGCGCCATCTCATATCGTGCTGGAGAACGTCACGTTTGGTCGTGATGGCCAGCCCGCGACCCTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACAACGCCATGTCGATACCATCCTGCTGGAAAAAGGGACGCTGAATCTCACCGACCAGACCGTGCCGCTACCGTTCAAAGCCGATCGTCTGc                                                                                                                                                                                                                          >  1:211957/1‑250 (MQ=255)                           aaTAGCGACTATCATCTGGCCTTCGGGGCGATGGATCACCGTTTTTCCGCGCCATCTCATATCGTGCTGGAGAACGTAACGTTTGGTCGTGATGGCCAGCCCGCGACCCTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAAAGGGACGCTGAATCTCACCGACCAGACCGTGCCGCTACCGTTCAAAGCAGATCGTCTGc                                                                                                                                                                                                                          >  1:376463/1‑250 (MQ=255)                              aGCGACTATCATCTGGCCTTCGGGGCGATGGATCACCGTTTTTCCGCGCCATCTCATATCGTGCTGGAGAACGTCACGTTTGGTCGTGATGGCCAGCCCGCGACCCTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGTGCCGCTACCGTTCAAAGCCGATCGTCTGCAACt                                                                                                                                                                                                                      >  1:299866/1‑251 (MQ=255)                                cGACTATCATCTGGCCTTCGGGGCGATGGATCACCGTTTTTCCGCGCCATCTCATATCGTGCTGGAGAACGTCACGTTTGGTCGTGATGGCCAGCCCGCGACCCTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGTGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGc                                                                                                                                                                                                                    >  1:332094/1‑251 (MQ=255)                                  aCTATCATCTGGCCTTCGGGGCGGTGGATCACCGCTTTTCCGCGCCATCTCATATCGTGCTGGAGAACGTCACGTTTGGTCGTGATGGCCAGCCCGCGACCCTGGTGGCCAAAAGAGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGACGATACCATCCTGCTGGAAAACAGGACGCTGAATCGCACCGACCAGACCGTGCCGCTCCCGGTGAAAGCCGATCGGCTACAACTGCgg                                                                                                                                                                                                                  >  2:290870/1‑250 (MQ=255)                                        atcTGGCCTTCGGGGCGATGGATCACCGTTTTTCCGCGCCATCTCATATCGTGCTGGAGAACGTCACGTTTGGTCGTGATGGCCAGCCCGCGACCCTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGTGCAGCTACAGTTCAAAGCCGAGCGTCTGCAACTGCGTGATATg                                                                                                                                                                                                            >  1:365591/1‑251 (MQ=255)                                           tGGCCTTCGGGGCGATGGATCACCGTTTTTCCGCGCCCTCTCATATCGTGCTGGAGAACGTCACGTTTGGTCGTGATGGCCAGCCCGCGACCCTGGTGGCCAAAAGTGTCGCCATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGCCGCTGAATCTCACCGACCAGACCGTCCCGCTACCGTTCAAAGCCGCTCGTCTGCAACTGCGTGATATGGCg                                                                                                                                                                                                         >  2:139892/1‑251 (MQ=255)                                                    ggtcGATGGATCCCCGTTGTTCCGCGCCATCTCATAGCGGGCTGGAGAACGTCACGTTTGGTCGTGATGGCCAGCCCGCGACCCTGGTGGCCAAAAGTCACGACATTGCGCGAAGCAGGCGGCAACTGACCGAACCACGCCATGTCGATAACATCCTGCTGGAAAACGGGACGCAGAATCTCACCGACCAGACCGTGCCGCTACCGTTCAAAGTCGATCGTCTGCAACTGCGTGATATGGCGGTTAATAGc                                                                                                                                                                                                <  2:285253/248‑1 (MQ=255)                                                               cACCGTTTTTCCGCGCCATCTCATATCGTGCTGGAGAACGTCACGTTTGGTCGTGATGGCCAGCCCGCGACCCTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGTGCCGCTACCGTTCAAAGCCGACAGTCTGCAACTGCGTGATATGGCGTTTAATAGCCAGAATAGCgc                                                                                                                                                                                     >  2:20031/1‑250 (MQ=255)                                                                       ttCCGCGCCATCTCATAGCGTGCTGGAGAACGTCACGTTTGGTCGTGATGGCCAGCCCGCGACCCTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCGGAATCTCACCGACCAGACCGTGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGaaa                                                                                                                                                                              <  1:253326/250‑1 (MQ=255)                                                                       gtCCGCGCCAGCTCATATCGTACTGGAGAACGGCTAGTTTGGTCTTGATGGCCAGCCCGCGACCCTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTTACCGAACCACGCCATGTCGATACCATCCTGCTAGAAAACGGGACGCTGAATCTCACCGACCAGACCGTGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGTCGTTTAATAGCCCGAATAGCGAATGGaaa                                                                                                                                                                              <  2:190455/249‑1 (MQ=255)                                                                         ccGCGCCATCTCATATCGTGCTGGAGAACGTCACGTTTGGTCGTGATGGACAGCCCGCGACCCTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGTGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACt                                                                                                                                                                            <  1:417999/250‑1 (MQ=255)                                                                              ccATCTCATATCGTGCTGGAGAACGTCACGTTTGGTCGTGATGGCCAGCCCGCGACCCTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGTGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAgcgc                                                                                                                                                                      >  1:42306/1‑251 (MQ=255)                                                                              ccATCTCATATCGTGCTGGAGAACGTCACGTTTGGTCGTGATGGCCAGCCCGCGACCCTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGTGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAgcgc                                                                                                                                                                      <  2:90347/251‑1 (MQ=255)                                                                                     atatCGTGCTGGAGAACGTCACGTTTGGTCGTGATGGCCAGCCCGCGACCCTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGTGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCgg                                                                                                                                                               >  1:262354/1‑251 (MQ=255)                                                                                          gTGCTGGCGAACGTCACGTTTGGACGTGATGGCCAGCCCGCGACCCTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCCGACCGTGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGTGCGCGCAGCGGGTaa                                                                                                                                                           >  2:415406/1‑250 (MQ=255)                                                                                          gTGCTGGAGAACGTCACGTTTGGTCGTGATGGCCAGCCCGCGACCCTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCTATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGTGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTaaa                                                                                                                                                          <  1:274914/251‑1 (MQ=255)                                                                                          gTGCTGGAGAACGTAACGTTTGGTCGTGATGGCCAGCCCGCGACCCTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGTGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCggggaa                                                                                                                                                           >  2:145925/1‑247 (MQ=255)                                                                                          gTGCGGGAGAACGTCACGTTTGGTCGTGATGGCCAGCCCGCGACCCTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGTGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTaaa                                                                                                                                                          <  1:335203/251‑1 (MQ=255)                                                                                            gCTGTAGAACGTCACGTTTGGTCGTGATGGCCAGACCGCGATCCTGGTGGCCAAAAGTGACGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGAGGCTGAATCTCACCGTCCAGACCGTGCCGCTACCGGTCAAAGCCGATCGTCTGCAACTGCATGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGATCGCGCAGCGGGTAAATg                                                                                                                                                        <  1:217065/251‑1 (MQ=255)                                                                                                   aaCGTCACGATTGCTCGTGCTGGCCCGCAAGAGACCATGGTGGCAAACAGTGTCGACATTGCGCTAAGCAGTCGGCACCTGACCGAACCACGCCATCTCGATACCATCCTGCTGGCAAACGGGACGCTGAATCTCACCGACCAGACCGAGCAGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGCTATGGAGGTTACCAGCCAGAATAACGCATGGAACCTACGCGCGCAGCGGGTAAATGacggcgg                                                                                                                                                 >  2:240595/1‑250 (MQ=255)                                                                                                           gTTTGGTCGTGATGGCCAGCCCGCGACCCTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGTGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGt                                                                                                                                         <  1:116555/251‑1 (MQ=255)                                                                                                              ttgTCGTGATGGCCAGCCCGCGACCCTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGTGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTgg                                                                                                                                       <  2:314878/248‑1 (MQ=255)                                                                                                              tGGTCGTGATGGCCAGCCCGCGACCCTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGATCCACGCCATGTCGATACCATCCTGCTTGAAAACGGGACGCTGAATCTCTCCTACCAGACCGTTCCGCTACCGTTCAAAGCCGATCGGCTGCAACTGCGTTATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTgg                                                                                                                                       <  1:182437/250‑1 (MQ=255)                                                                                                               ggTCATGATGGCCAGCCCGCGACCCTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGTGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTc                                                                                                                                     >  1:2680/1‑251 (MQ=255)                                                                                                                  cGTGATGGCCAGCCCGCGACCCTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGTGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCAc                                                                                                                                   <  1:417109/250‑1 (MQ=255)                                                                                                                    tGATGGCCAGCCCGCGACCCTGCTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCAACCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGTGCCGCTACCGATCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTACGTGGTCACCa                                                                                                                                 >  2:196528/1‑250 (MQ=255)                                                                                                                      aTGGCCAGCCCGCGACCCTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGTGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCaaaa                                                                                                                              <  2:82351/251‑1 (MQ=255)                                                                                                                                   gACCCTGGTGGCCAAAAGTGTCGACATTTCGCTAAGCAGTCGGCAACTGACCGACCCACTCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCACGCCGACCAGACCGTGCCGCTACCGTTCAAAGCCGATCGTCTGCACCTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGTTCACCAAAAGCCGGTAAAG‑TGc                                                                                                                <  1:353266/251‑1 (MQ=255)                                                                                                                                   gACCCTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGTGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAG‑TGc                                                                                                                >  1:91109/1‑251 (MQ=255)                                                                                                                                     cccTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCACCTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGTGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGACTGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAG‑TGCTg                                                                                                              >  2:218480/1‑251 (MQ=255)                                                                                                                                          gtggCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGTGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAG‑TGCTGGGTAc                                                                                                         >  1:235973/1‑251 (MQ=255)                                                                                                                                              ccAAAAGTGGCGACATTGCGCTAAGCAGGCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGTGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAG‑TGCTGGGTACGAAg                                                                                                     <  1:186630/251‑1 (MQ=255)                                                                                                                                               cAAAAGTTTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGTGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTACTAGCCCGAATAGCGAATGGAAACTGCGCGCGCAGCGGGAAAAGGGCGGCGTAGGTCCGTGGTCACCAAAAGCCGGTAAAG‑TGCTGGGTACGCAgg                                                                                                    >  2:229012/1‑251 (MQ=255)                                                                                                                                                 aaaGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGATCCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGTGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCTGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAG‑TGCTGGGTACGAAGgcg                                                                                                  <  2:244604/251‑1 (MQ=255)                                                                                                                                                 aaaGTGTAGACATTGCGATAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACACTGAATCTCACCAACCAGCCCGTGCCGCTACCGCTCAAAGCCGATCGCCTGCACCTGCGTGATATGGAGTTTACTAGACCGAGTAGCGAATGGAAACTGAGCGCTCAGCGGGTAAATGGCGGCGTAGTTCCGGGGTCCCCAAAATCCGTTCAAG‑TGCTGGGTACGCAGgcg                                                                                                  >  2:29126/1‑251 (MQ=255)                                                                                                                                                    gtgccGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGTGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCACGAATAGCGAATGGAACCTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAG‑TACTGGGTACGAAGgcgc                                                                                                 >  2:254090/5‑249 (MQ=255)                                                                                                                                                      gtCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGTGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAG‑TGCTGGGTACGAAGgcgc                                                                                                 >  1:223184/1‑247 (MQ=255)                                                                                                                                                                 gcTAAGCAGTAGGCAACTGCCCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGTGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAG‑TGCTGGGTACGAAGGCGCAGATTCAGTTTAGTg                                                                                  >  1:324396/1‑251 (MQ=255)                                                                                                                                                                    aaGCAGGCGGCAACTGACCGATCCACGCCATGTCGATTCCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGTTGCGATACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATTGAAGCTGAGCGCGCAGCGGTTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAG‑TGCTGGGTACGAAGCCGCAGATTCAGTTTAGTGCCg                                                                               <  2:282693/251‑1 (MQ=255)                                                                                                                                                                         gTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGTGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGAGATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGGTCCGTGGTCACCAAAAGCCGGTAAAG‑TGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATCg                                                                          <  2:360796/251‑1 (MQ=255)                                                                                                                                                                            ggCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGTGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAG‑TGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCtt                                                                       <  1:332140/251‑1 (MQ=255)                                                                                                                                                                             gCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGTGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAG‑TGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCttt                                                                      >  1:31426/1‑251 (MQ=255)                                                                                                                                                                             gCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGTGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAG‑TGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGACGGATCGCttt                                                                      >  1:206914/1‑251 (MQ=255)                                                                                                                                                                             gCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGGGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAG‑TGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCttt                                                                      <  1:139892/251‑1 (MQ=255)                                                                                                                                                                               aaCTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGTGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAG‑TGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGAGCGCTTTCg                                                                    <  1:352719/251‑1 (MQ=255)                                                                                                                                                                                  tGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGTGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAACAGCCGGTAAAG‑TGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTc                                                                 >  1:215130/1‑251 (MQ=255)                                                                                                                                                                                   gACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGTGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAACGCCGGTAAAG‑TGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTc                                                                 >  1:96236/1‑250 (MQ=255)                                                                                                                                                                                   gACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGTGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAG‑TGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCa                                                                <  1:392963/251‑1 (MQ=255)                                                                                                                                                                                     ccGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGTGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTAACCAAAAGACGGTAAAG‑TGCTGGGTACGAAGGCGCAGATTCAGTGTAGTGCCGGATCGCTTTCGCTCAat                                                              >  2:153596/1‑251 (MQ=255)                                                                                                                                                                                      cGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGTGCCGCTACCGTTCAAAGCCGAGCGGCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAAGGGCGGCGTAGTGCCGTGGTCACCAAAAGCCGGTAAAG‑TGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAatg                                                             <  1:205048/251‑1 (MQ=255)                                                                                                                                                                                        aaCCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGTGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCCGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAG‑TGCTGGGTACGAAGGCGCAGATTCCGTTTAGTGCCGGATCGCTTTCGCTCAatgat                                                           >  1:26861/1‑251 (MQ=255)                                                                                                                                                                                           cACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGTGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAG‑TGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGtt                                                        >  1:407221/1‑251 (MQ=255)                                                                                                                                                                                            aCGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGTGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAG‑TGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGtgc                                                       <  2:203902/251‑3 (MQ=255)                                                                                                                                                                                             cGCCATGTAGCTAACATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGTGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATCGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAACGCCGGTAAAG‑TGCTGGGTCCGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTcc                                                      >  2:358744/1‑251 (MQ=255)                                                                                                                                                                                              cccATGTCGATACCATCCTTCTCGAAATCGTGACTCTGAATCTCACCGACCAGCCCGTGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTTCGTGATATGGCGGTTAATTGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCTGGTAAATGGCGGCGTGGTTCCGTGGTCTCCAAAAGCCTTTAAAG‑TTCTGGGTACGAATGCGCAGATTCAGTTTATTGCCGGATCGCTTTCGCTCactggtgtgcct                                                     <  2:365176/250‑12 (MQ=255)                                                                                                                                                                                                 acgTCGATACCATCCTGCTGGAAAACGATACTCTGTATCTCACCGACCAGCCCTTGCCGCTACCGTTCAAAGCCGGTCGGCTGCACCTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAATTGAGCGCGCAGGGGGTAAACGGCGGCGGAGGGCCGTGGTCTCCAATAGCCGGTAAAGGTG‑TGGGTACGAAGGCGCAGATTCAGGTCAGAGCCGGATCGCTTTCTCGCAATGATTTTCCTGcc                                                  <  2:257614/249‑1 (MQ=255)                                                                                                                                                                                                  tGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGTGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAG‑TGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGcc                                                  >  1:70786/1‑250 (MQ=255)                                                                                                                                                                                                      gATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGTGCCGCTACCGTTCAAAGCCGATCGTCTGCAGCTGCGTGATATGACGTTTAATAGCCCGAATAGCGAATGGAAAATGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAG‑TGATGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTACTGccacca                                              >  1:139341/1‑250 (MQ=255)                                                                                                                                                                                                           cATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGTGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAG‑TGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTAc                                        >  2:205452/1‑251 (MQ=255)                                                                                                                                                                                                            aTCCTGCTGGAAAACGGGCCGCTGAATCTCACCGACCAGACCGTGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGAGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAAATGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAG‑TGCTGGGAACGAAGGCGCAGATTCAGTTTAGTGCCGGAACGCTTTCGCTCAATGATGTTCCTGCCACCAATGtcat                                       >  1:74916/1‑248 (MQ=255)                                                                                                                                                                                                                   tgGAAAACGGGACGCTGAATCTCACCGACCAGACCGTGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAG‑TGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGa                                 >  1:287088/1‑250 (MQ=255)                                                                                                                                                                                                                    ggAAAACGGGACGCTTAATCTCACCGACCAGACCGGGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGGTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGCGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAG‑TGCTGGGTACGAATGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGa                                 <  1:35986/249‑1 (MQ=255)                                                                                                                                                                                                                        aaCGGGCCGCTGAATCTCACAGACCAGACCGTGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAG‑TGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGCTTGAAGGc                             >  1:336101/1‑249 (MQ=255)                                                                                                                                                                                                                          ctggACGCTGAATCTCTCCGACCAGACCGTGCCGATCCCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGCATGGAAACTGAGCGCGCCGCGGGAAAATGGCGGCGTAGTTCCGGGGTCACCAAAAGCCGGTAAAG‑AGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCGTTCGCTCCATGATGTTACTACAACCAATGTACTGATTGCAGGCAGt                          >  2:27751/3‑250 (MQ=255)                                                                                                                                                                                                                           gggACGCTGAATCTCACCGACCAGACCGTGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAG‑TGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTAt                        <  1:118397/251‑1 (MQ=255)                                                                                                                                                                                                                           gggACGCTGAATCTCACCGACCAGACCGTGCCGCTACCGTTCAAAGCCGATCGTCAGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGCGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAG‑TGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGTTTGAAGGCAGTAt                        >  1:269337/1‑251 (MQ=255)                                                                                                                                                                                                                                 cTGAATCTCACCGACCAGACCGTGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGACTAGCGCATGGAAACTGAGCGCGCCGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAG‑TGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGTTGTTCCTGCCACCATTGTACTGATTGAGGGCAGTATTgtaa                   >  2:42989/1‑247 (MQ=255)                                                                                                                                                                                                                                 cTGAATCTCACCGACCAGACCGTGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAG‑TGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATa                   >  2:128532/1‑250 (MQ=255)                                                                                                                                                                                                                                   gaacctccccGACCAGACCGTGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCTTGATATGGCGTTTAATATCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGTGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGTTAAAG‑TGCTGGGTACTAAGGCTCAGATTCAGTATAGTGCCGGATCGCTTTCGCTCATTGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTAATGATAAc                 <  2:61432/242‑1 (MQ=255)                                                                                                                                                                                                                                    aaTCTCACCGACCAGACCGTGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAG‑TGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACg                >  1:122222/1‑250 (MQ=255)                                                                                                                                                                                                                                          aCCGACCAGACCGTGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAG‑TGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGt         >  1:386630/1‑251 (MQ=255)                                                                                                                                                                                                                                          aCCGACCAGACCGTGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACCGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAG‑TGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGt         >  1:31719/1‑251 (MQ=255)                                                                                                                                                                                                                                            cGACCAGACCGTGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAG‑TGCTGGGTACGAAGGCGCAGATTCAGTTTATTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGtt        >  1:250561/1‑250 (MQ=255)                                                                                                                                                                                                                                             gACCAGACCGTGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAG‑TGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTAc      <  2:324396/251‑1 (MQ=255)                                                                                                                                                                                                                                             gACCAGACCGTGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAG‑TGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTAc      <  2:299866/251‑1 (MQ=255)                                                                                                                                                                                                                                                 aGACCGTGCCGCTTCCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAG‑TGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTg  <  2:91109/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      AACATATCAGCGCATGGGTTTCCGAGAATAGCGACTATCATCTGGCCTTCGGGGCGATGGATCACCGTTTTTCCGCGCCATCTCATATCGTGCTGGAGAACGTCACGTTTGGTCGTGATGGCCAGCCCGCGACCCTGGTGGCCAAAAGTGTCGACATTGCGCTAAGCAGTCGGCAACTGACCGAACCACGCCATGTCGATACCATCCTGCTGGAAAACGGGACGCTGAATCTCACCGACCAGACCGCGCCGCTACCGTTCAAAGCCGATCGTCTGCAACTGCGTGATATGGCGTTTAATAGCCCGAATAGCGAATGGAAACTGAGCGCGCAGCGGGTAAATGGCGGCGTAGTTCCGTGGTCACCAAAAGCCGGTAAAG‑TGCTGGGTACGAAGGCGCAGATTCAGTTTAGTGCCGGATCGCTTTCGCTCAATGATGTTCCTGCCACCAATGTACTGATTGAAGGCAGTATTGATAACGATCGCGTTACGCTG  >  NC_000913/3830554‑3831044

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