Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,904,907	C→T	E231K (GAA→AAA)	bglF ←	beta‑glucoside specific PTS enzyme II/BglG kinase/BglG phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,904,907	0	C	T	98.2%	176.1 / ‑3.7	55	E231K (GAA→AAA)	bglF	beta‑glucoside specific PTS enzyme II/BglG kinase/BglG phosphatase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (28/26);  total (28/27)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.68e-01

CACAGCGGCACCAGGCCCCAGTGCAGTCCGAACATGACGAAGATTTGCCAGAAGCCGCCCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGC  >  NC_000913/3904660‑3905150                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    cacaGCGGCACCAGGCCCCAGTGCAGTCCGAACATGACGAAGAGTTGCCAGAAGCCGCCCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAg                                                                                                                                                                                                                                                  <  2:224494/251‑1 (MQ=255)            cAGGCCCCAGTGCAGTCCGAACATGACGAAGATTTGCCAGAAGCCGCCCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGcacc                                                                                                                                                                                                                                       >  1:245339/1‑250 (MQ=255)             aGGCCCCAGTGCAGTCCGAACATGTCGAAGATTTGCCAGAAGCCGCCCATTTCCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGTTTGATAGCGGCCCCACCATCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATTGTGTGAAGAAATTTTTGATTGCCGATGGTAACCGCGCATTAAGTCGGCGTTTCAGAATGCTGCACaa                                                                                                                                                                                                                                      <  1:432349/251‑1 (MQ=255)             aGGCCCCAGTGCAGTCCGAACATGCCGAAGATTTGCCAGAAGCCGCCCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACAACGAATTAAGTCGGCGTTTCAGATTGATGcaca                                                                                                                                                                                                                                       >  2:13169/1‑250 (MQ=255)                ccccAGTGCAGTCCGAACATGACGAAGATTTGCCAGAAGCCGCCCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCa                                                                                                                                                                                                                                   <  2:547457/251‑1 (MQ=255)                ccccAGTGCAGTCCGAACATGACGAAGATGTGCCAGAAGCCGCCCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCa                                                                                                                                                                                                                                   <  1:408158/251‑1 (MQ=255)                      tGCAGTCCGAACATGACGAAGCTTTGCCAGAAGCCGCCCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAATTCGGCGTTTCAGAATGCTGCACAACCAGGc                                                                                                                                                                                                                                >  1:457034/1‑248 (MQ=255)                         aGTCCGAACATGACGAAGATTTGCCAGAAGCCGCCCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGaac                                                                                                                                                                                                                           >  2:183083/1‑249 (MQ=255)                             cGAACATGACGAAGATTTGCCAGAAGCCGCCCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAaat                                                                                                                                                                                                                       >  1:358770/1‑250 (MQ=255)                                aCATGACGAAGATTTGCCAGAAGCCGCCCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAAACGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATg                                                                                                                                                                                                                   >  1:377730/1‑251 (MQ=255)                                 aaTGACGAAGATTTGCCAGAAGCCGCCCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCTACGGTAACCACGCATTAAGTCGGCGGTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATgg                                                                                                                                                                                                                  <  2:87115/250‑1 (MQ=255)                                    gACGAAGATTTGCCAGAAGCCGCCCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGaa                                                                                                                                                                                                               <  1:2407/251‑1 (MQ=255)                                       gAAGATTTGCCAGAAGCCGCCCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATa                                                                                                                                                                                                             >  2:163716/1‑250 (MQ=255)                                        aaGATTTGCCAGAAGCCGCCCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAAc                                                                                                                                                                                                           <  1:349095/251‑1 (MQ=255)                                           aTTTGCCAGAAGCCGCCCAGTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTGCGCTTATCCAGGTTGATAGCGGCCCCTCCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGa                                                                                                                                                                                                        <  2:46997/251‑1 (MQ=255)                                                           ccATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTCTCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTcaac                                                                                                                                                                                        >  1:254993/1‑251 (MQ=255)                                                                     gcCCGCAAATGCAGGAACCGCCTTATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGCGACCg                                                                                                                                                                              <  1:99745/251‑1 (MQ=255)                                                                          cAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATa                                                                                                                                                                         <  1:373061/251‑1 (MQ=255)                                                                                          cTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGcc                                                                                                                                                         >  2:229616/1‑251 (MQ=255)                                                                                           tGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGccc                                                                                                                                                        >  2:499799/1‑251 (MQ=255)                                                                                                aaGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAgcg                                                                                                                                                   <  1:496083/251‑1 (MQ=255)                                                                                                   ccAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCAt                                                                                                                                                >  1:191486/1‑251 (MQ=255)                                                                                                       agaTAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGc                                                                                                                                            >  1:278043/1‑251 (MQ=255)                                                                                                                 cggcAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCg                                                                                                                                  <  2:215103/251‑1 (MQ=255)                                                                                                                         aGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGAGTGAAGAAATTTTTGATTGCCGACGGTAACCAAGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGa                                                                                                                           >  1:388389/1‑250 (MQ=255)                                                                                                                              gCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGc                                                                                                                       >  1:63572/1‑249 (MQ=255)                                                                                                                              gCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAg                                                                                                                     <  1:521498/251‑1 (MQ=255)                                                                                                                                 tATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGt                                                                                                                    <  1:238775/249‑1 (MQ=255)                                                                                                                                 tATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTc                                                                                                                   <  1:222015/250‑1 (MQ=255)                                                                                                                                     cAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATACCCATCAGACATAGCAATGGTGCGAAGAACTTTTTGATTGCCGACGGTAAACACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAACAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACACAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAAt                                                                                                              >  1:532798/1‑251 (MQ=255)                                                                                                                                     cAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAAt                                                                                                              >  1:331645/1‑251 (MQ=255)                                                                                                                                        gTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATGGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTaa                                                                                                           <  2:526015/251‑1 (MQ=255)                                                                                                                                            atggCGGCCCCACCAGCAGCAAGGTGACGGGGGTGATAACCATCAGAAATAGCAATGGTATGAAGAAATTCTTGATTGCCGACGGGAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGACGAAAAACTAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCACGTTCTCCCAAGCAGTCAGAATTAATga                                                                                                        >  2:80930/4‑249 (MQ=255)                                                                                                                                                 ggCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCAc                                                                                                  <  2:35838/251‑1 (MQ=255)                                                                                                                                                    cccACCAGCAGAACGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGATCTCGAAAGCAGTCAGAATTAATGGATGCACTac                                                                                               >  2:267970/1‑250 (MQ=255)                                                                                                                                                                            gtgATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAg                                                                        <  2:388389/250‑1 (MQ=255)                                                                                                                                                                            gtgATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGt                                                                       <  2:367074/251‑1 (MQ=255)                                                                                                                                                                                     aTCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGAtt                                                              >  1:368193/1‑251 (MQ=255)                                                                                                                                                                                          aCATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGAGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTAccgc                                                         >  2:563477/1‑251 (MQ=255)                                                                                                                                                                                           cATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTAccgcc                                                        >  2:172682/1‑251 (MQ=255)                                                                                                                                                                                              aGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAGCCGATGACGAGTAATCCAACAATGTGACCGGAATTCCCAGGAATTCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGTGGCGCTAACGCTCCACCAATCCCCATGGCAGTAAATGGATTACAGCCGa                                                      <  2:307067/250‑1 (MQ=255)                                                                                                                                                                                                   tGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGACAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGCAAGCCGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGtt                                                >  2:462768/1‑251 (MQ=255)                                                                                                                                                                                                    ggTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATGACCGCCGAAGCGttt                                               <  1:163716/251‑1 (MQ=255)                                                                                                                                                                                                          aagaaATTTTTGGTTGCCGACGGCAACCACGCATTAAGTCGGCTTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATATCCGATGACGAGTAATGCAACAATTTGTCCGGAATACCCAGGAATGCCAGCCCCAGCGCATCCGTTTTTTGCCCTTTCTCGAAAGCATTCAGAATTAATGGATGTTCTTACTCTCCGCCAGCCACCATGGCAAAAAAGGGATTGCCGCGGAAGCGGTTCCCcgc                                         <  2:66322/251‑1 (MQ=255)                                                                                                                                                                                                          aagaaATTTTTGATTGCCGACGGTAACAACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGTTTCCCcgc                                         >  2:565516/1‑251 (MQ=255)                                                                                                                                                                                                                          ccGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGTTTCCCCGCGGTGTATCCCAGGAtt                         >  1:156365/1‑251 (MQ=255)                                                                                                                                                                                                                             aCGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCCCCATGGCAGTAAATGGATTACCGCCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATc                      <  1:229616/251‑1 (MQ=255)                                                                                                                                                                                                                             aCGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGACGAAGCGGTTCCCCGCGGTGTATCCCAGGATTATc                      >  2:152041/1‑251 (MQ=255)                                                                                                                                                                                                                                gTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCggg                   >  2:532237/1‑251 (MQ=255)                                                                                                                                                                                                                                gTAACCACGCATTAAGTCGGCGTTTCAGAAGGCTGCAAAACCAGGCAAAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCGGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCGCCACCAATCACCATGGCAGGAAAGGGATTACCGCCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCggg                   <  1:267970/251‑1 (MQ=255)                                                                                                                                                                                                                                    ccACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGga                  >  2:90799/1‑248 (MQ=255)                                                                                                                                                                                                                                       cGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAAcc            <  1:517190/251‑1 (MQ=255)                                                                                                                                                                                                                                          aTTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTGTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGTACTAACGCTCCACCAATCACCATTGCAGTAAATGGATTACCGCTGAAGCGTTTCCCCGCGGTGTATCCCAGTATTATCGGGAAGAACCaaa         <  1:299021/251‑1 (MQ=255)                                                                                                                                                                                                                                                 cGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGGAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCCCCAATCACCATGGCAGTAAATGGAGTACCGCCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGc  <  2:245339/251‑1 (MQ=255)                                                                                                                                                                                                                                                  ggCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGc  >  1:121893/1‑250 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    CACAGCGGCACCAGGCCCCAGTGCAGTCCGAACATGACGAAGATTTGCCAGAAGCCGCCCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGGC  >  NC_000913/3904660‑3905150

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