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

CCAGATTCAGCAGCTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGAC‑GCTCCTTACGCTGCCGTGGC‑TGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGCCGCAACGTGAC  >  NC_000913/4258827‑4259122                                                                                                                                                    |                                                                                                                                                       ccAGATTCAGCAGCTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGgttt                                                                                                                                                     >  2:970678/1‑151 (MQ=255)       tCAGCAGCTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAGCGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCt                                                                                                                                                <  1:95002/150‑1 (MQ=255)            agcTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGGGGCCGCCGCCGGGTTTTCGCCGAc                                                                                                                                             >  1:1095024/1‑148 (MQ=255)                ggctgACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGAC‑GCTCCt                                                                                                                                      <  2:13482/150‑1 (MQ=255)                         cTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGGGGGCCCCGCAGGGTTTGGGCTGAC‑GCTCCTTACgccgcc                                                                                                                             >  1:208266/1‑150 (MQ=255)                             cttatctGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCCCCGCCGGGTTTGCGCTGAC‑GCTCCTTTCGCTTCCGTg                                                                                                                          >  2:407519/1‑149 (MQ=255)                               tatctGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGTCGCAGGTTTTGCGCTGAC‑GCTCCTTACGCTGCCGTgg                                                                                                                         <  1:86944/148‑1 (MQ=255)                               tatctGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGAC‑GCTCCTTACGCTGCCGTg                                                                                                                          <  2:684298/147‑1 (MQ=255)                                 tctGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCCGGTTTTTCGCTGGC‑GCTCCCTTCGCTGCCGTGGC‑Tggg                                                                                                                   >  2:684873/1‑151 (MQ=255)                                   tGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCCCAGGGTTTTCGCCGAC‑GCTCCTTTCGCCGCCGCGGC‑Tggg                                                                                                                   >  1:955111/1‑149 (MQ=255)                                    ggCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGAC‑GCTCCTTACGCTGCCGTGGC‑TGGGTa                                                                                                                 <  2:227402/150‑1 (MQ=255)                                         gTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGGCGCCGCCGGGTTTTTGCTGGA‑GCTCCCTTCGCTGCCGGGGG‑TGGGGaa                                                                                                                >  1:933742/1‑146 (MQ=255)                                           gATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGAC‑GCTCCTTACGCTGCCGTGGC‑TGGGTAATCATGCt                                                                                                         <  1:646960/151‑1 (MQ=255)                                             ttttGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGCTGCCGCCGCAGGGTTTTGGCCGGC‑GCGCCCTTCCCTGCCCTGGCGTGGGCAAtatgcctt                                                                                                        >  1:553249/1‑144 (MQ=255)                                              tttGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGAC‑GCTCCTTACGCTGCCGTGGC‑TGGGTAATCATGCttt                                                                                                       <  2:549352/150‑1 (MQ=255)                                              tttGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGAC‑GCTCCTTACGCTGCCGTGGC‑TGGGTAATCATGCttt                                                                                                       <  2:968629/150‑1 (MQ=255)                                                  ccGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCCGGGTTTGCGCTGAC‑GCCCCTTACGCCGCCGTGGC‑TGGGTAATCATGCTTTGtgg                                                                                                   >  2:308486/1‑150 (MQ=255)                                                       ggtTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGTCGCCGCAGGTTTTGCGCTGAC‑GCTCCTTACGCTGCCGTGGC‑TGGGTAATCATGCTTTGtggctggc                                                                                              <  2:280722/150‑1 (MQ=255)                                                        gtTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCCCCGCCGGGTTTGCGCTGCC‑GCCCCCTTCGCTGCCGTGGC‑TGGGGAATCATGCttttt                                                                                                     >  1:851679/1‑140 (MQ=255)                                                              gTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCCGGGTTTGCGCTGAC‑GCCCCTTACGCTGCCGTGGC‑TGGGTAATCATGCTTTTTGGCTGGCAtttaac                                                                                       >  2:1183450/1‑146 (MQ=255)                                                                  ggTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGAC‑GCTCCTTACGCTGCCGTGGC‑TGGGTAATCATGCTTTGTGGCTGGCATTAACCGTc                                                                                    <  2:213500/149‑1 (MQ=255)                                                                  ggTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGAC‑GCTCCTTACGCTGCCGTGGC‑TGGGTAATCATGCTTTGTGGCTGGCATTAACCGTc                                                                                    <  2:383220/149‑1 (MQ=255)                                                                  ggTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGAC‑GCTCCTTACGCTGCCGTGGC‑TGGGTAATCATGCTTTGTGGCTGGCATTAACCGTc                                                                                    <  2:781038/149‑1 (MQ=255)                                                                  ggTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTTCGCTGAC‑GCCCCTTACGCTGCCGGGGC‑TGGGGAATCATGCTTTGTGGCCGGGATTAAACGTCtt                                                                                  >  2:1058669/1‑151 (MQ=255)                                                                          ctgctgGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGTCGCCGCCGGCTTTCCGCGGGC‑GCTCCCTTCTCTGCCGTGGC‑ATGGTAATCACTGTTTTTGGCTGGCCTTTACAGGTTTTCt                                                                               >  1:765831/1‑146 (MQ=255)                                                                            gctgGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGGCGCCGCAGGGGTTTCGCTGAC‑GCTTCTTACGCTCCCGGGGG‑TGGGGAATAATGCCTTGTGGCACGCCTTAAACGTCTTTCtcggg                                                                           >  2:72341/1‑144 (MQ=255)                                                                             ctgGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCCCAGGGTTTGCGCTGAC‑GCTCCTTACGCTGCCGTGGC‑TGGGTAATCATGCTTTGTGGGTGGCATTAAACGTCttttttgggtt                                                                         >  2:1007391/1‑141 (MQ=255)                                                                                 aCGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTTCGCTGGC‑GCTCCTTACGCTGCCGTGGG‑GGGGGAATTATGCTTTTTGGGTGGCCTTTACCCTTTTTTTGGGGTTTGGCgg                                                                   >  1:796099/1‑151 (MQ=255)                                                                                     cATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGAC‑GCTCCTTACGCTGCCGTGGC‑TGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCtt                                                               <  1:754592/151‑1 (MQ=255)                                                                                      aTGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGGCCCCGCAGGGTTTTCGCTGCC‑GCCCCCTTCGCTGCCGTGGC‑GGGGGTAACATTCTTTTTTGCCGGTGTTAAAACCCTTTCTTGCGTTTTGCggg                                                                  >  1:793460/1‑147 (MQ=255)                                                                                      aTGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTTTGCTGAC‑GCTCCTTTCGCGGCCGTGGG‑TGGGGAATAATGCTTTGTGGCTGGCATTAAACCGCTTTCTGGCGTTGCGCggg                                                                  >  2:292088/1‑147 (MQ=255)                                                                                      aTGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGAC‑GCTCCTTACGCTGCCGTGGC‑TGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCttt                                                              <  1:32318/151‑1 (MQ=255)                                                                                      aTGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGAC‑GCTCCTTACGCTGCCGTGGC‑TGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCttt                                                              <  1:25911/151‑1 (MQ=255)                                                                                      aTGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCCCAGGGTTTGCGCTGAC‑GCTCCTTACGCTTCCGTGGC‑TGGGGAATCATGCTTTGTGGCTGGCATTTACCCTCTTTTTTGCGTTGCGCGGGCttt                                                              >  2:979209/1‑151 (MQ=255)                                                                                        gTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGAC‑GCTCCTTACGCTGCCGTGGC‑TGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTtc                                                             <  1:1043405/150‑1 (MQ=255)                                                                                        gTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTGCGCTGAC‑GCTCCTTACGCTGCCGTGGC‑TGGGGAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCttt                                                              >  2:393073/1‑149 (MQ=255)                                                                                            ataGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGAC‑GCTCCTTACGCTGCCGTGGC‑TGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTg                                                         <  2:966751/150‑1 (MQ=255)                                                                                              aGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGAC‑GCTCCTTACGCTGCCGTGGC‑TGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTggc                                                       <  2:85369/150‑1 (MQ=255)                                                                                                gcgcAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTGCGCTGAC‑GCTCCCTACGCTGCCGTGGC‑TGGGTAATCATGCTTTTTGGCTGGCCTTAACCGTCTTTTTGGCGTTGCGCGGGCTTTCTCTggcg                                                      >  2:306640/1‑149 (MQ=255)                                                                                                gcgcAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGACCGCCGCAGGCTTTGTGCTGAC‑GCTCCTTACGGTGCCGGAGG‑CGGGGAATACTTGTTTGTGGGCGGGCTTAACCGTCTTTCCGGCCGTGGGCggggatttcctt                                                         >  1:939925/1‑137 (MQ=255)                                                                                                  gcAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCGGCCGCCGGGTTTTCGCTGAC‑GCCCCCTACGCTGCCGGGGC‑TGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTTGCGGTGCGCGGGCCTTCTCTggcgg                                                     >  1:376645/1‑148 (MQ=255)                                                                                                  gcAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCCGGGTTTGCGCTGAC‑GCTCCTTACGCTGCCGGGGC‑TGGGGAATCATGCTTTGTGGGTGGCATTAACCGTCTTTCTGGCGGTGCGCGGGGTTTCTCTggcgg                                                     >  2:451625/1‑148 (MQ=255)                                                                                                   cAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGAC‑GCTCCTTACGCTGCCGTGGC‑TGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTAt                                                 <  2:1145832/151‑1 (MQ=255)                                                                                                   cAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCCCCGGGTTTTTGCTGAC‑GCGCCCTACCCCGCCGTGGG‑TGGGGAAACATGCCTTTTGGGTGGCGTTAACCCTCTTTTTTGGGTTGCGCGGGGTTTTTCTGGCGGGTAt                                                 >  1:1008672/1‑151 (MQ=255)                                                                                                     aCGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGAC‑GCTCCTTACGCTGCCGTGGC‑TGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCt                                               <  1:828371/151‑1 (MQ=255)                                                                                                         gTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTGCGCTGAC‑GCTCCTTACGCTGCCGTGGC‑TGGGGAATCATGCTTTTTGGCTGGCATTTAACGTCTTTTTTGCGTTGCGCGGGGTTTCTCTGGCGGCTATCTgg                                             >  2:728836/1‑149 (MQ=255)                                                                                                          tGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCCGGGTTTTCGCTGAC‑GCTCCTTACGCTGCCGGGGG‑TGGGGAATCATGGTTTGTGGCTGGGATTAACCGGCTTTTTTGCGGTGCGCGGGCTTTTTCTGGGGGGTTTCTtgcgg                                          >  2:619565/1‑151 (MQ=255)                                                                                                          tGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTGCGCTGAC‑GCCCCTTACCCTGCCGGGGC‑TGGGGAATCAAGCCTTTTGGGTGGGGCTTAACGTCTTTTTGGCGGTTTCCGGGGGTTCTCTggcggccg                                                  >  1:1067761/1‑141 (MQ=255)                                                                                                           gCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGAC‑GCTCCTTACGCTGCCGTGGC‑TGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTggcggc                                         <  1:215610/151‑1 (MQ=255)                                                                                                              aCCGAAATGCGTAACAGTATGGCGGCGGGCGCCGCACGGTTTTCGCTGGCAGCCCCTTTC‑CTGCCGTGGC‑GGGGTAATCAAGCTTTTTGGCTGGGATTTAACGTCTTTTTTGGGTTGCGGGGGGTTTCTTTGGAGGCTAttttgtggggt                                       >  2:79027/1‑140 (MQ=255)                                                                                                                 gAAATGCGTAACAGTATGGCGGTGGCCGCCCCAGGGTTTTCGCTGAC‑GCCCCTTACGCCGCCGGGGC‑TGGGGAATCATGCTTTGTGGCTGGGCTTAACCCTCTTTTTTGGGGTGGGCGGGGTTTCTCTGGGGGGCTTCTggcgggggccc                                    >  1:89859/1‑144 (MQ=255)                                                                                                                  aaaTGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGAC‑GCTCCTTACGCTGCCGTGGC‑TGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTg                                  <  2:315098/151‑1 (MQ=255)                                                                                                                          aaCAGTATGGCGGTGGCCGCCGCAGGGTTTGCGCTGAC‑GCTCCTTACGCTGCCGTGGG‑TGGGGAATCATGCTTTGTGGCTGGCATTAACCGTCTTTTTTGCGTTGCGCGGGGTTTTTCTGGCGGCTATTTGGCGGCGGCACTGGCGCAATg                          >  1:1153277/1‑151 (MQ=255)                                                                                                                                aTGGCGGTGGGCGCCGCAGGTTTTGCGCTGAC‑GCTCCTTACGCTGCCGTGGC‑TGGGTAATCATGCTTTGTGGCTGGCATTAAACGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTAcc                     >  1:488640/1‑150 (MQ=255)                                                                                                                                aTGGCGGTGGCCGCCGCAGGTTTTGCGCTGAC‑GCTCCTTACGCTGCCGTGGC‑TGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCt                    <  2:77557/151‑1 (MQ=255)                                                                                                                                  ggCGGTGGCCGCCGCAGGTTTTGCGCTGAC‑GCTCCTTACGCTGCCGTGGC‑TGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTg                   <  2:369817/150‑1 (MQ=255)                                                                                                                                    cGGTGGCCGCCGCAGGTTTTGCGCTGAC‑GCTCCTTACGCTGCCGTGGC‑TGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGtt                <  2:1187290/151‑1 (MQ=255)                                                                                                                                     ggTGGCCGCCGCAGGTTTTGCGCTGAC‑GCTCCTTACGCTGCCGTGGC‑TGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGttt               >  1:134668/1‑151 (MQ=255)                                                                                                                                        ggCCGCCGCAGGTTTTGCGCTGAC‑GCTCCTTACGCTGCCGTGGC‑TGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGcc            >  2:939092/1‑151 (MQ=255)                                                                                                                                           cgccgcAGGTTTTGCGCTGAC‑GCTCCTTACGCTGCCGTGGC‑TGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGCCGCa         <  2:249501/151‑1 (MQ=255)                                                                                                                                               gcAGGTTTTGCGCTGAC‑GCTCCTTACGCTGCCGTGGC‑TGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGCCGCAACGt     <  1:336928/151‑1 (MQ=255)                                                                                                                                               gcAGGTTTTGCGCTGAC‑GCTCCTTACGCTGCCGTGGC‑TGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGCCGCAACGt     <  1:717470/151‑1 (MQ=255)                                                                                                                                                  ggTTTTGCGCTGAC‑GCTCCTTACGCTGCCGTGGC‑TGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGCCGCAACGTGAc  <  1:7466/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                       CCAGATTCAGCAGCTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGAC‑GCTCCTTACGCTGCCGTGGC‑TGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGCCGCAACGTGAC  >  NC_000913/4258827‑4259122

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: ‑