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

ACCCAGATTCAGCAGCTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTG‑CGCTGACGCT‑CCTTACG‑CTGCCGTGG‑CTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGCCGCAACGTG  >  NC_000913/4258825‑4259120                                                                                                                                                       |                                                                                                                                                       aCCCAGATTCAGCAGCTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGtt                                                                                                                                                       <  1:1156214/151‑1 (MQ=255)  cccAGATTCAGCAGCTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGttt                                                                                                                                                      <  2:1089908/151‑1 (MQ=255)     aGATTCAGCAGCTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGGCGCCGCAGGtttt                                                                                                                                                     >  1:444127/1‑149 (MQ=255)        ttCAGCAGCTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGGTTTt‑cg                                                                                                                                                 >  1:757427/1‑149 (MQ=255)         tCAGCAGCTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGgttt                                                                                                                                                     >  1:424461/1‑145 (MQ=255)          cagcagCTGGCTGACCGCTATCTTATCTGGCGGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTG‑CGCTGa                                                                                                                                             <  1:145005/151‑1 (MQ=255)          cagcagCTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTG‑CGCTGa                                                                                                                                             <  2:1465871/151‑1 (MQ=255)           agcagcTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTG‑CGCTGAc                                                                                                                                            <  1:1135587/151‑1 (MQ=255)            gcagcTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGGCCCCGCAGGGTTTT‑CGCTggc                                                                                                                                            >  1:127290/1‑148 (MQ=255)             cagcTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTG‑CGCTGACGc                                                                                                                                          <  2:1199485/151‑1 (MQ=255)                 tggctgACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGGCGCCGCCGGGTTTT‑CGCCGACGcc‑cc                                                                                                                                      >  2:588786/1‑150 (MQ=255)                 tggctgACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGGTTTG‑CGCTGACGCT‑cc                                                                                                                                      >  2:995775/1‑150 (MQ=255)                  ggctgACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTG‑CGCTGACGCT‑CCt                                                                                                                                     <  1:1177740/150‑1 (MQ=255)                    ctgACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGGTTTG‑CGCTGGCGCT‑CCtt                                                                                                                                    >  1:348796/1‑149 (MQ=255)                           cTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCCGGGTTTT‑CGCTGGCGCT‑CCCTTCg‑ccgcc                                                                                                                           >  2:488738/1‑150 (MQ=255)                            tatcttatctGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCCGGGTTTT‑CGCTGACGCC‑CCCTACG‑CTGcc                                                                                                                           >  2:641670/1‑149 (MQ=255)                             atcttatctGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTG‑CGCTGACGCT‑CCTTACG‑CTGCCGt                                                                                                                         <  2:568722/150‑1 (MQ=255)                               cttatctGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCCCCGGGTTTT‑CGCTGACGCT‑CCTTTCG‑CTTCCGCgg                                                                                                                       >  1:1356181/1‑150 (MQ=255)                                 tatctGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTG‑CGCTGACGCT‑CCTTACG‑CTGCCGTg                                                                                                                        <  1:993086/147‑1 (MQ=255)                                  atctGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGGTTTT‑CGCTGACGCT‑CCTTACG‑CTGCCGTGG‑CTgg                                                                                                                  >  1:470729/1‑151 (MQ=255)                                     tGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGGCGCCGCCGGGTTTT‑TGCTGACCCG‑CCTTACC‑CTGCCGTGG‑CTgg                                                                                                                  >  1:216130/1‑148 (MQ=255)                                      ggCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTG‑CGCTGACGCT‑CCTTACG‑CTGCCGTGG‑CTGGGt                                                                                                                <  2:757349/149‑1 (MQ=255)                                        cAGGTGATTTTGCCGTTGGTTGGCGTCTGGTTTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTG‑CGCTGACGCT‑CCTTACG‑CTGCCGTGG‑CTGGGTAATc                                                                                                            <  2:797940/151‑1 (MQ=255)                                          ggTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGGTTTT‑CGCTGACGCT‑CCTTACG‑CCGCCGTGG‑CTGGGGAATc                                                                                                            >  1:1061312/1‑149 (MQ=255)                                             gATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCCGGGTTTG‑GGCTGACCCT‑CCTTACG‑CTGCCGTGG‑CTGGGTaa                                                                                                              >  1:909194/1‑144 (MQ=255)                                                 ttGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCCGGGTTTG‑CGCTGACGCT‑CCTTACG‑CTGCCGGGG‑CTGGGGAATCATGCTTtgt                                                                                                   >  1:820948/1‑151 (MQ=255)                                                     cGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTG‑CGCTGACGCT‑CCTTACG‑CTGCCGTGG‑CTGGGTAATCATGCTTTGtggct                                                                                               <  2:246462/151‑1 (MQ=255)                                                        tggtTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCCGGGTTTT‑CGCTGACGCT‑CCTTACG‑CTGCCGGGG‑CCGGGGAATCATGCttt                                                                                                     >  1:1404061/1‑142 (MQ=255)                                                         ggtTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTG‑CGCTGACGCT‑CCTTACG‑CTGCCGTGG‑CTGGGTAATCATGCTTTGTGGCTGGCa                                                                                           <  1:119964/151‑1 (MQ=255)                                                             ggCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGGCGCCGCAGGGTTTT‑CGCTGACGCT‑CCTTACG‑CCGCCGGGG‑GCGGGGAATTAATCTTTGTGGGTGGGAtt                                                                                         >  2:514191/1‑149 (MQ=255)                                                             ggCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTG‑CGCTGACGCT‑CCTTACG‑CTGCCGTGG‑CTGGGTAATCATGCTTTGTGGCTGGCATTa                                                                                        <  1:1014668/150‑1 (MQ=255)                                                                gTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGGCGCCCCCCGCTTTGCCGCCGA‑GCTCCCTTAC‑‑CTGCCGTGG‑CCGGGGAATCATGCtttttagcgggggctaa                                                                                       >  2:1256306/1‑134 (MQ=255)                                                                gTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCCCCGGGTTTT‑CGCTGGCGCC‑CCTTTCG‑CGGCCGTGG‑CTGGGGGATTATTCTTTTTGGCTGGCGTTAAcc                                                                                     >  1:327758/1‑150 (MQ=255)                                                                   tGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCCTACCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGGTTTT‑CGCTGGCGCT‑CCTTTTT‑CTGCCGTGG‑CTGGGGAATCATGCTTTGTGGGTGGGAtttaacggct                                                                                 >  1:891505/1‑143 (MQ=255)                                                                    ggTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTG‑CGCTGACGCT‑CCTTACG‑CTGCCGTGG‑CTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCtt                                                                                <  1:778372/151‑1 (MQ=255)                                                                     gtgtTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTG‑CGCTGACGCT‑CCTTACG‑CTGCCGTGG‑CTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCttt                                                                               <  1:995870/151‑1 (MQ=255)                                                                        ttATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGGCGCCCCCGGTGTTG‑AGCTGACGCT‑CCTTACC‑ACGCCGTGG‑CGGGGTAAATAAAATTTGTGGGGGGCAATTACCGTCCTTCt                                                                             >  2:783536/1‑150 (MQ=255)                                                                         tATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGGCGCCGCAGGGTTTG‑CGCTGGCGCT‑CCTTACG‑CTGCCGTGG‑GTGGGGAATCATGCTTTGTGGGTGGGGTTTACCGTCttttt                                                                             >  2:1198282/1‑147 (MQ=255)                                                                          aTCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGCGGCCCCCGCCGGTTTTG‑CGCTGACGCT‑CCTTACG‑CCGCCCGGG‑CCGGGTAAACCTGCTTTGTGGCTGGCATTAAACGTCTTTCt                                                                             >  1:103329/1‑148 (MQ=255)                                                                           tCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCCGGGTTTT‑CGCTGACGCT‑CCTTACC‑CTGCCCGGG‑CTGGGGAATCGTGCCTTTTGGGCGGCATTAAACGTCTTTTTgg                                                                           >  1:285748/1‑149 (MQ=255)                                                                              gctgGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTG‑CGCTGACGCT‑CCTTACG‑CTGCCGTGG‑CTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGtt                                                                       <  2:991262/150‑1 (MQ=255)                                                                                tgGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCCCCCCAGGGTTTT‑CGCTGACCCT‑CCCTACG‑CTGCCGTGG‑CTGGGGCATCATGTTTTTGGGGTGGAATTAACCGTCTTTCTGGCg                                                                         >  1:1268031/1‑146 (MQ=255)                                                                                   aCGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTG‑CGCTGACGCT‑CCTTACG‑CTGCCGTGG‑CTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGcgcg                                                                  <  1:1461847/150‑1 (MQ=255)                                                                                       cATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGGCGCCGCAGGGTTTT‑CGCTGACGCC‑CCTTACG‑CTGCCGTGG‑CTGGGGAAACATGCTTTTTGGCTGGCATTAAACGTCTTTTTTGCGTTGCGCggg                                                                >  2:668170/1‑148 (MQ=255)                                                                                       cATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTG‑CGCTGACGCT‑CCTTACG‑CTGCCGTGG‑CTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCtt                                                             <  2:1303007/151‑1 (MQ=255)                                                                                        aTGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGGTTTT‑CGCTGACGCC‑CCTTACG‑CTGCCGTGG‑CTGGGGAATCATGCTTTTTGGCTGGCATTAAACGGCTTTTTGGCGTTGCGCGGGCgtt                                                            >  1:889318/1‑151 (MQ=255)                                                                                        aTGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGAGGCCCCCCCAGGTTTTT‑AGCCGAGACA‑CCCTACA‑CTGCCGCGG‑GTGGGGGAATATGTTTTTTGGGTCGGATTTAACGTCTTTTTTGTGCTGGGCggggctt                                                            >  2:1212319/1‑147 (MQ=255)                                                                                          gTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGGTTTG‑CGCTGAAGCT‑CCTTACG‑CTGCCGTGG‑CTGGGGAATCATGCTTTGTGGCTGGCATTAAACGTCTTTCTTGCGGTGCGCGGGGttt                                                            >  1:275753/1‑149 (MQ=255)                                                                                          gTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCCCCGGGTTTT‑CGGTGTCGCC‑CCCTACC‑CTGCCGCGG‑CTGGGGGATCATTCTTTGTGGCTGGGGTTAACCCTCTTTCTTGCGTTGCGCGGGCgtt                                                            >  2:1111750/1‑149 (MQ=255)                                                                                             tataGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGGCGCCCCCGGTTTTG‑CGGTGTCCCC‑CCTTTCG‑CTGCCGAGG‑CTGGGGGATTATGCTTTTTGGGGGGCATTAACAGCATTTCTGGCGTTGGGCCGGCCTTttctt                                                       >  2:1222345/1‑150 (MQ=255)                                                                                              ataGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTG‑CGCTGACGCT‑CCTTACG‑CTGCCGTGG‑CTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTg                                                       <  2:110166/150‑1 (MQ=255)                                                                                               taGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCCCCCCCGGGTTTG‑GGCTGACGGG‑CCTTACG‑CTGCCGTGG‑GTGGGGAAACATGCTTTTTGGGTGGCATTAACCGTCTTTCTCGCGTTTCGCGGGCGTTttct                                                        >  1:1069633/1‑148 (MQ=255)                                                                                                aGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTG‑CGCTGACGCT‑CCTTACG‑CTGCCGTGG‑CTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTggc                                                     <  1:164977/150‑1 (MQ=255)                                                                                                    gcAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGGTTTG‑CGCCGACGCT‑CCCTACG‑CTGCCGTGG‑CTGGGGAATCATGCTTTGTGGCTGGCCTTTACCGTCTTTCTTGCGGTGCGCGGGCTTTTTCTtgcg                                                    >  2:198059/1‑147 (MQ=255)                                                                                                     cAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTG‑CGCTGACGCT‑CCTTACG‑CTGCCGTGG‑CTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTCTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTAt                                               <  2:1061216/151‑1 (MQ=255)                                                                                                     cAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGGTTTT‑CGCTGACGCT‑CCTTACG‑CTGCCGTGG‑CTGGGGAATCATGCTTTGTGGCTGGCATTTAACCTCTTTCTTGCGTTGCGCGGGCTTTTTCTtgcgg                                                   >  1:625402/1‑147 (MQ=255)                                                                                                     cAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCCCCGCAGGGTTTT‑CGCCGCCGCT‑CCTTACG‑CCGCCGTGG‑CTGGGGAATCATGCTTTGTGGCTGGCCTTAAACGCCTTTCTGGCGGTGCGCGGGGTTTTTCTtgcgg                                                   >  2:1292692/1‑147 (MQ=255)                                                                                                      aaCGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCCGGGTTTG‑CGCTGACGCT‑CCTTACG‑CTGCCGTGG‑CTGGGTAATCATGCTTTTTGGCTGGCATTAAACGTCTTTCTGGCGGTGCGCGGGCTTTTTCTggcgg                                                   >  1:151045/1‑146 (MQ=255)                                                                                                        cgcgTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTG‑CGCTGACGCT‑CCTTACG‑CTGCCGTGG‑CTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTg                                            <  1:147644/151‑1 (MQ=255)                                                                                                            tGCCACCGAAATGCGTAACAGTATGGCGGGGGGCACCCGCGGGTTTT‑‑AGCTGACCCT‑CCCTAAA‑CTCCCGCGGCCGGGGGAATAAT‑ATTTGTGGGTGGGAATCACCGTCCTTTTTGCGTTGCGCCGGCTTTTTCTCGGGGCt                                                 >  2:1040827/1‑142 (MQ=255)                                                                                                             gCCACCGAAATGCGTAACAGTATGGC‑GGTGGGCGCCGCAGGGTTTT‑GGCTGACGCC‑CCTTACG‑CTGCCGGGG‑CTGGGGAAAAATGCTTTTTGGGTGGCATTTAACGGCTTTCTGGGGGTGCGCGGGGTTTTTCTGGCGGCTATCttgggg                                        >  2:312956/1‑145 (MQ=255)                                                                                                             gCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTG‑CGCTGACGCT‑CCTTACG‑CTGCCGTGG‑CTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTggcgg                                        <  2:30115/150‑1 (MQ=255)                                                                                                             gCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTG‑CGCTGACGCT‑CCTTACG‑CTGCCGTGG‑CTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTggcgg                                        <  1:233640/150‑1 (MQ=255)                                                                                                                  cGAAATGCGTAACAGTATGGC‑GGTGGCCGCCCCAAGGTTTG‑CGCTGGCCCC‑CCTTATGTCTCCCG‑GG‑CGGGGGAACATTCCTTTTTGGGCGGCAATAACAGCTTTTTTGGCGTTGCTTGGGGTTTTTTTGTGGGCTAttgttcggcggg                                    >  1:936280/1‑138 (MQ=255)                                                                                                                   gAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTT‑CGCTGACGCT‑CCTTACG‑CTGCCGTGG‑CTGGGGAATCATGCTTTGTGGCTGGCATTAACCGTCTTTTTGGCGTTGCGCGGGGTTTTTCTGGCGGCTATCTGGCGGCGTCAc                                  >  1:89241/1‑150 (MQ=255)                                                                                                                     aaTGCGTAACAGTATGGC‑GGTGGCCGCCGCCGGGTTTG‑CGCTGACGCT‑CCTTACG‑CTGCCGGGG‑GTGGGGGATCATGCTTTTTGGGTGGGATTTAACCTCTTTTTTGCGGTTCGCGGGGTTTTTCCTGCGGCTATTTTGCGGGGGCACCgg                               >  2:809283/1‑151 (MQ=255)                                                                                                                        gCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTG‑CGCTGACGCT‑CCTTACG‑CTGCCGTGG‑CTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGgcg                             <  1:992516/150‑1 (MQ=255)                                                                                                                          gTAACAGTATGGC‑GGTGGCCGCCGCAGGGTTTG‑CGCCGACGCC‑CCTTACG‑CTGCCGTGG‑CTGGGGAATCATGCTTTGTGGCGGGCCATAACCCCCTTTCTGGCCGTGCGCGGGCCTTCTCTCGCGGCTATTTGGCGGCGTccc                                  >  2:975232/1‑141 (MQ=255)                                                                                                                              cAGTATGGC‑GGTGGCCGCCGCCGGTTTTT‑CGCTGACGCT‑CCTTACG‑CTGCCGTGG‑CTGGGTAATCATGCTTTGTGGCTGGCATTAAACGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATgg                       >  2:301179/1‑150 (MQ=255)                                                                                                                                gTATGGC‑GGTGGCCGCCGCAGGTTTTG‑CGCTGACGCT‑CCTTACG‑CTGCCGTGG‑CTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTAc                    <  1:965742/151‑1 (MQ=255)                                                                                                                                gTATGGC‑GGTGGCCGCCGCAGGTTTTG‑CGCTGACGCT‑CCTTACG‑CTGCCGTGG‑CTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTAc                    <  1:34226/151‑1 (MQ=255)                                                                                                                                  aTGGC‑GGTGGCCGCCGCAGGTTTTG‑CGCTGACGCT‑CCTTACG‑CTGCCGTGG‑CTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCt                  <  1:460912/151‑1 (MQ=255)                                                                                                                                    ggC‑GGTGGCCGCCGCAGGTTTTG‑CGCTGACGCT‑CCTTACG‑CTGCCGTGG‑CTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTgg                <  2:444074/151‑1 (MQ=255)                                                                                                                                        ggTGGCCGCCGCAGGTTTTG‑CGCTGACGCT‑CCTTACG‑CTGCCGTGG‑CTGGGTAATCATTCTTTGTGGCTGGCATTACCCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGttt             <  2:865240/151‑1 (MQ=255)                                                                                                                                             ccgccgCAGGTTTTG‑CGCTGACGCT‑CCTTACG‑CTGCCGTGG‑CTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGCCGc        <  2:1069532/151‑1 (MQ=255)                                                                                                                                               gccgcAGGTTTTG‑CGCTGACGCT‑CCTTACG‑CTGCCGTGG‑CTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGCCGCa       >  1:538063/1‑150 (MQ=255)                                                                                                                                                ccgcAGGTTTTG‑CGCTGACGCT‑CCTTACG‑CTGCCGTGG‑CTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGCCGCAAc     >  1:867448/1‑151 (MQ=255)                                                                                                                                                ccgcAGGTTTTG‑CGCTGACGCT‑CCTTACG‑CTGCCGTGG‑CTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGCCGCAAc     >  2:348460/1‑151 (MQ=255)                                                                                                                                                   cAGGTTTTG‑CGCTGACGCT‑CCTTACG‑CTGCCGTTG‑CTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGGGTTTCGCGGGCTTTCTCTGTCGTTTATCTGGCGGCGTCAATTGCGCAATGGTACCTGGTTTGCCGCAACGTg  <  1:1369765/151‑1 (MQ=255)                                                                                                                                                   cAGGTTTTG‑CGCTGACGCT‑CCTTACG‑CTGCCGTGG‑CTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGCCGCAACGTg  >  2:1234961/1‑151 (MQ=255)                                                                                                                                                       |                                                                                                                                                       ACCCAGATTCAGCAGCTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTG‑CGCTGACGCT‑CCTTACG‑CTGCCGTGG‑CTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGCCGCAACGTG  >  NC_000913/4258825‑4259120

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