Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,307,151	G→A	P217S (CCT→TCT)	alsK ←	D‑allose kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,307,151	0	G	A	93.6%	125.1 / ‑3.3	47	P217S (CCT→TCT)	alsK	D‑allose kinase
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base A (21/23);  total (22/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCGGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCA‑CCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTT‑TCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAG‑CGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCC  >  NC_000913/4306904‑4307377                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                     cGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATACAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCA‑CCGCATCGGTATCGAGCAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGGTTTCAAGCAGACTCTGGACGGAAGAGGc                                                                                                                                                                                                                                   >  2:166870/1‑251 (MQ=255)  gTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCA‑CCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGAGGCg                                                                                                                                                                                                                                  >  1:196924/1‑251 (MQ=255)     gatacgcCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCA‑CCGCATAGGGATCGAACAGATTAGTGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAATCAGACTCTGGACGAAAGAGGCGttt                                                                                                                                                                                                                               >  1:481542/5‑251 (MQ=255)               tatTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCA‑CCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGAGGCGTTT‑TCCGCATgg                                                                                                                                                                                                                     >  1:200806/1‑250 (MQ=255)                atTGCTGCACCCTGAGCGCCATGAAAGTCAGAAGATGATGCAGCAATAAAGCGCTCGACCTGATGCGGTAGTGGACGGCGCAGGTTCTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCA‑CCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATTGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACTATAGATGCGTTT‑TCTGCATGGa                                                                                                                                                                                                                    <  2:410051/250‑1 (MQ=255)                     ggcACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCGGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCA‑CCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGAGGCGTTT‑TCCGCATGGACGaaaa                                                                                                                                                                                                              <  2:196924/250‑1 (MQ=255)                      gcACCCTGAGCGCCATTAAAGGCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGTCGGGCATATCCATCACGCCACCGCCCAGGATCA‑CCGCATCGGGATCGAACAGATTAATTCTGGTGGCAATGTCCCGTGCCGCTTTTTCAAGCAGTCGCTGGACGAAAGAGGCGTTT‑GCCTCATGGACGaaat                                                                                                                                                                                                              <  2:240642/250‑2 (MQ=255)                       cACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCAAGCCACCGCCCAGGATCA‑CCGCATCGGGATCGAACAGATTACTGCTGGTGGCAATGGCCCGAGCCGCGTTTTCAAGCAGACTCTGGGCGAAAGAGGCGTTT‑TCCGCATGGACgcaa                                                                                                                                                                                                               >  2:247397/1‑245 (MQ=255)                        aCCCTGAGCGCCATTAAAGTCAGAAGATGAGGCAGCACTACAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCA‑CCGCATCGGGATCGAACCGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGAGGCGTTT‑TCCGCCTGGCCGAAAAGa                                                                                                                                                                                                            >  2:410608/1‑250 (MQ=255)                                cgcCATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCA‑CCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGAGGCGTTT‑TCCGCATGGACGAAAAGATCGCGCAAt                                                                                                                                                                                                   <  1:247397/251‑1 (MQ=255)                                 cccATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGTAGTGGACGGCGCCGGTACTTTTGGGTCATGGCAACGAGAGTCCCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCA‑CCGCATCGGGATCGAACAGATTAATGCTGGTGGCATTGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGAGGCGTTT‑TCCGCATGGACGAAAAGATCGCGCAATg                                                                                                                                                                                                  <  2:416941/250‑1 (MQ=255)                                  ccATTAAAGTCAGAAGATGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACGGTTGGGTCATGGCAACGAGAGGCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCA‑CCGCATCGGGATCGAACAGATTAAGGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGATTCTGGACGAAAGAGGCGTTT‑TCCGCATGGACGAAAAGATCGCGCAATgg                                                                                                                                                                                                 <  1:170458/251‑1 (MQ=255)                                    aTTAAAGTCAGAAGATGAGGCAGCAATGAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCA‑CCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGAGGCGTTT‑TCCGCATGGACGAAAAGATCGCGCAATGGGt                                                                                                                                                                                               <  2:326975/251‑1 (MQ=255)                                             agaagaTGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTTGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCA‑CCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGAGGCGTTT‑TCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCggg                                                                                                                                                                                      >  1:182598/1‑251 (MQ=255)                                                  aTGAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCATGTCCTTTTGGGTCATGGCAACGAGAGTCGCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCA‑CCGCATCGGGATCGAACAGATTAAGGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGAGGCGTTT‑TCCGCATGGACTAAAAGATCGCGCAATGGGTAATTTCGGGGCtgt                                                                                                                                                                                 <  1:341984/251‑1 (MQ=255)                                                   ggAGGCAGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCA‑CCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGAGGCGTTT‑TCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCtgtt                                                                                                                                                                                <  2:404852/250‑1 (MQ=255)                                                      ggcagcAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCA‑CCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCCAGCAGACTCTGGACGAAAGAGGCGTTT‑TCCGCATGGACGAAAAGATCGCGCAATGGGGAATTTCGGGGCTgttgtt                                                                                                                                                                             >  2:474298/5‑251 (MQ=255)                                                                 aGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCA‑CCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGAGGCGTTT‑TCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGt                                                                                                                                                                  <  1:315245/251‑1 (MQ=255)                                                                                  ggCAGTGGACGGCGCAGGGACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCA‑CCGCCTCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGAGGCGTTT‑TCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGGTCGTCCCAGCGTCTTAG‑GGCCCTTCCagg                                                                                                                                                >  2:90450/1‑250 (MQ=255)                                                                                    cAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCA‑CCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGAGGCGTTT‑TCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCCGCGTCTTAG‑CGCCATTCCAGAGc                                                                                                                                              >  1:176886/1‑251 (MQ=255)                                                                                        ggACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCA‑CCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGAGGCGTTT‑TCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAG‑CGCCATTCCAGAGCAAt                                                                                                                                           <  1:359353/250‑1 (MQ=255)                                                                                          aCGGCGCAGGTAATTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCAGCACGCCACCGCCCAGGATCA‑CCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGAGGCGTTT‑TCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAG‑CGCCATTCCAGAGCAATTgg                                                                                                                                        <  2:422568/251‑1 (MQ=255)                                                                                                                                 gcgTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCA‑CCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGAGGCGTTT‑TCCGCATGGACGaaaa                                                                                                                                                                                                              >  1:3231/1‑143 (MQ=255)                                                                                                                                 gcgTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCA‑CCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGAGGCGTTT‑TCCGCATGGACGaaaa                                                                                                                                                                                                              <  2:3231/143‑1 (MQ=255)                                                                                                                                   gTGGGAAGGGGGGCATCTCCATCCCGCCACCGCCCAGGA‑CAGCCGCATCGGGATCGAACAGATTAATGAGGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTATGGACGAAAGAGGCGTTT‑TCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAG‑CGCCATTCCAGAGCAAGTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTc                                                                                               <  2:178254/251‑1 (MQ=255)                                                                                                                                   gTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCA‑CCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGAGGCGTTT‑TCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAG‑CGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTc                                                                                               >  1:440572/1‑251 (MQ=255)                                                                                                                                   gTGGGAAGGCGGGCATAGCCATCTCGCCACCGCCCAGGATCA‑CCGCAGCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGAGGCGTTT‑TCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGTGCTGTTGTTCGTACCAGCGTCTTAG‑CGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTc                                                                                               <  2:178240/251‑1 (MQ=255)                                                                                                                                    tGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCA‑CCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGAGGCGTTT‑TCCGCATGGACGAAAAGATCGCGCAATg                                                                                                                                                                                                  >  1:97236/1‑152 (MQ=255)                                                                                                                                    tGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCA‑CCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGAGGCGTTT‑TCCGCATGGACGAAAAGATCGCGCAATg                                                                                                                                                                                                  <  2:97236/152‑1 (MQ=255)                                                                                                                                           gCGGGCATATCCATCACGCCACCGCCCAGGATCA‑CCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGAGGCGTTT‑TCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAG‑CGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTcc                                                                                       <  1:222645/251‑1 (MQ=255)                                                                                                                                                              ccccGCCCAGGATCA‑CCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGAGGCGTTT‑TCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAG‑CGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGt                                                                    <  2:11182/249‑1 (MQ=255)                                                                                                                                                               cccGCCCAGGATCA‑CCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGAGGCGTTT‑TCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAG‑CGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGtt                                                                   <  1:166870/250‑1 (MQ=255)                                                                                                                                                                 cGCCCAGGATCA‑CCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGTACGAAAGAGGCGTTT‑TCCGCATGGACGAAAAGATCGCGCAATGGGTAAGTTCGGGGCTGTTGTTCGTACCAGCGTCTTAG‑CGCCATTCCAGAGCAATTGGTTTCCAGGCTGCCAGGAGTGCCACACGCGCAGTGGTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCg                                                                 <  2:240272/251‑1 (MQ=255)                                                                                                                                                                     cAGGATCA‑CCGCATCGGGATCGAACATATTAAGGCTGGTGGCAAGGGCCCGTGCCGCGGTTTCAAGCAGACTCTGGACGAAAGAGGCGTTT‑TCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAT‑CGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTg                                                             <  2:84631/251‑1 (MQ=255)                                                                                                                                                                          tCA‑CCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTT‑TCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAG‑CGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCCGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCcac                                                         >  2:310114/1‑250 (MQ=255)                                                                                                                                                                            a‑CCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGAGGCGTTT‑TCCGCATGGCCGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAG‑CGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCAACAACCGCAGTGTTGGTTCATATCTCCCAAGGGGATATGGCACAGTTCCCCTGCCACAcc                                                      >  2:207069/1‑251 (MQ=255)                                                                                                                                                                            a‑CCGCATCGAGATCGAACAGATTAATGCTGGTGGCAATGGCCCGCGCAGCGGTTTAAAGCAGACTCTGGCCGAAAGAGGCGTTT‑TCCGCATGGACGAAACGATCGCGCAATGGGGAATTTCGGGACTGTTGTTAGGACCAGAGTATTAG‑CGCCATTCCAGAACAGTTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGGCTTGGGACATCTCTCCCAAGGGGATATAACCCAGTTAGCCTGCCAAAcc                                                      >  2:32127/1‑251 (MQ=255)                                                                                                                                                                                    cGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGAGGCGTTT‑TCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAG‑CGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACCCCGTGTGcc                                               >  2:124909/1‑250 (MQ=255)                                                                                                                                                                                              cAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGAGGCGTTT‑TCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAG‑CGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACg                                     >  1:174786/1‑251 (MQ=255)                                                                                                                                                                                                               cAATGGCCCGTTCCGCTTTTCCAACCAGCCTCGGGAGAAAAGGGGCGTTT‑TCCGCATGGAGAAAAAGATCGCGCAATGGTTAATTTCGGGGCTGTTGTTCGTACCAGGGTTTTGG‑CGCCATTCCAGAGCAATTGGTTCCCAGGTACCCAGTTTTGCCCCACGCTCAGTGTTGGTTCTTATCTCCCAGGGGGTTATGGCCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCcac                    <  1:65299/251‑1 (MQ=255)                                                                                                                                                                                                               cAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGAGGCGTTT‑TCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAG‑CGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCCGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCCCCGTTCATCcac                    >  2:161574/1‑251 (MQ=255)                                                                                                                                                                                                                agtgacccGTGCCGCGTTTTAAAGCAGACCTTGGACGAAATAGGTGTTT‑TCCGAATGGACTAAAAGACCGCGGAATGGTTTATTTCGGGGCTTTTTTTCGTCCCAGCGTCTTAG‑CTCCATTCCATGGCAATGGTTTTCCAGGCCCCGATGTTTGCCAGACGCGCTGTTTTGGGTAATATCTCCCAGGGGGACTTGCCCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCcaca                   <  1:317591/246‑1 (MQ=255)                                                                                                                                                                                                                     cccGTGCCGTGTTTTCAAGCAGACT‑TGGGAGAAAAGGGGGTTTGTCCGCGTGAACAAAAAGTCCGCGAAATGGTTATTTTGGGGGTTGTTGTCCGACCCAGCG‑CTTAGGCTCCTTTCCAGAGCTATTGGTTCCCAGGCCCCCAGTATTGCCACCCGCGCGGTGTTGGTTAATACCCCCCGGGGGGATATGACCCATTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCg              <  2:174786/251‑1 (MQ=255)                                                                                                                                                                                                                            cgcgTTTTCAAGCAGACTCTGGACGAAAGAGGCGTTT‑TCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAG‑CGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGAcc                                                                        <  1:448908/186‑1 (MQ=255)                                                                                                                                                                                                                            cgcgTTTTCAAGCAGACTCTGGACGAAAGAGGCGTTT‑TCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAG‑CGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGAcc                                                                        >  2:448908/1‑186 (MQ=255)                                                                                                                                                                                                                               gTTTTCAAGCAGACTCTGGACGAAAGAGGCGTTT‑TCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAG‑CGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGACCCCCATcc    >  1:305261/1‑251 (MQ=255)                                                                                                                                                                                                                                 tttCAAGCAGACTCTGGACGAAAGAGGCGTTT‑TCCGCATGGACGAAAAGATCGCGCAATGGGAAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAG‑CGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATcccc  >  1:103320/1‑251 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                     CGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCGGCAATAAAGCGCACGACCTGATGCGGCAGTGGACGGCGCAGGTACTTTTGGGTCATGGCAACGAGAGTCTCGCGTGGGAAGGCGGGCATATCCATCACGCCACCGCCCAGGATCA‑CCGCATCGGGATCGAACAGATTAATGCTGGTGGCAATGGCCCGTGCCGCGTTTTCAAGCAGACTCTGGACGAAAGGGGCGTTT‑TCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAG‑CGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCC  >  NC_000913/4306904‑4307377

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