Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,492,834	0	T	G	44.8%	149.5 / 24.3	87	G89G (GGT→GGG)	tsgA	putative transporter
Reads supporting (aligned to +/- strand):  ref base T (11/37);  new base G (39/0);  total (50/37)
Fisher's exact test for biased strand distribution p-value = 6.48e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCCTGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTAC‑ATTCCTGGTAACACAAATGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTACCG  >  NC_000913/3492686‑3492979                                                                                                                                                      |                                                                                                                                                   tGCCTGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGGCGGgtt                                                                                                                                                  >  1:517191/1‑151 (MQ=255)     tGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGGCGGgtt                                                                                                                                                  >  1:1341002/1‑147 (MQ=255)        ttCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCCGGGTTGATGtt                                                                                                                                            >  1:862381/1‑150 (MQ=255)         tCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCATCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCCAGTTTTGTGtt                                                                                                                                            >  2:1071777/1‑149 (MQ=255)          ccAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCCGGGTTGGTGtt                                                                                                                                            >  1:1708900/1‑148 (MQ=255)            aGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCCGGTTTGATGTTCAGcc                                                                                                                                       <  1:372119/151‑1 (MQ=255)            aGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCCGGGTTGTTGTTTAGcc                                                                                                                                       >  2:970821/1‑151 (MQ=255)             gTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCCGGTTTGATGTTCAGCca                                                                                                                                      <  1:552255/151‑1 (MQ=255)               aTGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCCGGTTTGATGTTCAGCcaca                                                                                                                                    <  2:33410/151‑1 (MQ=255)               aTGAGAAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCCGGGTTGTTTTTCAGcc                                                                                                                                       >  1:1150352/1‑148 (MQ=255)                 gAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCCGGGTTTATGTTCAGCCACAGc                                                                                                                                  >  2:1501684/1‑151 (MQ=255)                   gTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCCGGTTTGATGTTCAGCCACAGCCt                                                                                                                                <  1:724928/151‑1 (MQ=255)                   gTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCCGGTTTGATGTTCAGCCACAGCCt                                                                                                                                <  1:773678/151‑1 (MQ=255)                   gTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCCGGTTTGATGTTCAGCCACAGCCt                                                                                                                                <  2:693742/151‑1 (MQ=255)                   gTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCCGGTTTGATGTTCAGCCACAGCCt                                                                                                                                <  2:916193/151‑1 (MQ=255)                        accttcaccttcCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCCGGGTTTATGTTCAGCCACAACCTcgcg                                                                                                                            >  2:1018451/1‑150 (MQ=255)                         ccttcaccttcCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCCGGTTTGATGTTCAGCCACAGCCTGgcgc                                                                                                                           <  1:979143/150‑1 (MQ=255)                          cttcaccttcCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCCGGTTTGATGTTTAGCCACCACCCGgcgc                                                                                                                           >  1:89910/1‑149 (MQ=255)                            tcaccttcCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTg                                                                                                                         <  2:882940/149‑1 (MQ=255)                             caccttcCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCCCGGTTTTTGTTTAGCCACAGCCTGGCGCTGtt                                                                                                                       >  2:614635/1‑150 (MQ=255)                              accttcCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTtc                                                                                                                      <  2:237556/150‑1 (MQ=255)                                 ttcCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCCGGGTTGGTGTTCAGCCACAACCCGGCGCTGTTTTcgg                                                                                                                  >  2:545867/1‑151 (MQ=255)                                  tcCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCCGGTTTGATGTTCAGCCACAGCCCGGCGCTGTtttccg                                                                                                                  >  2:233604/1‑148 (MQ=255)                                      cAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTcggcgg                                                                                                               <  1:1554183/149‑1 (MQ=255)                                      cAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCggcggc                                                                                                              <  2:963972/150‑1 (MQ=255)                                      cAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGgcggcg                                                                                                             <  2:851427/151‑1 (MQ=255)                                       aaCGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCCGGGTTTATGTTCAGCCACCACCTGGCGCCGTTTTCGgcggcgg                                                                                                            >  2:1315320/1‑150 (MQ=255)                                       aaCGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCCGGGTTTATGGTTAGGCACAGACCGGGGGTGTTTTCGgcggcgg                                                                                                            >  1:983439/1‑150 (MQ=255)                                        aCGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCCGGGTTTATTTTTATCCCCAACCCGGCGCTGTtttccgggggggg                                                                                                           >  1:1641506/1‑142 (MQ=255)                                         cGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATg                                                                                                          <  2:112932/151‑1 (MQ=255)                                         cGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCCGGGTTGTTGGTCAGCCCCAGCCTGGCGCTGTTCTccgcgggggt                                                                                                           >  2:386087/1‑146 (MQ=255)                                           ccGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTTCCGGTTTGATGTTTAGCCACAGCCTGGCGCTGTTTTCGgcggcggggtt                                                                                                        >  1:822441/1‑146 (MQ=255)                                               cATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGGCGGGTTGATGTTTAGCCACAGCCTGGCGCTGTTTTCGGCGGGGAAGTTCAtt                                                                                                    >  2:102078/1‑151 (MQ=255)                                                aTTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTTGCGGGTTTG‑TGTTTAGCCCCACCCTGGGGCGGTTCTTGGCGGCGTTGTTTAtt                                                                                                    >  1:1177255/1‑150 (MQ=255)                                                  tttAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCCGGGTTTGTGTTTAGCCACAACCCGGCGCTGTTCTCCGCGGCGGTggttctttt                                                                                                  >  2:462536/1‑142 (MQ=255)                                                   ttAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATtctc                                                                                                 <  2:1370703/150‑1 (MQ=255)                                                     aaTCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTTTCGGGTTGCTGGTCAACCACAGCCCTGTGCTGTTGTTGGCGGGGATGCTCTTTCTCggg                                                                                              >  2:1192364/1‑151 (MQ=255)                                                     aaTCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTCTCggg                                                                                              <  1:1047976/151‑1 (MQ=255)                                                           tATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCCGGGTTTTTGTTTAGCCACACCCCGGCGCTGTTCTCGGCGGGGGTGGTTATTCTCGggggggg                                                                                         >  1:366685/1‑149 (MQ=255)                                                            aTCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTCTCGGGGTGGTc                                                                                        <  1:442866/150‑1 (MQ=255)                                                             tCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCCGGGTTTTTGTTTCGCCACAACCTGGCGCTGTTTTCGGCGGGGGGGTTTCTTTTCGgggggggcgg                                                                                      >  1:787763/1‑147 (MQ=255)                                                                 ccTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTCTCGGGGTGGTCAGCg                                                                                    <  2:235285/149‑1 (MQ=255)                                                                 ccTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGc                                                                                  <  1:196599/151‑1 (MQ=255)                                                                 ccTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCCGGGTTGATGTTTAGCCACAGCCCGGCGCTGTTTTCGGCGGCGGTGTTTATTCTCTGGGGGGGCAGCgg                                                                                   >  2:910476/1‑150 (MQ=255)                                                                   tCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCCGGGTTTATGTTTAGCCACAGCCTCGCGCTGTTTTCGGCGGCGGTGTTCCTTTCCGGGGGGGGTAGCGGCAt                                                                                >  2:694486/1‑151 (MQ=255)                                                                    cAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTTCCGGGTTGATGTTTAGCCACAGCCTGGCGCTGTTTTCGGCGGGGGTGTTTTTTCTCGGGGTGGTCAGCGGCAt                                                                                >  2:843063/1‑150 (MQ=255)                                                                    cAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCCGGGTTGATGTTTAGCCCCAACCTGGCGCTGTTCTCGGCGGGGGTGTTTATTCTCGGGGGGGGCAGGGGGAt                                                                                >  1:1590885/1‑150 (MQ=255)                                                                       cGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGGTTTGCCGGTTT‑ATGTTTAGCCACAGCCCGGCGCTGTTTTCGGCGGGGGTGTTTATTCTCGGGGGGGGCAGCGGGATCAcc                                                                            >  2:1458678/1‑151 (MQ=255)                                                                       cGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTTCCGGGTTGGTGTTCAGCCCCAACCTGGCGCTGTTCTCGGCGGGGGTGTTTTTTTTCGGGGGGGGCAGCGGCATCAcc                                                                            >  2:1631076/1‑151 (MQ=255)                                                                        gCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCa                                                                           <  2:473192/151‑1 (MQ=255)                                                                           tGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCCGGGTTGGTGTTTAGCCACAACCCGGCGCTGTTCTCGGCGGCGGTGTTTATTTTCGGGGGGGGCAGCGGCCTcccc                                                                            >  1:393269/1‑147 (MQ=255)                                                                           tGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCCGGGTTGATTTTTAGCCACCACCCGGCGCTGTTTTCGGCGGGGGCGTTCTTTTTCGGGGGGGGCAGCgggcgcacc                                                                            >  1:1631710/1‑140 (MQ=255)                                                                            ggCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATg                                                                         <  1:1071691/149‑1 (MQ=255)                                                                            ggCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGt                                                                        <  2:1248904/150‑1 (MQ=255)                                                                            ggCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTc                                                                       <  1:297043/151‑1 (MQ=255)                                                                             gCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCg                                                                      <  2:755170/151‑1 (MQ=255)                                                                             gCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCg                                                                      <  2:392019/151‑1 (MQ=255)                                                                               tGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCCGGGTTGATGTTTAGCCCCAGCCTGGCGCTGTTCTCGGCGGGGGTGTTCATTCTCGGGGGGGGCCGCGGCATCACCCTGTCggt                                                                    >  2:1371979/1‑149 (MQ=255)                                                                                gATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGAt                                                                    <  2:1153340/150‑1 (MQ=255)                                                                                   ggAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTTCCGGGTTTATGTTTAGCCACAGCCCGGGGCTGTTTTCGGCGGGGGTGTTCCTTCTCGGGGGGGGCAGCGGCGACACCATGTCGGtt                                                                   >  2:565824/1‑148 (MQ=255)                                                                                   ggAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑GTGCCGGGTTGGAGTTTAGCCACAGCCCTGCGCTGCTTTCGGCGGGGGTGTTTCTTCTCGGGGGGGGCCGCGGGATCCCCCTGttggttgg                                                                 >  1:305156/1‑150 (MQ=255)                                                                                    gAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTa                                                               <  1:1117999/151‑1 (MQ=255)                                                                                        tCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCCGGGTTGATGTTTAGCCCCAGCCTGGCGCTGTTCTCGGCGGGGGTGTTTATTTTTGGGGGGGTCAGCGGCCTCACCATGTCGATTTGTTC‑Att                                                          >  1:898175/1‑151 (MQ=255)                                                                                          gTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGGCGGGTTGATGTTTAGCCCCAGCCCGGCGCTGTTTTCGGCGGCGGGGTTTATTTTCTGGGGGGGCAGGGGCGTCACCATTTCGGTTGGGAC‑Att                                                          >  1:1339142/1‑149 (MQ=255)                                                                                          gTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTAC‑ATTcc                                                        <  2:850161/151‑1 (MQ=255)                                                                                          gTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCCGGGTTGATGTTTAGCCACAGCCTGGCGCTGTTTTCGGCGGGGGTGTTTATTTTCGGGGGGGTCAGCGGGATCACCCTGTCGGTTGGTAC‑ATTcc                                                        >  2:1381637/1‑151 (MQ=255)                                                                                            cccGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCCGGGTTTATGTTCAGCCACAACCTGGCGCCGGTCTCTGCGGCGGGGGTTATTTTCGGGGGGGGGAGCGGCCCCAACATGTCGATTGGGAC‑Att                                                          >  2:219860/1‑147 (MQ=255)                                                                                                  gAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTAC‑ATTCCTGGTAaca                                                <  1:380127/151‑1 (MQ=255)                                                                                                        gCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGTCGGGTTGATTTTTTGCCCCAACCTGGCGCCGGTTTCGGGGGGGCCGTTCTTTTTCGGGGGGGGCAGCGGGATCACCAAGTCTATTTGGAC‑ATTTCTTGTaaaa                                                >  2:1414031/1‑143 (MQ=255)                                                                                                          aGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGTCGGGTTGATGTTCAGCCCCAGCCTGGCGCTGGTCTCGGCGGGGGTGTTTTTTTTCGGGGGGGTCCGGGGCCTCTCCACGTCGTTTTGTACAATTCCTGTTaaacca                                              >  1:1556777/1‑142 (MQ=255)                                                                                                           gTTACGTTTTGGCTTTCTCCTGAGGGTGCTGGCTG‑TTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTAC‑ATTCCTGCTAACACAAATGTa                                        <  1:1101789/150‑1 (MQ=255)                                                                                                            ttACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCCGGGTTGATGTTTAGCCACAGCCTGGCGCTGTTCTCGGCGGCGGTGTTTATTTTCGGGGGGGTCAGCGGCATCACCATGTCGGTTGGTAC‑ATTTCTGGGAAAACaaa                                            >  2:292987/1‑145 (MQ=255)                                                                                                                gTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTAC‑ATTCCTGGTAACACAAATGTATGAAgg                                  <  2:1709016/151‑1 (MQ=255)                                                                                                                   ttGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTAC‑ATTCCTGGTAACACAAATGTATGAAGGGc                                <  1:1388240/150‑1 (MQ=255)                                                                                                                     ggCTTTCTCCTGATGGTGCTGGCGG‑TTGCCGGGTTGATGTTTAGCCACACCCTGGCGCTGTTTTCGGCGGCGGTGTTCATTCTCGGGGGGGGCAGCGGCATCACCATGTCGATTGGGAC‑ATTCCTTGGAAAACAAAAGTATGAAggg                                 >  2:960411/1‑147 (MQ=255)                                                                                                                        tttCTCCTGATGGTGCTGGCGG‑GTGCCGGGTTTGTGTTCAGCCACAGCCCTGGGGTGTTCTCGGCGGCGGTTTTTCTTTTCGGGGGGGGCCGCGGCCTCACCCTGTCGGTTTGGTC‑ATTCCTTGGAACACAAAAGTATGAAGGGGGTCAg                           >  1:254484/1‑150 (MQ=255)                                                                                                                         ttCTCCTGATGGTGCTGGCGG‑TTGCCGGGTTTATGTTTAGCCACAACCTGGCGCTGTTCTCGGCGGCGGTGTTCATTTTCGGGGGGGGCAGCGGCATCACCAAGTCGATTGGGAC‑ATTTCTGGGAAAACAAAATTATGAAagggggca                            >  2:53474/1‑143 (MQ=255)                                                                                                                          tctcCTGATGGTGCTGGCGG‑TTGGCGGGTTGATGTTCAGCCACAGCCCGGCGCTGTTCTCGGGGGGGGTGTTTATTTTCGGGGGGGGGAGCGGCATCACCATGTCGGGTGGGAC‑ATTTCTGGGGAAACAAATGTATGAAggggg                               >  1:662583/1‑142 (MQ=255)                                                                                                                             ccTGATGGTGCTGGCGG‑TTGCCGGGTTTATGTTCAGCCACAGCCCGGCGCTGTTTTCGGCGGGGGTGGTTATTCTCTGGGGGGGGAGCGGCATCACCCTGTCGGTTGGGAC‑ATTCCTGGTAACACAAAAGTATGAAGGGGGGCCGCGCgg                      >  1:1031915/1‑150 (MQ=255)                                                                                                                                gATGGTGCTGGCGG‑TTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTAC‑ATTCCTGGTAACACAAATGTATGAAGGGCGTCAGCGCGGTTcc                  <  2:1419397/151‑1 (MQ=255)                                                                                                                                 aTGGTGCTGGCGG‑TTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTAC‑ATTCCTGGTAACACAAATGTATGAAGGGCGTCAGCGCGGTTccc                 <  2:1684756/151‑1 (MQ=255)                                                                                                                                 aTGGTGCTGGCGG‑TTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTAC‑ATTCCTGGTAACACAAATGTATGAAGGGCGTCAGCGCGGTTccc                 <  1:577546/151‑1 (MQ=255)                                                                                                                                         ggCGG‑TTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTAC‑ATTCCTGGTAACACAAATGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTa         <  2:912840/151‑1 (MQ=255)                                                                                                                                          gCGG‑TTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTAC‑ATTCCTGGTAACACAAATGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTAt        >  2:138924/1‑151 (MQ=255)                                                                                                                                               ttGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTAC‑ATTCCTGGTAACACAAATGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTAc    >  1:1380244/1‑151 (MQ=255)                                                                                                                                                tGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTAC‑ATTCCTGGTAACACAAATGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTAcc   >  2:545913/1‑151 (MQ=255)                                                                                                                                                tGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTAC‑ATTCCTGGTAACACAAATGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTAcc   >  1:962071/1‑151 (MQ=255)                                                                                                                                                tGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTAC‑ATTCCTGGTAACACAAATGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTAc    <  2:393297/150‑1 (MQ=255)                                                                                                                                                  ccGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTAC‑ATTCCTGGTAACACAAATGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTACCg  >  1:864024/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                   TGCCTGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTGCTGGCGG‑TTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTCGGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTAC‑ATTCCTGGTAACACAAATGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTACCG  >  NC_000913/3492686‑3492979

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