Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,163,371	0	C	T	46.9%	73.3 / 35.3	64	L38F (CTC→TTC)	thiK	thiamine kinase
Reads supporting (aligned to +/- strand):  ref base C (0/34);  new base T (30/0);  total (30/34)
Fisher's exact test for biased strand distribution p-value = 6.17e-19
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGGTGCAGACGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGAGTTTTCTCATTGAACATCAGGGCCAG‑CG‑TTTTGTTGTGCGTCAGCCGCACGATCCTGATGCGCCGCAGTCCGCGTTCTTGCGCCAGTATCGGGCTTTATCACAACTACCCGCATGCATTGCACCGAAGCCGCATTTATATCTCCGTGACTGGA  >  NC_000913/1163222‑1163518                                                                                                                                                       |                                                                                                                                                      gCTGGTGCAGACGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGAGTTTtct                                                                                                                                                      <  2:32613/151‑1 (MQ=255)   tGGTGCAGACGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGGGGTTtttct                                                                                                                                                    >  1:288061/1‑150 (MQ=255)       gCAGACGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGAGTTTTCTCATTg                                                                                                                                                 <  1:840004/150‑1 (MQ=255)        cAGACGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGAGTTTTCTCATTGaa                                                                                                                                               <  1:1429016/151‑1 (MQ=255)         aGACGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGGGGTTTTTCttttta                                                                                                                                               >  2:864164/1‑147 (MQ=255)         aGACGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGAGTTTTCTCATTGAAc                                                                                                                                              <  1:857081/151‑1 (MQ=255)         aGACGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGAGTTTTCTCATTGAAc                                                                                                                                              <  2:1078064/151‑1 (MQ=255)          gACGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGAGTTTTCTCATTGAACa                                                                                                                                             <  1:384006/151‑1 (MQ=255)             gggCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGGGGTTTTTCCTTTAAAATc                                                                                                                                           >  2:135808/1‑150 (MQ=255)             gggCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGAGTTTTCTCATTGAACATCa                                                                                                                                          <  2:292986/151‑1 (MQ=255)                       aTCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGGGGTTTTTCGTTTATCAACACGGCCAgagg                                                                                                                                >  2:161713/1‑148 (MQ=255)                         cTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGGGGTTTTTCTTTTAAAATCAGGGCCGG‑CG‑ttt                                                                                                                            >  2:1280082/1‑151 (MQ=255)                              cAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGGGGTTTTTTTTTTTACAACACGGCCCG‑CG‑TTttgttt                                                                                                                       >  1:651037/1‑150 (MQ=255)                               aGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGAGTTTTCTCATTGAACATCAGGGCCAG‑CG‑TTTtgttgt                                                                                                                      <  2:74339/151‑1 (MQ=255)                               aGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGAGTTTTCTCATTGAACATCAGGGCCAG‑CG‑TTTtgttgt                                                                                                                      <  1:1176355/151‑1 (MQ=255)                                ggTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGCGGGGGGTTTTTTCTTTTAA‑ATCACGGCCAG‑CG‑GTTTTTtgtg                                                                                                                     >  1:194241/1‑151 (MQ=255)                                      ggTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGGGGTTTTTCATTTAACATCAAGGCCCG‑CG‑GTTTTTTTTGGGTc                                                                                                                 >  2:737633/1‑149 (MQ=255)                                      ggTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGAGTTTTCTCATTGAACATCAGGGCCAG‑CG‑TTTTGTTGTGCGTCAg                                                                                                               <  1:1323885/151‑1 (MQ=255)                                           cGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGAGTTTTCTCATTGAACATCAGGGCCAG‑CG‑TTTTGTTGTGCGTCAGCCGCa                                                                                                          <  1:705734/151‑1 (MQ=255)                                            gTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGAGTTTTCTCATTGAACATCAGGGCCAG‑CG‑TTTTGTTGTGCGTCAGCCGCAc                                                                                                         <  2:351839/151‑1 (MQ=255)                                            gTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGAGTTTTCTCATTGAACATCAGGGCCAG‑CG‑TTTTGTTGTGCGTCAGCCGCAc                                                                                                         <  1:711191/151‑1 (MQ=255)                                             tttCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGAGGTTTTTTCTTTTACATCAAGGGCAG‑CG‑GTTTTTTTTGCGTCAGCCGCCCg                                                                                                        >  2:182097/1‑151 (MQ=255)                                              ttCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGAGTTTTCTCATTGAACATCAGGGCCAG‑CG‑TTTTGTTGTGCGTCAGCCGCACg                                                                                                        <  2:1256894/150‑1 (MQ=255)                                                cGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGGGGTTTTTCCTTTAACATCAGGGCCCG‑CG‑GTTTTTTGTGCGTCAGCCGCCCGaac                                                                                                     >  2:399362/1‑149 (MQ=255)                                                      aataatCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGAGGTTTTTTTTTTAACATCAAGGCCCG‑CG‑TTTTTTTTTTGGTGCGCCGCACGCACCCGATg                                                                                               >  2:971552/1‑151 (MQ=255)                                                       ataatCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGAGTTTTCTCATTGAACATCAGGGCCAG‑CG‑TTTTGTTGTGCGTCAGCCGCACGATCCTGATgc                                                                                              <  1:1429055/151‑1 (MQ=255)                                                         aatCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGAGTTTTCTCATTGAACATCAGGGCCAG‑CG‑TTTTGTTGTGCGTCAGCCGCACGATCCTGATgcg                                                                                             <  2:425490/150‑1 (MQ=255)                                                                tCACGCGCGACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGGGTTTTTTCATTTAACATCAAGGCCAG‑CG‑TTTTGTTTTGTGTCAGCCGCCCGCACCTGATGCGCCGCAggc                                                                                     >  2:820790/1‑149 (MQ=255)                                                                 cACGCGCGACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGGGTTTTTTTTTTTTAAATAAAGGCCAG‑CG‑GTTTTTTTTGGGTCCGCCGCACGATCCCTATGCTGCGCCGTcc                                                                                    >  1:898423/1‑151 (MQ=255)                                                                 cACGCGCGACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGGGGTTTTTCTTTTTAAATCAAGGCCCG‑CG‑GTTTGTTTTGGGTCAGCCGCACGATCCTCCTGCGCCGCCGTcc                                                                                    >  2:1353007/1‑151 (MQ=255)                                                                 cACGCGCGACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGAGTTTTCTCATTGAACATCAGGGCCAG‑CG‑TTTTGTTGTGCGTCAGCCGCACGATCCTGATGCGCCGCAGTc                                                                                     <  1:793054/150‑1 (MQ=255)                                                                    gcgcgACGAATTGCTGTCGCGC‑TTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGAGGTTTTTCCTTTTACATCAAGGCCAG‑CG‑GTTTTTTGTGGGGCAGCCGCCCCATCCTGATGCGCCGCCGTCcgcgc                                                                                >  2:66825/1‑150 (MQ=255)                                                                     cgcgACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGGGGTTTTTTTTTTTAAAACAAGGCCAC‑CG‑GTTTTTTTTTGGTGCGCCGCACCCTCCCGGTGcgccgccggcccccg                                                                                >  1:1067858/1‑142 (MQ=255)                                                                     cgcgACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGAGTTTTCTCATTGAACATCAGGGCCAG‑CG‑TTTTGTTGTGCGTCAGCCGCACGATCCTGATGCGCCGCAGTCCGCGt                                                                                <  1:1091283/151‑1 (MQ=255)                                                                     cgcgACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGAGTTTTCTCATTGAACATCAGGGCCAG‑CG‑TTTTGTTGTGCGTCAGCCGCACGATCCTGATGCGCCGCAGTCCGCGt                                                                                <  1:536355/151‑1 (MQ=255)                                                                         acgaATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGGGGTTTTTTTTTTTAAATCAAGGCCAG‑CG‑CTTTTTTTTGTGGCAGCCCCCCGAGCCCCATTCGCCCCCGTCCGCGGTCtt                                                                            >  2:36833/1‑151 (MQ=255)                                                                         acgaATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGAGTTTTCTCATTGAACATCAGGGCCAG‑CG‑TTTTGTTGTGCGTCAGCCGCACGATCCTGATGCGCCGCAGTCCGCGTTCtt                                                                            <  1:1311329/151‑1 (MQ=255)                                                                           gaATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGGGGTTTTTTTTTTTAAAACAAGGCCCG‑CG‑TTTTTTTTTGTGTCAGCCGCCCCAACCTGATGCGCCGCCGTCCCCGGTCtttc                                                                          >  2:966346/1‑149 (MQ=255)                                                                           gaATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGAGGTTTTTTTTTGTACATCAAGGCCAG‑CG‑GTTTTTTTTGCGTCAGCCGCACGATCCCGATGCGCCGCAGTCCGCGTTCTTgc                                                                          >  2:1457890/1‑151 (MQ=255)                                                                            aaTTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGAGTTTTCTCATTGAACATCAGGGCCAG‑CG‑TTTTGTTGTGCGTCAGCCGCACGATCCTGATGCGCCGCAGTCCGCGTTCTTgcg                                                                         <  2:700371/151‑1 (MQ=255)                                                                             aTTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGAGGTTTTTTTTTTTACATCAAGGCCAGACG‑TTTT‑TTTTGGGTTTGGCGCCCGCACCCGATGCGGCGCCGCCCGCGGTcctttgc                                                                        >  1:889561/1‑145 (MQ=255)                                                                               tGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGGGGTTTTTTCTTTAAAAACAAGGCCAG‑CG‑GTTTGTTGTGGGTCAGGCCCCCGATCCTGATGCGCCGCCGGCCCCGTTCTTTCGccc                                                                      >  2:322724/1‑150 (MQ=255)                                                                                   gTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGAGTTTTCTCATTGAACATCAGGGCCAG‑CG‑TTTTGTTGTGCGTCAGCCGCACGATCCTGATGCGCCGCAGTCCGCGTTCTTGCGCCAGTAt                                                                  <  1:342275/151‑1 (MQ=255)                                                                                        gcTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGAGTTTTCTCATTGAACATCAGGGCCAG‑CG‑TTTTGTTGTGCGTCAGCCGCACGATCCTGATGCGCCGCAGTCCGCGTTCTTGCGCCAGTATCGGGc                                                             <  2:870913/151‑1 (MQ=255)                                                                                        gcTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGAGGTTTTTTATTTGACAACAGGGGCAG‑CG‑TTTTTTTTTGGGGCAGCCGCACGATCCCGGTGCGCCGCCGTCCCCGTTTTTTCGCCCGTATCgggg                                                             >  2:239951/1‑150 (MQ=255)                                                                                                 cccattATCACCCCGTCACGACGTTTAATAGTGGGTTTAGTGCC‑GGGAGTTTTCTCATTGACCATCA‑GGCCAG‑CG‑TTTTGTTGTGCGTCAGCCGCACGATCCTGATGCGCCGCAGTCCGCGTTCTTGCGCCAGTATCGGGCTTTATCACaa                                                   <  1:1279869/146‑1 (MQ=255)                                                                                                 cGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGGGGTTTTTTATTTTAAATCACGGGCCG‑CG‑GTTTTTTTTGGGTCGGCCGCCCGACCCTCATTCGCCCCCGGCCCCGTTCTTTCGCCCGCATCggggctttttac                                                     >  2:896804/1‑141 (MQ=255)                                                                                                  gCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGAGTTTTCTCATTGAACATCAGGGCCAG‑CG‑TTTTGTTGTGCGTCAGCCGCACGATCCTGATGCGCCGCAGTCCGCGTTCTTGCGCCAGTATCGGGCTTTATCACaa                                                   <  1:300318/151‑1 (MQ=255)                                                                                                   cAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGAGGTTTTTTTTTTAACATCAAGGCCAG‑CGCTTTTTTTTTGGGT‑TGGCGCCCGCGACCGATGCGCCGCAGGCCCCGTGCTTTCGCCCCGATCGGGCcttttttcaa                                                   >  1:877022/1‑143 (MQ=255)                                                                                                    aGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGGGGTTTTTTCTTTAAAATCAAGGCCCG‑CG‑GTTTTTTTTGGGTCAGCCGCCCGATCCCTATGCGCCGCCGTCCGCGTTCTTTCGCCCGTATCGGGGTTTTTCACaaa                                                  >  2:684182/1‑149 (MQ=255)                                                                                                        tcatcCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGAGGTTTTTTATTTAACATCAAGGCCCG‑CG‑GTTTTTTGTGGGGCAGCCGCCCGATCCCGATGCGCCGCCGTCCCCGGTCTTTCGCCCGTATCGGGCTTTTTTACAACTAccc                                             >  2:1596071/1‑151 (MQ=255)                                                                                                            cccGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGAGTTTTTTCATTTTAAAACAAGGCCAG‑CG‑CTTTGTTTTGGGGCAGCCGCCCGATCCTCATGCGCCGCCGGCCGCGTTCTTGCGCCCGTATAGGGCCTTTTTCCAACTAcccccct                                         >  1:1454325/1‑147 (MQ=255)                                                                                                                tCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGAGGTTTTTTATTTTAAAACAAGGCCAG‑CG‑GTTTTTTTTTGGGCAGCCCCCCCATCCTGATGCGCCGCAGGCCCCGGTCTTTCGCCCGTATCGGGGTTTTTTACAACTAcccccccgcat                                     >  2:824801/1‑143 (MQ=255)                                                                                                                 cACGACGTTTAATAGTGGGCTTAGTGGC‑GGGAGTTTTCTCATTGAACATCAGGGCCAG‑CG‑TTTTGTTGTGCGTCAGCCGCACGATCCTGATGCGCCGCAGTCCGCGTTCTTGCGCCAGTATCGGGCTTTATCACAACTACCCGCATGCAtt                                    <  1:1439921/151‑1 (MQ=255)                                                                                                                   cgacgTTTAATAGTGGGCTTAGTGGCGGGGGGTTTTTTCTTTTAA‑ATCAAGGGCCG‑GG‑TTTTTTTTTGGGGGGGGCGCCCGATCCCGATGCGGCGCCGGCCGCGGTCTTTCGCCAGCATCGGGGGTTTTTACAACCAcccccccgcctttt                                  >  1:388178/1‑140 (MQ=255)                                                                                                                     acgTTTAATAGTGGGCTTAGTGGC‑GGGAGGTTTTTTTTTTTAAATCAGAGCCAG‑CG‑TTTTTTTTGGGGGGCGCCGCACGATCCTGTTGCGCCGCCGGCCGCGTTCTTTCGCCCGCATCGGGCTTTTTTACAACTACCCCCCTGCATTGCAc                                >  1:773866/1‑151 (MQ=255)                                                                                                                            aTAGTGGGCTTAGTGGC‑GGGGGTTTTTTTTTTTAAAAACAGGGGCCC‑CG‑TTTTTTTGTTCGGCTGCCGCGCGCTCCTTGTGGGCCGCAGTCCGCGGTTTTTCGCGCGTATCGGGGGTTTTCATAACTAACCCCATGCCCTGCACCgcaacc                         >  2:737845/1‑146 (MQ=255)                                                                                                                              aGTGGGCTTAGTGGC‑GGGAGTTTTCTCATTGAACATCAGGGCCAG‑CG‑TTTTGTTGTGCGTCAGCCGCACGATCCTGATGCGCCGCAGTCCGCGTTCTTGCGCCAGTATCGGGCTTTATCACAACTACCCGCATGCATTGCACCGAAGCCg                        <  1:18587/150‑1 (MQ=255)                                                                                                                                      tAGTGGC‑GGGAGTTTTCTCATTGAACATCAGGGCCAG‑CG‑TTTTGTTGTGCGTCAGCCGCACGATCCTGATGCGCCGCAGTCCGCGTTCTTGCGCCAGTATCGGGCTTTATCACAACTACCCGCATGCATTGCACCGAAGCCGCATTtata                <  1:300602/150‑1 (MQ=255)                                                                                                                                      tAGTGGC‑GGGAGTTTTCTCATTGAACATCAGGGCCAG‑CG‑TTTTGTTGTGCGTCAGCCGCACGATCCTGATGCGCCGCAGTCCGCGTTCTTGCGCCAGTATCGGGCTTTATCACAACTACCCGCATGCATTGCACCGAAGCCGCATTtat                 <  1:69542/149‑1 (MQ=255)                                                                                                                                       aGTGGC‑GGGAGTTTTCTCATTGAACATCAGGGCCAG‑CG‑TTTTGTTGTGCGTCAGCCGCACGATCCTGATGCGCCGCAGTCCGCGTTCTTGCGCCAGTATCGGGCTTTATCACAACTACCCGCATGCATTGCACCGAAGCCGCATTTatat               <  1:1206066/150‑1 (MQ=255)                                                                                                                                        gTGGC‑GGGAGTTTTCTCATTGAACATCAGGGCCAG‑CG‑TTTTGTTGTGCGTCAGCCGCACGATCCTGATGCGCCGCAGTCCGCGTTCTTGCGCCAGTATCGGGCTTTATCACAACTACCCGCATGCATTGCACCGAAGCCGCATTTATAtct             <  1:347366/151‑1 (MQ=255)                                                                                                                                               ggAGTTTTCTCATTGAACATCAGGGCCAG‑CG‑TTTTGTTGTGCGTCAGCCGCACGATCCTGATGCGCCGCAGTCCGCGTTCTTGCGCCAGTATCGGGCTTTATCACAACTACCCGCATGCATTGCACCGAAGCCGCATTTATATCTCCGTGa       <  2:760043/151‑1 (MQ=255)                                                                                                                                                  gTTTTCTCATTGAACATCAGGGCCAG‑CG‑TTTTGTTGTGCGTCAGCCGCACGATCCTGATGCGCCGCAGTCCGCGTTCTTGCGCCAGTATCGGGCTTTATCACAACTACCCGCATGCATTGCACCGAAGCCGCATTTATATCTCCGTGACTg    <  2:1463197/151‑1 (MQ=255)                                                                                                                                                    tttCTCATTGAACATCAGGGCCAG‑CG‑TTTTGTTGTGCGTCAGCCGCACGATCCTGATGCGCCGCAGTCCGCGTTCTTGCGCCAGTATCGGGCTTTATCACAACTACCCGCATGCATTGCACCGAAGCCGCATTTATATCTCCGTGACTgg   <  1:805451/150‑1 (MQ=255)                                                                                                                                                     ttCTCATTGAACATCAGGGCCAG‑CG‑TTTTGTTGTGCGTCAGCCGCACGATCCTGATGCGCCGCAGTCCGCGTTCTTGCGCCAGTATCGGGCTTTATCACAACTACCCGCATGCATTGCACCGAAGCCGCATTTATATCTCCGTGACTGGa  <  1:1413204/150‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                      GCTGGTGCAGACGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GGGAGTTTTCTCATTGAACATCAGGGCCAG‑CG‑TTTTGTTGTGCGTCAGCCGCACGATCCTGATGCGCCGCAGTCCGCGTTCTTGCGCCAGTATCGGGCTTTATCACAACTACCCGCATGCATTGCACCGAAGCCGCATTTATATCTCCGTGACTGGA  >  NC_000913/1163222‑1163518

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