Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,258,9740TG42.5% 125.4 / 10.6 73G411G (GGT→GGGdinFoxidative stress resistance protein; putative MATE family efflux pump; UV and mitomycin C inducible protein
Reads supporting (aligned to +/- strand):  ref base T (10/32);  new base G (31/0);  total (41/32)
Fisher's exact test for biased strand distribution p-value = 3.90e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAGATTCAGCAGCTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGCCGCAACGTG  >  NC_000913/4258827‑4259120
                                                                                                                                                   |                                                                                                                                                  
ccAGATTCAGCAGCTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGgttt                                                                                                                                                 >  2:287225/1‑151 (MQ=255)
    aTTCAGCAGCTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCCGGGTTTgc                                                                                                                                               >  1:286991/1‑149 (MQ=255)
      tCAGCAGCTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCt                                                                                                                                            <  2:842552/150‑1 (MQ=255)
       cagcagCTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTtcgc                                                                                                                                             >  2:422022/1‑148 (MQ=255)
        agcagcTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGgttt                                                                                                                                                 >  1:263374/1‑143 (MQ=255)
         gcagcTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCCGGGTTTTCGCTGACg                                                                                                                                        >  1:946338/1‑151 (MQ=255)
           agcTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGGCGCCGCCCCAGGTTTTGCGATg                                                                                                                                           >  1:453021/1‑146 (MQ=255)
            gctggctgACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTc                                                                                                                                     <  2:252728/151‑1 (MQ=255)
            gctggctgACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTc                                                                                                                                     <  2:623349/151‑1 (MQ=255)
               ggctgACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCCGGGTTTGGGCCGACGCTccct                                                                                                                                  >  1:867622/1‑149 (MQ=255)
               ggctgACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCtt                                                                                                                                  <  2:847657/151‑1 (MQ=255)
               ggctgACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCt                                                                                                                                   <  1:443194/150‑1 (MQ=255)
                     ccGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTTCGCTGACGCCCCTTACGc                                                                                                                              >  2:679560/1‑149 (MQ=255)
                     ccGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCCCCGCAGGGTTTGCGCTGACGCCCCCTACGc                                                                                                                              >  2:913195/1‑149 (MQ=255)
                      cGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGGCGCCGCAGGGTTTGCGCTGACGCTCCTTACGCTg                                                                                                                            >  2:333405/1‑150 (MQ=255)
                       gCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGCCGCTCCTTACGCTGcc                                                                                                                          <  1:804403/151‑1 (MQ=255)
                        cTATCTTATCTGGCGGGTGATTGTGCCGGTGGTTGGGGTCTGCTGTTATCTGCGGGGCGGCGTGTTTATGGGCGCGACGCGTGCCGCGGAAATGCTTAACAGGATGGCGGTGGCAGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGc                                                                                                                           <  1:936648/149‑1 (MQ=255)
                          atcttatctGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGt                                                                                                                        <  1:618972/150‑1 (MQ=255)
                          atcttatctGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGt                                                                                                                        <  1:778092/150‑1 (MQ=255)
                              tatctGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCCGGGTTTGCGCTGACGCTCCCTACGCCGCCGGGGCt                                                                                                                    >  2:440246/1‑150 (MQ=255)
                              tatctGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTgg                                                                                                                      <  2:892353/148‑1 (MQ=255)
                              tatctGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCt                                                                                                                    <  1:1003120/150‑1 (MQ=255)
                               atctGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTGCGCTGACGCTCCTTTCGCCGCCgtgggtgg                                                                                                                  >  1:504460/1‑151 (MQ=255)
                               atctGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCCCCCCAGGGTTTGCGCGGGCGCTCCCTTTGCTGcccggggtgg                                                                                                                  >  1:870687/1‑143 (MQ=255)
                                      aGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTGCGCTGACGCCCCTTACGCTGCCGCGGCTGGGGAAt                                                                                                             >  2:902061/1‑149 (MQ=255)
                                          gATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGc                                                                                                        <  2:218482/150‑1 (MQ=255)
                                          gATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGGAATCATGc                                                                                                        >  2:299920/1‑150 (MQ=255)
                                              ttGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTtg                                                                                                    <  1:299270/150‑1 (MQ=255)
                                               tGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGGCGCCGCCGGGTTTGCGCTGGCGCCCCTTTCGCTGCCGGGGCTGGGGAATCATGGttt                                                                                                     >  2:666484/1‑148 (MQ=255)
                                               tGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTtg                                                                                                    <  1:83194/149‑1 (MQ=255)
                                               tGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCCCCGCAGGGTTTTCGCTGACGCTCCTTACGCCGCCGCGGCTGGGGAATCATGCTTtgt                                                                                                   >  2:837221/1‑150 (MQ=255)
                                                gCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGtgg                                                                                                 <  2:576878/151‑1 (MQ=255)
                                                  cGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGtggct                                                                                               <  1:59712/151‑1 (MQ=255)
                                                   ggtggtTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGtggctg                                                                                              <  2:825142/151‑1 (MQ=255)
                                                    gtggtTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTGCGCTGACGCTCCTTACGCTGCCGTGGGTGGGGAATCATGCTTTGtggctgg                                                                                             >  2:175119/1‑151 (MQ=255)
                                                         tGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTGCGCTGACGCCCCTTACGCTGCCGTGGCTGGGGAATCATGATTTGTGGCTGGCAtt                                                                                         >  1:560061/1‑150 (MQ=255)
                                                          ggCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTaa                                                                                       <  2:843914/151‑1 (MQ=255)
                                                              tCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGGGGCCGCCGCCGGGTTTGCGCTGACGCCCCCTTCGCTGCCGCGGCTGGGGAATCATGCTTTTtggctgggaattaac                                                                                     >  2:840663/1‑141 (MQ=255)
                                                               cTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGAGGCCGCCGCAGGGTTTTCGCTGACGCTCCTTACGCTGCCGTGGCTGGGGAATCATGTTTTTTGGCCGGGAtttaac                                                                                     >  1:210045/1‑144 (MQ=255)
                                                                   tgtTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGGCGCCGCCGGTTTTTCGCTGAAGCTCCTTAAGCTGCCGGGGCGGGGGAATCATGCTTTGTGGCTGGCATTTAACCTCttt                                                                               >  1:196776/1‑150 (MQ=255)
                                                                     ttATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCt                                                                             <  1:320938/150‑1 (MQ=255)
                                                                           gctgGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCg                                                                         <  2:324311/148‑1 (MQ=255)
                                                                           gctgGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGtt                                                                       <  2:380006/150‑1 (MQ=255)
                                                                           gctgGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCCCAGGGTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGGAATCATGCTTTGTGGCTGGCATTAACCCTCTTTCTGGCggt                                                                       >  1:625847/1‑148 (MQ=255)
                                                                              ggACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTgcgc                                                                   <  1:469377/151‑1 (MQ=255)
                                                                                  ggCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTTCGCTGACGCTCCTTACGCTGCCGTGGCTGGGGAATCATGCTTTGTGGCTGGCATTAAACGTCTTTTTGGCGTTGCGCggg                                                                >  2:884648/1‑150 (MQ=255)
                                                                                   gCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCt                                                              <  2:325439/151‑1 (MQ=255)
                                                                                     aTGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCCTTTACCCTCTTTTTGGCGTTGCGCGGGCttt                                                            >  2:809741/1‑151 (MQ=255)
                                                                                        tttATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGGCGCCCCCGGGTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATAATGCTTTTTGGGTGGCATTAACCCTCTTTCTTGCGTTTCGCGGGCttt                                                            >  1:1015039/1‑148 (MQ=255)
                                                                                           ataGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTTCGCTGACGCCCCTTTAGCTGCCGTGGGTGGGTAATCATGCTTTGTGGCTGGGATTAAACGTCTTTTTGGCGGTGCGCGGGGTTTCTCTgg                                                      >  2:811259/1‑151 (MQ=255)
                                                                                               gcgcAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCCGGGTTTTCGCTGACGCCCCTTACGCTGCCCGGGCTGGGGAAACATGCTTTTTGGCTGGCCTTTACCCTCTTTTTTGGGTTGCGCGGGCTTTTTCTggcgg                                                   >  2:506616/1‑150 (MQ=255)
                                                                                               gcgcAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTggcggc                                                  <  2:946511/151‑1 (MQ=255)
                                                                                                   aaCGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGGAATCATGCTTTTTGGCTGGCATTAAACGTCTTTCTGGCGGTTCGCGGGGTTTCTCTTGCGGCTATc                                              >  1:298990/1‑151 (MQ=255)
                                                                                                    aCGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCt                                             <  2:691171/151‑1 (MQ=255)
                                                                                                     cgcgTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTGCGCTGACGCCCCTTTCGCTGCCGTGGGTGGGTAAACATTCTTTGTGGCTGGGATTTAACGTCTTTTTTGCGTTTCGCGGGGTTTTTTTGGGGGCTATCt                                             >  2:718782/1‑150 (MQ=255)
                                                                                                       cgTGCCACCGAAATGCGTAACTGTATGGCGGGGGCCGCTGCAGGGTTTGGGCTTGAGCTCCCTAAGCCGCCGTGGCTTGGGAATCATGTTTTGTGGGTGGTAGTTAACGGCTTTCTTGCGTAGGCGGGGGTTTTTCTAGCGGGTTTCtttg                                          >  1:207696/1‑148 (MQ=255)
                                                                                                        gTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTggcg                                         <  2:286995/151‑1 (MQ=255)
                                                                                                        gTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTggc                                          <  1:566902/150‑1 (MQ=255)
                                                                                                         tGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTTCGCTGACGCTCCTTACGCTGCCGTGGCTGGGGAATCATGCTTTTTGGCTGGCATTAAACGTCTTTCTGGCGTTGCGCGGGGTTTCTCTGGCGGCTATCttg                                           >  2:468583/1‑146 (MQ=255)
                                                                                                           ccACCGAAATGCGTAACAGTATGGCGGCGGCCGCCGCAGGGTTTGCGCCGACGGACCCTACGCTGGCGGGGGTGGGGAAAACTGCTTTTTGGCTGGCCCTTACAGGCTTTCTGGCGTTTGGCGGGGCTTTTCTGGGGGGTATCTtgcgg                                        >  1:224292/1‑149 (MQ=255)
                                                                                                             aCCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTc                                    <  1:472288/151‑1 (MQ=255)
                                                                                                              ccGAAATGCGTAACAGTATGGCGGTGGGCGCCGCCGGGTTTTCGCTGACGCCCCTTAAGCTGGCCGGGCTGGGGGATAATGCCTTGTGGGTGGGCTTTACCCTCTTTTTTGCGGTGCGCGGGGCTTTTCTGGCGGGTGTCttg                                           >  1:715164/1‑141 (MQ=255)
                                                                                                                  aaTGCGTAACAGTATGGCGGTGGCCGCCGCCGGGTTTTCGCCGACGCCCCTTACGCTGCCGGGGCTGGGGAAAAATGCTTTGTGGGTGGCATTTACCGGCTTTTTGGGGTTTCGCGGGGCTTTTCTGGGGGGTATCTTGGGGGGTCCCTg                                >  1:908371/1‑150 (MQ=255)
                                                                                                                   aTGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGgc                              <  2:335391/151‑1 (MQ=255)
                                                                                                                       gTAACAGTATGGCGGTGGCCGCCGCAGGGTTTGCGCTGACGCCCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGGCCCTGgcg                             >  2:642126/1‑148 (MQ=255)
                                                                                                                         aaCAGTATGGCGGTGGCCGCCCCCGGGTTTGCGCCGACGCTCCTTACGCTTCCGGGGGGGGGGAATCCTGCTTTGTGGCTGGGCTTAAACGTCTTTCTGGGGTTGCGCGGGCCTTTTCTGGGGGCTATCTTGCGGGGTCACTGGCGCAAt                         >  1:419612/1‑150 (MQ=255)
                                                                                                                             gTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTAc                    <  1:474278/151‑1 (MQ=255)
                                                                                                                               aTGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCt                  <  2:676308/151‑1 (MQ=255)
                                                                                                                                 ggCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTg                 <  2:370572/150‑1 (MQ=255)
                                                                                                                                       ggCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCTGTGGCTGGGGAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGGGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGcc          >  1:127115/1‑151 (MQ=255)
                                                                                                                                       ggCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGcc          >  2:421638/1‑151 (MQ=255)
                                                                                                                                       ggCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGcc          >  2:679542/1‑151 (MQ=255)
                                                                                                                                              gcAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGCCGCAACGt   >  1:341174/1‑151 (MQ=255)
                                                                                                                                               cAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCTCTGGCGCAATGGTCCCTGGTTTGCCGCAACGTg  >  1:726930/1‑151 (MQ=255)
                                                                                                                                                aGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGCCGCAACGTg  >  1:525309/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CCAGATTCAGCAGCTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGCCGCAACGTG  >  NC_000913/4258827‑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: