Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,875	A→G	V513V (GTA→GTG)	thrA →	fused aspartate kinase/homoserine dehydrogenase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,875	0	A	G	95.7%	209.4 / ‑1.4	70	V513V (GTA→GTG)	thrA	fused aspartate kinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base G (27/40);  total (29/41)
Fisher's exact test for biased strand distribution p-value = 5.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.80e-01

CTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAA‑TGTACATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCAGTTGCGTTATG  >  NC_000913/1628‑2119                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                     cataGGGATCTTCGGAAAGCTCAATCTCTGTCGGGGTAAATAACGAAGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCAt                                                                                                                                                                                                                                                   <  2:125033/248‑1 (MQ=255)     gggATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGCTCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGGAGCGTCCACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACCCGAAGGCAATACCCACCAA‑TGTGCATGGcc                                                                                                                                                                                                                                               >  2:463057/1‑251 (MQ=255)      ggATCTTCTGAAAGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGATATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCt                                                                                                                                                                                                                                              >  1:112668/1‑251 (MQ=255)            tctGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATc                                                                                                                                                                                                                                         <  1:253589/250‑1 (MQ=255)               gAACGCTCAATCGCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGTCTAAATCTg                                                                                                                                                                                                                                       <  1:107304/249‑1 (MQ=255)                aaCGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGa                                                                                                                                                                                                                                     >  1:666637/1‑250 (MQ=255)                      cAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACt                                                                                                                                                                                                                                <  1:227325/249‑1 (MQ=255)                      cAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTgg                                                                                                                                                                                                                              <  1:256849/251‑1 (MQ=255)                      cAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTg                                                                                                                                                                                                                               <  1:51793/250‑1 (MQ=255)                      cAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCGGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTTTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACt                                                                                                                                                                                                                                <  2:374563/249‑1 (MQ=255)                              gTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAAAATAAACATATCGACTTACGTGTCTGCGGGGTTTACAACTCGAAGGCACTGCTCACACA‑TGTGCATGGCCTAAATCTGGAAAACTGCCAGga                                                                                                                                                                                                                         >  2:576554/1‑248 (MQ=255)                                   ggTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCTACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTgg                                                                                                                                                                                                                 <  1:626377/251‑1 (MQ=255)                                        aTAACGATGAAGCGACCACTGGCATGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGGGGCGTTGGCGGTGCGATGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTATCTGCGGTGTTGCCAACCCGAAGTACCTGCTCACCAAATGGAC‑TGGACCAAATCTGGAAAACTGGCAGgaaaaaacggcccc                                                                                                                                                                                                             >  2:444113/1‑239 (MQ=255)                                          aaCGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGGACTGGCGCAAGc                                                                                                                                                                                                          <  1:501096/251‑1 (MQ=255)                                            cGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCa                                                                                                                                                                                                        <  1:290569/251‑1 (MQ=255)                                                 atgCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAAt                                                                                                                                                                                                                                          <  1:3764/212‑1 (MQ=255)                                                 atgCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAAt                                                                                                                                                                                                                                          >  2:3764/1‑212 (MQ=255)                                                 atgCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAagag                                                                                                                                                                                                   <  1:268555/251‑1 (MQ=255)                                                      aCCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATTGCCTAAATCTGGAAAACTGGCCGGAAGAACTGGCGCAAGCCAAAGAGCCGtt                                                                                                                                                                                              >  2:20377/1‑251 (MQ=255)                                                         aCTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTaa                                                                                                                                                                                           <  1:20377/251‑1 (MQ=255)                                                            gtcgttctcgTTACTCTGCAGATGCTGTTCAAGTCCGATCAGGTTATCGACGTGTTTGTGATTGGCGTCGGTGTCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACTTGTCTTCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCGAAATCTTGAAAACTGGCAGGAAGAACTTGCGTAAGCCAAATAGCCGTTTAAt                                                                                                                                                                                          <  2:34004/241‑1 (MQ=255)                                                               gTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCgg                                                                                                                                                                                     <  1:249673/251‑1 (MQ=255)                                                                     gTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTAGCGTCGGTGGCGTTGGCGGTGCGATGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACAAATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAAACTGGAAAACTGGCAGGAAGAACTGGACCAAGCCAAAGAACCGTTTAATATCGGGCGCtt                                                                                                                                                                               >  2:477947/1‑251 (MQ=255)                                                                         cTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATCTCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGACAAAGAGCCGTTTAATCTAGGGCGCTTAAtt                                                                                                                                                                           >  2:163745/1‑251 (MQ=255)                                                                           catctGATGCTGGTCAATTCCGATCAGGTTATCGAAGTGTTTGGGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAATAATAAACATGTCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCGTACCAA‑TGTGCATGGCCTAAATCTGGAAAACGGGCAGGAAGAACTGGCGCAAGCTATAGAGCCGTTTAATCTCGGGCGCTTAATTc                                                                                                                                                                          <  2:655808/250‑1 (MQ=255)                                                                           catcaGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGGCGGTGGCGTTGGCGGTGCGTTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTc                                                                                                                                                                          <  1:388932/250‑1 (MQ=255)                                                                           catcaGATGCTGTTCAATACCGATCAGGTTATCGAAGTGGTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCg                                                                                                                                                                         <  1:247855/251‑1 (MQ=255)                                                                           catcaGAGGCTGTTCAATACCGATCAGGTTATCGAAGTGGTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCg                                                                                                                                                                         <  1:617089/251‑1 (MQ=255)                                                                               aGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTc                                                                                                                                                                     <  1:609608/251‑1 (MQ=255)                                                                                ggtgtggTTCAATACCGAGCAGGTTTTCGAAGTGGTTGTGATTGGCTGCGGTGGCGTGGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACGGGCAGGAAGAACTGGCGCAAGCCAAAGAGCTGTTTAATCTCTGGCGCTTAATTCGCCTCg                                                                                                                                                                    <  2:400836/245‑1 (MQ=255)                                                                                             aCCGATCAGGTTCTCGTAGTGTTTGTGATTCACGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAGAGCTGGCTGGGGAATAACCATATCGGCTTACGTGTCTGCGGTGTTGCCAACTCGCAGGCACTGCTCCCCAC‑TGAGCATGGCCTAAATCTGGaa                                                                                                                                                                                                                                    >  2:73224/1‑174 (MQ=255)                                                                                             aCCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGaa                                                                                                                                                                                                                                    <  1:73224/174‑1 (MQ=255)                                                                                                    aGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGAGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAA‑TGTACATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATctgctg                                                                                                                                                >  1:41477/1‑251 (MQ=255)                                                                                                       ttATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGCATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACAAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTAGAGCAAGCAAAAGCGCAGTTTTATATCGGGCGCTTTATTCGCCTCGTGAAAGAATATCATCTGATGAAc                                                                                                                                             >  2:29525/1‑251 (MQ=255)                                                                                                       ttATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAAc                                                                                                                                             >  1:439176/1‑251 (MQ=255)                                                                                                        tATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAAt                                                                                                                                                                                          >  2:125305/1‑205 (MQ=255)                                                                                                        tATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAAt                                                                                                                                                                                          <  1:125305/205‑1 (MQ=255)                                                                                                         aTCGAAGTGTTTGGGATTGGCGGCGGGGGAGTAGGCGGTGCGCTTCTGGTGCAACTGAATCGTCAACAAAGCGTGCTGATGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCATTGCTCACCAA‑TGTGCATGGCCTAAATTTGGAAAACTGGCAGGAAGAACTGGCGCAAGACAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAAcc                                                                                                                                            <  2:375278/250‑1 (MQ=255)                                                                                                           cGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCg                                                                                                                                          <  1:156721/250‑1 (MQ=255)                                                                                                                tgtTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCCAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGCAAGAATATCATCTGCTGAACCCGGTCAtt                                                                                                                                    >  2:563405/1‑251 (MQ=255)                                                                                                                 gtTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCAttg                                                                                                                                   <  1:545447/251‑1 (MQ=255)                                                                                                                       gATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACt                                                                                                                             >  2:176529/1‑251 (MQ=255)                                                                                                                       gATTGGAGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCCCCCCGTGAAAAAATATCATCTTCTGACCCCGGTCATTGTTGACt                                                                                                                             >  2:280962/1‑251 (MQ=255)                                                                                                                        aTTGGTGTCTGTGGCGTTGGGGGTTCGTTGATGGAGCAACTGAGGCGTCAATAACGCTGGCTGAAGAAGAAACATATCGTCTGACGTGTCTGCTGTGTTGCCAACTCGAAGGCACAGCACACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACGGGCGCAAGCCAAATAGCCGTTTAATCACGGGCGCGTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGGCATTGTTGACt                                                                                                                             <  2:543366/250‑1 (MQ=255)                                                                                                                                         tGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGACTAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGCACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGccg                                                                                                           >  2:261741/1‑249 (MQ=255)                                                                                                                                                          aGCAACTGAAGCGTCAACAAAGCTGTCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGcc                                                                                          <  1:611862/251‑1 (MQ=255)                                                                                                                                                          aGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGcc                                                                                          <  1:344392/251‑1 (MQ=255)                                                                                                                                                            cAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATgacgg                                                                                        >  1:89258/1‑247 (MQ=255)                                                                                                                                                                tGAAGCGTCAACAAAGCTGTCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGACCAATATGCCgccctc                                                                                    >  2:673986/1‑246 (MQ=255)                                                                                                                                                                   aGCGTCAACAAATATGGCTGAAGAATAATCATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACGGGCGCAAGCCAACGAGACGTTTAATCTCGGGCGCTTAATTAGCCTCGTGAACGAATATCATCTGCTGAACCCGGTCCTTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGGTCAATATGCCGACtttctt                                                                                 >  2:378470/1‑247 (MQ=255)                                                                                                                                                                              aGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGttt                                                                      <  1:373944/251‑1 (MQ=255)                                                                                                                                                                                        agaaTAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATTAGGCGAATTCTTGCGCGAAGGTTTCCACGTTGTc                                                            >  1:476267/1‑251 (MQ=255)                                                                                                                                                                                             aaaCATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGcc                                                       <  1:463057/251‑1 (MQ=255)                                                                                                                                                                                             aaaCATAACGACTTACGTGTCGGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGcc                                                       <  2:163702/251‑1 (MQ=255)                                                                                                                                                                                                     cgattgaCGTGTCGTCGGTGGTGCCAACTAGAAGGAACTTCGAACCAA‑TGTGCCTGGCCTAAATCTGGAAAACTTGCAGGAAGAACTGGCGTAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTTAGTCGCCTCGTGAAAGAAGATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCTGAGCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGGCTCGCCGAACaaa                                                <  2:289873/244‑1 (MQ=255)                                                                                                                                                                                                      gACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGGTGTCACGCCGAACcaaa                                               >  1:378867/1‑250 (MQ=255)                                                                                                                                                                                                       aCTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGGTTCCACGTTGTCACGCCGAACaaaa                                               >  1:329781/1‑249 (MQ=255)                                                                                                                                                                                                         ttACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGCCTTCTTGCGCGAAGGGTTCCACGTTGTCACGCCGAACAAAAAGGc                                           >  1:113085/1‑251 (MQ=255)                                                                                                                                                                                                           aCGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGcc                                          <  1:110532/250‑1 (MQ=255)                                                                                                                                                                                                                  tGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCAACCCGCACAAAAAAGCCAACAccac                                  >  2:599376/1‑249 (MQ=255)                                                                                                                                                                                                                        gtTGCCAACTCGAAGGCTCTGCTCACCAA‑TGTACATGGCCTTAATCGGGAAAACTGGCAGGAAGAACTGGCGAAAGCCAAAGTGCCGCTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCAGCTGCCGAACCCGGTCATTGTTGACTGCACTTCCAGCCATGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGTCCAACACCTCGTCGAt                            <  2:41477/251‑1 (MQ=255)                                                                                                                                                                                                                        gtTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGAt                            <  1:176529/251‑1 (MQ=255)                                                                                                                                                                                                                                cTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTcca                                                                                                                    <  2:68169/155‑1 (MQ=255)                                                                                                                                                                                                                                cTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTcca                                                                                                                    >  1:68169/1‑155 (MQ=255)                                                                                                                                                                                                                                  cGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCTCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATtactac                  <  1:132325/251‑1 (MQ=255)                                                                                                                                                                                                                                            gCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTg         <  1:217660/250‑1 (MQ=255)                                                                                                                                                                                                                                            gCTCAACAA‑TGTGCATGGCCTAAAACTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTg         <  1:17525/250‑1 (MQ=255)                                                                                                                                                                                                                                                aCCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCCTCTGTTGCGtc     >  1:274721/1‑249 (MQ=255)                                                                                                                                                                                                                                                aCCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTAGTCGATGGATTACTACCATCTGTTGCGtc     >  1:125735/1‑249 (MQ=255)                                                                                                                                                                                                                                                  cAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATg  <  1:599376/251‑1 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                     CTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAA‑TGTACATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCAGTTGCGTTATG  >  NC_000913/1628‑2119

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