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→GGGdinFoxidative 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

AACCCAGATTCAGCAGCTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGCCGCAACGTG  >  NC_000913/4258824‑4259120
                                                                                                                                                       |                                                                                                                                                    
aaCCCAGATTCAGCAGCTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGt                                                                                                                                                      <  1:302810/151‑1 (MQ=255)
         tCAGCAGCTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGGGGCCGCCGCAGGGTTTg                                                                                                                                                  >  2:1189919/1‑146 (MQ=255)
          cagcagCTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGGTGGCCCCCCAGGGTTTgcg                                                                                                                                                >  2:969719/1‑147 (MQ=255)
          cagcagCTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTg                                                                                                                                             <  1:846755/150‑1 (MQ=255)
          cagcagCTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGGTGGCGCCGCAGGGTTTgcg                                                                                                                                                >  2:1217334/1‑147 (MQ=255)
           agcagcTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGAc                                                                                                                                           <  1:656347/151‑1 (MQ=255)
              agcTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGGCGCCCCAGGGTTTGCGCTGAc                                                                                                                                           >  1:762382/1‑148 (MQ=255)
               gctggctgACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTTCGCTTACGCTc                                                                                                                                       >  1:382863/1‑151 (MQ=255)
                 tggctgACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCCCCCCAGGGTTTGTGGTGACGCTcc                                                                                                                                      >  2:165044/1‑150 (MQ=255)
                  ggctgACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCtt                                                                                                                                    <  1:1232382/151‑1 (MQ=255)
                   gctgACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTa                                                                                                                                   <  1:231431/151‑1 (MQ=255)
                    ctgACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCCCAGAGTGTGCGCGGTCGCTCCtttc                                                                                                                                  >  2:1196511/1‑149 (MQ=255)
                      gACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGc                                                                                                                                <  1:762770/151‑1 (MQ=255)
                       aCCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCt                                                                                                                               <  1:217942/151‑1 (MQ=255)
                                 tatctGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCCGGGTTTGCGGTGCAGCCCCTTACGCCTCCGTGGGTg                                                                                                                     >  2:1105151/1‑151 (MQ=255)
                                 tatctGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTgg                                                                                                                        <  2:716297/148‑1 (MQ=255)
                                 tatctGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTgg                                                                                                                        <  2:978048/148‑1 (MQ=255)
                                      ggCGGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTaa                                                                                                                <  2:435766/151‑1 (MQ=255)
                                      ggCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGt                                                                                                                  <  2:520231/149‑1 (MQ=255)
                                      ggCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTaa                                                                                                                <  1:134753/151‑1 (MQ=255)
                                        cAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCCGGGTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGGaa                                                                                                                >  1:418508/1‑149 (MQ=255)
                                         aGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCa                                                                                                             <  2:148426/151‑1 (MQ=255)
                                           gTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATg                                                                                                           <  1:109654/151‑1 (MQ=255)
                                              aTTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCtt                                                                                                        <  2:457830/151‑1 (MQ=255)
                                              aTTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTTCGCTGACGCTCCTTACGCCGCCGGGGCTGGGGAATCaatctt                                                                                                        >  1:1023171/1‑146 (MQ=255)
                                               ttttGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGGAATCATGCttt                                                                                                       >  1:827165/1‑151 (MQ=255)
                                                tttGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCttt                                                                                                       <  1:133204/150‑1 (MQ=255)
                                                   gCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTtgtg                                                                                                    <  1:848752/150‑1 (MQ=255)
                                                    ccGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCGGGGTTTTGGCTGACGCTCCCTACGCTCGCGTGGGTGGGGAAACATGCttt                                                                                                       >  2:807992/1‑146 (MQ=255)
                                                      ggtggtTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGtggctg                                                                                              <  1:1194624/151‑1 (MQ=255)
                                                      ggtggtTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGtggctg                                                                                              <  1:922542/151‑1 (MQ=255)
                                                       gtggtTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCGGGGTTTGCGCTGACGCTCCTTACGCTTCCGTGGGTGGGGAATCATGCTTTGtggctgg                                                                                             >  1:644899/1‑151 (MQ=255)
                                                         ggtTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGtggctggc                                                                                            <  2:115987/150‑1 (MQ=255)
                                                         ggtTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGtggctggaa                                                                                           <  2:507386/151‑3 (MQ=255)
                                                          gtTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGGCGCCGCCGGGTTTGCGCCCACGGCCCTTACCCCGCCCTGGCTGGGGAATCATTATTTTtgggcg                                                                                                >  2:920413/1‑144 (MQ=255)
                                                            tGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCCGGGTTTGCGCTGACGCTCCTTACGCTGCCGGGGCTGGGTAATCATGCTTTGTGGCTGGGAtt                                                                                         >  2:801032/1‑150 (MQ=255)
                                                            tGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCAtt                                                                                         <  1:195545/150‑1 (MQ=255)
                                                                    ggTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCt                                                                                 <  1:653523/150‑1 (MQ=255)
                                                                    ggTGTTATCTGCTGGACGGCATGTTTATAGGAGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCtt                                                                                <  2:814184/151‑1 (MQ=255)
                                                                     gtgtTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCttt                                                                               <  2:864510/151‑1 (MQ=255)
                                                                      tgtTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGGCGCCCCAGGGTTTGCGCTGACGCCCCTTACGCTGTCGTGGGTGGGGAATCAAGCTTTTTGGGTGGGATTTAccccccttt                                                                              >  2:1003783/1‑144 (MQ=255)
                                                                         tATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTTCGCTGACGCTCCTTACGCTGCCGGGGCTGGGGAATCATGCTTTTTGGCTGGGCTTTACCGGCTTTCttg                                                                           >  1:403943/1‑149 (MQ=255)
                                                                             tgctgGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTTCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTTTGGCTGGCATTAACCGTCTTTCTgg                                                                           >  1:561541/1‑147 (MQ=255)
                                                                              gctgGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGtt                                                                       <  1:469259/150‑1 (MQ=255)
                                                                                tgGACGGCATGTTTATAGGCGCAACACGTGCCACCGAAATGCGTAACAGTATGGCGGTTGCCGACGCGGGGTTTACGCTGATTCTCCTTACGCCGGCTTTGCGGGGTGATCATACGTTTTTGCTGGCATTTAACATCttttttacgtt                                                                       >  2:284119/1‑140 (MQ=255)
                                                                                 ggACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTgcg                                                                    <  1:66023/150‑1 (MQ=255)
                                                                                    cGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCggg                                                                <  1:178080/151‑1 (MQ=255)
                                                                                    aggCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGGCGCCCCAGGGTTTGCGGCGGAGCTCCTTTCGCTGCCGGGGGCGGGGAAATATGCTTTTTTGGGGGGATTAACCGGCTTTCttggggttggcggg                                                                >  2:737795/2‑138 (MQ=255)
                                                                                       cATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGGCGCCGCCGGGTTTTTGCTGACGCTCCCTACGCTGCCGGGGCTGGGGAATTATGGTTTTTGGTTGGGATTTAACGTCTTTCTTGGGttttgggggggt                                                             >  2:399138/1‑141 (MQ=255)
                                                                                        aTGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCCCAGGGTTTGCGCTGTCGCTCCTTTCGCTGCCGCTGCTGGGTAATCATGCTTTTTGGCTGGCATTAAACGTCTTTTTTGCGTTTTGCGGGCttt                                                            >  1:214929/1‑151 (MQ=255)
                                                                                         tGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTGTGCTGGCGCCCCTTACGCTGCCCAGGGTGGGTAAACCATCTTTTTTGGGGGCATTTACCGGCCTTTTTGGGTTGCGCGGGCttt                                                            >  2:1237368/1‑150 (MQ=255)
                                                                                            ttATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGGCGCCGCAGGGTTTGCCGTGACGCTCCCTACGCTGCCGGGGGCGGGGAATCATGCTTTTTGGCTGGGCTTTAACGCTCTTTTTGCGTGTCGCggggatt                                                            >  2:767135/1‑143 (MQ=255)
                                                                                                aGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTggc                                                     <  2:1054446/150‑1 (MQ=255)
                                                                                                aGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCCCCGCTGGGTTTGCGCTGGCCCGCCCTACGCGGTCGTGGGTGGGGAAGCCTTGTTTTTGGCATGGCCTTACCCTCTTTTCGGCGttttccggggtttttc                                                         >  1:685075/1‑132 (MQ=255)
                                                                                                      aaCGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTAt                                               <  2:465316/150‑1 (MQ=255)
                                                                                                      aaCGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATc                                              <  2:570798/151‑1 (MQ=255)
                                                                                                       aCGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCt                                             <  1:485469/151‑1 (MQ=255)
                                                                                                        cgcgTGCCACCGAAATGCGTAACAGTATGGCGGAGGCCGCCGCAGGGGTTTCGCCGACACCCCTTACACTGCCGTGGGGGGGGGATCAGGCCTTGTGGCTGGCATTTACAGCCCTTTTGGCGGTGTGCGGGGTTTCTCTggcgggtg                                                >  2:740322/1‑144 (MQ=255)
                                                                                                            tGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCCCCGGGTTTGCCCTGACGCTCCTTTCGCCGCCCTGGCTGGGGAAACATGCTTTTTGGGTGGGATTTACCCTCTTTCTTGCGGTGCGCGGGGCTTTTCTGGCGGCGATTTtgcgg                                        >  2:274935/1‑151 (MQ=255)
                                                                                                             gCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTggcgg                                        <  2:119910/150‑1 (MQ=255)
                                                                                                             gCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTggcgg                                        <  1:200092/150‑1 (MQ=255)
                                                                                                             gCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCCGACGCCCCTTACGCCGCCGTTGCCCGGCAATCATGCCTTTTGGCTGGTATTTAACGTCTTTCTGGCGGTTCCTGGGCCTTCTTTGGCGGGTATCTagcgg                                        >  1:603114/1‑150 (MQ=255)
                                                                                                              ccACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGgcggcg                                      <  2:1109940/151‑1 (MQ=255)
                                                                                                               cACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTggcgg                                        <  2:704797/148‑1 (MQ=255)
                                                                                                                 ccGAAATGCGTAACAGTATGGCGGTGGCCGCCGCCGGGTTTTCGCTGACGCTCCTTACGCTGCCGTGGCTGGGGAATCATGCTTTTTGGGTGGCATTAACCGTCTTTCTTGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTTgcggcg                                      >  2:730825/1‑148 (MQ=255)
                                                                                                                    aaaTGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTg                                <  1:515348/151‑1 (MQ=255)
                                                                                                                     aaTGCGTAACAGTATGGCGGTGGGCCCCGCAGGGTTTTCGCTGACGCTCCTTACGCTGCCGTGGGTGGGGAAACATGGTTTTTGGCTGGGATTAAACGTCTTTTTGGCGGTGCGCGGGGTTTTTCTGGCGGCTCTCTTGCGGCGTCCCttg                               >  2:582714/1‑149 (MQ=255)
                                                                                                                       tGCGTAACAGTATGGCGGTGGCCGCCCCCGGTTTTGCGCCGGCGCTCCTTTCCCTGCCGGGGGTGGGGTATCATGCTTTGTTGCGGGCATTTACCCGCTTTCTGGCGTTGCGCGGGGCTTCTTTTGCGGCTTTTTtgcgggggcc                                   >  2:814591/1‑140 (MQ=255)
                                                                                                                        gCGTAACAGTATGGCGGTGGGCGCCCCCGGGTTTTCGCCTGCGCCCCCTACGCCGCCGTGGGTGGGGAATCATGCTTTTTGGCCGGCCTTAAACGCCTTTTTGGCGGTGCGCGGGGCTTTTCTGGCGGCTATCTtgcgg                                        >  2:153970/1‑139 (MQ=255)
                                                                                                                        gCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGgcg                             <  1:807999/150‑1 (MQ=255)
                                                                                                                            aaCAGTATGGCGGTGGCCGCCGCAGGCTTTGTGCTGACGCCCCCTACGCAGCCGGGGGCGGGGAATCATGCTTTGTGGCTGGCATTAACCCTCTTTTTGGCGGTGGGCGGGCCTTTTCTGGTGGGTTTCTTGCGGGGGCACTg                                >  2:160684/1‑143 (MQ=255)
                                                                                                                                        ggTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCCATGGTACCTGGttt             >  2:63050/1‑151 (MQ=255)
                                                                                                                                        ggTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGttt             <  1:525476/151‑1 (MQ=255)
                                                                                                                                            gccgccGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGCCg         <  1:935219/151‑1 (MQ=255)
                                                                                                                                             ccgccgCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGCCGc        <  1:1217319/151‑1 (MQ=255)
                                                                                                                                                    aGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGCCGCAACGTg  >  1:323390/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                    
AACCCAGATTCAGCAGCTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGCCGCAACGTG  >  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: