Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,553,671	C→T	F257F (TTC→TTT)	uxuB →	D‑mannonate oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,553,671	0	C	T	97.0%	191.4 / NA	67	F257F (TTC→TTT)	uxuB	D‑mannonate oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (38/27);  minor base G (0/2);  total (38/29)
Fisher's exact test for biased strand distribution p-value = 1.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.56e-01

GGTTGAAAGCGTTTACGGTGATGTCCTGCGATAACGTGCGTGAAAACGGTCATGTGGCGAAGGTCGCGGTACTGGGGCTGGCTCAGGCGCGTGACCCGCAGCTGGCGGCATGGATTGAAGAAAATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCTTACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTCCGTCAGTGGGTGATTGAAGATAACTTCGTTAATGGTCGCCCGGATTGGGATAAAGTGGGCGCACAGTTCGTTGCAGACGTTGTGCCGTTCGAAATGATGAAGCTGCGTATGCTGAACGGCAGCCACTCTTTTCTGGCGTACCTCGGTTACCTCGGCGGCTATGAAACCATTGCCGACACCGTGACTAACCCGGCTTATCGCAAAGCGGCCTTTGCCCTGATGATGCAGGAACAAGCG  >  NC_000913/4553427‑4553908                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                              ggTTGAAAGCGTTTACGGTGATGTCCTGCGATAACGTGCGTGAAAACGGTCATGTGGCGAAGGTCGCGGTACTGGGGCTGGCTCAGGCGCGTGACCCGCAGCTGGCGGCATGGATTGAAGAAAATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCTTACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTTCGTCAg                                                                                                                                                                                                                                         <  1:112581/251‑1 (MQ=255)      aaagggtctACGGTGATGGCCTGCGATAACGTGCGTGAAAACGGTCATGTGGCGTAGGGCGAGGTACTGGGGCTTGATCAGGCGCGTGACCCGCAGCTGGCGGCATGGATTGAAGAAAATGTCACCGTCCCGGGCACCAGGGTGGACAGCATCGGTCCGGCGGCGACGCCAGAAACCTTACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTTCGTCAGTGGGt                                                                                                                                                                                                                                    <  1:401080/243‑1 (MQ=255)         gCGTTTACGGTGATGTCCTGCGATAACGTGCGTGAAAACGGTCATGTGGCGAAGGTCGCGGTACTGGGGCTGGCTCAGGCGCGTGACCCGCAGCTGGCGGCATGGATTGAAGAAAATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCTTACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTTCGTCAGTGGGTGAt                                                                                                                                                                                                                                 >  1:313599/1‑251 (MQ=255)          cGTTTACGGTGATGTCCTGCGATAACGTGCGTGAAAACGGTCATGTGGCGAAGGTCGCGGTACTGGGGCTGGCTCAGGCGCGTGACCCGCAGCTGGCGGCATGGATTGAAGAAAATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCTTACAGGCAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTTCGCCAGTGGGCGAtt                                                                                                                                                                                                                                >  2:172313/1‑251 (MQ=255)           ggttACGGTGATGTCCTGCGATAACTGGCGTGAAAACGGTCATGTGGCGAAGGTGGCGGTACTGGGGCTGGCTCAGGCGCGTGAGCCGCAGCTGGCGGCATGGATTGAAGAAAATGTCACCGTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCTTACAGGAAATTGCTGCCCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTTCGTCAGTGGGTGATTg                                                                                                                                                                                                                               <  2:298649/249‑1 (MQ=255)              tACGGTGATGTGCTGGGATAACGGGCGTGAAAACGGTCAGGTGGCGAAGGTCTCGGGACTGGGGCTGGCGCAGGCGCGTGACCCGCAGCTGGCGGCATGGATTGAAGTAAATGTCGCCTTCCCGTGCACCAGGGTGGACCGCATCGTTCCGGCGGCGACGCCAGAAACCTTGCAGGAGATGGCTGACCGGCGGTGGGTTTACGTGCCGTGCTCCAGTGCCTGCGGGCCGGTTCGTCAGTGTGTGGTTGaa                                                                                                                                                                                                                             <  2:180415/250‑1 (MQ=255)               aCGGTGATGTCCTGCGATAACGTGCGTGAAAACGGTCATGTGGCGAAGGTCGCGGTACTGGGGCTGGCTCAGGCGCGTGACCCGCAGCTGGCGGCATGGATTGAAGAAAATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCTTACAGGAAATTGCGGACCAGCTGGGAGTTTACGACCCGTGCGCCATTACCTGCGAACCGTTTCGTCAGTGGGTGATTGa                                                                                                                                                                                                                              >  2:305331/1‑248 (MQ=255)               aCGGTGATGTCCTGCGATAAAGTGCGTGAAAACGGTCATGTGGCGAAGGTCGAGGTACTGGGGCTGGCTCCGGCGCGTGACACGCAGCTGGCGGCATGGATTGAAGCAAATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCTTACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTTCGTCAGTGGGTGATTGa                                                                                                                                                                                                                              >  2:179716/1‑248 (MQ=255)                cGGTGATGTCCTGCGATAACGTGCGTGAAAACGGTCATGTGGCGAAGGTCGCGGTACTGGGGCTGGCTCAGGCGCGTGACCCGCAGCTGGCGGCATGGATTGAAGAAAATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCTTACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGTCTGCGAACCGTTGCGTCAGTGGGTGATTGAAGAt                                                                                                                                                                                                                          <  1:326227/251‑1 (MQ=255)                cGGTGATGTCCTGCGATAACGTGCGTGAAAACGGTCATGTGGCGAAGGTCGCGGTACTGGGGCTGGCTCAGGCGCGTGACCCGCAGCTGGCGGCATGGATTGAAGAAAATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCTTACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTTCGTCAGTGGGTGATTGAAGAt                                                                                                                                                                                                                          <  2:145076/251‑1 (MQ=255)                    gATGTCCTGCGATAACGTGCGTGAAAACGGTCATGTGGCGAAGGTCGCGGTACTGGGGCTGGCTCAGGCGCGTGACCCGCAGCTGGCGGCATGGATTGAAGAAAATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCTTACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTTCGTCAGTGGGTGATTGAAGATAAc                                                                                                                                                                                                                       <  1:4362/250‑1 (MQ=255)                    gATGTCCTGCGATAACGTGCGTGAAAACGGTCATGTGGCGAAGGTCGCGGTACTGGGGCTGGCTCAGGCGCGTGACCCGCAGCTGGCGGCATGGATTGAAGAAAATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCTTACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTTCGTCAGTGGGTGATTGAAGATAACt                                                                                                                                                                                                                      <  2:111443/251‑1 (MQ=255)                             cGATAACGTGCGTGAAAACGGTCATGTGGCGAAGGTCGCGGTACTGGGGCTGGCTCAGGCGCGTGACCCGCAGCTGGCGGCATGGATTGAAGAAAATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCTTACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTTCGCCAGTGGGGGATTGAAGATAACTTCGTTAATg                                                                                                                                                                                                             >  2:310703/1‑251 (MQ=255)                               ataaagggcgggAAAACGGTCATGTGGCGAAGGTCGCGGTACGGGGGCAGGCTCAGGCGCGTGACCCGCAGCTGGCGGCATGGATTGAAGAAAATGGCACCTTCCCGTGCACCATGGTTGACCGCAACGTTCCGGCGGCGACGCCAGAAACCTTACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTGCGTCAGTGGGTGATTGAAGATAACTTCGTTAATggg                                                                                                                                                                                                           <  2:108928/240‑2 (MQ=255)                               aTAACGTGCGTGAAAACGGTCATGTGGCGAAGGTCGCGGTACTGGGGCTGGCTCAGGCGCGTGACCCGCAGCTGGCGGCATGGATTGAAGAAAATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCTTACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGATCCGTTTCGTCAGTGGGTGATTGAAGATAACTTCGTTAATgg                                                                                                                                                                                                            <  2:338883/250‑1 (MQ=255)                               aTAACGTGCGTGAAAACGGTCATGTGGCGAAGGTCGCGGTACTGGGGCTGGCTCAGGCGCGTGACCCGCAGCTGGCGGCATGGATTGAAGAAAATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCTTACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTTCGTCAGTGGGTGATTGAAGATAACTTCGTTAATgg                                                                                                                                                                                                            >  1:338883/1‑250 (MQ=255)                               aTAACGTGCGTGAAAACGGTCATGTGGCGAAGGTCGCGGTACTGGGGCTGGCTCAGGCGCGTGACCCGCAGCTGGCGGCATGGATTGAAGAAAATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCTTACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTTCGTCAGTGGGTGATTGAAGATAACTGCGTTAATGGt                                                                                                                                                                                                           <  1:342670/251‑1 (MQ=255)                                 aaCGTGCGTGAAAACGGTCATGTGGCGAAGGTCGCGGTACTGGGGCTGGCTCAGGCGCGTGTCCCGCAGCTGGCGGCATGGATTGAAGAAAATGTCACCTTCCCGTGCACCATGGTTTACCGCATCGTTCCGGCGGCGACGCCAGAAACCTTACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTTCGTCAGTGGGTGATTGAAGATAACTTCGTGAATGGTc                                                                                                                                                                                                          <  1:345245/250‑1 (MQ=255)                                 aaCGTGCGGGAAAACGGGCATGTGGCGAAGGTCGCGGTACTGGGGCTGGCTCAGGCGCGAGACCCGCAGCTGGCGGCCTGGATTGAAGAAAATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACATTACAGGCAATTGCTGAACAGCTGGGAGTTTACGACCAGTGCGCCAATGCCTGCGAACCGTTTCGTCAGTGGGTGCGTGAAGGTAACTTCGTTAATGATCg                                                                                                                                                                                                         >  2:11813/1‑251 (MQ=255)                                       cgtgAAAACGGTCATGTGGCGAAGGTCGCGGTACTGGGGCTGGCTCAGGCGCGTGACCCGCAGCTGGCGGCATGGATTGAAGAAAATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCTTACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTTCGTCAGTGGGTGATTGAAGATAACTTCGTTAATGGTCGCCCgg                                                                                                                                                                                                    >  1:54847/1‑250 (MQ=255)                                               cGGTCATGTGGCGAAGGTCGCGGTACTGGGGCTGGCTCAGGCGCGTGACCAGCAGCTGGCGGCATGGATTGAAGAAACTGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGCCGCCAGAAACCTTACAGGAAATTGCTGACCAGCTGGGAGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTTCGTCAGTGGGTGATTGAAGATAACTTCGTTAATGGTCGCCCGGATTGGGAt                                                                                                                                                                                            >  1:180340/1‑250 (MQ=255)                                                 gTCATGTGGCGAAGGTCGCGGTACTGGGGCTGGCTCAGGCGCGTGACCCGCAGCTGGCGGCCTGGATTGAAGAAAATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCTTACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTTCGTCAGGGGGTGATTGCAGATAACTTCGTTAATGGTCGCCCGGATTGGGATaa                                                                                                                                                                                          >  2:126848/1‑250 (MQ=255)                                                        ggCGAAGGTCGCGGTACTGGGGCTGGCTCAGGCGCGTGACCCGCAGCTGGCGGCATGGATTGAAGAAAATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCTTACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTTCGTCAGTGGGTGATTGAAGATAACTTCGTTAATGGTCGCCCGGATTGGGATAAAGTGGgcg                                                                                                                                                                                  <  2:103990/251‑1 (MQ=255)                                                           gAAGGTCGCGGTACTGGGGCTGGCTCAGGCGCGTGACCCGCAGCTGGCGGCATGGATTGAAGAAAATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCTTACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTTCGTCAGTGGGTGATTGAAGATAACTTCGTTAATGGTCGCCCGGATTGGGATAAAGTGGGCGcac                                                                                                                                                                               >  1:108745/1‑251 (MQ=255)                                                            aaGGTCGCGGAACTGGGGCTGGCTCAGGCGCGTGAACCGCAGCTGGCGACATGGATTGAAGAAAATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCTTACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTTCGTCAGTGGGTGATTGACGATAACTTCGTTAATGGTCGCCCGGATTGGGATAAAGTGGGCGcac                                                                                                                                                                               >  1:241708/1‑250 (MQ=255)                                                                 cgcgGTACTGGGGCTGGCTCAGGCGCGTGACCCGCAGCTGGCGGCATGGATTGAAGAAAATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCTTACAGGAACTTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCAGTGCCTGCGAACCGTTTCGCCAGTGGGTGATTGAAGATAACTTCGTTAATGGTCGCCCGGATTGGGATAAAGTGGGCGCACAGTTCg                                                                                                                                                                         >  2:339010/1‑251 (MQ=255)                                                                 cgcgGTACTGGGGCTGGCTCAGGCGCGTGACCCGCAGCTGGCGGCATGGATTGAAGAAAATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCTTACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTTCGTCAGTGGGTGATTGAAGATAACGTCGTTAATGGTCGCCCGGATTGGGATAAAGTGGTCGCACAGTTCg                                                                                                                                                                         <  1:336353/251‑1 (MQ=255)                                                                    ggTACTGGGGCTGGCTCAGGCGCGTGACCCGCAGCTGGCGGCATGGATTGAAGAAAATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCTTACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTTTGTCAGTGGGTGATTGAAGATAACTTCGTTAATGGTCGCCCGGATTGGGATAAAGTGGGCGCACAGTTCGTTg                                                                                                                                                                      >  1:180882/1‑251 (MQ=255)                                                                       aCTGGGGCTGGCTCAGGCGCGTGACCCGCAGCTGGCGGCATGGATTGAAGAAAATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCTTACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTTCGTCAGTGGGTGATTGAAGATAACTTCGTTAATGGTCGCCCGGATTGGGATAAAGTGGGCGCACAGTTCGTTGCAg                                                                                                                                                                   <  1:322893/251‑1 (MQ=255)                                                                            ggctggctCAGGCGCGTGACCCGCAGCTGGCGGCATGGATTGAAGAAAATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCTTACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTTCGTCAGTGGGTGATTGAAGATAACTTCGTTAATGGTCGCCCGGATTGGGATAAAGTGGGCGCACAGTTCGTTGCAGACGtt                                                                                                                                                              <  1:333260/251‑1 (MQ=255)                                                                                          cgTGACCAGCAGCTGGCGGCATGGATTGAAGAAAATGTAACCTTCCCGGGCACCATGGATGACCGCATCGTTCCGGCGGCGGCGCCAGAAACCTTACAGGAACTTGCTGAACAGCTGGGTGTTTACGCCCCGTGCGCAATTGCCTGCGAACCGTTTCGTCAGTGGTAGTTTGAAGATAACTTAGTTATTGGTCGCCCGGATTGGGATAAAGTGGGAGCACAGTTCGTTGCAGACGtagggacgttggcaat                                                                                                                                                >  2:246557/1‑236 (MQ=255)                                                                                           gTGACCCGCAGCTGGCGGCATGGATTGAAGAACATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGAGACGCCAGCAACCTTACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTTCGTCAGTGGGTGCTTGAAGATCACTTCGTTAATGGTCGCCCGGATTGGGATAAAGTGGGCGCACAGTTCGTTGCAGACGTTGTGCCGTTCGAAatg                                                                                                                                               >  2:48584/1‑251 (MQ=255)                                                                                           gTGACCCGCAGCTGGCGGCATGGATTGAAGAAAATGCCACCTTCCCGTGCACCATGGTTGCCCGCATCGCTCCGGCGGCGACGCCAGAAACCTTACAGGAAACTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTTCGTCAGTGGGGGATTGACGATAACTTCGTTAATAGACGCCCGGATTGGGATAAAGTGGGCGCACAGTTCGTTGCAGACGTTGTGCCGTTCGAAatg                                                                                                                                               >  2:384284/1‑251 (MQ=255)                                                                                             gACCCGCAGCTGGCGGCATGGATTGAAGACAATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCTTACAGGAAATGGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTTCGTCAGTGGGTGATTGAAGATAACTTCGTTAATGGTCGCCCGGGTTGGGATAAAGTGGGCGCACAGTTCGTTGCAGACGTTGGGCCGTTCGAAatgat                                                                                                                                             >  1:82432/1‑251 (MQ=255)                                                                                                          cggcATGGATTGAAGAAAATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCTTACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTTCGTCAGTGGGTGATTGAAGATAACTTCGTTAATGGTCGCCCGGATTGGGATAAAGTGGGCGCACAGTTCGTTGCAGACGTTGTGCCGTTCGCAATGATGAAGCTGCGTATg                                                                                                                                >  2:152693/1‑251 (MQ=255)                                                                                                          cggcATGGAGTGAAGAAAATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCTTACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTTCGTCAGTGGGTGATTGAAGATAACTTCGTTAATGGTCGCCCGGATTGGGATAAAGTGGGCGCACAGTTCGTTGCAGACGTTGTGCCGTTCGAAATGATGAAGCTGCGTATg                                                                                                                                <  2:290374/251‑1 (MQ=255)                                                                                                           ggcATGGATTGAAGAAAATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCTTACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTTCGTCAGTGGGGGATTGAAGATAACTTCGTTAATGGTCGCCCGGATTGGGATAAAGTGGGCGCACAGTTCGTTGCAGACGTTGTGCCGTTCGAAATGATGAAGCTGCGTATGc                                                                                                                               >  2:89528/1‑251 (MQ=255)                                                                                                              aTGGATTGAAGAAAATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCTTACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTTCGTCAGTGGGTGATTGAAGATAACTTCGTTAATGGTCGCCCGGATTGGGATAAAGTGGGAGCACAGTTCGTTGCAGACGTTGTGCCGTTCGAAATGATGAAGCTGCGTATGCTg                                                                                                                             >  1:18109/1‑250 (MQ=255)                                                                                                                  aTTGAAGCAAATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCTTACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTTCGTCAGTGGGTGATTGAAGATAACTTCGTTAATGGTCGACCGGATTGGGATAAAGTGGGCGCACAGTTCGTTGCAGACGTTGTGACGTACGAAATGATGAAGCTGCGTATGCTGAACgg                                                                                                                        >  1:295656/1‑251 (MQ=255)                                                                                                                     gaagaaAATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCTTACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTTCGTCAGTGGGTGATTGAAGATAACTTCGTTAATGGTCGCCCGGATTGGGATAAAGTGGGCGCACAGTTCGTTGCAGACGTTGTGCCGTTCGAAATGATGAAGCTGCGTATGCTGAACGGc                                                                                                                       >  1:56339/1‑249 (MQ=255)                                                                                                                          aaaTGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGAGACGCCAGAAACCTTACAGGAAATTGCTGACCAGTTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTTCGTCAGTGGGTGATTGAAGATAACTTCGTTAATGGTCGCCCGGATTGGGATAAAGTGGGCGCACAGTTCGTTGCAGACGTTGTGCCGTTCGAAATGATGAAGCTGCGTATGCTGCACGGCAGCCAct                                                                                                                <  2:7365/251‑1 (MQ=255)                                                                                                                           aaTGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCTTACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTTCGTCAGTGGGTGATTGACGATAACTTCGTTAATGGTCGCCCGGATTGGGATAAAGTGGGCGCACAGTTCGTTGCAGACGTTGTGCCGTTCGAAATGATGAAGCTGCGTATGCTGAACGGCAGCCCctc                                                                                                               >  2:148847/1‑251 (MQ=255)                                                                                                                                cACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCTTACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTTCGTCAGTGGGTGATTGAAGATAACTTCGTTAATGGTCGCCCGGATTGGGATAAAGTGGGCGCACAGTTCGTTGCAGACGTTGTGCCGTTCGAAATGATGAAGCTGCGTATGCTGAACGGCAGCCACTCTTTTc                                                                                                          >  1:259491/1‑251 (MQ=255)                                                                                                                                      cccGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCTTACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGACATTGCCTGCGAACCGTTTCGTCCGTGGGTGATTGAAGATAACTTCGTTAATGGTCGCCCGGATTGGGATAAAGTGGGCGCACAGTTCGTTGCAGCCGTTGTGCCGTTCGAAATGATGAAGCTGCGTATGCTGAACGGCAGCCACTCTTTTCTGGCGt                                                                                                    >  1:20001/1‑251 (MQ=255)                                                                                                                                       ccGTGCACCATGGTTGACCGCAGCGTTCCGGCGGCGACGCCAGAAACCTTACAGGAAATTGCTGACCAGCTGGGTGTTTACTACCCGTGCGCCATTGCCTGCGAACCGTTTCGTCAGTGGGTGATTGAAGATAACTGCGTTAATGGTCGCCCGGATTGGGATAAAGTGGGCGCACAGTTCGTTGCAGACGTTGTGCCGTTCGAAATGATGAAGCTGCGTATGCTGAACGGCAGCCACTCTTTTCTGGCGTa                                                                                                   <  2:83591/251‑1 (MQ=255)                                                                                                                                             aCCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCTTCCAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTTCGTCAGGGGGTGGTTGCAGCTAACTTCGTTAATGGTCGCCCGGATTGGGATAATGTGGGCGCACAGTTCGTTGCAGACGTTGTGCCGTTCGAAATGATGCAGCTGCGTATGCTGAACGGCAGCCACTCTTTTCTGGCGtcgctcgg                                                                                             >  2:83203/1‑244 (MQ=255)                                                                                                                                                   gTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCTTACAGGAAATTGCTGACCAGCTGGGCGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTTCGTCAGGGGGTGATTGCAGATAACTTCGTTAATGGTCGCCCGGATTGGGATAAAGTGGGCGCACAGTTCGTTGCAGACGATGTGCCGTTCGAAATGATGAAGCTGCGTATGCTGACCGGCAGCCACTCTTTTCTGGCGGACCTCGGTTACCt                                                                                       >  2:307662/1‑251 (MQ=255)                                                                                                                                                           cATCGTTCCGGCGGCGACGCCAGAAACCTTACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTTCGTCAGTGGGTGATTGAAGATAACTTCGTTAATGGTCGCCCGGATTGGGATAAAGTGGGCGCACAGTTCGTTGCAGACGTTGTGCCGTTCGAAATGATGAAGCTGCGTATGCTGAACGGCAGCCACTCTTTTCTGGCGTACCTCGGTTACCTCGGCGGCt                                                                               >  1:51808/1‑251 (MQ=255)                                                                                                                                                              cGTTCCGGCGGCGACGCCAGAAACCTTACAGGAAAGTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTTCGTCAGTGGGTGATTGAAGATAACTTCGTTAATGGTCGCCCGGATTGGGATAAAGTGGGCGCACAGTTCGTTGCAGACGTTGTGCCGTTCGAAATGATGAAGCTGCGTATGCTGAACGGCAGCCACTCTTTTCTGGCGTACCTCGGTTACCTCGGCGGCTAt                                                                             <  2:56339/250‑1 (MQ=255)                                                                                                                                                                   cggcggCGACGCCAGAAACCTTACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTTCGTCAGTGGGTGATTGAAGATAACTTCGTTAATGGTCGCCCGGATTGGGATAAAGTGGGCGCACAGTTCGTTGCAGACGTTGTGCCGTTCGAAATGATGAAGCTGCGTATGCTGAACGGCAGCCACTCTTTTCTGGCGTCCCTCGGTTACCTCGGCGGCTATGAAAcc                                                                       >  2:331104/1‑251 (MQ=255)                                                                                                                                                                         cgACGCAAGAAACCTTACAGGAAATTGCTGACCAGCTGGGTGTTTACGACTCGTGCGCCATTGCCTGCGAACCGTTTCGTCAGTGGGTGATTGAAGATAACTTCGTTAATGGTCGCCCGGATTGGGATAAAGTGGGCGCACAGTTCGTTGCAGACGTTGTGCCGTTCGAAATGATGAAGCTGCGTATGCTGAACGGCAGCCACTCTTTTCTGGCGTACCTCGGTTACCTCGGCGGCTATGAAACCATTGcc                                                                 <  2:31496/251‑1 (MQ=255)                                                                                                                                                                                  aaaCCTTACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTTCGTCAGTGGGTGATTGAAGCTAACTTCGTTAATGGTCGCCCGGATTGGGATAAAGTGGGCGCACAGTTCGTTGCAGACGTTGTGCCGTTCGAAATGATGAAGCTGCGTATGCTGAACGGCAGCCACTCTTTTCTGGCGTACCTCGGTTACCTCGGCGGCTATGAAACCATTGCCGACACCGTg                                                        >  2:107091/1‑251 (MQ=255)                                                                                                                                                                                              taaTTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTTCGTCAGTGGGTGATTGAAGATAACTTCGTTAATGGTCGCCCGGATTGGGATAAAGTGGGCGCACAGTTCGTTGCAGACGGTGTGCCGTTCGAAATGATGAAGCTGCGTATGCTGAACGGCAGCCACGCTTTTCTGGCGTACCTCGGTTACCTCGGCGGATCTGACACCATTGCCGGCACCGAGACTAACCCGGCt                                            >  2:224089/2‑251 (MQ=255)                                                                                                                                                                                              aaaTTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTTCGTCAGTGGGTGATTGAAGATAACTTCGTTAATGGTCGCCCGGATTGGGATAAAGTGGGCGCACAGTTCGTTGCAGACGTTGTGCCGTTCGAAATGATGAAGCTGCGTATGCTGAACGGCAGCCACTCTTTTCTGGCGTACCTCGGTTACCTCGGCGGCTATGAAACCATTGCCGACACCGTGACTAACCCGGCt                                            <  2:221162/251‑1 (MQ=255)                                                                                                                                                                                               aaTTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTTCGTCAGTGGGTGATTGAAGATAACTTCGTTAATGGTCGCCCGGATTGGGATAAAGTGGGCGCACAGTTCGTTGCAGACGTTGTGCCGTTCGAAATGATGAAGCTGCGTATGCTGAACGGCAGCCACTCTTTTCTGGCGTACCTCGGTTACCTCGGCGGCTATGAAACCATTGCCGACACCGTGACTAACCCGGCtt                                           <  1:100760/251‑1 (MQ=255)                                                                                                                                                                                                  tgcggACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTGTCGTCAGTGGGTGATGGAAGATATCTTCGTTAATGGTCGCCCGGATTGGGATAAAGTGGGCGCACAGTTCGTTGCAGACGTTGTGCCGTTCGAAATGATGAAGCTGCGTATGCTGAACGTCAGCCACTCTTTTCTGGCGTACCTCGGTTACCTCGGCGGCTATGAAACCATTGCCTACTTCGTGACTAATCCGGCTTAt                                         <  2:239692/246‑1 (MQ=255)                                                                                                                                                                                                       aCCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTTCGTCAGTGGGTGATTGAAGATAACTTCGTTAATGGTCGCCCGGATTGGGATAAAGTGGGCGCACAGTTCGTTGCAGACGTTGTGCCGTTCGAAATGATGAAGCTGCGTATGCTGAACGGCAGCCACTCTTTTCTGGCGTACCTCGGTTACCTCGGCGGCTATGAAACCATTGCCGACACCGTGACTAACCCGGCTTATCGCaa                                    >  1:312248/1‑250 (MQ=255)                                                                                                                                                                                                       aCCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGATTCGTCAGTGGGTGATTGAAGATAACTTCGTTAATGGTCGCCCGGATTGGGATAAAGTGGGCGCACAGTTCGTTGCAGACGTTGTGCCGTTCGAAATGATGAAGCTGCGTCTGCTGAACGGCAGCCACTCTTTTCTGGCGTACCTCGGTTACCTCGGCGGCTATGAAACCATTGCCGACACCGTGACTAACCCGGCTTATCGCaa                                    >  1:330833/1‑250 (MQ=255)                                                                                                                                                                                                          aGCTGGGCGTTAACGACCCGTGCGCCATTGCCTGCCAACCGCTTCGTCAGTAGGTGATTGAAGATAACTTCGTTAATGGACGCCCGGATTGGGATCAAGTGGGCGCACAGTTCGTTGCAGACGTTGTGCCGTTCGACATGATGAAGCTGCGTATGCTGAACGGCAGCCACTCTTTTCTGGCGTACCTCGGTTACCTAGGCGGCTATGAAACCATTGCCGAAACCGGGACTACCCCGGCTCATCGCAAAGCg                                >  2:30289/1‑251 (MQ=255)                                                                                                                                                                                                                  gtTTACGACCCGTGCGCCATTGCCTGCGAACCGTTTCGTCAGTGGGTGATTGAAGATAACTTCGTTAATGGTCGCCCGGATTGGGATAAAGTGGGCGCACAGTTCGTTGCAGACGTTGTGCCGTTCGAAATGATGAAGCTGCGTATGCTGAACGGCAGCCACTCTTTTCTGGCGTACCTCGGTTACCTCGGCGGCTATGAAACCATTGCCGACACCGTGACTAACCCGGCTTATCGCAAATCGGCCTTTGc                        <  1:305331/251‑1 (MQ=255)                                                                                                                                                                                                                         aCCCGTGCGCCATTGCCTGCGAACCGTTTCGTCAGTGGGTGATTGAAGATAACTTCGTTAATGGTCGCCCGGCTTGGGATAAAGCGGGCGCACAGTTCGTTGCAGACGTTGTGCCGTTCGCAATGCTGAAGCTGCGTATGCTGAACGGCAGCCACTCTTTTCTGGCGTACCTCGGTTACCTCGGCGGCTATGAAACCATTGCCGACACCGTGACTAACCCGGCTTATCGCAAAGCGGCCTTTGCCCtgatg                 >  1:400954/1‑251 (MQ=255)                                                                                                                                                                                                                                cgcCATTGCCTGCGAACCGTTTCGTCAGTGGGTGATTGAAGATAACTTCGTTAATGGTCGCCCGGATTGGGATAAAGTGGGCGCACAGTTCGTTGCAGACGTTGTGCCGTTCGAAATGATGAAGCTGCGTATGCTGAACGGCAGCCACTCTTTTCTGGCGTACCTCGGTTACCTCGGCGGCTATGAAACCATTGCCGACACCGTGACTAACCCGGCTTATCGCAAAGCGGCCTTTGCCCTGATGATGCAgg          <  1:152693/251‑1 (MQ=255)                                                                                                                                                                                                                                 gcCATTGCCTGCGAACCGTTTCGTCAGTGGGTGATTGAAGATAACTTCGTTAATGGTCGCCCGGATTGGGATAAAGTGGGCGCACAGTTCGTTGCAGACGTTGTGCCGTTCGAAATGATGAAGCTGCGTATGCTGAACGGCAGCCACTCTTTTCTGGCGTACCTCGGTTACCTCGGCGGCTATGAAACCATTGCCGACACCGTGACTAACCCGGCTTATCGCAAAGCTGCCTTTGCCCAGATGATGCAGGa         <  2:330833/251‑1 (MQ=255)                                                                                                                                                                                                                                    aTTGCCTGCGAACCGTTTCGTCAGTGGGTGATTGAAGATAACTTCGTTAATGGTCGCCCGGATTGGGATAAAGTGGGCGCACAGTTCGTTGCAGACGTTGTGCCGTTCGAAATGATGAAGCTGCGTATGCTGAACGGCAGCCACTCTTTTCTGGCGTACCTCGGTTACCTCGGCGGCTATGAAACCATTGCCGACACCGTGACTAACCCGGCTTATCGCAAAGCGGCCTTTGCCCTGATGCTGCAGGAAc       >  2:222162/1‑250 (MQ=255)                                                                                                                                                                                                                                      tgacggCGAACCGTTTCGTCAGGGGGTGATTGAAGATAACTTCGTTAATGGTCGCCCGGATTGGGATAAAGTGGGCGCACAGTTCGTTGCAGACGTTGTGCCGTTCGAAATGATGAAGCTGCGTATGCTGAACGGCAGCCACTCTTTTCTGGCGTACCTCGGTTACCTCGGCGGCTATGAAACCATTGCCGACACCGTGACTAACCCGGCTTATCGCAAAGCGGCCTTTGCCCTGATGATGCAGGAACaa     <  2:312618/245‑1 (MQ=255)                                                                                                                                                                                                                                       gCCTGCGAACCGTTTCGTCAGTGGGTGATTGAAGATAACTTCGTTAATGGTCGCCCGGATTGGGATAAAGTGGGCGCACAGTTCGTTGCAGACGTTGTGCCGTTCGAAATGATGAAGCTGCGTATGCTGAACGGCAGCCACTCTTTTCTGGCGTACCTCGGTTACCTCGGCGGCTATGAAACCATTGCCGACACCGTGACTAACCCGGCTTATCGCAAAGCGGCCTTTGCCCTGATGATGCAGGAACAAg    <  2:313599/250‑1 (MQ=255)                                                                                                                                                                                                                                        ccTGCGAACCGTTTCGTCAGTGGGTGATTGAAGATAACTTCGTTAATGGTCGCCCGGATTGGGATAAAGTGGGCGCACAGTTCGTTGCAGACGTTGTGCCGTTCGAAATGATGAAGCTGCGTATGCTGAACGGCAGCCACTCTTTTCTGGCGTACCTCGGTTACCTCGGAGGCTATGAAACCATTGCCGACACCGTGACTAACCCGGCTTATCGCAAAGCGGCCTTTGCCCTGATGATGCAGGAACAAgcg  >  1:4232/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                              GGTTGAAAGCGTTTACGGTGATGTCCTGCGATAACGTGCGTGAAAACGGTCATGTGGCGAAGGTCGCGGTACTGGGGCTGGCTCAGGCGCGTGACCCGCAGCTGGCGGCATGGATTGAAGAAAATGTCACCTTCCCGTGCACCATGGTTGACCGCATCGTTCCGGCGGCGACGCCAGAAACCTTACAGGAAATTGCTGACCAGCTGGGTGTTTACGACCCGTGCGCCATTGCCTGCGAACCGTTCCGTCAGTGGGTGATTGAAGATAACTTCGTTAATGGTCGCCCGGATTGGGATAAAGTGGGCGCACAGTTCGTTGCAGACGTTGTGCCGTTCGAAATGATGAAGCTGCGTATGCTGAACGGCAGCCACTCTTTTCTGGCGTACCTCGGTTACCTCGGCGGCTATGAAACCATTGCCGACACCGTGACTAACCCGGCTTATCGCAAAGCGGCCTTTGCCCTGATGATGCAGGAACAAGCG  >  NC_000913/4553427‑4553908

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