Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,581,678	A→C	6.4%	intergenic (+39/+56)	tktB → / ← ypfG	transketolase 2, thiamine triphosphate‑binding/DUF1176 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,581,678	0	A	C	6.4%	253.1 / 4.3	94	intergenic (+39/+56)	tktB/ypfG	transketolase 2, thiamine triphosphate‑binding/DUF1176 family protein
Reads supporting (aligned to +/- strand):  ref base A (49/38);  new base C (3/4);  total (52/42)
Fisher's exact test for biased strand distribution p-value = 6.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGA  >  NC_000913/2581530‑2581827                                                                                                                                                     |                                                                                                                                                      cGGGTTACGGGGAATCTGCTCCGGCGGATATGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGAc                                                                                                                                                      >  1:1523753/1‑150 (MQ=255)    gTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGAcgcg                                                                                                                                                   >  1:1451658/1‑150 (MQ=255)       aCGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCt                                                                                                                                                <  2:1632719/150‑1 (MQ=255)       aCGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTCGGGCAAAAGCGCATAAGGTGCGGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCt                                                                                                                                                <  1:104241/150‑1 (MQ=255)         ggggAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTa                                                                                                                                              >  2:140140/1‑150 (MQ=255)          gggAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTAt                                                                                                                                             >  1:91423/1‑150 (MQ=255)           ggAAACTTCTCCGGCGGAAAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGGGCTGGGAGTGAAAGGGGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATc                                                                                                                                            <  1:655344/150‑1 (MQ=255)              aTCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGGGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCgg                                                                                                                                         <  2:1721671/150‑1 (MQ=255)               tCTGCTCCGGCTGATAAGCTCTTAACGTTCTTTGGCTTTACCTCCGTGAATATCGTGGCTAACGCGCATAAGGTGCTGGGAGGCAAAGGAGCCTGATGGTGATTGCCGTATGCTGATTGCCGGATAAGACGCGTCAGCGTCGCATCCGGc                                                                                                                                        <  2:288001/150‑1 (MQ=255)               tCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGc                                                                                                                                        <  1:382351/150‑1 (MQ=255)                cTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCTCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGTTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGcc                                                                                                                                       >  2:1466308/1‑150 (MQ=255)                cTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGAGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGcc                                                                                                                                       <  2:213486/150‑1 (MQ=255)                cTGCTCCGGCGGATAAGCTGTTCCCGTGCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCAT‑AGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCt                                                                                                                                      <  1:40372/150‑1 (MQ=255)                 tGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTCCCGCCGAGAATATTGGGGCAAACGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGACGGTGATTGCCGGATGCTGATTGCCGGATAAGACGCGTCAGCGTCGCATCCGGcaa                                                                                                                                      <  2:827838/150‑3 (MQ=255)                 tGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCt                                                                                                                                      <  2:1231608/150‑1 (MQ=255)                     ccGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAg                                                                                                                                  >  2:868855/1‑150 (MQ=255)                     ccGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAg                                                                                                                                  >  2:811067/1‑150 (MQ=255)                       ggcggATAAGCTGTTCCCGTTCTTTGGCGTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCg                                                                                                                                <  1:8843/150‑1 (MQ=255)                        gcggATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGt                                                                                                                               >  1:765894/1‑150 (MQ=255)                        gcggATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGt                                                                                                                               <  1:1004142/150‑1 (MQ=255)                        gcggATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGt                                                                                                                               <  2:1691391/150‑1 (MQ=255)                        gcggATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGt                                                                                                                               >  1:1076647/1‑150 (MQ=255)                             tAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTAt                                                                                                                          >  2:1709684/1‑150 (MQ=255)                               aGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATcc                                                                                                                        <  2:641075/150‑1 (MQ=255)                               aGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATcc                                                                                                                        >  2:26670/1‑150 (MQ=255)                                gCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCATATCCGGCCTACAGCGTCTTATCCg                                                                                                                       >  1:898814/1‑150 (MQ=255)                                 cTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCgg                                                                                                                      >  1:1427105/1‑150 (MQ=255)                                 cTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCgg                                                                                                                      >  2:888948/1‑150 (MQ=255)                                   gTTCCCGTTCTTAGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGCCGCGTCTTATCCCGCCTACAGCGTCTTATCCGGcc                                                                                                                    >  2:441500/1‑150 (MQ=255)                                            cTTTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTcc                                                                                                           >  1:145996/1‑150 (MQ=255)                                             tttGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATCGCCGGATAAGACGCGTCAGCGTCGCATCCGGcaatc                                                                                                                                    >  1:167460/1‑120 (MQ=37)                                             tttGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATCGCCGGATAAGACGCGTCAGCGTCGCATCCGGcaatc                                                                                                                                    <  2:167460/124‑5 (MQ=37)                                              ttGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCg                                                                                                         <  1:470171/150‑1 (MQ=255)                                              ttGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCg                                                                                                         >  1:105834/1‑150 (MQ=255)                                                ggCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGcc                                                                                                       >  2:1561008/1‑150 (MQ=255)                                                ggCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGcc                                                                                                       >  1:1444257/1‑150 (MQ=255)                                                  cTTTTCCGCCGAGAATATTGTGGCAAAAGCGGATAAGGGGCTGGGAGTGAAACGTGCCTGATGGTGATTGCCGGATGCGGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCAt                                                                                                     <  1:842055/150‑1 (MQ=255)                                                  cTTTACAGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGACTACATGTCCCGCCAt                                                                                                     >  1:1512834/1‑150 (MQ=255)                                                      aCCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTg                                                                                                 <  1:1570284/150‑1 (MQ=255)                                                         gccgAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGttt                                                                                              >  1:1010851/1‑150 (MQ=255)                                                         gccgAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGttt                                                                                              <  1:94550/150‑1 (MQ=255)                                                           cgAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTCCAGGGCCCGCCATTTTGTTTac                                                                                            >  2:942287/1‑149 (MQ=255)                                                            gagaTTATTGTGGCAACAGCGCATTAGGGGCTGGGAGTGAAACGGGCCTGTTGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAAc                                                                                           <  1:85149/150‑1 (MQ=255)                                                                atatTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGTAAGGTGCCTGATGGTGATTGCCGGAGGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGt                                                                                       <  2:733273/150‑1 (MQ=255)                                                                     gtgGCAAAAGCGCATAAGGTGCGGGGAGTGAAAGGGGCCTGATGGTGATTGCCGGAT‑CTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCca                                                                                 <  1:856117/150‑1 (MQ=255)                                                                       ggCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCcac                                                                                <  1:685923/150‑1 (MQ=255)                                                                       ggCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCcac                                                                                >  2:848731/1‑150 (MQ=255)                                                                        gCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCcaca                                                                               <  2:1623747/150‑1 (MQ=255)                                                                         cAAAAGCGCATAAGGGGCGGGGAGTGAAAGGGGGCGGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGGCTTATGCGGCCTACCGCGTCTTTTACCGGCGACATGTGCCGCCATTTTGTTTTACGGGTGATacacca                                                                              >  2:1186/1‑148 (MQ=255)                                                                          aaaaGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAAc                                                                             <  1:707436/150‑1 (MQ=255)                                                                           aaaGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACg                                                                            <  2:988863/150‑1 (MQ=255)                                                                           aaaGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACg                                                                            >  1:1330473/1‑150 (MQ=255)                                                                              gcgcATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTgg                                                                         <  1:694013/150‑1 (MQ=255)                                                                                gcATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGc                                                                       >  2:177769/1‑150 (MQ=255)                                                                                 cATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGcc                                                                      <  2:1715418/150‑1 (MQ=255)                                                                                     aGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGc                                                                  >  1:9231/1‑150 (MQ=255)                                                                                      ggTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCg                                                                 >  2:1411741/1‑150 (MQ=255)                                                                                         gCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCt                                                              <  2:1427105/150‑1 (MQ=255)                                                                                           tGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTgg                                                            >  1:212182/1‑150 (MQ=255)                                                                                            gggAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGc                                                           >  2:1098667/1‑150 (MQ=255)                                                                                             ggAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGcc                                                          >  1:108015/1‑150 (MQ=255)                                                                                             ggAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGcc                                                          >  2:446533/1‑150 (MQ=255)                                                                                             ggAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGcc                                                          >  1:913145/1‑150 (MQ=255)                                                                                             ggAGTGAAAGGTGCCTGATGGGGATTGCCGTATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCGACAGCGTATTATCCGGCCTAGATGTCCCGCCATTTTGTTTAACGGGTTATCCACAACGTGGGCCAGGCGTCTGGcc                                                          <  1:738651/150‑1 (MQ=255)                                                                                              gAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACTGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGccc                                                         >  2:1027246/1‑150 (MQ=255)                                                                                              gAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGccc                                                         >  1:1010503/1‑150 (MQ=255)                                                                                               aGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCa                                                        >  2:686815/1‑150 (MQ=255)                                                                                                   aaaGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGACGCCTGGCCCacgtt                                                    <  2:1610209/150‑5 (MQ=255)                                                                                                     tggTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAg                                                  <  2:377911/149‑1 (MQ=255)                                                                                                        tGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTa                                               >  2:491605/1‑150 (MQ=255)                                                                                                        tGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTa                                               >  2:461416/1‑150 (MQ=255)                                                                                                         gCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTAt                                              >  1:1255938/1‑150 (MQ=255)                                                                                                            tGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACACCGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGc                                           >  1:219740/1‑150 (MQ=255)                                                                                                                     ttgccggatgctgattgccggatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGtt                                  >  1:592360/1‑150 (MQ=255)                                                                                                                       gccggatgctgattgccggatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCt                                <  2:1266519/150‑1 (MQ=255)                                                                                                                       gccggatgctgattgccggatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCt                                >  2:195069/1‑150 (MQ=255)                                                                                                                        ccggatgctgattgccggatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTg                               >  1:651534/1‑150 (MQ=255)                                                                                                                           gatgctgattgccggatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTgcgg                            <  1:177769/150‑1 (MQ=255)                                                                                                                              gctgattgccggatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGt                         <  1:978124/150‑1 (MQ=255)                                                                                                                               ctgattgccggatgcGACGCTGCCGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTa                        <  2:1403975/150‑1 (MQ=255)                                                                                                                                 gattgccggatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAAc                      <  1:140140/150‑1 (MQ=255)                                                                                                                                 gattgccggatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAAc                      >  2:1322683/1‑150 (MQ=255)                                                                                                                                     gccggatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCAc                  <  2:1010851/150‑1 (MQ=255)                                                                                                                                      ccggatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCAcc                 >  2:848462/1‑150 (MQ=255)                                                                                                                                      ccggatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCAcc                 >  2:1398265/1‑150 (MQ=255)                                                                                                                                      ccggatgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCAcc                 >  2:1405192/1‑150 (MQ=255)                                                                                                                                          atgcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGc             <  2:662279/150‑1 (MQ=255)                                                                                                                                           tgcGACGCTGACGCGTCTTATCCGGCCTACAGCGGCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCg            <  2:530941/150‑1 (MQ=255)                                                                                                                                            gcGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTAACCGGCATACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGa           >  1:531615/1‑150 (MQ=255)                                                                                                                                             cGACGCTGCCGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCAAAACATGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCTGAGTCACGCACCAGGCGaa          >  2:324774/1‑150 (MQ=255)                                                                                                                                             cGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGaa          >  2:645538/1‑150 (MQ=255)                                                                                                                                             cGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGaa          >  1:1294637/1‑150 (MQ=255)                                                                                                                                                cGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAAc       <  1:686815/150‑1 (MQ=255)                                                                                                                                                 gCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACg      <  2:212182/150‑1 (MQ=255)                                                                                                                                                     aCGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGa  <  2:983881/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                      CGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACATGTCCCGCCATTTTGTTTAACGGGTGATCCACAACGTGGGCCAGGCGTCTGGCCCATGCCAGTTATCGCAGGTGGGTTCTGCGGCGTAACGCACCAGGCGAAAACGCTGA  >  NC_000913/2581530‑2581827

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