Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,243,185	C→T	T156T (ACG→ACA)	folE ←	GTP cyclohydrolase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,243,185	0	C	T	98.4%	222.3 / NA	64	T156T (ACG→ACA)	folE	GTP cyclohydrolase 1
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (31/32);  minor base A (1/0);  total (32/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 = 9.69e-01

AACAAAATCGAGCGTGACGTTGCGCTCCATGGTTCCTGCCTTTTAATCAGTTGTGATGACGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGCGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATTTCATCGACCTTCATTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAAT  >  NC_000913/2242938‑2243430                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      cgcataaTCGAGCGTGACGTTGCGCTCCATTGTTCCTGCCTTTTAATCAGTTGTGATGACGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCCCCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTc                                                                                                                                                                                                                                                    >  2:282438‑M2/8‑251 (MQ=255)   cataaTCGAGCGTGACGTTGCGCTCCATTGTTCCTGCCTTTTAATCAGTTGTGATGACGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTc                                                                                                                                                                                                                                                    >  2:238998‑M2/6‑249 (MQ=255)         cGAGCGTGCCGTTGCGCTCCATTGTTCCTGCCTTTTAATAAATTATGACGACCCACAGCGAGCAGAACCTCGTCGCCCGTATTCTGACTGGATTTGACCTATCCCCCAACAGAGGACGTCGTCATGGCACTGGTTGCATCGCGGCCGAAACGCGCATTCACGCAGTAATGCACCGCGTCGATCGAGACAGCACCGTTATTGGTGCCCAGCAGCGTTTGTAGAGCAATAAGAATTTGCTGTGTCAGTCGTTc                                                                                                                                                                                                                                            >  1:660247/1‑251 (MQ=255)           aGCGTGACGTTGCGCTCCATTGTTCCTGCCTTTTAATCAGTTGTGATGACGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGACACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCt                                                                                                                                                                                                                                          >  1:570127/1‑251 (MQ=255)             cGTGACGTTGCGCTCCATTGTTCCTGCCTTTTAATCAGTTGTGATGACGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGc                                                                                                                                                                                                                                        >  2:55231/1‑251 (MQ=255)               tgacgcttcccccaacccttCCTCCCTTTTAACCAGATGAGATGACGAACAGCACCAAGAACCTCGTAGCGCGTATTCTGACTCGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCACTAAGAATTTGCTGTGTCAGACGTTCCTGCAc                                                                                                                                                                                                                                      >  2:394579/19‑251 (MQ=255)                                gTTCCTGCCTTTTAATCAGTTGTGATGACGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCa                                                                                                                                                                                                                      >  2:154795/1‑250 (MQ=255)                                   ccttccTTTTAATCAGTTGTGATGACGCACAGTGCGCAGAAACTCGCGGTGCGTATTCTGACTGGATTTGAACAATCCCCCAAACGAGGTCGTTGTAGACCCACTGGTTGCATCGCGGATGCCAAGCGCCTTCACGCAGCAATGCACCGCGTCGATCGAGACAGCCACGCTATTGGTGACCAGCAGCGTTTGTAGCGCAATAAGAAGTTGATGTGTCAGACGTTCCTGCACCTGCGGAACCTGGGCaa                                                                                                                                                                                                                     >  2:583328/5‑248 (MQ=255)                                      gCCTTTTAATCAGTTGTGATGACGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAAGCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACt                                                                                                                                                                                                               <  1:82670/251‑1 (MQ=255)                                       ccTTTTAATCAGTTGTGATGACGCAACGCCCGCAGACACTCGAGGAGCGAATACTGACTGGATTAGAAAAATCCACCCAGCGAGGGCGTTGTCGCGGCTCTGGTTGCATCACGAATGCCACGCGCCTTCACGCACTAATGCACCGCGTCGATCGCGACCGCCCCGTTATTGGTGCCCAGCAGCGTTTGTAGCGCATTCAGAATTTGCTGTGTCAGACGTTCCTGCACCTGAGGACGCTGGGCAAAGAACTg                                                                                                                                                                                                              >  1:316475/1‑251 (MQ=255)                                            tAATCAGTTGAAATGACGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGAGGCACTGGTTCCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGGTAGTGGTGCCCAGCAGCGTTTGTAGCGCAATAATCATTTGCTTTGTCAGACGTTCCTGCACCTGCGGACGCGGGGCACAGAACTGcaca                                                                                                                                                                                                          >  2:466928/1‑250 (MQ=255)                                             aaTCAGTTGTGATGACGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACaa                                                                                                                                                                                                         <  1:238998/250‑1 (MQ=255)                                                 aGTTGTGATGACGCACAACGCGCAGAAACTCGTGGCGCGAACTCTGCCTGGCTTTGAACAATCCACCAAGAAAGTTCGCTGTCGTGGCACTGGTCGCATCGCAGATGCCCCGCGCCTTAACGCAGTAATGCACCCCGTCGACCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCgg                                                                                                                                                                                                    >  1:266918/1‑251 (MQ=255)                                                        atgaCGCACAGAGCGCAGTAACTCGTGGCGCGTCTTCTGACTGGAGTTGAACAATCCACCAAGAGCGGCCGTTGTCGTGCAACTGGTTGCATCGCGGATGCCGCGCGCCTTCACGCAGTAATGCACCGCGTCGATCGCGACAGCCACGTTATTGGAGCCCAGCAGCGTTTGTAGCGCGATAAGAATTTGCTGTGTCAGACGTTCCTGAACCTGCGGCCGCTCGGCCACGAACTGCCCAATGAGGTTAAttt                                                                                                                                                                                             >  2:465593/1‑251 (MQ=255)                                                           aCGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTg                                                                                                                                                                                          >  1:580906/1‑251 (MQ=255)                                                               acaGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAAt                                                                                                                                                                                                        <  1:231131/233‑1 (MQ=255)                                                               acaGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAAt                                                                                                                                                                                                        >  2:231131/1‑233 (MQ=255)                                                                 agcgcgcagaacttagttgcgtgcaTTCTGACAGGATTTGCACCATCAAACCCGAGAGGCCGTTGTCGCGGGACTGGCTCCTTAGCGGATGCCACGCGCCTTCACGCAGTAATGCCCCGCGTCGCTCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTCTAGCGCAATAAGAATTTGCTGTGCCAGACGTTCCTGCAACTGCGGTCGCTGGGCAAAGAACTGCACACTGCGGTTAATTTTTGACAGAc                                                                                                                                                                                    >  2:672619/25‑251 (MQ=255)                                                                   cgcgcCGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAATAGAGGTCGTTGTCGTGGCACTGGTGGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCg                                                                                                                                                                                  <  2:13481/251‑1 (MQ=255)                                                                   cgcgcAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCg                                                                                                                                                                                  <  1:397338/251‑1 (MQ=255)                                                                      gcAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTTCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATc                                                                                                                                                                               <  1:121085/251‑1 (MQ=255)                                                                      gcAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATc                                                                                                                                                                               <  1:335903/251‑1 (MQ=255)                                                                                         aTTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACGGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGatat                                                                                                                                                            <  2:243732/251‑1 (MQ=255)                                                                                          ttCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGatat                                                                                                                                                            >  1:668205/1‑250 (MQ=255)                                                                                               aCTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCa                                                                                                                                                      <  2:580906/251‑1 (MQ=255)                                                                                                  ggATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTTGCACGGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGAGAGACCGATCACCGAATCTTTCGGGATATAGGCCAc                                                                                                                                                     <  2:239486/249‑1 (MQ=255)                                                                                                      atGAACAATCCACCAAGAGCGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGc                                                                                                                                               >  2:55783/2‑251 (MQ=255)                                                                                                       tGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTGCACGCAGTAATGCACCGCGTCGATCGAGGCATCGACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAATAATTTGCTGTGGCAGACGTTCCTGCTCCTGCGGTCGGTGGGCAAAGAACGGCACAATGCGGTTAATTTGTGTCAGGCCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGc                                                                                                                                               <  1:300357/250‑1 (MQ=255)                                                                                                            aaTCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTGTTGACAGACCGATCACCGAATCTGTCGTGATATAGGCCACCGTCGCTTTGcc                                                                                                                                         <  2:662306/251‑1 (MQ=255)                                                                                                            aaTCCACCAAGAGAGGTCGGTGTCGTGGCACTGGGTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAAGTTTTGACAGACCGATCGCCGAAGCTTTCGGGATATAGGCCACCGTCGCTTTGcc                                                                                                                                         <  2:444762/251‑1 (MQ=255)                                                                                                                         agGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCttt                                                                                                                                            >  1:651364/1‑235 (MQ=255)                                                                                                                         agGTCGTTGTCGTGGCACTGGTTGCAGCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGGCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAAGCTTTCGGTATATAGGCCACCGTCGCttt                                                                                                                                            <  2:651364/235‑1 (MQ=255)                                                                                                                                 gTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGgtg                                                                                                                    <  1:55231/251‑1 (MQ=255)                                                                                                                                        cACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAg                                                                                                             <  2:315828/251‑1 (MQ=255)                                                                                                                                                 cATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTGTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTc                                                                                                    <  1:330740/251‑1 (MQ=255)                                                                                                                                                   tCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCCTCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTc                                                                                                    >  2:506241/1‑249 (MQ=255)                                                                                                                                                    cgcgGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTTCTGGTCaga                                                                                                 <  1:484312/251‑1 (MQ=255)                                                                                                                                                      cgGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGt                                                                                               <  1:587766/251‑1 (MQ=255)                                                                                                                                                       ggATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTg                                                                                              <  1:282438/251‑1 (MQ=255)                                                                                                                                                           gCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGatat                                                                                          >  1:281197/1‑251 (MQ=255)                                                                                                                                                            ccACGCGCCTTCACTCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGTCGTTCCTTCCCCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAAATTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCGTTGCCATCGGGGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATc                                                                                         <  1:465593/251‑1 (MQ=255)                                                                                                                                                             cACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATcg                                                                                        <  1:162437/251‑1 (MQ=255)                                                                                                                                                              aCGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCTCGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATcgc                                                                                       >  1:243799/1‑251 (MQ=255)                                                                                                                                                                 cgcCTTCACGCAGTACTGCACCGCGTCGCTCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACACAATGGTGTTCCCAGGTGCTGGTCAGAGTGATATCgcgcc                                                                                    >  2:85351/1‑250 (MQ=255)                                                                                                                                                                  gcCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAAt                                                                                                                        <  2:375409/214‑1 (MQ=255)                                                                                                                                                                  gcCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAAt                                                                                                                        >  1:375409/1‑214 (MQ=255)                                                                                                                                                                  gcCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCAc                                                                                   <  1:154795/251‑1 (MQ=255)                                                                                                                                                                    cTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGAGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACgg                                                                                 >  1:49826/1‑251 (MQ=255)                                                                                                                                                                            aGTAATGCACCGCGTCGATCGCGACAGCCACGTTATTGGTACCCAGCAGCGTTTGAAGCGCAATACGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGCGGGCTATCGCGCACGGTGCCCAtt                                                                         >  2:649452/1‑251 (MQ=255)                                                                                                                                                                                      cgcgTCAAACGAGAAAGCCAAGATATTGGTGCACAGCAGCGTATGTAGCGCAGTAAGAATTTGCTGAGTCAGACGTTCCTGCACCTGCGGACGCTGGGAAACGAACTGCAAACTGCGGTTCATTTTTGACAGACCGCTCAACGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGCGTGATATCGCGCACGGTGACCATTTCATCGACCt                                                               >  2:636070/1‑251 (MQ=255)                                                                                                                                                                                       gcgTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAAGTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAg                                                                                                                                                                                                                   <  1:662626/102‑1 (MQ=255)                                                                                                                                                                                       gcgTCGATCGACAAAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAg                                                                                                                                                                                                                   >  2:662626/1‑102 (MQ=255)                                                                                                                                                                                        cgTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATTTCATCGACCTTc                                                             <  2:188678/251‑1 (MQ=255)                                                                                                                                                                                         gTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATTTCATCGACCTTc                                                             >  1:212316/1‑250 (MQ=255)                                                                                                                                                                                          ccgatcgaGACAGTCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGTCCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATTTCATCGACCTTCAt                                                           <  1:466928/250‑1 (MQ=255)                                                                                                                                                                                                  gaCAGCCACGATATTGGACCCCAGCAGCGATTGTAGCGCAGTAAGACTTTGCTGTGTCCGACGTTCCTGACCCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGCCCGCTCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGCGTTCACTGGTGCTGGTCAGCGTGATATCGCGCACGGTGACCATTTCATCGACCTTCgttttgtttt                                                   >  2:616678/1‑251 (MQ=255)                                                                                                                                                                                                       ccACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCTCCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATCTCATCGACCTTCATTTTGTTTTCAATg                                              <  2:281197/251‑1 (MQ=255)                                                                                                                                                                                                          cGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATGTGCGGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCGTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTTGTCAGAGTGATATCGCGCACGGTGACCATTTCATCGACCTTCATTTTGTTTTCAATGAgg                                           <  2:570127/251‑1 (MQ=255)                                                                                                                                                                                                              aTTGGGGCCCAGCAGCGTTTGTAGCGCAATAGGAATTTGCTGTGTCAGACGGTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCGCCGTCGCGTTGCCATCGATGGGATCAAAATGGTGTTCACATGTGCTGGTCAGGGGGATAGCGTGCACGTTGACCATGTCATCTACCTACAGTTTGTTTTCAATGAGGGTGa                                       <  2:242739/251‑1 (MQ=255)                                                                                                                                                                                                                 ggTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATTTCATCGACCTTCATTTTGTTTTCAATGAGGGTGAttt                                    >  2:478008/1‑251 (MQ=255)                                                                                                                                                                                                                 ggTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATTTCATCGACCTTCATTTTGTTTTCAATGAGGGTGAttt                                    >  2:182380/1‑251 (MQ=255)                                                                                                                                                                                                                      ccAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATTTCATCGACCTTCAGTTTGTTTTCAATGAGGGTGATTTTCggg                               <  2:212316/251‑1 (MQ=255)                                                                                                                                                                                                                                     gcgcAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATTTCATCGACCTTCATTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATcc                <  1:171094/251‑1 (MQ=255)                                                                                                                                                                                                                                           tAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATTTCATCGACCTTCATTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCg          >  2:114786/1‑251 (MQ=255)                                                                                                                                                                                                                                                    gctgcGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAATAACTGCACAACACTGTGAATGTTTGACAGGCCGATCACCTAATCTTGCGGGATATAGGCCACCGTCGCGTTGCCATCGATGGTAACAAAATGTTGTTCACTGGTTCTGGTCAGAGTGATATCTCGCACGGGTACCATGTCATCTCCCGTCATTTTGTTCTCAATTAAGGTGATTTTCGGGAAAGTGTCGTAATCCATACCAGAGAAAAt  <  2:295712/250‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      AACAAAATCGAGCGTGACGTTGCGCTCCATGGTTCCTGCCTTTTAATCAGTTGTGATGACGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGCGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATTTCATCGACCTTCATTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGGAGAAAAT  >  NC_000913/2242938‑2243430

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