Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,258,974	0	T	G	48.3%	87.0 / 15.3	58	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 (4/26);  new base G (28/0);  total (32/26)
Fisher's exact test for biased strand distribution p-value = 1.25e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATTCAGCAGCTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGCCGCAACGTGACG  >  NC_000913/4258830‑4259123                                                                                                                                                 |                                                                                                                                                      gATTCAGCAGCTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGGCGCCGCTGGgttt                                                                                                                                                    >  1:393207/1‑148 (MQ=255) gATTCAGCAGCTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCCCCCCAGGGTTTtcg                                                                                                                                                 >  2:841571/1‑151 (MQ=255)    tCAGCAGCTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCt                                                                                                                                               <  1:669258/150‑1 (MQ=255)       gcagcTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGGCGCCGCCGGGTTTgcggcgac                                                                                                                                            >  2:973828/1‑145 (MQ=255)         agcTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCt                                                                                                                                         <  1:176892/151‑1 (MQ=255)             ggctgACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCtt                                                                                                                                     <  2:429866/151‑1 (MQ=255)               ctgACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTTCGCTGACGCCCCtt                                                                                                                                     >  1:337268/1‑149 (MQ=255)                tgACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCACGGTTTGCGCTGACGCTCCTTAc                                                                                                                                   >  2:1014985/1‑150 (MQ=255)                tgACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGGGGGCGCCGCCGGGTTTGCGCCGGCGCTCCTTAc                                                                                                                                   >  2:545427/1‑150 (MQ=255)                      cTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGGCGCCGCCGGGTTTTCGCTGACGCTCCTTACgccgcc                                                                                                                             >  2:1125319/1‑150 (MQ=255)                        atcttatctGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCg                                                                                                                            <  2:214942/149‑1 (MQ=255)                             atctGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTg                                                                                                                      <  1:208442/150‑1 (MQ=255)                                tGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCCCCGCAGGGTTTGCGCTGACGCCCCTTACGCTGCCGTGGCTggg                                                                                                                    >  1:1320684/1‑149 (MQ=255)                                  gCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAAt                                                                                                                <  1:1057893/151‑1 (MQ=255)                                     ggTGATTTTGCCGGTGGTTGGGGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGACGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCAt                                                                                                             <  2:832848/151‑1 (MQ=255)                                     ggTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCAt                                                                                                             <  1:244630/151‑1 (MQ=255)                                     ggTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGGAAACa                                                                                                              >  1:991065/1‑150 (MQ=255)                                      gTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGGCGCCGCAGGGTTTGCGCTGACGCTCCCTACGCTGCCGGGGGTGGGGAATTAt                                                                                                             >  1:1219000/1‑150 (MQ=255)                                      gTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATg                                                                                                            <  1:1172401/151‑1 (MQ=255)                                      gTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTGCGCTGACGCTCCCTACGCTGCCGTGGGTGGGTAATCa                                                                                                              >  1:1134260/1‑149 (MQ=255)                                        gATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGGTTGCGCTGGCGCTCCTTACGCCGCCGTGGCCggg                                                                                                                    >  2:381205/1‑141 (MQ=255)                                         aTTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGGCGCCGCCGGGTTTGCGCTGGCGCCCCCTACGCTGCCCGGGGTGGGGGATCATGc                                                                                                           >  2:164925/1‑149 (MQ=255)                                           tttGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTGCGCTGACGCTCCTTACGCCGCCGGGGGTGGGGAATCATGCttt                                                                                                        >  1:459193/1‑150 (MQ=255)                                              gCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGtgg                                                                                                    <  1:276179/151‑1 (MQ=255)                                                            tCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTTCGCTGACGCTCCCTACGCTGCCGTGGCTGGGGAATCATGCTTTTTGGCTGGGATTaaac                                                                                        >  1:440052/1‑147 (MQ=255)                                                             cTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTTCGCTGACGCTCCTTACGCCGCCGGGGCTGGGGAATCATGCTTTGTGGCTGGCATTAAcc                                                                                        >  1:405020/1‑148 (MQ=255)                                                                  gtTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGGAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCt                                                                                >  1:389451/1‑151 (MQ=255)                                                                     aTCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCCGGGTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGGAATCATGCTTTTTGGCTGGCCTTAACCCTCttttttg                                                                              >  2:370471/1‑146 (MQ=255)                                                                      tCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTTCGCTGACGCCCCTTACGCTGCCGTGGCTGGGGAATCATGCTTTTTGGGTGGCATTAACCGTCTTTCTTGCg                                                                            >  1:1193230/1‑151 (MQ=255)                                                                       ctgctgGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGGCGCCGCAGGGTTTGCGCTGACGCCCCTTACGCTGCCGTGGGTGGGGAATCATGCTTTGTGGCTGGCATTAACCGTCTTTTTGGCg                                                                            >  1:431414/1‑150 (MQ=255)                                                                       ctgctgGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTGCGCTGACGCTCCCTTCGCTGCCGTGGCTGGGGAATCATGCTTTTTGGGTGGCATTAACCGTCTTTTTGGCg                                                                            >  2:891507/1‑150 (MQ=255)                                                                         gctgGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCTGCCGCAGGTTTTGCGCTGACGATCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGtt                                                                          <  2:299659/150‑1 (MQ=255)                                                                         gctgGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTg                                                                         <  2:431611/151‑1 (MQ=255)                                                                          ctgGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGGCGCCGCAGGGGTTGCGCTGACGCTCCTTACGCTGCCGGGGGTGGGGAAACATGCTTTGTGGGGGGGGTTAAACGTCTTTTTGGGGtt                                                                          >  2:232847/1‑149 (MQ=255)                                                                              aCGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGACGCCGCCCCCGGGTTTGCGCTGACGCTCCCTACCCTGGCGTGGCTTGGGAATCCTGCTTTGTGGCCGGCATTTACCGCCTTTCTGGCGGTGGGCgg                                                                    >  1:1077592/1‑151 (MQ=255)                                                                                ggCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCCCCGGGTTTGCGCTGACGCTCCTTTCTCTGCCGCGGCTTGGGAAACCTGATTTTTGGCTGGCATTTAACCCCTTTCTTGGGTTGCGCggg                                                                   >  2:683545/1‑150 (MQ=255)                                                                                   aTGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGCGGCCGCCGCAGGGTTTGCGCCGTCCCCCCCTACGCCGCCGTGGCTGGGGAAACAAGCTTTGTGGCTGGGATTAAACGCCTTTCTTGCGTTGCGCGGGGttt                                                               >  1:873938/1‑151 (MQ=255)                                                                                     gTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTtct                                                             <  1:232772/151‑1 (MQ=255)                                                                                           aGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTggc                                                        <  1:758494/150‑1 (MQ=255)                                                                                               gcAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCCGGGTTTTCGCTGACGCTCCTTACGCTGCCGTGGGTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTTGCGTTGCGCGGGCTTTCTCTggcgg                                                      >  1:881699/1‑148 (MQ=255)                                                                                                cAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTAt                                                  <  1:1326699/151‑1 (MQ=255)                                                                                                cAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTGCGCTGACGCCCCTTACGCTGCCGCGGGTGGGTAATCCTTCTTTTTGGCTGGGCTTTACCCTCTTTTTTGCGGTGGGCGGGCTTTCTTTggcgg                                                      >  2:1261549/1‑147 (MQ=255)                                                                                                   cgcgTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTg                                               <  1:983631/151‑1 (MQ=255)                                                                                                       tGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCCCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAAACGTCTTTTTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTggcgg                                           >  1:483217/1‑151 (MQ=255)                                                                                                               aaaTGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTg                                   <  1:1234081/151‑1 (MQ=255)                                                                                                                 aTGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGgc                                 <  1:590940/151‑1 (MQ=255)                                                                                                                   gCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTTCGCTGACGCCCCTTACGCTGCCGTGGCTGGGGAATCATGCTTTTTGGCTGGCATTAACCGTCTTTTTGGCGGTGCGCGGGGTTTTTCTGGCGGCTATCTTGCGGCGTCACttg                                  >  1:540046/1‑146 (MQ=255)                                                                                                                     gTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCaa                             <  1:996489/151‑1 (MQ=255)                                                                                                                          aGTATGGCGGTGGCCGCCGCAGGGTTTTCGCCGACGCTCCTTACGCTGCCGTGGGTGGGGAATCATGCTTTGTGGCTGGCCTTAACCGTCTTTTTGGCGTTGGGCGGGCTTTCTCTGGCGGGTATTTCGGGGCGTCACCGgggcca                             >  2:927570/1‑144 (MQ=255)                                                                                                                           gTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTAc                       >  1:122950/1‑151 (MQ=255)                                                                                                                             aTGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCt                     <  1:896431/151‑1 (MQ=255)                                                                                                                                 cGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGtt                 <  2:731189/151‑1 (MQ=255)                                                                                                                                  ggTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGATTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGttt                <  2:522576/151‑1 (MQ=255)                                                                                                                                       ccgccgCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGCCGc           <  2:803738/151‑1 (MQ=255)                                                                                                                                        cgccgcAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGCCGc           >  1:1269171/1‑150 (MQ=255)                                                                                                                                          ccgcAGGTTTTGCGCTGACGCTCCTTTCGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGCCGCAAc        <  2:522001/151‑1 (MQ=255)                                                                                                                                           cgcAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGCCGCAACg       <  2:577428/151‑1 (MQ=255)                                                                                                                                               ggTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGCCGCAACGTGAc   >  2:1173788/1‑151 (MQ=255)                                                                                                                                                gTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGCCGCAACGTGACg  <  2:337442/151‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                      GATTCAGCAGCTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGCCGCAACGTGACG  >  NC_000913/4258830‑4259123

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