Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,258,974	0	T	G	49.4%	137.7 / 12.5	85	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 (6/37);  new base G (42/0);  total (48/37)
Fisher's exact test for biased strand distribution p-value = 4.21e-18
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACCCAGATTCAGCAGCTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCT‑GGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGCCGCAACGTGACG  >  NC_000913/4258824‑4259123                                                                                                                                                       |                                                                                                                                                       aaCCCAGATTCAGCAGCTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATATTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGt                                                                                                                                                        <  1:1132136/151‑1 (MQ=255)    ccAGATTCAGCAGCTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGtttt                                                                                                                                                     <  1:1271588/151‑1 (MQ=255)     cAGATTCAGCAGCTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCCGGgtttt                                                                                                                                                    >  1:91034/1‑150 (MQ=255)      aGATTCAGCAGCTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTgc                                                                                                                                                   >  1:356227/1‑151 (MQ=255)       gATTCAGCAGCTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTgcg                                                                                                                                                  <  1:1363379/151‑1 (MQ=255)       gATTCAGCAGCTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTgcg                                                                                                                                                  <  1:944628/151‑1 (MQ=255)         ttCAGCAGCTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTTCGCt                                                                                                                                                >  1:600093/1‑151 (MQ=255)            agcagcTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCCGGTTTTGCGCTGAc                                                                                                                                             >  1:283993/1‑151 (MQ=255)             gcagcTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCCGGGTTTGCGCTGACg                                                                                                                                            >  1:949739/1‑151 (MQ=255)             gcagcTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACg                                                                                                                                            <  2:456795/151‑1 (MQ=255)              cagcTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTGCGCTGCCGc                                                                                                                                           >  1:207733/1‑151 (MQ=255)               agcTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCCGGGTTTGCGCTGACGcc                                                                                                                                          >  1:978912/1‑150 (MQ=255)               agcTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTTCGCTGACGcc                                                                                                                                          >  1:1346914/1‑150 (MQ=255)                gctggctgACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTc                                                                                                                                         <  1:419552/151‑1 (MQ=255)                   ggctgACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGGCGCCGCCGGGTTTGCGCTGACGCTCCtt                                                                                                                                      >  2:909287/1‑151 (MQ=255)                   ggctgACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCt                                                                                                                                       <  1:507125/150‑1 (MQ=255)                   ggctgACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTTCGCTGACGCTCCtt                                                                                                                                      >  1:954708/1‑151 (MQ=255)                    gctgACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTa                                                                                                                                     <  2:900338/151‑1 (MQ=255)                    gctgACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTTCGCTGACGCTCCtt                                                                                                                                      >  1:632851/1‑150 (MQ=255)                     ctgACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCCGGGTTTGCGCTGACGCTCCTTAc                                                                                                                                    >  1:928580/1‑151 (MQ=255)                        aCCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTTCGCTGACGCTCCTTACGCt                                                                                                                                 >  1:679983/1‑151 (MQ=255)                          cGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGGCGCCGCCGGGTTTTCGCTGACGCTCCCTTCGCTGc                                                                                                                               >  2:1026563/1‑151 (MQ=255)                          cGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGGCGCCCCCGGGTTTTGGCTGACGCCCCCTAACCTGc                                                                                                                               >  2:907960/1‑151 (MQ=255)                          cGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCCCAGGGTTTTGGCTGACGCCCCTTTcgccgc                                                                                                                               >  1:821126/1‑151 (MQ=255)                                  tatctGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTgg                                                                                                                          <  1:779045/148‑1 (MQ=255)                                  tatctGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTgg                                                                                                                          <  2:828413/148‑1 (MQ=255)                                  tatctGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCt                                                                                                                        <  2:708297/150‑1 (MQ=255)                                   atctGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGGCGCCGCCGGGTTTTTGCTGGCGCCCCTTACGCCGCCGGGGCggg                                                                                                                      >  2:714501/1‑151 (MQ=255)                                   atctGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTg                                                                                                                       <  1:133416/150‑1 (MQ=255)                                    tctGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCCGGGTTTTCGCTGACGCTCCCTACGCTGCCGGGGCTggg                                                                                                                     >  1:220177/1‑151 (MQ=255)                                     ctGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGGCGCCCCCGGGTTTGCGCTGACGCCCCTTACCCTGCCGTGGCTgggg                                                                                                                    >  2:1220268/1‑150 (MQ=255)                                       ggCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTa                                                                                                                   <  1:577211/150‑1 (MQ=255)                                       ggCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGGGGGCGCCCCCGGGTTTTCGCTGTCGCTCCTTAAGCTGCCgggggggggga                                                                                                                   >  2:194453/1‑143 (MQ=255)                                          aGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCCGGGTTTGCGCTGACGCCCCTTTCGCCGCCGGGGGCGGGGAATc                                                                                                                >  2:1249456/1‑150 (MQ=255)                                          aGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGCGGCCCCCGCAGGTTTTTCGCTGAAGCTCCCTACGCCGCCGTGGCTGGGGAATc                                                                                                                >  1:441540/1‑150 (MQ=255)                                            gTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATg                                                                                                             <  2:753421/151‑1 (MQ=255)                                             tGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTTCGCTGACGCCCCTTACGCTGCCGTGGCTGGGGAATCATGc                                                                                                            >  1:1330903/1‑151 (MQ=255)                                               aTTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCtt                                                                                                          <  1:289989/151‑1 (MQ=255)                                               aTTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCCCAGGGTTTTCGCTGACGCCCCTTACGCTGCCGTGGCTGGGGAATCATGCtt                                                                                                          >  1:777426/1‑151 (MQ=255)                                                ttttGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCCCCGGGTTTGCGCTTACGCTCCTTTCGCCGCCGGGGGCGGGGAATTATTCttt                                                                                                         >  2:421654/1‑151 (MQ=255)                                                  ttGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTtg                                                                                                        <  2:1323931/150‑1 (MQ=255)                                                        gtggtTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGtggct‑gg                                                                                                <  2:201405/151‑1 (MQ=255)                                                          ggtTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGtggct‑ggc                                                                                               <  1:8620/150‑1 (MQ=255)                                                           gtTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTTTGCTGGCGCTCCTTACGCCGCCGTGGGTGGGGAATCATGCTTTTTGGGT‑GGCAt                                                                                             >  1:1349031/1‑151 (MQ=255)                                                           gtTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTGCGCTGACGCTCCTTACGCTGCCGGGGCTGGGGAATCATGCTTTGTGGCT‑GGCAt                                                                                             >  1:761347/1‑151 (MQ=255)                                                                cGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCT‑GGCATTAAcc                                                                                        <  2:196754/151‑1 (MQ=255)                                                                 gTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCT‑GGCATTAACCg                                                                                       <  2:914895/151‑1 (MQ=255)                                                                  tCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGGCGCCCCCGGGTTTTGGCTGACGCTCCTTTCCCTGCCGGGGGTGGGGAATCATGCTTTTTGGCT‑GGGATTAAACGt                                                                                      >  2:311675/1‑151 (MQ=255)                                                                     ggTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCT‑GGCATTAACCGt                                                                                      <  2:920920/148‑1 (MQ=255)                                                                       tgtTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGGAATCATGCTTTTTGGCT‑GGCATTAACCGTCtttt                                                                                 >  1:846289/1‑150 (MQ=255)                                                                        gtTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCT‑GGCATTAACCGTCTTTca                                                                                <  2:974333/151‑2 (MQ=255)                                                                        gtTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCT‑GGCATTAACCGTCTTTCt                                                                                <  2:413392/151‑1 (MQ=255)                                                                        gtTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCCCCGGGTTTGCGCTGACGCTCCCTACGCCGCCGTGGCTGGGGAATCATGCTTTGTGGCTGGGCTTTAACCG‑CTTTCt                                                                                >  2:136527/1‑151 (MQ=255)                                                                          tATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCT‑GGCATTAACCGTCTTTCTg                                                                               <  1:1269881/150‑1 (MQ=255)                                                                            tCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCCGGGTTTGCGCTGACGCTCCTTTCGCTGCCGTGGGTGGGGAATTATGCTTTGTGGCT‑GGGATTTAACGTCTTTTTGGCg                                                                            >  1:525111/1‑151 (MQ=255)                                                                               gctgGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCT‑GGCATTAACCGTCTTTCTGGCGtt                                                                          <  2:1108820/150‑1 (MQ=255)                                                                                ctgGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCT‑GGCATTAAACCTCTTTCTGGCGTTgc                                                                        >  1:686223/1‑151 (MQ=255)                                                                                  ggACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCT‑GGCATTAACCGTCTTTCTGGCGTTgcgc                                                                      <  1:770949/151‑1 (MQ=255)                                                                                    aCGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCT‑GGCATTAACCGTCTTTCTGGCGTTGCGCgg                                                                    <  2:1000153/151‑1 (MQ=255)                                                                                            tttATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTGCGCTGAAGCTCCTTAAGCTGCCGTGGGTGGGGAATAAAGCTTTGTGGCT‑GGCATTAACCCGCTTTTTGGCGGTGCGCGGGCTTtctc                                                            >  2:279160/1‑151 (MQ=255)                                                                                                taGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGGAATCATGCTTTGTGGCT‑GGCATTAACCGGCTTTTTTGCGTTGCGCGGGCTTTCTCTggc                                                        >  1:1103183/1‑151 (MQ=255)                                                                                                 aGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCT‑GGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTggc                                                        <  1:1170115/150‑1 (MQ=255)                                                                                                 aGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCT‑GGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTggc                                                        <  1:1199783/150‑1 (MQ=255)                                                                                                   gcgcAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCT‑GGCATTAAACGTCTTTTTGGCGGTGCGCGGGGTTTTTCTggcggc                                                     >  1:1357047/1‑151 (MQ=255)                                                                                                      cAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCCGGGTTTTCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCT‑GGCATTAACCGTCTTTTTGGCGTTTCGCGGGGTTTTTCTGGCGGCTAt                                                  >  1:816399/1‑151 (MQ=255)                                                                                                      cAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCCGGGTTTTCGCTGACGCTCCTTACGCCGCCGTGGCTGGGTAATCATGCTTTTTGGGT‑GGCATTAACCGGCTTTCTGGCGGTGGGCGGGGTTTTTCTGGGGGCttt                                                  >  1:1255888/1‑149 (MQ=255)                                                                                                      cAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCT‑GGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTAt                                                  <  2:791402/151‑1 (MQ=255)                                                                                                       aaCGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTGCGCTGACGCTCCTTACGCTGCCGCGGCTGGGGAATCATGCTTTTTGGCT‑GGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTTGCGGCTATc                                                 >  1:635830/1‑151 (MQ=255)                                                                                                        aCGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCT‑GGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCt                                                <  2:1346590/151‑1 (MQ=255)                                                                                                         cgcgTGCCACCGAAATGCGTAACAGTATGGCGGTGGGCCCCGCCGGGTTTTCGCCGACGCCCCTTACGCCGCCGTGGCTGGGGGATCATGCTTTTTGGGT‑GGGATTAAACGTCTTTTTGGCGTTTCGCGGGGGTTTTCTTGGGGCCATCtt                                               >  2:1233001/1‑150 (MQ=255)                                                                                                         cgcgTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCT‑GGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTg                                               <  1:243508/151‑1 (MQ=255)                                                                                                         cgcgTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTTCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTTTGGCT‑GGCATTTACCGGCTTTTTTGCGTTGGGCGGGGTTTTTCTGGCGGCTAtttt                                               >  1:833098/1‑148 (MQ=255)                                                                                                          gcgTGCCACCGAAATGCGTAACAGTATGGCGGTGGGCGCCGCAGGGTTTTCGCTGACGCTCCTTAAGCTGCCGGGGGTGGGGAAACATGCTTTGTGGGG‑GGGATTAAACGGCTTTTTTGGGGTGGGGGGGGGTTTTCTGGGGGGTATCTgg                                              >  2:535177/1‑151 (MQ=255)                                                                                                             tGCCACCGAAATGCGTAACAGTATGGCGGTGGGCGCCGCCGGGTTTTCGCTGACGCTCCCTACGCTGCCGGGGGTGGGGAATAATGCTTTTTTGGT‑GGGCTTAAACGGCTTTTTTGCGGTGCGCGGGGGTTTTCTGGCGGCGATCTtgcgg                                           >  2:770345/1‑151 (MQ=255)                                                                                                               ccACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTGCGCTGACGCTCCTTACGCTGCCGGGGCTGGGGAATCATGCTTTGTGGCT‑GGGATTAACCCTCTTTTTTGCGTTGCGCGGGGTTTCTCTGGCGGCTATCTGgcggcg                                         >  1:1096272/1‑151 (MQ=255)                                                                                                                      aaTGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCT‑GGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTgg                                  <  2:1326988/151‑1 (MQ=255)                                                                                                                        tGCGTAACAGTATGGCGGTGGCCCCCGCCGGGTTTGCGCTGGCGCTCCTTACGCTGCCGGGGCTGGGGGATCATGCTTTTTGGCT‑GGCATTAAACGTCTTTTTTGCGGTGCGCGGGGTTTTTCTGGCGGGTTTCTggcgggggccccggcg                                >  2:457211/1‑140 (MQ=255)                                                                                                                         gCGTAACAGTATGGCGGTGGGCCCCGCAGGGTTTTGGCTGACGCTCCCTACGCTGCCGTGGGTGGGGAATCATGCTTTGGGGGT‑GGCATTAAACGTCTTTTTGGCGGTGCGGGGGGCTTTTCTGGCGGCGATCTGGCGGCGTCCCTGgcgc                               >  2:286052/1‑151 (MQ=255)                                                                                                                              aCAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCT‑GGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATgg                          <  1:1047958/151‑1 (MQ=255)                                                                                                                               cAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCT‑GGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGt                         <  1:25769/151‑1 (MQ=255)                                                                                                                               cAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCT‑GGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGt                         <  2:224067/151‑1 (MQ=255)                                                                                                                                aGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCT‑GGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTa                        <  2:1180526/151‑1 (MQ=255)                                                                                                                                aGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCT‑GGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTa                        <  1:691966/151‑1 (MQ=255)                                                                                                                                  tATGGCGGTGGCCGCCGCAGGGTTTTCGCTGACGCTCCTTTCGCTGCCGGGGCGGGGGAATCATGCTTTGTGGCT‑GGCATTAAACGTCTTTCTTGCGTTGGGCGGGCTTTCTCTGGCGGGCATTTTGCGGGGGCCCTGGCGCAATGGTAcc                      >  2:436629/1‑151 (MQ=255)                                                                                                                                         gTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCT‑GGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTg               >  1:1305224/1‑151 (MQ=255)                                                                                                                                                 cgcAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCT‑GGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGCCGCAACg       >  2:581276/1‑151 (MQ=255)                                                                                                                                                    aGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCT‑GGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGCCGCAACGTg     >  1:568441/1‑150 (MQ=255)                                                                                                                                                      gTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCT‑GGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGCCGCAACGTGACg  <  1:136532/151‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                       AACCCAGATTCAGCAGCTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCT‑GGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGCCGCAACGTGACG  >  NC_000913/4258824‑4259123

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