Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,788,382	A→G	L94L (CTA→CTG)	yibQ →	divergent polysaccharide deacetylase domain‑containing protein YibQ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,788,382	0	A	G	98.5%	215.3 / NA	67	L94L (CTA→CTG)	yibQ	divergent polysaccharide deacetylase domain‑containing protein YibQ
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (27/39);  minor base T (0/1);  total (27/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CGCTCTGTTGGCGCTCTCCTCCCCCGTACTTGCTGGCAAACTTGCCATCGTCATTGATGATTTTGGGTATCGCCCGCACAACGAAAACCAGGTGCTGGCGATGCCTTCCGCTATCTCCGTCGCTGTATTACCCGATTCACCGCACGCCAGAGAGATGGCGACCAAAGCGCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCGTTGAGTAAACAGCCGCTGGAGAAAAATACGCTACGCCCGGAGATGAGCAGCGACGAAATTGAGC‑GCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGCACTGGCGTGAAGGTGATCAAACGGAAGGTA  >  NC_000913/3788136‑3788628                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                        cGCTCTGTTGGCGCTCTCCTCCCCCGTACTTGCTGGCAAACTTGCCATCGTCATTGATGATTTTGGGTATCGCCCGCACAACGAAAACCAGGTGCTGGCGATGCCCTCCGCTATCTCCGTCGCTGTATTACCCGATTCACCGCACGCCAGAGAGATGGCGACCAAAGCGCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGcc                                                                                                                                                                                                                                                     >  2:198229/1‑251 (MQ=255)         tGGCGCTCTCCTCCCCCGTACTTGCTGGCAAACTTGCCATCGTCATTGATGATTTTGGGTATCGCCCGCACAACGAAAACCAGGTGCTGGCGATGCCCTCCGCTATCTCCGTCGCTGTATTACCCGATTCACCGCACGCCAGAGAGATGGCGACCAAAGCGCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGAt                                                                                                                                                                                                                                              <  1:364671/250‑1 (MQ=255)          ggCGCTCTCCTCCCCCGTACTTGCTGGCAAACTTGCCATCGTCATTGATGATTTTGGGTATCGCCCGCACAACGAAAACCAGGTGCTGGCGATGCCCTCCGCTATCTCCGTCGCTGTATTACCCGATTCACCGCACGCCAGAGAGATGGCGACCAAAGCGCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATg                                                                                                                                                                                                                                             <  2:243873/250‑1 (MQ=255)                ctcctcCCCCGTACTTGCTGGCAAACTTGCCATCGTCATTGATGATTTTGGGTATCGCCCGCACAACGAAAACCAGGTGCTGGCGATGCCCTCCGCTATCTCCGTCGCTGTATTACCCGATTCACCGCACGCCAGAGAGATGGCGACCAAAGCGCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGagcag                                                                                                                                                                                                                                        <  1:425942/249‑1 (MQ=255)                            ccTTGCTGGCAAACTTGCCATCGTCATTGATGATTTTGGGTATCGCCCGCACAACGAAAACCAGGTGCTGGCGATGCCCTCCGCTATCTCCGTCGCTGTATTACCCGATTCACCGCACGCCAGAGAGATGGCGACCAAAGCGCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAgc                                                                                                                                                                                                                          <  2:46036/250‑1 (MQ=255)                               ttcTGGCAAACTTGCCATCGTCATTGATGATTTTGGGGATCGCCCGCACAACGAAAACCAGGTGCTGGCTATGCCCTCCGCTATCTCCGTCGCTGTATTACCCGATTCACCGCACGCCAGAGAGATTGCGACCAAAGCGCATAACAGCGGGCACGAAGTGTTGAGTCATCTCCCGATGGCACCACTGAGTAAACAGTCGCTTGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAgc‑g                                                                                                                                                                                                                        <  1:276699/247‑1 (MQ=255)                                           tGCCATCGTCATTGATGATTTTGGGTATCGCCCGCACAACGAAAACCAGGTGCTGGCGATGCCCTCCGCTATCTCCGTCGCTGTATTACCCGATTCACCGCACGCCAGAGAGATGGCGACCAAAGCGCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGC‑GCATTATTCGTAGt                                                                                                                                                                                                           <  1:46219/250‑1 (MQ=255)                                           agCCATCGTCATTGATGATTTTGGGTATCGCCCGCACAACGAAAACCAGGGGCTGGCGATGCCCTCCGCTATCTCCGTCGCTGTATTACCCGATTCACCGCACGCCAGATAGATGGCGACCAAAGCGCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCTCTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGATATTGAGC‑GCATTATTCGTAg                                                                                                                                                                                                            <  1:140410/248‑1 (MQ=255)                                            tccATCGTCATTGCTGATTTTGGGGATCGCCCGCACAACGAAAACCAGGTGCTGGCGATGCCCTCCGCTATCTCCGTCGCTGTATTACCCGATTCACCGCACGCCAGAGAGATGGCGACCAAAGCGCATAACAGCGGGCACGAAGTGGTGATTCATCTCCCGATGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAAGTGAGC‑GCATTATTCGTAGTg                                                                                                                                                                                                          <  2:296613/249‑1 (MQ=255)                                            gCCATCGTCATTGATGATTTTGGGTATCGCCCGCACAACGAAAACCAGGTGCTGGCGATGCCCTCCGCTATCTCCGTCGCTGTATTACCCGATTCACCGCGCGCCAGAGAGATGGCGACCAAAGCGCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGC‑GCATTATTCGTAGTg                                                                                                                                                                                                          <  1:86124/250‑1 (MQ=255)                                            gCCATCGTCATTGATGATTTTGGGTATCGCCCGCACAACGAAAACCAGGTGCTGGCGATGCCCTCCGCTATCTCCGTCGCTGTATTACCCGATTCACCGCACGCCAGAGAGATGGCGACCAAAGCGCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGC‑GCATTATTCGTAGTGc                                                                                                                                                                                                         >  1:285170/1‑251 (MQ=255)                                                   tCATTGATGATTTTGGGTATCGCCCGCACAACGAAAACCAGGTGCTGGCGATGCCCTCCGCTATCTCCGTCGCTGTATTACCCGATTCACCGCACGCCAGAGAGATGGCGACCAAAGCGCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGC‑GCATTATTCGTAGTGCGGTCAAt                                                                                                                                                                                                  >  1:242862/1‑251 (MQ=255)                                                     aTTGATGATTTTGGGTATCGCCCGCACAACGAAAACCAGGTGCTGGCGATGCCCTCCGCTATCTCCGTCGCTGTATTACCCGATTCACCGCACGCCAGAGAGATGGCGACCAAAGCGCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGC‑GCATTATTCGTAGTGCGGTCAATaa                                                                                                                                                                                                <  1:45908/251‑1 (MQ=255)                                                        gctgatGTTGGGTATCTCACTCACAACGAAAACCAGGTGGTGGCGATGCCCTCCGCGATCTCCGTCGCTGTATTACCCGATTCACCGCACGTCAGATAGATGGCGACCAAAGCGCATAACAGCGGGCACGAAGTGGTGATTCATCTCCCTATGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTTAGC‑GCATTATTCGTAGTGCGGTCAATAACg                                                                                                                                                                                              <  2:138582/248‑1 (MQ=255)                                                        gatgatTTTGGGTATCGCCCGCACAACGAAAACCAGGTGCTGGCGATGCCCTCCGCTATCTCCGTCGCTGTATTACCCGATTCACCGCACGCCAGAGAGATGGCGACCAAAGCGCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGC‑GCATTATTCGTAGTGCGGTCAATAACGt                                                                                                                                                                                             >  2:152713/1‑251 (MQ=255)                                                          ggattttggggtgcGCCCGCAAAACGGAAACCAGGTGCGGGGGAGGGCCTACGCCATCTCCGTCGCTGTAGTACCCGATGCACCGGTCGCCAGAGAGATGGCGACCAAAGCGCATAACAGCTGGCACGAAGTGTTGAGTGATCTACCGATGGCACGACTGAGTAAATAGCCGCTGGAGAAAAATACGGTTCGCCGGGATATGAGCAGCGACGAAATTGAGC‑TCATGAGTGGTAGTGCGGTCAAGAACGTGc                                                                                                                                                                                           <  2:438730/238‑1 (MQ=255)                                                           gatTTTGGGTATCGCCCGCACAACGAAAACCAGGTGCTGGCGATGCCCTCCGCTATCTCCGTCGCTGTATTACCCGATTCACCGCACGCCAGCGAGATGGCGACCAAAGCGCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGC‑GCATTATTCGTAGTGCGGTCAATAACGTGcc                                                                                                                                                                                          >  1:45050/1‑251 (MQ=255)                                                           gatTTTGGGTATCGCCCGCACAACGAAAACCAGGTGCTGGCGATGCCCTCCGCTATCTCCGTCGCTGTATTACCCGATTCACCGCACGCCAGAGAGATGGCGACCAACGCGCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGC‑GCATTATTCGTAGTGAGGTCAATAACGTGcc                                                                                                                                                                                          >  2:438860/1‑251 (MQ=255)                                                               ttGGGTGTCTCCCCCACCACGAAAACCAGGTCCTGCCGATGCCCTCCGCTATCTCCGTCGCTGTATTACCCGACTCACCGCACGCCAGACAGATGGCGCCCAAAGCGCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGCTGGCACCACTGCGTAAACAGCCGCTGGAGAACAATACGCTGCGCCCGGAGCTGCGCAGCGACCAAATTG‑GCAGCATTATTCGTAGTGCGGTCCATAACCTGCCCTAt                                                                                                                                                                                      >  2:56351/1‑251 (MQ=255)                                                                 gggTATCGCCCGCACAACGAAAACCAGGTGCTGGCGATGCCCTCCGCTATCTCCGTCGCTGTATTACCCGATTCACCGCACGCCAGAGAGATGGCGACCAAAGCGCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCACTGAGTAAACAGCCGCTGGATAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAAGTGAGC‑GCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGc                                                                                                                                                                                    <  1:171066/251‑1 (MQ=255)                                                                            gCACAACGAAAACCAGGTGCTGGCGATGCCCTCCTCTATCTCCGTCGCTGTATTACCCGATTCACCGCACGCCAGAGAGATGGCGACCAAAGCGCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGC‑GCATTATTCGTAGTGCGGGCAATAACGTGCCCTATGCCGTGGGGATca                                                                                                                                                                         <  1:305710/251‑1 (MQ=255)                                                                                aaCGAAAACCAGGTGCTGGCGATGCCCTCCGCTATCTCCGTCGCTGTATTACCCGATTCACCGCACGCCAGAGAGATGGCGACCAAAGCGCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGCAATTGAGC‑GCATTATTCGTAGTGCGGTCACTAACGTGCCCTATGCCGTGGGGATcaaca                                                                                                                                                                      >  2:357701/1‑250 (MQ=255)                                                                                   gAAAACCAGGTGCTGGCGATGCCCTCCGCTATCTCCGTCGCTGTATTACCCGATTCACCGCACGCCAGAGAGATGGCGACCAAAGCGCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGC‑GCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAAcc                                                                                                                                                                   >  1:197604/1‑250 (MQ=255)                                                                                    aaaaCCAGGTGCTGGCGATGCCCTCCGCTATCTCCGTCGCTGTATTACCCGATTCACCGCACGCCAGAGAGATGGCGACCAAAGCGCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGC‑GCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGCTCAACAACcac                                                                                                                                                                 >  1:384996/1‑251 (MQ=255)                                                                                     aaaCCAGGTTCTGGCGATGCCCTCCGCTATCTCCGTCGCTGTATTACCCGATTCACCGCACGCCAGAGAGATGGCGACCAAAGCGCATAACAGCGGGCACGAATTGTTGATTCATCTCCCGATGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGC‑GCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACcaca                                                                                                                                                                <  1:423312/251‑1 (MQ=255)                                                                                         cAGGTGCTGGCGATGCCCTCCGCTATCTCCGTCGCTGTATTACCCGATTCACCGCACGCCAGAGAGATGGCGACCAAAGCGCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGC‑GCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATggg                                                                                                                                                            >  1:400136/1‑251 (MQ=255)                                                                                               cTGGCGATGCCCTCCGCTATCTCCGTCGCTGTATTACCCGATTCACCGCACGCCAGAGAGATGGCGACCAAAGCGCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGC‑GCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGgtagca                                                                                                                                                       >  2:274655/1‑250 (MQ=255)                                                                                                   cGATGCCCTCCGCTATCTCCGTCGCTGTATTACCCGATTCACCGCACGCCAGAGAGATGGCGACCAAAGCGCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGC‑GCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATg                                                                                                                                                  <  2:329874/251‑1 (MQ=255)                                                                                                    gATGCCCTCCGCTATCTCCGTCGCTGTATTACCCGATTCACCGCACGCCAGAGAGATGGCGACCAAAGCGCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGC‑GCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGTATCAACAAGCACATGGGTAGCAAGATGa                                                                                                                                                 <  2:269375/251‑1 (MQ=255)                                                                                                     aTGCCCTCCGCTATCTCCGGCGCTGTATTACCCGATTCACCGCACGCCAGAGAGATGGCGACCAAAGCGCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGC‑GCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGAc                                                                                                                                                <  2:153591/251‑1 (MQ=255)                                                                                                        ccctCCGCTATCTCCGTCGCTGTATTACCCGATTCACCGCACGCCAGAGAGATGGCGACCAAAGCGCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCACTGAGAAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGC‑GCATTATTCGTAGTGCGGTCAATAACGTGCCATATGCCGTGGGGATCAACAACCACATGGGAAGCAAGATGACctc                                                                                                                                             >  1:96977/4‑251 (MQ=255)                                                                                                         cctCCGCTATCTCCGTCGCTGTATTACCCGATTCACCGCACGCCAGAGAGATGGCGACCAAAGCGCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGC‑GCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACctct                                                                                                                                            <  1:178606/249‑1 (MQ=255)                                                                                                            ccGCTATCTCCGTCGCTGTATTACCCGATTCACCGCACGCCAGAGAGATGGCGACCAAAGCGCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGC‑GCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTaa                                                                                                                                          <  2:171976/250‑1 (MQ=255)                                                                                                                   ctcCGTCGCTGTATTACCCGATTCACCGCACGCCAGAGAGATGGCGACCAAAGCGCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGC‑GCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGCCCTCTAACCTGTTTg                                                                                                                                  >  1:25908/1‑251 (MQ=255)                                                                                                                     ccGTCGCTGTATTACCCGATTCACCGCACGCCAGAGAGATGGCGACCAAAGCGCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGC‑GCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGttt                                                                                                                                   <  1:39872/248‑1 (MQ=255)                                                                                                                               aTTACCCGATTCACCGCACGCCAGAGAGATGGCGACCAAAGCGCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCACTTAGTAAACAGCCGCTGGAGAAAAAGACGCTGCGCCCGGAGAt                                                                                                                                                                                                                                              <  1:145288/132‑1 (MQ=255)                                                                                                                               aTTACCCCATTCACCGCACGCCAGAGAGATGGAGACCAAAGCGCATAACAGCGGGCACGACGTGTTGGTTCATCTCCCGATGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGAt                                                                                                                                                                                                                                              >  2:145288/1‑132 (MQ=255)                                                                                                                                    ccGATTCACCGCACGCCAGAGAGATGGCGACCAAAGCGCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGC‑GCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATg                                                                                                                 <  1:275195/251‑1 (MQ=255)                                                                                                                                     cGATTCACCGCACGCCAGAGAGATGGCGACCAAAGCGCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGC‑GCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGc                                                                                                                >  1:284931/1‑251 (MQ=255)                                                                                                                                          cACCGCACGCCAGAGGGAGGGCGACCAAAGCGCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCATCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCTTCCGGAGATGAGCAGCGACGAAATTGAGC‑GCATTATTCGGAGTGCGGTCAAGAACGTGCCTTATGCCGTGGGAATCACCAGCCACATGGGTAGCAAGATGGCCTCTAGCCTGTTTGGTATTCAGAAAGTGATACAGgc                                                                                                            <  2:413293/250‑1 (MQ=255)                                                                                                                                                    cAGAGAGATGGCGACCAAAGCGCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGC‑GCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGtt                                                                                                 <  2:240877/251‑1 (MQ=255)                                                                                                                                                    cAGAGAGATGGCGACCAAAGCGCATAACAGCGGGCACGAAGTGTTGAGTCATCTCCCGAGGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGC‑GCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGtt                                                                                                 <  2:375220/251‑1 (MQ=255)                                                                                                                                                    cAGAGAGATGGCGACCAAAGCGCATAACAGCGGGCACGAAGTGGTGATTCATCTCCCGATGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGC‑GCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGAGCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGtt                                                                                                 <  2:58687/251‑1 (MQ=255)                                                                                                                                                     agagagATGGCGACCAAAGCGCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGC‑GCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTa                                                                                                <  2:400136/251‑1 (MQ=255)                                                                                                                                                     agagagATGGCGACCAAAGCGCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGC‑GCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTa                                                                                                <  1:32234/251‑1 (MQ=255)                                                                                                                                                      gagagaTGGCGACCAAAGCGCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGC‑GCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGTGTTAc                                                                                               <  2:284931/251‑1 (MQ=255)                                                                                                                                                             ggCGACCAAAGCGCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGC‑GCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCttt                                                                                        <  1:14222/251‑1 (MQ=255)                                                                                                                                                               cGACCAAAGCGCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCACTGAGTAAACAGCAGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGC‑GCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTAc                                                                                      >  1:101628/1‑251 (MQ=255)                                                                                                                                                                  ccAAAGCGCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGC‑GCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTc                                                                                   >  2:16452/1‑251 (MQ=255)                                                                                                                                                                    aaaGCTCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGC‑GCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCt                                                                                 <  2:25908/251‑1 (MQ=255)                                                                                                                                                                    aaaGCGCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATGCGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGC‑GCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTGTACTTCCt                                                                                 <  2:144163/251‑1 (MQ=255)                                                                                                                                                                     aaGCGCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGC‑GCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTc                                                                                <  1:409862/251‑1 (MQ=255)                                                                                                                                                                         gcATAACCGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCACGGAGATGAGCAGCGACGAAATTGAGC‑GCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAACGTGATGCCGGCACTGGAGCGTTACAATCTTTACTTCCTCGAc                                                                             >  2:266222/1‑250 (MQ=255)                                                                                                                                                                               cAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGC‑GCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTaa                                                                      <  1:289611/251‑1 (MQ=255)                                                                                                                                                                                          aaGTGTTGCTTCATCTCCAGATGGCACCACTGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGC‑GCATTATACGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGCCCTCTACCCTGTTTTGTCTGCAGAAAGTGATGCAGGCACTGTAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAAt                                                           >  1:271234/1‑251 (MQ=255)                                                                                                                                                                                              gtTGATTCCTCTCCCGCCGGCAACACTGAGTAAAAAGCCGCTGGAGAAAAATACGCTGCGCCCGGCGATGAGCAGCGACGAAATTGAGC‑GCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGCAAGTGATGCAGGCACTGGAGCGTTACAATCTTTCCTTCCTCGACAGCGTAACCATCGGTAATAccc                                                       >  1:299311/1‑251 (MQ=255)                                                                                                                                                                                                             gATGGCACCACTGAGTAAACAGCCGCTGGAGAAAACTACGCTGCGCCCGGAGATGAGCAGCGACGCAATTGAGC‑GCATTATTCGTAGCGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGcgcgccg                                        >  1:20887/1‑251 (MQ=255)                                                                                                                                                                                                                   accactGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGC‑GCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTc                                                                                   >  2:152494/6‑202 (MQ=255)                                                                                                                                                                                                                   accactGAGTAAACAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGC‑GCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTc                                                                                   <  1:152494/197‑1 (MQ=255)                                                                                                                                                                                                                                cAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGC‑GCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGa                      >  1:397569/1‑250 (MQ=255)                                                                                                                                                                                                                                 aGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGC‑GCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGa                      >  2:171257/1‑249 (MQ=255)                                                                                                                                                                                                                                        ggAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGC‑GCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGAGCACCAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCGTTACTTCCTCGACAGCGGAACCATCGGTAATACCCAGGCGATGAGCGCCGCGCTAGGTACTGGCGTGaattggatca             <  2:385750/251‑9 (MQ=255)                                                                                                                                                                                                                                          agaAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGC‑GCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCaa            >  1:55652/1‑250 (MQ=255)                                                                                                                                                                                                                                           gaAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGC‑GCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAAc          <  2:335228/251‑1 (MQ=255)                                                                                                                                                                                                                                             aaaaTACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGC‑GCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACgg        <  2:63326/251‑1 (MQ=255)                                                                                                                                                                                                                                                aTACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGC‑GCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAg     <  1:274655/251‑1 (MQ=255)                                                                                                                                                                                                                                                aTACGCTGCGCACGGCGATGAGCAGCGACGAAATTGCGC‑GCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGACGGTGATCAAACGGa       >  2:275531/1‑249 (MQ=255)                                                                                                                                                                                                                                                   cgctgcgcCCGGAGATGAGCAGCGACGAAATTGAGC‑GCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGtg  >  2:82789/1‑250 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                        CGCTCTGTTGGCGCTCTCCTCCCCCGTACTTGCTGGCAAACTTGCCATCGTCATTGATGATTTTGGGTATCGCCCGCACAACGAAAACCAGGTGCTGGCGATGCCTTCCGCTATCTCCGTCGCTGTATTACCCGATTCACCGCACGCCAGAGAGATGGCGACCAAAGCGCATAACAGCGGGCACGAAGTGTTGATTCATCTCCCGATGGCACCGTTGAGTAAACAGCCGCTGGAGAAAAATACGCTACGCCCGGAGATGAGCAGCGACGAAATTGAGC‑GCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGCACTGGCGTGAAGGTGATCAAACGGAAGGTA  >  NC_000913/3788136‑3788628

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: ‑