Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,351	A→G	L206L (CTA→CTG)	thrC →	threonine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,351	0	A	G	98.4%	190.7 / ‑4.7	63	L206L (CTA→CTG)	thrC	threonine synthase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (31/31);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.45e-01

AAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAACTGAAAGTGGCGCTAGGGTTAAACTCGGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCTGTTGCGCAGCTGCCGCAGGA‑GACGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACGTTTCCTGCACGACGGTCAG  >  NC_000913/4106‑4594                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     aaGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTACCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGtt                                                                                                                                                                                                                                                 <  2:305792/251‑1 (MQ=255)    ccAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTaa                                                                                                                                                                                                                                               >  1:295962/1‑250 (MQ=255)                  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)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     AAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAACTGAAAGTGGCGCTAGGGTTAAACTCGGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCTGTTGCGCAGCTGCCGCAGGA‑GACGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACGTTTCCTGCACGACGGTCAG  >  NC_000913/4106‑4594

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