Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,185,199	0	T	G	43.8%	116.9 / 31.0	90	T35P (ACC→CCC)	thiM	hydoxyethylthiazole kinase
Reads supporting (aligned to +/- strand):  ref base T (10/40);  new base G (39/0);  total (50/40)
Fisher's exact test for biased strand distribution p-value = 3.73e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑GGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGT  >  NC_000913/2185052‑2185348                                                                                                                                                      |                                                                                                                                                      caggCAGCGCCTGGGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGCCACCGGCGATAGCCGCAAACTGACTTGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GTTAAA‑GGt                                                                                                                                                       <  2:503151/145‑1 (MQ=255)  aCGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑GGTTTGc                                                                                                                                                   <  2:525359/151‑1 (MQ=255)   cGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑GGgtt                                                                                                                                                     >  2:249648/1‑148 (MQ=255)    gCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑GGGTTGcac                                                                                                                                                 >  2:496521/1‑151 (MQ=255)    gCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGACAGCAAGGTATTGGCGGGGAAA‑GGtt                                                                                                                                                      >  1:679266/1‑147 (MQ=255)     cATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGGAAA‑GGGTTGcccc                                                                                                                                                >  2:595673/1‑151 (MQ=255)        cGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑GGTTTGcacc                                                                                                                                                >  1:572988/1‑148 (MQ=255)         gCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑GGTTTGCACCAcat                                                                                                                                            <  1:44952/151‑1 (MQ=255)          ccTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑GGTTTGCACCAcatc                                                                                                                                           <  1:1130296/151‑1 (MQ=255)          ccTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑AGGTTGCACCAcatc                                                                                                                                           >  1:477157/1‑151 (MQ=255)             gAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGAAAA‑GGTTTGGACCACATCAtt                                                                                                                                        >  1:299150/1‑151 (MQ=255)               gcgcGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑GGTTTGCACCACATCATTgg                                                                                                                                      <  1:784429/151‑1 (MQ=255)               gcgcGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGCAAA‑GGGTTTGACCACATCCtt                                                                                                                                        >  1:968574/1‑149 (MQ=255)                   gTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTTAA‑AGATTGcccccccacccttgt                                                                                                                                     >  2:586277/1‑137 (MQ=255)                   gTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGGAAA‑GGGTTGAACCCCATCATTgg                                                                                                                                      >  1:764270/1‑147 (MQ=255)                    tGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGTCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGCGCGCCAGCAAGGTATTGGC‑GGTAAA‑GGTTTGCACCACATCATTGGTCa                                                                                                                                   <  2:1196770/149‑1 (MQ=255)                     ggCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑AGGTTGCAACAcgtctgttg                                                                                                                                      >  2:1080825/1‑140 (MQ=255)                       cTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑AGGTTGCACCCCACCAtttg                                                                                                                                      >  2:1092876/1‑141 (MQ=255)                        tGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGGAAAGGGTTGGCACACCATCTTTggcat                                                                                                                                   >  1:840155/1‑143 (MQ=255)                         gCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑GGTTTGCACCACATCATTGGTCATGCAGTg                                                                                                                            <  2:298363/151‑1 (MQ=255)                            tCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTGAA‑AGGATGCACCCcatcactttttattta                                                                                                                               >  2:350700/1‑134 (MQ=255)                             cAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑GGTTTGCACCACATCATTGGTCATGCAGTGcaca                                                                                                                        <  1:641626/151‑1 (MQ=255)                             cAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑GGTTTGCACCACATCATTGGTCATGCAGTGcac                                                                                                                         <  1:879805/150‑1 (MQ=255)                             cAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑GGGTTTCACCCCATCCTTTGTCAAGCAggggac                                                                                                                         >  1:1117188/1‑145 (MQ=255)                              aGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTGAA‑GGGTTGCACCACATCATTGGGCATGccagg                                                                                                                            >  2:683915/1‑142 (MQ=255)                              aGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑GGTTTGCACCACATCATTGGTCATGCAGTGCACaa                                                                                                                       <  2:282685/151‑1 (MQ=255)                                 gtgCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑GGTTTGCACCACATCATTGGTCATGCAGTGCACAagag                                                                                                                    <  2:730968/151‑1 (MQ=255)                                 gtgCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑GGGTTGCACCACATCATTGGTCATGCAGTGCACAag                                                                                                                      >  1:359775/1‑149 (MQ=255)                                  tgCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑GGGTTGGACCACATCATTGGTCATGCAGGGCACAAGAgg                                                                                                                   >  1:882791/1‑151 (MQ=255)                                   gCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑GGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAggg                                                                                                                  <  2:1172353/151‑1 (MQ=255)                                   gCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑GGGTTTCACCACATCATTTGTCATGCAGGGCACAAGAggg                                                                                                                  >  2:248511/1‑151 (MQ=255)                                    ccAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑GGTTTTGTCCACATCATTGGTCCTGCCGTGCACAAGAgggg                                                                                                                 >  2:415835/1‑151 (MQ=255)                                    ccAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑GGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAgggg                                                                                                                 <  1:594964/151‑1 (MQ=255)                                        cGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑GGGTTGCACCACATCCTTGGTCATGCAGTGCACAAGAgggg                                                                                                                 >  2:589179/1‑147 (MQ=255)                                        cGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑GGGTTGCAACACCTCATTGGTCATGCAGTGCACAAGAggggg                                                                                                                >  2:486209/1‑147 (MQ=255)                                         gTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑GGGTTGGAACACATCATTGGGCATGCAGAGCACAAGAgggg                                                                                                                 >  1:340014/1‑146 (MQ=255)                                         gTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑GGGTTGCACCACATTCTTGGTCATGCAGGGCAAAAGAggggg                                                                                                                >  1:645993/1‑146 (MQ=255)                                            aaTCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑GGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTg                                                                                                          <  1:1127463/150‑1 (MQ=255)                                             aTCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑GGGTTGCACCCCATCATTTGTCATGCAGTGCACAAGAGGGGGATGtt                                                                                                           >  1:174189/1‑148 (MQ=255)                                                aacaaGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑GGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGaa                                                                                                     <  1:577795/151‑1 (MQ=255)                                                  caaGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑GGGTTTCACCACATCATTTGTCATGCAGTGGAAAAGAGGGGGATGTTgggggaa                                                                                                    >  2:970135/1‑145 (MQ=255)                                                  caaGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGGATTGGC‑GGGAAA‑GGTTTGGGCCACAACATTGGGCATGCCAGGCAAAAAAGGGGAGAGTTGGTg                                                                                                       >  2:692279/1‑147 (MQ=255)                                                   aaGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑GGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGaaaaa                                                                                                  <  1:707657/151‑1 (MQ=255)                                                      gCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑GGGTTGCACCACATCATTTGTCATGCCGTGCACAAGAGGGGGAAGTTGGTGaaaaa                                                                                                  >  2:660315/1‑148 (MQ=255)                                                      gCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑GGGTTGCACCACATCATTTGTCATGCAGGGCACAAGAGGGGGATGTTGGTGGAAAAg                                                                                                 >  1:1125548/1‑149 (MQ=255)                                                      gCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑GGGTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGGATGTTGGTGAAAAAgg                                                                                                >  1:1200634/1‑150 (MQ=255)                                                       cACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑GGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGgtg                                                                                              <  2:981417/151‑1 (MQ=255)                                                        aCTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑GGGTTTCACCACATCATTGGTTATGCAGGGCCCAAGAGGGGGAAGGTGGGGgaaaa                                                                                                  >  2:1071809/1‑146 (MQ=255)                                                        aCTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑AGGTTGCACCACATCCTTTGGCATGCAGTGGACAAGAGGGGGATGTTGGTGGAAAAgggg                                                                                              >  2:26589/1‑148 (MQ=255)                                                           ggCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑GGGTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGGATGTTGGGGAAAAAGGGGTaa                                                                                           >  1:497167/1‑150 (MQ=255)                                                              gATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑GGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAAcgcg                                                                                       <  2:174340/151‑1 (MQ=255)                                                              gATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑GGGTTGCACCACATCATTGGTCATGCCGTGCACAAGAGGGGGATGTTGGGGAAAAAgggggaa                                                                                           >  1:741649/1‑142 (MQ=255)                                                               aTAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑AGGTTGCACCCCCTCATTTGTCATGCAGGGCACAAGAGGGGGATGTTGGGGAAAAAAGGGTa                                                                                            >  2:837789/1‑145 (MQ=255)                                                                 aGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑GGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTgc                                                                                    <  1:239226/151‑1 (MQ=255)                                                                     gCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGGATTGGC‑GGCAAA‑AGGTTGCTCCAACAAATTGGCTTTGTGGGGCTCAAGTGGGGGATTTTGGGGGAAAAGGCGTAACTCGTgcgccg                                                                                >  1:772384/1‑149 (MQ=255)                                                                       aaaCTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑GGTGTGCACCACATCCTTGGTCATGCCGGGCACAAGAGGGGGATGGTGGTGAAAAAGGGGGAACGCGTgcg                                                                                   >  2:30346/1‑146 (MQ=255)                                                                            gactgGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑GGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTgcgc                                                                         <  1:360653/151‑1 (MQ=255)                                                                            gactgGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑GGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTgcg                                                                          <  1:953472/150‑1 (MQ=255)                                                                            gactgGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑GGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTgc                                                                           <  1:858617/149‑1 (MQ=255)                                                                            gactgGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑GGGTTGCACCACAACATTGGTCATGCAGTGGACCAGAGGGGGATGTTGGTGGAAAAGGTGTAACGCGTGGGCAGGtt                                                                             >  1:413298/1‑147 (MQ=255)                                                                                ggCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑GGGTTGCACCACATCATTTGTCATGCAGTGCACAAGAGGGGGATGTTGGTGAAAAAGGGGTAACGCGTGCGCAGATTGCGCTGa                                                                      >  1:593718/1‑150 (MQ=255)                                                                                 gCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑GGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAAc                                                                    <  2:845117/151‑1 (MQ=255)                                                                                       tCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑GGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCca                                                                 <  2:142637/148‑1 (MQ=255)                                                                                        cGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑GGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCcagca                                                              <  1:582761/150‑1 (MQ=255)                                                                                        cGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑GGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCcagca                                                              <  1:779804/150‑1 (MQ=255)                                                                                           tttCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑GGGTTGCACCACATCATTTGTCATGCAGTGGACAAAAGGGGGATGTTGGTGGAAAAAGGGTAAAGCGGGCGCAGGTTGCGCTGAACCCAGCAgg                                                            >  2:415195/1‑149 (MQ=255)                                                                                           tttCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑GGGTTGCACCACATCATTTGTCATGCAGGGCACAAGAGGGGAATGTTGGGGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAgg                                                            >  2:814447/1‑149 (MQ=255)                                                                                            ttCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑GGTTTGTACCCCATTATTGGTCATGCAGTGCACAAGAGGGGGAAGGTGGGGAAAAAGGTGGAACGCGGGCGCAGAATGCGCTTAACCCCGAAGGGCg                                                         >  2:1003127/1‑151 (MQ=255)                                                                                                aTAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑GGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACtt                                                     <  2:299138/151‑1 (MQ=255)                                                                                                  aaCCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑GGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTg                                                    <  2:670057/150‑1 (MQ=255)                                                                                                     cATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑AGGTTGCACCACAACATTGGGCATGCAGTGGACAAGAGGGGGATGTTGGGGAAAAAAGTGTAAAGCGTGCGCAGAGTGCGCTGAAACCAGCAGGGCGACTTGCAt                                                 >  2:565547/1‑150 (MQ=255)                                                                                                      aTCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑GGGTTGCACCACATCATTTGTCATGCAGTGCACAAGAGGGGGATGTTGGTGAAAAAGGGGGAAAGCGCGCGCAGATTGCGCTGAAACCAGCAGGGCGACTTGCAt                                                 >  1:395504/1‑149 (MQ=255)                                                                                                        cGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑GGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGt                                              <  2:205995/150‑1 (MQ=255)                                                                                                         gCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑GGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGttt                                            <  1:574323/151‑1 (MQ=255)                                                                                                           tGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑GGGTTTCACCCCATCATTTGGCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGGGGAAAGCGGGCGCAGATTGCGCTGAACCCAGCAGGGCGACTTGCATAGttt                                            >  1:259690/1‑149 (MQ=255)                                                                                                             gCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑GGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTc                                        <  1:401049/151‑1 (MQ=255)                                                                                                              cGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑GGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTcc                                       <  2:590897/151‑1 (MQ=255)                                                                                                              cGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑GGGTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGGATGTTGGTGGAAAAGGGGTAAAGCGTGCGCAGATTGCGCTGAACCCCGCAGGTCGACTTGCATAGTTTGCTcc                                       >  1:916876/1‑151 (MQ=255)                                                                                                                aTGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑GGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTg                                     <  2:572962/151‑1 (MQ=255)                                                                                                                aTGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑GGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTg                                     <  2:661919/151‑1 (MQ=255)                                                                                                                aTGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑GGGTTGCACCCCATCAATGGTCCTGCAGGGCACAAAAGGGGGGAGTTGGGGAAAAAAGTGGTAAGCGGGCGCAGAATGCGGTGAAACCAACAGGGTGGCTTGCACAAttttttcc                                       >  2:809493/1‑144 (MQ=255)                                                                                                                   cACCGAGCGCCAGCAAGGTATTGGC‑GGGAAA‑GGTTTGCACCACAACATTTGGCATGCCGGGGACAAGAGGGGGATGGTGGTGAAAAAGGGGTAAAGCGGGCGCAGATTTGGCTGAAACCAGCAGGGCGACTTGCATAGTTTTCTCCCGcc                                   >  2:713546/1‑150 (MQ=255)                                                                                                                    aCCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑GGGTTGCACCACATCATTGGGCATGCAGTGGACAAGAGGGGGATGTTGGGGGAAAAGGGGTTAAGCGTGCGCAGATTGCGCTGAACCCCGCAGGGCGACTTGCATAGTTTGCTCCCGcc                                   >  2:854243/1‑149 (MQ=255)                                                                                                                       gAGCGCCAGCAAGGTATTGGC‑GGTAAA‑TGGTTTTACCACACCATTGGTCCAGGCGTGTACAAAAGGGGGAAGTTGGGGAAAAAAGGGTAACCAGTACGGAGAATTAGCGTGACCCAGGAGGGCGACTTTCCTAATTTGGTcccggccta                                >  2:554866/1‑142 (MQ=255)                                                                                                                         gcgcCAGCAAGGTATTGGC‑GGTAAA‑GGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGt                            <  2:1034500/151‑1 (MQ=255)                                                                                                                               gCAAGGTATTGGC‑GGTAAA‑GGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTgaaga                       <  1:649292/150‑1 (MQ=255)                                                                                                                               gCAAGGTATTGGC‑GGTAAA‑GGTTTGCACCACATCATTGGTCATGCAGGGCACAAGAGGGGGATGTTGGTGAAAAAGGGGTAACGCGTGCGCAGATTTCGCTGAACCCAGCAGGTCGACTTTCATAGTTTGCTCCTGCCATAACGCgaaga                       >  2:760065/1‑150 (MQ=255)                                                                                                                                 aaGGTATTGGC‑GGTAAA‑GGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGc                    <  1:1149137/151‑1 (MQ=255)                                                                                                                                    gTATTGGC‑GGTAAA‑GGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGTGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAAt                 <  2:477134/151‑1 (MQ=255)                                                                                                                                      aTTGGC‑GGTAAA‑GGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATg                <  2:549514/150‑1 (MQ=255)                                                                                                                                      aTTGGC‑GGTAAA‑GGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATg                <  2:385628/150‑1 (MQ=255)                                                                                                                                             ggTAAA‑GGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCtggt         <  2:71970/151‑1 (MQ=255)                                                                                                                                                     gTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGt  >  2:545289/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                                      CACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGC‑GGTAAA‑GGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGT  >  NC_000913/2185052‑2185348

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