Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,010,883	C→T	G228G (GGC→GGT)	yigL →	phosphosugar phosphatase YigL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,010,883	0	C	T	91.3%	116.4 / ‑5.4	46	G228G (GGC→GGT)	yigL	phosphosugar phosphatase YigL
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (18/24);  minor base G (4/0);  total (22/24)
Fisher's exact test for biased strand distribution p-value = 4.48e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCACCTGCGATTCCCATGAACAACTGCTGCCGCTGGAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGTCAGTTGTCAACCTGATACTTCGCTACAATGGATACCCGTTAATCAAAGAGTTTTCCATTGTGGCGCTACTTATCATCACCACGATTCTGTGGGCCTTCTCCTTTAGC  >  NC_000913/4010645‑4011123                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                 cccccTGCGATTCCCATCAACAACTGCTGCCGCTTGAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTCCCTTAACCGGTCTGGAAGTGATGGCGGGCGGCGGTGCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGTATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGTTGCAGTATGGGc                                                                                                                                                                                                                                      <  1:307687/248‑1 (MQ=255)   aCCTGCGATTCCCATGAACAACCGCTGCCGCTGGAGCAGGCGGTTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCTGCGTTTCAAAAAGCCATGCGCTGGAAGCGGTGGAGAAGAGACTGGGCTCCAACCTGAAGGGTTGGATTGCGTTTGGGGACGGGGTGACCGAAGCCGAAATGCTGGCGATGGCGGGGAAAGGGTGCAATATGGGcgg                                                                                                                                                                                                                                    >  1:109727/1‑249 (MQ=255)         gATTCCCATGAACAACTGCTGCCGCTGGAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGGAAGGTTGCCTTATGGGCAGTgcgc                                                                                                                                                                                                                               >  1:308608/1‑250 (MQ=255)             cccATGAACAACTGCTGCCGCTGGAGCAGGCGCTTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGTTGCCTTATGGGCAGTGCGCAcc                                                                                                                                                                                                                            >  1:294795/1‑249 (MQ=255)                   aacaccTGCTGCCGCTGGAGCAGGCGATTAACGCTCGGTGGGGCGATCGCGTCAACGTCAGTTTCTCTTGCTTAACCTGTCTGTAAGTTAGGGCGGTCGGCGTTTCAAAAGTCCATGCGCTTGAAGCGGTGGCGAATAAACTGTGCTGCAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTTGATGGCGGGGTGAGGTTGCATTATTGTCAGTGCGCTCCAGCGTCTg                                                                                                                                                                                                                    <  2:171614/251‑1 (MQ=255)                             gcCGCTGGAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTTTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGTTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTc                                                                                                                                                                                                          <  1:68572/251‑1 (MQ=255)                                          gCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGGAAGGTTGCATTATTGGCAGTTCGCACCAGCGTCTGAAAGACCTTCCTCCCGCGCTgg                                                                                                                                                                                              >  1:103136/1‑250 (MQ=255)                                                aaCGCGCGTTGGGGCGAGCGCGTCAACGTCAGTTTCTCTTCCTTAACCTGGCTGGAAGTGATGGCGGGCTGCGGTTCAAAAGGCCATGCGCTGGAAGCGGTGTCGAAGAAACTGGGCGACAGCCTGAAGGATTGTATTGCGTGTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGTTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCAACCCGAGCTGGAAGTGAt                                                                                                                                                                                       <  2:89675/251‑1 (MQ=255)                                                           gggCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGTTGCATTATGGGCAGTGCGCACCAGCGTCTGGAAGACCTTCATCCCGAGCTGGGAGTGATTGGTACTAATg                                                                                                                                                                            >  1:19866/1‑251 (MQ=255)                                                            ggCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGGAAGGGTGCATTTTGGGCAGTGCGCACCAGCGTCTGAAAGGCCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGc                                                                                                                                                                           >  1:112597/1‑251 (MQ=255)                                                                  cgcgTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGTTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCcgacga                                                                                                                                                                     <  1:5834/251‑1 (MQ=255)                                                                   gcgTCAACGTCAGGATCTCTCCCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGTGGAAAGGTTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCgacgac                                                                                                                                                                    <  1:53845/251‑1 (MQ=255)                                                                    cgTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCCTGCGCTGGGAGCGGTGGCGAAGGCACTGGGGTACAGCATGAAGGATTGTATTGCGTTTGGGGACGGGATGAACGACGCCGAAATGCTGTCGGTGGCGGGGGAAGGTTGACTTTTGGGCAGTGCGCACCAGCGGCTGAAAGACCTTCATCCCGAGCTGGAAGTGCTTGGTACTAATGCCGccgacg                                                                                                                                                                   >  1:33554/1‑251 (MQ=255)                                                                             cAGTTTCTCTACCTTAACCTGTCTTGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGTTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGc                                                                                                                                                          <  1:252984/251‑1 (MQ=255)                                                                             cAGTTTCTCTACCTTAAACTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGTTGCAGTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCg                                                                                                                                                           <  1:137482/250‑1 (MQ=255)                                                                                        ccTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGCAGGATTGTATTGCGTTTGGTGACGGGGTGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGTTGCCTTATGGGCAGTGCGCACCCGCGGCTGAAAGACCTTCATCCCGCGCTGGAAGTGATTGGGACTAATGCCGACGACGCGGTGCCGCATTATCTGCGt                                                                                                                                               >  1:4945/1‑251 (MQ=255)                                                                                         cTTAACCTGTCTGGAAGTGATGGAGGGCGGCGTTTAAAACCGCCATGAGCTAGAAGCGGAGGCGCAGAAACTGGGCTACAGCCTCGAGGATTGTATTGCGTTTGGAGACGGGATGAACGACGCCGAAATGCTGTCGATGGAGGGGAAAGGGTGCATTTTGGGCAGTGCGCCCCCGCGCCTGAAAGACCTTCATTCCGCGCTGGAAGTGATTGGTACTAATGCCGGAGCCGCGGTGCCGAATTATATGCGt                                                                                                                                               >  2:414117/1‑250 (MQ=255)                                                                                                gtgcTGGAAGTGATGGCGGGCGTCGTTGCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGAGTGTATTGCGTGTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCTGGGAAAGGTTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAATTGATGGGGACTAATGCCGACGACGCGGTGCCGCATTAGCTTCGTAAACTCTa                                                                                                                                       <  2:44450/248‑1 (MQ=255)                                                                                                         gTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGTTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTa                                                                                                                              <  1:296420/251‑1 (MQ=255)                                                                                                          ggATGGCTGGCGTCGTTTCAAAAGTCCATGCGCTGGACGCGGTGGCGTATAATCTGGGCTACAGCCGGAAGGATTGTATAGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGGCGATGTCGTGGAAAGGTTGCATTATGGGCAGTGGACACCAGCGTCTGAAAGACCGTCATCCCGAGCTGGAAGTTACTGGTACTAATGCCGACGAGGCGGTGCCGCAGTATCTGCGTAAACTCTATGTATCGTaa                                                                                                                             <  2:267748/250‑1 (MQ=255)                                                                                                                cgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGTTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTc                                                                                                                       <  1:319993/251‑1 (MQ=255)                                                                                                                        gTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCAGAAATGCTGTCGATGGCGGGGAAAGGTTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTg                                                                                                               >  1:196116/1‑251 (MQ=255)                                                                                                                        gTTTCAAAAGGCCATGCGCTGGAAGCGGCGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGGAAGGTTGCATTATGGGCCGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGGGATTGGTACTAATGCCCACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCCTTATTTg                                                                                                               >  2:69243/1‑251 (MQ=255)                                                                                                                            cAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAATCTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGTTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGTCa                                                                                                           <  1:236721/251‑1 (MQ=255)                                                                                                                                 ggCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTTAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGTTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGTCAGTTGt                                                                                                      <  1:134430/251‑1 (MQ=255)                                                                                                                                                     ggCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGTTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGTCAGTTGTCAACCTGATACTTCGCTAc                                                                                   >  1:227041/1‑250 (MQ=255)                                                                                                                                                          agaaACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGTTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGTCAGTTGTCAACCTGATACTTCGCTACAATGGa                                                                             <  1:100902/251‑1 (MQ=255)                                                                                                                                                                           cTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGCCGCCGAAATGCTGTCGATGGCGGGGAAAGGTTGCATTATGGGCAGTGCGCACCCGCGACTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGAGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGTCCGTTGTCAACCTGATACTTCGCTACAATGGATACCCGTTAATCaa                                                               >  2:376279/1‑248 (MQ=255)                                                                                                                                                                              aaGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGTTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGTCAGTTGTCAACCTGATACTTCGCTACAATGGATACCCGTTAATCAAAGAGtt                                                         <  1:98713/251‑1 (MQ=255)                                                                                                                                                                               aGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGTTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGTCAGTTGTCAACCTGATACTTCGCTACAATGGATACCCGTTAATCAAAGAGttt                                                        >  1:130217/1‑251 (MQ=255)                                                                                                                                                                                ggATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGTTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGTCAGTTGTCAACCTGATACTTCGCTACAATGGATACCCGTTAATCAAAGAGtttt                                                       <  1:184513/251‑1 (MQ=255)                                                                                                                                                                                    tGTATTGCGGTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGTTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGTCAGTTGTCAACCTGATACTTCGCAACAATGGATACCCGTGAATCAAAGAGTTTGCCa                                                    <  1:96418/250‑1 (MQ=255)                                                                                                                                                                                     gTATTGCGTGTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGTTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTAGTTATCGTAATCGTTCTTTATTTGGTCAGTTGTCAACCTGATTCTTCGCTACAATGGATACCCGTTAATCAAAGATTTTTCCAtt                                                  <  2:294795/251‑1 (MQ=255)                                                                                                                                                                                        ttGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGTTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGTCAGTTGTCAACCTGATACTTCGCTACAATGGATACCCGTTAATCAAAGAGTTTTCCATtgtg                                               >  1:327735/1‑251 (MQ=255)                                                                                                                                                                                          gCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGCAAGGTTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGTCAGTTGTCAACCTGATACTTCGCTACAATGGATAc                                                                          >  1:188695/1‑222 (MQ=255)                                                                                                                                                                                          gCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGTTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGTCAGTTGTCAACCTGATACTTCGCTACAATGGATAc                                                                          <  2:188695/222‑1 (MQ=255)                                                                                                                                                                                               gggggACGGGATGAACGAAGCCGAAATGCTGGCGATGGCGGGGAAAGGTTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCAACCCGAGCTGGAAGTGATTGGTACTAATTGCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGTCAGTTGTCAACCTGATACTTCGCTACAATGGATGCCCGTTAATCAAAGAGTTTTCCATTGTGGCGCTa                                         <  2:115127/246‑1 (MQ=255)                                                                                                                                                                                                 gTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGCAAGGTTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGTCAGTTGTCAACCTGATACTTCGCTACAATGGATACCCGTTAATCAAAGAGTTTTCCATTGTGGCGCTACtt                                      >  1:283467/1‑251 (MQ=255)                                                                                                                                                                                                 gTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGTTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCTTAAACTCTATTTATCGTAATCGTTCTTTATTTGGTCAGTTGTCAACCTGATACTTCGCTACAATGGATACCCGTTAATCAAAGAGTTTTCCATTGTGGCGCTACtt                                      <  1:383195/251‑1 (MQ=255)                                                                                                                                                                                                     cGGGATGAACGACGCCGAAATGCTGTCGATGCCGGGGAAAGGGTGCATTATGGGCAGTGCGCACCAGCGTCTGACAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTCATCGTTCTTTATTTGGTCAGTTGTCAACCTGATACTTCGCTACAATGGATACCCGTTAATCAAAGAGTTTTCCATTGTGGCGCTACTTatc                                   >  1:354367/1‑250 (MQ=255)                                                                                                                                                                                                      gggATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGTTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGTCAGTTGTCAACCTGATACTTCGCTACAATGGATACCCGTTAATCAAAGAGTTTTCCATTGTGGCGCTACTTatcat                                 >  1:86375/1‑251 (MQ=255)                                                                                                                                                                                                               gacgCCGAAATGCTGTCGATGGCGGGGAAAGGTTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGTCAGTTGTCAACCTGATACTTCGCTACAATGGATACCCGTTAATCAAAGAGTTTTCCATTGTGGCGCTACTTATCATCACCACGAt                        >  1:318506/1‑251 (MQ=255)                                                                                                                                                                                                                     gAAATGCTGACGATGGCGGGGAAAGGTTGCCTTTTGGGCAGTGCGAACCCGCGTCTGAAAGACCTTCCTCCCGAGATGGACGAGATTGGTACTAATGACGACGACGCGGGGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGTCAGTTGTCAACCTGATACTTCGCTACAATGGATACCCGTTAATCAATGAGTTTTCCCTTGTGGCGCTACTTATCATCACCACGATACtgtg                  >  2:28487/1‑251 (MQ=255)                                                                                                                                                                                                                                   ggCGGGGAAAGGTTGCATTATGGTCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAATCTCTATTTATCGTAATCTTTCTTTATTTGGTCAGTTGTCAACCTGATACTTCGCTACAATGGATTCCCGTTAATCAAAGAGTTTTCCATTGTGGCGCTACTTATCATCACCACGATTCTGTGGGCCTTCTCCTTTa    <  1:266492/251‑1 (MQ=255)                                                                                                                                                                                                                                   ggCGGGGAAAGGTTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGTCAGTTGTCAACCTGATACTTCGCTACAATGGATACCCGTTAATCAAAGAGTTTTCCATTGTGGCGCTACTTATCATCACCACGTTTCTGTGGGCCTTCTCCttt     <  2:408681/250‑1 (MQ=255)                                                                                                                                                                                                                                     cGGGGAAAGGTTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGTCAGTTGTCAACCTGATACTTCGCTACAATGGATACCCGTTAATCAAAGAGTTTTCCATTGTGGCGCTACTTATCATCACCACGATTCTGTGGGCCTTCTCCTTTAGc  >  1:261895/1‑251 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                 TCACCTGCGATTCCCATGAACAACTGCTGCCGCTGGAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTAATGCCGACGACGCGGTGCCGCATTATCTGCGTAAACTCTATTTATCGTAATCGTTCTTTATTTGGTCAGTTGTCAACCTGATACTTCGCTACAATGGATACCCGTTAATCAAAGAGTTTTCCATTGTGGCGCTACTTATCATCACCACGATTCTGTGGGCCTTCTCCTTTAGC  >  NC_000913/4010645‑4011123

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