Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,363	G→A	S210S (TCG→TCA)	thrC →	threonine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,363	0	G	A	98.4%	185.6 / ‑4.4	62	S210S (TCG→TCA)	thrC	threonine synthase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (31/30);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.43e-01

GACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAACTGAAAGTGGCGCTAGGGTTAAACTCGGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCTGTTGCGCAGCTGCCGCAGGA‑GACGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACGTTTCCTGCACGACGGTCAGTGGTCACCC  >  NC_000913/4123‑4603                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                  gACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGTGCTGGGGTTAAACTCAGCTAACt                                                                                                                                                                                                                                            <  1:129398/248‑1 (MQ=255) gACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTc                                                                                                                                                                                                                                           <  1:457014/249‑1 (MQ=255)    cgcgACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTa                                                                                                                                                                                                                                      <  2:236895/251‑1 (MQ=255)      ctaCCTCCGGTGATACCGGAGCGGCAGTGGCTCCTGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAAc                                                                                                                                                                                                                                    >  2:238573/3‑251 (MQ=255)       gACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGCATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAACTCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGACATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAACGAGGCGCTGGGGTTAAACTCAGCTAACTCGATTAAc                                                                                                                                                                                                                                    >  1:303860/1‑250 (MQ=255)       gACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAAc                                                                                                                                                                                                                                    >  2:337115/1‑250 (MQ=255)       gACCTCCGGAGATACCGGAACGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGAGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGAGCAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAAATCAGCTAACTCGATTAAc                                                                                                                                                                                                                                    >  2:370704/1‑250 (MQ=255)            ccGGTGGTACCGTCGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGTGTGTTATGATGAAGAGCTGAAAGTTGCGCGGGTGTTAAACTCAGCTAACTCGATTAACATCAgt                                                                                                                                                                                                                              <  2:355499/251‑2 (MQ=255)            ccGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCa                                                                                                                                                                                                                                <  1:232986/249‑1 (MQ=255)             cGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTTCACTGCAAGAAAAACTGTTCTGTACATTTGGCGGCAATATCGAAACTGTTGCCATCGACGGCTATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGcc                                                                                                                                                                                                                             <  1:388515/251‑1 (MQ=255)                       ggAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGgc                                                                                                                                                                                                                   <  2:160289/251‑1 (MQ=255)                        gAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGGGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGAGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACCTCAGCCGTTTGCTGgcg                                                                                                                                                                                                                  >  2:290014/1‑251 (MQ=255)                           cGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAAATCGATTAACATCAGCCGTTTGCTGgcgc                                                                                                                                                                                                                 >  1:7306/1‑249 (MQ=255)                                 tGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGc                                                                                                                                                                                                         >  1:627517/1‑251 (MQ=255)                                                aCGGTTAACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACAGTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAg                                                                                                                                                                                            <  2:69375/249‑1 (MQ=255)                                                  ggTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGcagt                                                                                                                                                                                        <  1:372365/251‑4 (MQ=255)                                                       aCCGAATGTGAAAGTGGTTATTCTTTATTCACGAGGCAAAATCAGTACACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAATTGTCGCCATCGACGGCGATTTCGATGCCTGTCAGGCGATGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTGGCTACTACTTTGAAGCAGTTgcg                                                                                                                                                                                    <  2:277547/250‑1 (MQ=255)                                                       aCCGAAGGTGAAAGTGGTTATCCTCTATCATCGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTgcgc                                                                                                                                                                                   <  2:8070/251‑1 (MQ=255)                                                        ccGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCa                                                                                                                                                                                  <  1:456669/251‑1 (MQ=255)                                                         cGAATGTGACAGTGGTTATCCTCTATCCCCGATGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGAAATATCGAAACTGTTGCCATCGACGGAGAATTCGATGCCTGTCAGGCGCTGGAGGCGCTAGCGGTTGATGATGAAGAGCTGAAAATGGCGCTGGGGTTAAACTTAGCTAACTCGATTAACATCAGCCGATTGCTCGCGCAGATTTGCaactacattaaggcagttgcaccg                                                                                                                                                                                 >  2:480703/1‑233 (MQ=255)                                                            aTGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTg                                                                                                                                                                              >  1:23299/1‑251 (MQ=255)                                                            aTGTGAAAGTGGTTATCATATATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTCCATTGGGCGGCAATATCGAAACTGTTGCCATCGCCGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAGCATCAGACGTTTGCTGGCGCCGATTTTCTACTACTTTGAAGAAGTTGCGCAGcgg                                                                                                                                                                              >  2:520112/1‑249 (MQ=255)                                                             tgtgAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGAAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGCAGCAGTTGCGCAGCTGc                                                                                                                                                                             >  2:412990/1‑251 (MQ=255)                                                                  aaGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAg                                                                                                                                                                        <  1:414830/251‑1 (MQ=255)                                                                   aGTGGTTAACCTCTATCCACGAGGCCAAACCAGTCGAATGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAACCTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGTTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGGGGCGATGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAGGCCGATGCGCAGCTACCGAAgg                                                                                                                                                                       >  2:379601/1‑251 (MQ=255)                                                                          aTCCTCTATCCACGAGGCTAAATAAGTTCACTTCAAGAAAAACTTTACTGTCCATTGGTAGTCAATATCGAAACTGTTGCCATCGACGGCGATTTCTATGCCTGTCAGTCGCTGGTGAAGCAGGCGTTTGATGATGAAGAACTGAAAGTGGCGCTAGGGTTAATCTCGGCTAACTCGATTAACATCAGCGGGATGCTGGCTCAGATTTCCTACTACTTTGAAGCTGTTTCGCAGCTGCCGCAGGT‑GActcg                                                                                                                                                               <  2:409211/251‑4 (MQ=255)                                                                               ctATCCACGAGGCAAAATCATTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCTATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTGGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCaa                                                                                                                                                            <  2:585389/249‑1 (MQ=255)                                                                               ctATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAAcc                                                                                                                                                          >  2:213540/1‑251 (MQ=255)                                                                                          gCAAAATCAGTCCACTGCAAGAAAAACTTTTCTGTACATTGGGCTGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGGAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAATCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGtc                                                                                                                                               <  2:357277/251‑1 (MQ=255)                                                                                          gCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACGTCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGtc                                                                                                                                               >  1:32305/1‑251 (MQ=255)                                                                                             aaagcAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCg                                                                                                                                            <  1:604965/247‑1 (MQ=255)                                                                                              aaTCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGACACTGATGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCgg                                                                                                                                           >  2:149560/1‑251 (MQ=255)                                                                                               aTCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCCGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCggg                                                                                                                                          >  2:431329/1‑250 (MQ=255)                                                                                                      cACTGCAAGAAAAACTGTTCTGTTCAGTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGAGGCCTGTCAGGCGCTGTTGAAGCAGGCGTTTGATGATGTAGAGCTGAAAGTGGCGCTGGTGTTAAACGCAGCTAACTCGATGAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAATCAGTTGCCCAGCTGTCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGc                                                                                                                                   <  2:277319/251‑1 (MQ=255)                                                                                                          gCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCATCCAGCTGGTTGTCTCGGTGCCAAGCGGaa                                                                                                                               <  2:526402/251‑1 (MQ=255)                                                                                                           cAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGaa                                                                                                                               >  1:419795/1‑250 (MQ=255)                                                                                                              gAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACtt                                                                                                                           >  2:330529/1‑251 (MQ=255)                                                                                                                 aaaCTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAATCGGAATCTTCgg                                                                                                                        <  1:499889/251‑1 (MQ=255)                                                                                                                 aaaCTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCgg                                                                                                                        >  2:105044/1‑251 (MQ=255)                                                                                                                       ttCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGAttt                                                                                                                  >  1:315340/1‑251 (MQ=255)                                                                                                                         cTGTCCATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTg                                                                                                                 >  2:129887/1‑250 (MQ=255)                                                                                                                                 tGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCgggg                                                                                                        >  2:511917/1‑250 (MQ=255)                                                                                                                                         aaTATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTggg                                                                                                >  1:15404/1‑250 (MQ=255)                                                                                                                                               gAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCggaggc                                                                                          >  2:174817/1‑246 (MQ=255)                                                                                                                                                     gtggcctacGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTGGATGATGAAGAACTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGAGTTGCTGCTACTGTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACGTCGGCGAGTTGACGGCGGGTCTGCTGGCGATGTCACTc                                                                                      <  2:547897/241‑1 (MQ=255)                                                                                                                                                     gTTGCCATCGACGGCTATTTAGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACGCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACGACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGGCTGCTGGCGAAGTCACTCg                                                                                     <  2:451502/250‑1 (MQ=255)                                                                                                                                                        gtcaacGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGGTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCACCCAGCTGGTTGTCTCGGTGCCAAGCGGAACCTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCt                                                                                 <  2:290957/246‑1 (MQ=255)                                                                                                                                                        gCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCgggggtgccggcggagtcccccgggct                                                                                 >  2:231244/1‑227 (MQ=255)                                                                                                                                                              gACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGt                                                                           <  1:171426/251‑1 (MQ=255)                                                                                                                                                                    gATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGGCACCCGGTCTGCCGGGGAAACg                                                                     >  1:172013/1‑251 (MQ=255)                                                                                                                                                                        tCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGCGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGCAGTCACTCCGTCTGCCGGTGaaaaggttt                                                                 >  2:14570/1‑245 (MQ=255)                                                                                                                                                                         cGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGtttt                                                                 <  1:529948/250‑1 (MQ=255)                                                                                                                                                                                           ggTGAAGCCGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTTTTTCTGCGGCCCACGTgg                                              >  1:200014/1‑250 (MQ=255)                                                                                                                                                                                              gAAGCAGGTGGTTGAAGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAATGCAGCTAACTCGATTAACCTCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGatcg                                           <  2:233944/251‑4 (MQ=255)                                                                                                                                                                                                            atgatgaAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCCGGCGAAAGCACACGGTCGGCCGGTGGAACGTTTTATTGCTGCGACCCAAGGGGACGCAACCGTGCCACGt                             >  1:404046/1‑251 (MQ=255)                                                                                                                                                                                                                tgaAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGCAGTCACCCGGTCTGCCGGGGTAACGTGTTTTTTCTGCGCCCAACGTGAACGAGAACGTGCCACggtttc                         >  2:506223/1‑246 (MQ=255)                                                                                                                                                                                                                  aagagCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAATTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACGTTTCCTg                       <  1:563956/251‑1 (MQ=255)                                                                                                                                                                                                                     aaCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACGCGATTAACATCAGCCATTTGATGGCGCAGATTTGCTACTATTTTTAAGCAGTTGCGCAGCTGCCGCAGGAAGA‑GCGCAACCAGCTTGTTGTCTCGGTGCCAAGCGGAAACTTCGGCAATTGGACGGCGGGTCTGCTGTCGACGGCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATGCCTTGCCGCGTTTCCTGcgc                    <  2:7764/251‑3 (MQ=255)                                                                                                                                                                                                                      gcTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACGTTTCCTGCacg                   <  1:14570/250‑1 (MQ=255)                                                                                                                                                                                                                         gAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGGTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACGTTTCCTGCacgac                 <  1:104435/250‑1 (MQ=255)                                                                                                                                                                                                                               ggCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACGTTTCCTGCACGACGGTCAg           <  1:129887/250‑1 (MQ=255)                                                                                                                                                                                                                                       ggTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGGCTGCTGGCGAAGGCACTCGGTCTGCCGGGGGAACGGTTTTTTGCTGCGACCAACGCGAAAGATACCGTGCCACGCTTTCTGCCAGACGGTCAGGGGGCAccc  >  2:430617/1‑251 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                  GACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAACTGAAAGTGGCGCTAGGGTTAAACTCGGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCTGTTGCGCAGCTGCCGCAGGA‑GACGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACGTTTCCTGCACGACGGTCAGTGGTCACCC  >  NC_000913/4123‑4603

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: ‑