Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,258,974	0	T	G	39.2%	128.3 / 14.8	75	G411G (GGT→GGG)	dinF	oxidative stress resistance protein; putative MATE family efflux pump; UV and mitomycin C inducible protein
Reads supporting (aligned to +/- strand):  ref base T (11/34);  new base G (29/0);  total (41/34)
Fisher's exact test for biased strand distribution p-value = 7.56e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAGATTCAGCAGCTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCT‑GGCATTAACCGTCTTTCTGGCGTTGCGCGGGC‑TTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGCCGCAACGTGACGG  >  NC_000913/4258827‑4259124                                                                                                                                                    |                                                                                                                                                         ccAGATTCAGCAGCTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGtttt                                                                                                                                                       <  2:744622/151‑1 (MQ=255) ccAGATTCAGCAGCTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGtttt                                                                                                                                                       <  1:851411/151‑1 (MQ=255) ccAGATTCAGCAGCTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGgttt                                                                                                                                                       >  1:392212/1‑151 (MQ=255) ccAGATTCAGCAGCTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTCGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGtttt                                                                                                                                                       <  1:728231/151‑1 (MQ=255)          gcagcTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACg                                                                                                                                              <  1:955834/151‑1 (MQ=255)           cagcTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGGCGCCCCCGGGTTTTCGCCGAc                                                                                                                                               >  2:1163743/1‑149 (MQ=255)            agcTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGGGGGCGCCGCCGGGTTTGCTTTGACGcc                                                                                                                                            >  2:1155289/1‑150 (MQ=255)             gctggctgACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGGCGCCGCCGGGTTTTCGCTGgcgccc                                                                                                                                           >  1:490012/1‑149 (MQ=255)                ggctgACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCtt                                                                                                                                        <  1:1057109/151‑1 (MQ=255)                  ctgACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGGCGCCGCAGGGTTTTCGCCGACGcccc                                                                                                                                          >  1:229327/1‑147 (MQ=255)                  ctgACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCCCCGGGTTTTCGCTGACGCAcc                                                                                                                                          >  1:229168/1‑147 (MQ=255)                      ccGTTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGTCGGCATGTTTATAGGCGCAACGCGTTCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTg                                                                                                                                  <  2:493602/151‑1 (MQ=255)                       cGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGGCGCCGCCGGGTTTTCGCTGACGCCCCCTACGc                                                                                                                                    >  1:1181813/1‑148 (MQ=255)                          tatcttatctGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGGCGCCGCCGGTTTTTGGCTTAAGCTCCTTACGCTg                                                                                                                                  >  1:1183474/1‑147 (MQ=255)                           atcttatctGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCg                                                                                                                               <  2:122848/149‑1 (MQ=255)                            tcttatctGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTTTGCTGACGGCCCTTACGCGGCCGTgg                                                                                                                            >  2:139781/1‑151 (MQ=255)                               tatctGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTgg                                                                                                                            <  2:306477/148‑1 (MQ=255)                               tatctGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTg                                                                                                                             <  2:844768/147‑1 (MQ=255)                                atctGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTg                                                                                                                         <  2:412696/150‑1 (MQ=255)                                    ggCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTggg                                                                                                                       <  2:108796/148‑1 (MQ=255)                                       aGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGGCGCCCCCGGGTTTTTGCTGACGCCCCTTACCCTGCCGTGGCTggg                                                                                                                       >  2:1189738/1‑145 (MQ=255)                                           gATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCt                                                                                                             <  2:1111880/151‑1 (MQ=255)                                                tGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTtgt                                                                                                         <  1:713855/150‑1 (MQ=255)                                                 gCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCCCCTTACGCTGCCCTGGCTGGGAAATCATTCttt                                                                                                           >  1:431656/1‑147 (MQ=255)                                                  ccGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCCCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTTGGCGCCGCCGGGTTTGGGCCGCAGCGCCTTCCCCCGCCCTGGGCGGGTAATCaggatctttggg                                                                                                      >  2:571488/1‑140 (MQ=255)                                                    ggtggtTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCCCAGGGTTTGCGCTGACGCCCCTTACGCTGCCGTGGGTGGGGGATCATGCTTTTtgg                                                                                                       >  2:941281/1‑148 (MQ=255)                                                     gtggtTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGGCGCCCCAGGGTTTTTGCTGAACGTCCCTACGCTGCCGGGGCTGGGTAATCAAGCTTTTtgggggg                                                                                                   >  2:955593/1‑147 (MQ=255)                                                      tggtTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCCGGTTTTTCGCCGACCCTCCTTACGCTTCCGTGGCTGGGGAATCATGCTTTgtgggt‑gg                                                                                                  >  1:942852/1‑150 (MQ=255)                                                        gtTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCT‑GGCAt                                                                                               <  1:294686/151‑1 (MQ=255)                                                           ggCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTGCGCTGACGCTCCCTACGCTGCCGGGGCTGGGGAATCATGCTTTTtgggt‑gggctt                                                                                              >  1:750927/1‑145 (MQ=255)                                                              gTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGGCGCCGCCGGGTTTGCGCTGACGCTCCTTACGCTGCCCTGGCTGGGGAATCATGCTTTTTGGCT‑GGCATTaaac                                                                                          >  2:647178/1‑148 (MQ=255)                                                              gTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGGCGCCGCAGGGTTTTCGCTGACGCTCCCTACGCTGCCGGGGGTGGGGGATCATTGTTTgtgggt‑ggggtttaacc                                                                                         >  2:375437/1‑142 (MQ=255)                                                                cTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTTGGCTGAGGCCCCTTAATCTGCCGAGGGTGGGGAAACCTTCCTTGtggct‑ggctct                                                                                              >  1:479525/1‑141 (MQ=255)                                                                  ggTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCT‑GGCATTAACCGt                                                                                        <  1:503044/148‑1 (MQ=255)                                                                  ggTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCT‑GGCATTAACCGTc                                                                                       <  1:708903/149‑1 (MQ=255)                                                                   gtgtTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGGCGCCGCACGGGTTGGGGTGACGATCCCTATGCTGCCCCGGCTGGGGAATCATTCTTTTTTGCT‑GGCCTTTAACGTCttt                                                                                    >  2:734235/1‑151 (MQ=255)                                                                     gtTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCT‑GGCATTAACCGTCTTTc                                                                                   <  1:10901/150‑1 (MQ=255)                                                                      ttATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGGGGGCGCCGCAGGGTTTGCGCTGACGCCCCTTACGCCGCCGTGGGTGGGGGATCATGCTTTTTTGCTGGGCTTTAACC‑TCttttt                                                                                  >  2:914858/1‑148 (MQ=255)                                                                        aTCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCT‑GGCATTAACCGTCTTTCTgg                                                                                <  1:420413/150‑1 (MQ=255)                                                                               ggACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCT‑GGCATTAACCGTCTTTCTGGCGTTgcg                                                                         <  1:616598/150‑1 (MQ=255)                                                                               ggACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCT‑GGCATTAACCGTCTTTCTGGCGTTgcg                                                                         <  2:527558/150‑1 (MQ=255)                                                                                gACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGGCGCCGCCGGGTTTGCGCTGACGCCCCTTACGCTGCCCGGGCTGGGGAAAAAAGCTTTGTGGCT‑GGCATTAAACGTCTTTTTTGCGGTGcgcg                                                                       >  1:888956/1‑151 (MQ=255)                                                                                    gCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGGCGCCGCCGGGTTTGCGCTGAAGCTCCCTACCCTGCCGGGGCTGGGGGATAAAGCCTTTTGGGT‑GGGATTTAACGTCTTTTTTGCGTTGTGCggg                                                                     >  2:1093781/1‑149 (MQ=255)                                                                                     cATGTTTATAGGCGCAACGCGTGCCACCGAGATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCT‑GGCATTAACCGTCTTTCTGGCGTTGCGCGGGC‑tt                                                                 <  1:648869/151‑1 (MQ=255)                                                                                     cATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCCGGGTTTGCGCTGACGCCCCTTACGCTGCCGGGGCTGGGGAAACATGCTTTGTGGGT‑GGCATTAACCCGCTTTTTGGGGGTGGGCGGGC‑tt                                                                 >  1:841653/1‑151 (MQ=255)                                                                                        gTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCT‑GGCATTAACCGTCTTTCTGGCGTTGCGCGGGC‑TTtc                                                               <  1:1180431/150‑1 (MQ=255)                                                                                          ttATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCCGGGTTTGTGCTGACGCTCCTTACGCTGCCGCGGCTGGGGAATCATGCTTTTTGGGTGGGCTTTAACC‑TCTTTTTTGCGGTTGGCGGGG‑ttttt                                                              >  2:174428/1‑147 (MQ=255)                                                                                          ttATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTTCGCCGACGCTCCCTACCCTGCCCTGGCTGGGGAATTATTCTTTGTGGCT‑GGGATTTAACCTCCTTTTGGCGGTGTGCGGGG‑ttt                                                                >  1:836599/1‑147 (MQ=255)                                                                                          ttATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTGTGCTGGCGCTCCTTAAGCTGCCGTGGGTGGGTAATCATGCTTTTTGGGT‑GGGATTAAACGTCTTTTTGGCGTTGGGCGGGC‑gtttttt                                                            >  1:1143338/1‑147 (MQ=255)                                                                                           tataGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCT‑GGCATTAACCGTCTTTCTGGCGTTGCGCGGGC‑TTTctct                                                            <  2:974945/150‑1 (MQ=255)                                                                                            ataGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCT‑GGCATTAACCGTCTTTCTGGCGTTGCGCGGGC‑TTTCTCTg                                                           <  2:383693/150‑1 (MQ=255)                                                                                                    aaCGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGGGCTGGCGCCCCTTAAGCTGCCCGGGCTGGGGAATCATGCTTTGTGGGT‑TGCATTTAACCTCTTTTTGGCGGTGCGCggggcttttcttgggg                                                        >  1:272625/1‑133 (MQ=255)                                                                                                    aaCGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTGCGCTGACGCTCCTTTCGCTGCCGTGGCGGGGGAATCATGCTTTTTGGGT‑GGCATTAACCGTCTTTTTGGGGTTGCGCGGGC‑TTTTTCTTGCGGGTAt                                                   >  2:71236/1‑150 (MQ=255)                                                                                                      cgcgTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCCGGGTTTTCGCTGACGCTCCTTACGCCGCCGGGGCTGGGGAATTATGCTTTTTGGCT‑GGGATTAACCGGCTTTTTGGGGTTGCGCGGGC‑TTTCTCTGGCGGCTATc                                                  >  1:949191/1‑149 (MQ=255)                                                                                                        cgTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTTCGCTGACGCCCCTTACGCTGCCGTGGCTGGGGAATAATGCTTTTTGGGT‑GGCATTTACCGTCTTTTTGGCGTTTCGCGGGC‑TTTTTCTTGCGGGTATCt                                                 >  1:408018/1‑148 (MQ=255)                                                                                                         gTGCCACCGAAATGCGTAACAGTATGGCGGTGGGCCCCGCCGGGTTTTCGCTGACGCCCCTTTCGCTGCCGTGGGTGGGGAATCATGCTTTTTGGGT‑GGGCTTTAACGTCTTTTTTGGGTTGCGCGGGG‑GTTTTCTGGGGGGTAttttg                                               >  2:1106506/1‑145 (MQ=255)                                                                                                           gCCACCGAAATGCGTAACAGTATGGCGGAGGCCCCCCCCGGATTTTCGCTGACGCCCCTTACCACGCCGGGGGCGGGGAAGCATTCTTTTTGGGT‑GGGATTAAACCTCTTTTTTGGGGTGCGGGGGCGTTTTTTCGGCGGC‑GTCTggcgg                                            >  2:88444/1‑150 (MQ=255)                                                                                                            cctccGAAATGCGTAACAGTATGGCGGTTGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCT‑GGCATTAACCGTCTTTCTGGCGTTGCGCGGGC‑TTTCTCTGGCGGCTATCTggcgg                                            <  1:667830/146‑1 (MQ=255)                                                                                                            ccACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCT‑GGCATTAACCGTCTTTCTGGCGTTGCGCGGGC‑TTTCTCTGGCGGCTATCTGgcggcg                                          <  1:808406/151‑1 (MQ=255)                                                                                                              aCCGAAATGCGTAACAGTATGGCGGTGGCCGCCCCAGGGTTTTCGCTGACGCTCCTTTCGCTGCCGTGGCCGGGTAATTATGCTTTTTGGCT‑GGCATTAACCGTCTTTTTGGGGTTTCGCGGGC‑TTTTTCTTGCGGCCTTCTGgcggcg                                          >  1:750713/1‑149 (MQ=255)                                                                                                                 gAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTTCGCTGACGCTCCTTACGCTGCCGTGGCTGGGGAATCATGCTTTGTGGCT‑GGCATTAAACGTCTTTTTTGCGTTGCGCGGGC‑TTTCTCTGGCGGCTATCTGGCGGCGTccc                                      >  1:1121043/1‑148 (MQ=255)                                                                                                                 gAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTGCGCTGGCGCTCCTTACGCTGCCGGGGCTGGGGAATCAAGCTTTTTGGCT‑GGCATTAAACGGCTTTCTGGCGTTGCGCGGGG‑TTTTTCTTGCGGCTATCTGggggcg                                          >  1:86581/1‑146 (MQ=255)                                                                                                                  aaaTGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCT‑GGCATTAACCGTCTTTCTGGCGTTGCGCGGGC‑TTTCTCTGGCGGCTATCTGGCGGCGTCACTg                                    <  1:215464/151‑1 (MQ=255)                                                                                                                    aTGCGTAACAGTATGGCGGTGGGCGCCCCAGGGTTTGCGCTGACGCTCCTTACGCTGCCCGGGCTGGGGAATTATGCTTTGTGGGT‑GGCGTTAACCGCTTTTTTGGCGTTGGGCGGGG‑GTTTTCTGGGGGGTATTTGGCGGGGTCACTgg                                   >  2:892151/1‑150 (MQ=255)                                                                                                                       cGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCT‑GGCATTAACCGTCTTTCTGGCGTTGCGCGGGC‑TTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCa                               <  2:1150496/151‑1 (MQ=255)                                                                                                                          aaCAGTATGGCGGTGGCCGCCCCAGGTTTTGCGCTTACGCCCCTTCCGCCGCCGTGGGTGGGGAAACATACCTTTGGGGT‑GGCCATAAACCTCTTTTCGGCGGTGGGCCGGG‑CTTCTCTTGGGGCTCTCTTGGGGGGGCACttg                                   >  2:287624/1‑142 (MQ=255)                                                                                                                                aTGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCT‑GGCATTAACCGTCTTTCTGGCGTTGCGCGGGC‑TTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTAcc                       <  2:1057238/150‑1 (MQ=255)                                                                                                                                aTGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCT‑GGCATTAACCGTCTTTCTGGCGTTGCGCGGGC‑TTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCt                      <  1:1029648/151‑1 (MQ=255)                                                                                                                                        ggCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCT‑GGCATTAACCGTCTTTCTGGCGTTGCGCGGGC‑TTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGc               <  2:351138/150‑1 (MQ=255)                                                                                                                                            gccgcAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCT‑GGCATTAACCGTCTTTCTGGCGTTGCGCGGGC‑TTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGCCGCa           >  1:706492/1‑150 (MQ=255)                                                                                                                                             ccgcAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCT‑GGCATTAACCGTCTTTCTGGCGTTGCGCGGGC‑TTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGCCGCAAc         <  1:171866/151‑1 (MQ=255)                                                                                                                                              cgcAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCT‑GGCATTAACCGTCTTTCTGGCGTTGCGCGGGC‑TTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGCCGCAACg        >  1:77822/1‑151 (MQ=255)                                                                                                                                              cgcAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCT‑GGCATTAACCGTCTTTCTGGCGTTGCGCGGGC‑TTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGCCGCAACg        <  1:564291/151‑1 (MQ=255)                                                                                                                                               gcAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCT‑GGCATTAACCGTCTTTCTGGCGTTGCGCGGGC‑TTTCTCTGGCGGCTATCTGGCGGCGTCGCTGGCGCAATGGTACCTGGTTTGCCGCAAc         <  1:247637/149‑1 (MQ=255)                                                                                                                                                cAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCT‑GGCATTAACCGTCTTTCTGGCGTTGCGCGGGC‑TTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGCCGCAACGTg      <  2:997827/151‑1 (MQ=255)                                                                                                                                                  ggTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCT‑GGCATTAACCGTCTTTCTGGCGTTGCGCGGGC‑TTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGCCGCAACGTGAc    >  1:959366/1‑151 (MQ=255)                                                                                                                                                   gTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCT‑GGCATTAACCGTCTTTCTGGCGTTGCGCGGGC‑TTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGCCGCAACGTGACg   <  1:301342/151‑1 (MQ=255)                                                                                                                                                   gTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCT‑GGCATTAACCGTCTTTCTGGCGTTGCGCGGGC‑TTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGCCGCAACGTGACg   >  1:698553/1‑151 (MQ=255)                                                                                                                                                   gTTTTGCGCTGACACTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCT‑GGCATTAACCGTCTTTCTGGCGTTGCGCGGGC‑TTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGCCGCAACGTGACg   <  1:592591/151‑1 (MQ=255)                                                                                                                                                    ttttGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCT‑GGCATTAACCGTCTTTCTGGCGTTGCGCGGGC‑TTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGCCGCAACGTGACgg  >  2:998872/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                         CCAGATTCAGCAGCTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCT‑GGCATTAACCGTCTTTCTGGCGTTGCGCGGGC‑TTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGCCGCAACGTGACGG  >  NC_000913/4258827‑4259124

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