Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,258,9740TG40.5% 146.4 / 17.2 75G411G (GGT→GGG) dinFoxidative stress resistance protein; putative MATE family efflux pump; UV and mitomycin C inducible protein
Reads supporting (aligned to +/- strand):  ref base T (4/40);  new base G (30/0);  total (35/40)
Fisher's exact test for biased strand distribution p-value = 9.89e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACCCAGATTCAGCAGCTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGC‑T‑GGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGCCGCAACGTG  >  NC_000913/4258824‑4259120
                                                                                                                                                       |                                                                                                                                                    
aaCCCAGATTCAGCAGCTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGt                                                                                                                                                      <  1:302810/151‑1 (MQ=255)
         tCAGCAGCTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGGGGCCGCCGCAGGGTTTg                                                                                                                                                  >  2:1189919/1‑146 (MQ=255)
          cagcagCTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGGTGGCC‑CCCCAGGGTTTgcg                                                                                                                                                >  2:969719/1‑147 (MQ=255)
          cagcagCTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCTg                                                                                                                                             <  1:846755/150‑1 (MQ=255)
          cagcagCTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGGTGGCGCCGCAGGGTTTgcg                                                                                                                                                >  2:1217334/1‑147 (MQ=255)
           agcagcTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCTGAc                                                                                                                                           <  1:656347/151‑1 (MQ=255)
              agcTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGGCGCCCCAGGGTTTGCGCTGAc                                                                                                                                           >  1:762382/1‑148 (MQ=255)
               gctggctgACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGGTTTTCGCTTACGCTc                                                                                                                                       >  1:382863/1‑151 (MQ=255)
                 tggctgACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCCCCCCAGGGTTTGTGGTGACGCTcc                                                                                                                                      >  2:165044/1‑150 (MQ=255)
                  ggctgACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCtt                                                                                                                                    <  1:1232382/151‑1 (MQ=255)
                   gctgACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTa                                                                                                                                   <  1:231431/151‑1 (MQ=255)
                    ctgACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCCCAGAGTGTGCGCGGTCGCTCCtttc                                                                                                                                  >  2:1196511/1‑149 (MQ=255)
                      gACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGc                                                                                                                                <  1:762770/151‑1 (MQ=255)
                       aCCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCt                                                                                                                               <  1:217942/151‑1 (MQ=255)
                                 tatctGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCCGGGTTTGCGGTGCAGCCCCTTACGCCTCCGTGGGTg                                                                                                                     >  2:1105151/1‑151 (MQ=255)
                                 tatctGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTgg                                                                                                                        <  2:716297/148‑1 (MQ=255)
                                 tatctGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTgg                                                                                                                        <  2:978048/148‑1 (MQ=255)
                                      ggCGGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTaa                                                                                                                <  2:435766/151‑1 (MQ=255)
                                      ggCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGt                                                                                                                  <  2:520231/149‑1 (MQ=255)
                                      ggCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTaa                                                                                                                <  1:134753/151‑1 (MQ=255)
                                        cAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCCGGGTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGGaa                                                                                                                >  1:418508/1‑149 (MQ=255)
                                         aGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCa                                                                                                             <  2:148426/151‑1 (MQ=255)
                                           gTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATg                                                                                                           <  1:109654/151‑1 (MQ=255)
                                              aTTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCtt                                                                                                        <  2:457830/151‑1 (MQ=255)
                                              aTTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGGTTTTCGCTGACGCTCCTTACGCCGCCGGGGCTGGGGAATCaatctt                                                                                                        >  1:1023171/1‑146 (MQ=255)
                                               ttttGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGGTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGGAATCATGCttt                                                                                                       >  1:827165/1‑151 (MQ=255)
                                                tttGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCttt                                                                                                       <  1:133204/150‑1 (MQ=255)
                                                   gCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTtgtg                                                                                                    <  1:848752/150‑1 (MQ=255)
                                                    ccGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCGGGGTTTTGGCTGACGCTCCCTACGCTCGCGTGGGTGGGGAAACATGCttt                                                                                                       >  2:807992/1‑146 (MQ=255)
                                                      ggtggtTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGtggc‑t‑g                                                                                              <  1:1194624/151‑1 (MQ=255)
                                                      ggtggtTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGtggc‑t‑g                                                                                              <  1:922542/151‑1 (MQ=255)
                                                       gtggtTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCGGGGTTTGCGCTGACGCTCCTTACGCTTCCGTGGGTGGGGAATCATGCTTTGtggc‑t‑gg                                                                                             >  1:644899/1‑151 (MQ=255)
                                                         ggtTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGtggc‑t‑ggc                                                                                            <  2:115987/150‑1 (MQ=255)
                                                         ggtTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGtggc‑t‑ggaa                                                                                           <  2:507386/151‑3 (MQ=255)
                                                          gtTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGGCGCCGCCGGGTTTGCGCCCACGGCCCTTACCCCGCCCTGGCTGGGGAATCATTATTTTtgggcg                                                                                                >  2:920413/1‑144 (MQ=255)
                                                            tGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCCGGGTTTGCGCTGACGCTCCTTACGCTGCCGGGGCTGGGTAATCATGCTTTGTGGC‑T‑GGGAtt                                                                                         >  2:801032/1‑150 (MQ=255)
                                                            tGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGC‑T‑GGCAtt                                                                                         <  1:195545/150‑1 (MQ=255)
                                                                    ggTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGC‑T‑GGCATTAACCGTCt                                                                                 <  1:653523/150‑1 (MQ=255)
                                                                    ggTGTTATCTGCTGGACGGCATGTTTATAGGAGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGC‑T‑GGCATTAACCGTCtt                                                                                <  2:814184/151‑1 (MQ=255)
                                                                     gtgtTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGC‑T‑GGCATTAACCGTCttt                                                                               <  2:864510/151‑1 (MQ=255)
                                                                      tgtTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGGCGCCCCAGGGTTTGCGCTGACGCCCCTTACGCTGTCGTGGGTGGGGAATCAAGCTTTTTGGG‑T‑GGGATTTAccccccttt                                                                              >  2:1003783/1‑144 (MQ=255)
                                                                         tATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGGTTTTCGCTGACGCTCCTTACGCTGCCGGGGCTGGGGAATCATGCTTTTTGGC‑T‑GGGCTTTACCGGCTTTCttg                                                                           >  1:403943/1‑149 (MQ=255)
                                                                             tgctgGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGGTTTTCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTTTGGC‑T‑GGCATTAACCGTCTTTCTgg                                                                           >  1:561541/1‑147 (MQ=255)
                                                                              gctgGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGC‑T‑GGCATTAACCGTCTTTCTGGCGtt                                                                       <  1:469259/150‑1 (MQ=255)
                                                                                tgGACGGCATGTTTATAGGCGCAACACGTGCCACCGAAATGCGTAACAGTATGGC‑GGTTGCCGACGCGGGGTTTACGCTGATTCTCCTTACGCCGGCTTTGCGGGGTGATCATACGTTTTTGC‑T‑GGCATTTAACATCttttttacgtt                                                                       >  2:284119/1‑140 (MQ=255)
                                                                                 ggACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGC‑T‑GGCATTAACCGTCTTTCTGGCGTTgcg                                                                    <  1:66023/150‑1 (MQ=255)
                                                                                    cGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGC‑T‑GGCATTAACCGTCTTTCTGGCGTTGCGCggg                                                                <  1:178080/151‑1 (MQ=255)
                                                                                    aggCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGGCGCCCCAGGGTTTGCGGCGGAGCTCCTTTCGCTGCCGGGGGCGGGGAAATATGCTTTTTTGG‑G‑GGGATTAACCGGCTTTCttggggttggcggg                                                                >  2:737795/2‑138 (MQ=255)
                                                                                       cATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGGCGCCGCCGGGTTTTTGCTGACGCTCCCTACGCTGCCGGGGCTGGGGAATTATGGTTTTTGGT‑T‑GGGATTTAACGTCTTTCTTGGGttttgggggggt                                                             >  2:399138/1‑141 (MQ=255)
                                                                                        aTGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCCCAGGGTTTGCGCTGTCGCTCCTTTCGCTGCCGCTGCTGGGTAATCATGCTTTTTGGC‑T‑GGCATTAAACGTCTTTTTTGCGTTTTGCGGGCttt                                                            >  1:214929/1‑151 (MQ=255)
                                                                                         tGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGGTTTGTGCTGGCGCCCCTTACGCTGCCCAGGGTGGGTAAACCATCTTTTTTGG‑G‑GGCATTTACCGGCCTTTTTGGGTTGCGCGGGCttt                                                            >  2:1237368/1‑150 (MQ=255)
                                                                                            ttATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGGCGCCGCAGGGTTTGCCGTGACGCTCCCTACGCTGCCGGGGGCGGGGAATCATGCTTTTTGGC‑TGGGCTTTAACGCTCTTT‑TTGCGTGTCGCggggatt                                                            >  2:767135/1‑143 (MQ=255)
                                                                                                aGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGC‑T‑GGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTggc                                                     <  2:1054446/150‑1 (MQ=255)
                                                                                                aGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCCCCGCTGGGTTTGCGCTGGCCCGCCCTACGCGGTCGTGGGTGGGGAAGCCTTGTTTTTGGCAT‑GGCCTT‑ACCCTCTTTTCGGCGttttccggggtttttc                                                         >  1:685075/1‑132 (MQ=255)
                                                                                                      aaCGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGC‑T‑GGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTAt                                               <  2:465316/150‑1 (MQ=255)
                                                                                                      aaCGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGC‑T‑GGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATc                                              <  2:570798/151‑1 (MQ=255)
                                                                                                       aCGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGC‑T‑GGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCt                                             <  1:485469/151‑1 (MQ=255)
                                                                                                        cgcgTGCCACCGAAATGCGTAACAGTATGGC‑GGAGGCCGCCGCAGGGGTTTCGCCGACACCCCTTACACTGCCGTGGGGGGGGGATCAGGCCTTGTGGC‑T‑GGCATTTACAGCCCTTTTGGCGGTGTGCGGGGTTTCTCTggcgggtg                                                >  2:740322/1‑144 (MQ=255)
                                                                                                            tGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCCCCGGGTTTGCCCTGACGCTCCTTTCGCCGCCCTGGCTGGGGAAACATGCTTTTTGGG‑T‑GGGATTTACCCTCTTTCTTGCGGTGCGCGGGGCTTTTCTGGCGGCGATTTtgcgg                                        >  2:274935/1‑151 (MQ=255)
                                                                                                             gCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGC‑T‑GGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTggcgg                                        <  2:119910/150‑1 (MQ=255)
                                                                                                             gCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGC‑T‑GGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTggcgg                                        <  1:200092/150‑1 (MQ=255)
                                                                                                             gCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCCGACGCCCCTTACGCCGCCGTTGCCCGGCAATCATGCCTTTTGGC‑T‑GGTATTTAACGTCTTTCTGGCGGTTCCTGGGCCTTCTTTGGCGGGTATCTagcgg                                        >  1:603114/1‑150 (MQ=255)
                                                                                                              ccACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGC‑T‑GGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGgcggcg                                      <  2:1109940/151‑1 (MQ=255)
                                                                                                               cACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGC‑T‑GGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTggcgg                                        <  2:704797/148‑1 (MQ=255)
                                                                                                                 ccGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCCGGGTTTTCGCTGACGCTCCTTACGCTGCCGTGGCTGGGGAATCATGCTTTTTGGG‑T‑GGCATTAACCGTCTTTCTTGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTTgcggcg                                      >  2:730825/1‑148 (MQ=255)
                                                                                                                    aaaTGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGC‑T‑GGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTg                                <  1:515348/151‑1 (MQ=255)
                                                                                                                     aaTGCGTAACAGTATGGC‑GGTGGGCCCCGCAGGGTTTTCGCTGACGCTCCTTACGCTGCCGTGGGTGGGGAAACATGGTTTTTGGC‑T‑GGGATTAAACGTCTTTTTGGCGGTGCGCGGGGTTTTTCTGGCGGCTCTCTTGCGGCGTCCCttg                               >  2:582714/1‑149 (MQ=255)
                                                                                                                       tGCGTAACAGTATGGC‑GGTGGCCGCCCCCGGTTTTGCGCCGGCGCTCCTTTCCCTGCCGGGGGTGGGGTATCATGCTTTGTTGC‑G‑GGCATTTACCCGCTTTCTGGCGTTGCGCGGGGCTTCTTTTGCGGCTTTTTtgcgggggcc                                   >  2:814591/1‑140 (MQ=255)
                                                                                                                        gCGTAACAGTATGGC‑GGTGGGCGCCCCCGGGTTTTCGCCTGCGCCCCCTACGCCGCCGTGGGTGGGGAATCATGCTTTTTGGC‑C‑GGCCTTAAACGCCTTTTTGGCGGTGCGCGGGGCTTTTCTGGCGGCTATCTtgcgg                                        >  2:153970/1‑139 (MQ=255)
                                                                                                                        gCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGC‑T‑GGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGgcg                             <  1:807999/150‑1 (MQ=255)
                                                                                                                            aaCAGTATGGC‑GGTGGCCGCCGCAGGCTTTGTGCTGACGCCCCCTACGCAGCCGGGGGCGGGGAATCATGCTTTGTGGC‑T‑GGCATTAACCCTCTTTTTGGCGGTGGGCGGGCCTTTTCTGGTGGGTTTCTTGCGGGGGCACTg                                >  2:160684/1‑143 (MQ=255)
                                                                                                                                        ggTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGC‑T‑GGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCCATGGTACCTGGttt             >  2:63050/1‑151 (MQ=255)
                                                                                                                                        ggTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGC‑T‑GGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGttt             <  1:525476/151‑1 (MQ=255)
                                                                                                                                            gccgccGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGC‑T‑GGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGCCg         <  1:935219/151‑1 (MQ=255)
                                                                                                                                             ccgccgCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGC‑T‑GGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGCCGc        <  1:1217319/151‑1 (MQ=255)
                                                                                                                                                    aGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGC‑T‑GGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGCCGCAACGTg  >  1:323390/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                    
AACCCAGATTCAGCAGCTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGC‑T‑GGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGCCGCAACGTG  >  NC_000913/4258824‑4259120

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