Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,938,783	G→T	V379F (GTT→TTT)	pykA →	pyruvate kinase 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,938,783	0	G	T	92.8%	177.1 / 5.2	69	V379F (GTT→TTT)	pykA	pyruvate kinase 2
Reads supporting (aligned to +/- strand):  ref base G (1/4);  new base T (27/37);  total (28/41)
Fisher's exact test for biased strand distribution p-value = 6.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.24e-02

ACGCGTGCAGAAGTCATGGACGTAGCAAACGCCGTTCTGGATGGTACTGACGCTGTGATGCTGTCTGCAGAAACTGCCGCTGGGCAGTATCCGTCAGAAACCGTTGCAGCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGACCTGGTGATTGTCACCC  >  NC_000913/1938537‑1939029                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       aCGCGTTCAGAAGTCATGGACGTACCAAACGCCGTTCTGGATGGTACTGACGCTGTGATGCTGTCTGCAGAAACTGCCGCTGGGCAGTATCCGTCAGAAACCGTTGCAGCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCCCCTGAAAGGCTTTAc                                                                                                                                                                                                                                                    <  2:127871/251‑1 (MQ=255) aCGCGTGCAGAAGTCATGGACGTAGCAAACGCCGTTCGGGATGGTACTGACGCTGTGATGCTGTCTGCAGAAACTGCCGCTGGGCAGTATCCGTCAGAAACCGTTGCAGCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCTTTAc                                                                                                                                                                                                                                                    <  1:148959/251‑1 (MQ=255)   gcgggCAGAAGGCATGGACGTAGCAAACGCCGGTCTGGATGGTACTGACGCTGTGATGCTGTCTGCAGAAACTGCCGCTGGGCAGTATCCGTCAGAAACCGTTGCAGCCATGGCGCGCGGTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCAGGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCTTTACgg                                                                                                                                                                                                                                                  <  2:35489/247‑1 (MQ=255)    cgTGCAGAAGTCATGGACGTAGCAAACGCCGTTCTGGATGGTACTGACGCTGTGATGCTGTCTGCAGAAACTGCCGCTGGGCAGTATCCGTCAGAAACCGTTGCAGCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAGGGCTTTACGGc                                                                                                                                                                                                                                                 >  2:106945/1‑251 (MQ=255)    agTGCAGAAGTCATGGACGTAGCAAACGCCGTTCTGGATGGTACTGACGCTGTGATGCTGTCTGCAGAAACTGCCGCTGGGCAGTATCCGTCAGAAACCGTTGCAGCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCTTTACGGc                                                                                                                                                                                                                                                 >  2:106958/2‑251 (MQ=255)        cAGAAGTCATGGACGTAGCAAACGCCGTTCTGGATGGTACTGACGCTGTGATGCTGTCTGCAGAAACTGCCGCTGGGCAGTATCCGTCAGAAACCGTTGCAGCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCTTTACGGCGatc                                                                                                                                                                                                                                             <  2:536471/251‑1 (MQ=255)        cAGAAGTCATGGACGTAGCAAACGACGTTCTGGATGGTACTGACGCTGTGATGCTGTCTGCAGAAACTGCCGCTGGGCAGTATCCGTCAGAAACCGTTGCAGCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGCCGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTCCGCAGCTAACCACCTGAAAGGCTTTACGGCGatc                                                                                                                                                                                                                                             >  2:123399/1‑251 (MQ=255)         aGAAGTCATGGACGTAGCAAACGCCGTTCTGGATGGTACTGACGCTGTGATGCTGTCTGCAGAAACTGCCGCTGGGCAGTATCCGTCAGAAACCGTTGCAGCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCTTTACGGCGatc                                                                                                                                                                                                                                             >  1:604640/1‑250 (MQ=255)                  ggACGTAGCAAACGCCGTTCTGGATGGTACTGACGCTGTGATGCTGTCTGCAGAAACTGCCGCTGGGCAGTATCCGTCAGAAACCGTTGCAGCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCTTTACGGCGATCATCACCATg                                                                                                                                                                                                                                    >  2:629315/1‑250 (MQ=255)                    aCGTAGCAAACGCCGTTCTGGATGGTACTGACGCTGTGATGCTGTCTGCAGAAGCTGCCGCTGGGCAGTTTCCGGCAGAAACCGTTGCAGCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCGAAACACGGTCTGGACGTGCAGTTCGGGAATGTGGAGGAATCTATTGCCATGTCAGTAATGTACGCAGCTAACCACCTGAAAGGCTTTACGGCGATCATCACCATGAcc                                                                                                                                                                                                                                 <  1:282841/251‑1 (MQ=255)                      gTAGCAAACGCCGTTCTGGATGGTACTGACTCTTTGATTCTGTCTGCAGAAACTGCCGCTTGGCAGGATCCGTCAGAAACCGTTGCAGCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTTGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACTCAGCTAACCACCTGAAAGGCTTTACGGCGATCACCACCATGTCCGa                                                                                                                                                                                                                               <  2:146124/251‑1 (MQ=255)                         gCAAACGCCGTTCTGGATGGTACTGACGCTGTGATGCTGTCTGCAGAAACTGCCGCTGGGCAGTATCCGTCAGAAACCGTTGCAGCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCTTTACGGCGATCATCACCATGACCGAATc                                                                                                                                                                                                                            >  2:155999/1‑251 (MQ=255)                            aaCGCCGTTCTGGATGGTACTGACGCTGTGCTGCTGTCTGCAGAAACTGCCGCTGGGCAGTATCCGTCAGAAAACGTTGCAGCCATGGCGAGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATATCAGCAATGTACGCAGCTAACCACCTGACAGGCTTTACGGCGATCATCACCATGACCGAATCggg                                                                                                                                                                                                                         >  2:302925/1‑251 (MQ=255)                             aCGCCGTTCTGGATGGTACTTACGCTGTGATGCTGTCTGCAGAAACTGCCGCTGGGCAGTATCCGTCAGAAACCGTTGCAGCCATGGCGCGCGTTTGCCTGGGTTCGGAAAAAATCCCGATCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTCCGCAGCTAACCACCTGAAAGGCTTTACGGCGATCATCACCATGACCGAATCGGGt                                                                                                                                                                                                                        <  1:392202/251‑1 (MQ=255)                                                 gACGCTGTGATGCTGTCTGCAGAAACTGCCGCTGGGCAGTATCCGTCAGAAACCGTTGCAGCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCTTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCt                                                                                                                                                                                                     <  1:153189/250‑1 (MQ=255)                                                 gACGCTGTGATGCTGTCTGCAGAAACTGCCGCTGGGCAGTATCCGTCAGAAACCGTTGCAGCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCTTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTc                                                                                                                                                                                                    <  2:239234/251‑1 (MQ=255)                                                         gATGCTGTCTGCAGAAACTGCCGCTGGGCAGTATCCGTCAGAAACCGTTGCAGCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCTTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCa                                                                                                                                                                                            <  1:182595/251‑1 (MQ=255)                                                               gTCTTCAGAAACTGCCGCGTGGCAGTATCCGTCAGAAACCGTTGCAGCCATGGCTCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTAGTGCCATGGCAGCAATGTACGCAGCTAACCACCTGAAAGGCTTTACGGCGAGCAGCACCATGACCGAATCGGGTCGTACCGCGGTGATGACCGCCCTTATCAGctc                                                                                                                                                                                        <  2:421655/249‑1 (MQ=255)                                                                   gCAGAAACTGCCGCTGGGCAGTATCCGTCAGAAACCGTTGCAGCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCTTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTgg                                                                                                                                                                                     <  1:261417/248‑1 (MQ=255)                                                                    cAGAAACTGCCGCTGGGCAGTATCCGTCAGAAACCGTTGCAGCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCTTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTg                                                                                                                                                                                 <  1:103915/251‑1 (MQ=255)                                                                        aaCTGCCGCTGGGCAGTATCCGTCAGAAACCGTTGCAGCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCTTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCaa                                                                                                                                                                             <  1:150610/251‑1 (MQ=255)                                                                             ccGCTGGGCAGTATCCGTCAGAAACCGTTGCAGCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCTTTACGGCGATCATCACCATGACCGAATCGGGTCGTACcg                                                                                                                                                                                                                 <  1:323014/210‑1 (MQ=255)                                                                             ccGCTGGGCAGTATCCGTCAGAAACCGTTGCAGCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCTTTACGGCGATCATCACCATGACCGAATCGGGTCGTACcg                                                                                                                                                                                                                 >  2:323014/1‑210 (MQ=255)                                                                             ccGCTGGGCAGTATCCGTCAGAAACCGTTGCAGCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTAGTGCCATTGCAGCAATGTACGCAGCTAACCACCTGAAAGGCTTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTc                                                                                                                                                                        <  1:289023/251‑1 (MQ=255)                                                                                         aTCCGTCAGAAAACGTTGCAGCCATGGCGCGCGTTTGCCTGGGAGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTCTTGCCATGTCAGCAATGTACGCCGCTAACCACCTGAACGGCTTTACGGAGATCATCACCATGACCGAATCGGTTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCgcgc                                                                                                                                                            >  1:355216/1‑251 (MQ=255)                                                                                           ccGTCAGAAACCGTTGCAGCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCTTTACGGCGATCAGCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTcgcgac                                                                                                                                                           <  1:106958/250‑3 (MQ=255)                                                                                           ccGTCAGAAACCGTTGCAGCCATGGCGCGCGTGTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCTTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGcc                                                                                                                                                           <  1:106945/250‑1 (MQ=255)                                                                                              tCAGAAACCGTTGCAGCCATGGCGCGCGTTTGCCTGGGTGCGGAACAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGCCGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCTTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTCTCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATg                                                                                                                                                        >  2:594221/1‑250 (MQ=255)                                                                                                      cGTTGCAGCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCTTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGc                                                                                                                                               >  1:339012/1‑251 (MQ=255)                                                                                                           cAGCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGCTCGACAATAGGGAAGAAGCTATTGACATGTCAGCAATGTACGCAGCTCACCACCTGAAAGGCTTTACGGCGATCATCAACATGACCGAATCGGGTCGTACAGCGCTGATGACCTCCCGTCTCAGCCCTGGTCTGCCAATTTGCGCCATGTCgcgcaaagaaagtacactgcac                                                                                                                                          >  2:342880/1‑233 (MQ=255)                                                                                                             gCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCt                                                                                                                                        <  1:188480/251‑1 (MQ=255)                                                                                                             gCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTATAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCTTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCACTTTTCGACATGTCGCGCCATGAACGTACGCTGAACCt                                                                                                                                        >  1:51777/1‑251 (MQ=255)                                                                                                              ccATGGCGCGAGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCTTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTg                                                                                                                                       >  1:406801/1‑251 (MQ=255)                                                                                                               cATGGCGCGCGTTTGCCTGTGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCTTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGa                                                                                                                                      <  1:497482/251‑1 (MQ=255)                                                                                                                       gcgTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCTTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGCTGACCTCCCGTATCCGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGCACGTACGCTGAACCTGACTGctctc                                                                                                                              >  2:218041/1‑251 (MQ=255)                                                                                                                           ttGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCTTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATc                                                                                                                          >  1:96445/1‑251 (MQ=255)                                                                                                                             gCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTGCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCAGCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCg                                                                                                                         <  2:262302/250‑1 (MQ=255)                                                                                                                              ccTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCTTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTg                                                                                                                       <  1:428741/251‑1 (MQ=255)                                                                                                                              ccTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGGCAGCAATGTACGCAGCTAACCCCCTGAAAGGCTTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTg                                                                                                                       <  2:547388/251‑1 (MQ=255)                                                                                                                                 gggTGCGGAAAAAATCCCGAGCAACAACGTTTCTAAACACCGGCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCTCCTGAAAGGCTTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACGGCTCTCTATCGTGGCg                                                                                                                    <  2:60992/251‑1 (MQ=255)                                                                                                                                    tGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGtt                                                                                                                  >  2:38466/1‑250 (MQ=255)                                                                                                                                      cGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCTTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTg                                                                                                                                             >  1:221739/1‑221 (MQ=255)                                                                                                                                      cGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCTTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTg                                                                                                                                             <  2:221739/221‑1 (MQ=255)                                                                                                                                                  cGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCTTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACttt                                                                                                   <  2:421236/251‑1 (MQ=255)                                                                                                                                                  cGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACttt                                                                                                   <  1:191766/251‑1 (MQ=255)                                                                                                                                                       aTCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCTTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCCTGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTgct                                                                                                >  2:345652/1‑247 (MQ=255)                                                                                                                                                        tCAACGTTTCTAAACACCGTCTGGACGGTCAGTTCGACAATGTGGACGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCTTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAgc                                                                                             >  1:76542/1‑251 (MQ=255)                                                                                                                                                           aCGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCt                                                                                          <  1:489187/251‑1 (MQ=255)                                                                                                                                                           aCGTTTCTAAACACCGTCTGGACGTGCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTGCGCAGCTAACCACCTGAAAGGCTGTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCt                                                                                          <  1:282684/251‑1 (MQ=255)                                                                                                                                                              tttCTAACCACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCTTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCCTGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAAt                                                                                       >  1:471878/1‑251 (MQ=255)                                                                                                                                                               ttCTAAACCCCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTCCGCAGCTAACCACCTGAAATGCTTTACGGCTATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGTCCTCCCGTATCAGCTTTGGTCTGCCAATTTTCTCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACTCCTGTGCACTTTGATAGCGCttat                                                                                       <  2:44498/250‑4 (MQ=255)                                                                                                                                                                   aaaCACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCTTGACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGa                                                                                     <  1:123399/248‑1 (MQ=255)                                                                                                                                                                            cTGGACGTTCAGTTAGACAATGTGGAAGAAGCTATTGCCATGACAGCAATGTACGCAGCTAACCACCTGAAAGGCTTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTATGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCCCTTTGATAGCGCTAATGACGGCGTagcagc                                                                         >  2:595536/1‑251 (MQ=255)                                                                                                                                                                             tGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCTTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGACCGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCGGGAGCACTTTGATATCGCTCATGCCGGCGTAGCAGCt                                                                        >  2:506078/1‑251 (MQ=255)                                                                                                                                                                                  gTTCAGTTCGACAATGTGGAAGAAGCTATTGCTACGTCAGCAATGTATGCAGATAACCACCTGAAAGGCTTTACGGCGCTCATCACCATGACTGAATCGGTTCGTACCGTGCTTATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAGCGCACGCTGAACCTGACTCATCTCTATCGTGGCGTTAAGCCGGTGCACTTGGATAGAGATAATGACGGCGGAGGAGCTGcc                                                                     >  2:133868/1‑249 (MQ=255)                                                                                                                                                                                            aCAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCTTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGtt                                                         <  1:261887/251‑1 (MQ=255)                                                                                                                                                                                               aTGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCTTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGAGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTAAAt                                                      >  1:298787/1‑251 (MQ=255)                                                                                                                                                                                                     aagaagCTATTGCCATGTCAGCAATGTACGCAGCTATTCTCCTGAAAGGCTTTACGTCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACTCCGGTGCACTTTTATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATctgctg                                                <  1:382410/251‑1 (MQ=255)                                                                                                                                                                                                     aagaagCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGCAAGGCTTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGCAGCGGTTAATctgctg                                                >  1:54374/1‑251 (MQ=255)                                                                                                                                                                                                                gCCATGTCAGCAATGTACGCATCTAACCACCTGAAAGGCTTGACGTCGATCATCCCCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCATCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACGTCTCTCTATCGTTGCGTTACGCCGGTGCACTTTGATAGCGCTAAGGACGGCGTAGCATCTGCCAGCGAGGCGGTTAATCTGCTGCTCGCTAAAg                                      <  2:618949/250‑1 (MQ=255)                                                                                                                                                                                                                  cATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCTTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGtt                                   >  1:172751/1‑251 (MQ=255)                                                                                                                                                                                                                          aaaTGTACGCAGCGAACCTCCTGAAAGGCTTTACGGCGATCATCACCATGACCGAATCGGGTCGTGCCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACGGCTCTCTATCGTGGCGTTACGCCGGTGCACGTTGAGAGCGCTAATGTCGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGTGCGATAAAGGTTACTTGATg                           <  2:302395/250‑1 (MQ=255)                                                                                                                                                                                                                                  gCAGCTAACCACCTGAAAGGCTTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGa                   <  1:629315/251‑1 (MQ=255)                                                                                                                                                                                                                                        aaCCATCTGAAAGGCTTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTGCGCCATGTCGCGCCAGGAACGTACTCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACGGTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGACCTGGt             <  1:620789/251‑1 (MQ=255)                                                                                                                                                                                                                                        aaCCACCTGAAAGGCTTTACGGCGATCATCACCATGACCGAAGCGTGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGACCTg               <  2:99541/249‑1 (MQ=255)                                                                                                                                                                                                                                        aaCCACCTGAAAGGCTTTACGGCCACCATCACCATGAACGAATAGGGTCGTACCGCGCTGATGACCTCCCGTATCAGATCTGGTCTGCCTAGTTTCGACATGTCGCGCCAAGAACGTACGATGAACCTGACTGCTCTCTATCGTGGCGTTAAGCCGGTGAACTTTGATAGCGCTAATGAAGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGAGATAAAGGTTACTTGAGGTCTGGTGACCTgg              <  2:355216/250‑1 (MQ=255)                                                                                                                                                                                                                                             ccTGAAAGGCTTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGACCTGGTGAt          <  2:604640/249‑1 (MQ=255)                                                                                                                                                                                                                                                gAAAGGCTTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGTTAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGACCTGGTGATTGt       <  1:351587/249‑1 (MQ=255)                                                                                                                                                                                                                                                   aGGCTTTACGGCGATCATAACAATCACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGACCGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAACGGTTACTTGATGTCTGGTGACCTGGTGATTGTCAccc  >  2:437515/1‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       ACGCGTGCAGAAGTCATGGACGTAGCAAACGCCGTTCTGGATGGTACTGACGCTGTGATGCTGTCTGCAGAAACTGCCGCTGGGCAGTATCCGTCAGAAACCGTTGCAGCCATGGCGCGCGTTTGCCTGGGTGCGGAAAAAATCCCGAGCATCAACGTTTCTAAACACCGTCTGGACGTTCAGTTCGACAATGTGGAAGAAGCTATTGCCATGTCAGCAATGTACGCAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGACCTGGTGATTGTCACCC  >  NC_000913/1938537‑1939029

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