Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,258,9740TG43.2% 169.5 / 22.1 95G411G (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 (10/44);  new base G (41/0);  total (51/44)
Fisher's exact test for biased strand distribution p-value = 1.06e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCAGCAGCTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGCCGCAACGTGACGG  >  NC_000913/4258833‑4259124
                                                                                                                                             |                                                                                                                                                      
tCAGCAGCTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCt                                                                                                                                                <  2:1218283/150‑1 (MQ=255)
 cagcagCTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGGCGCCGCAGGGTTTTCGCTg                                                                                                                                               >  1:1139998/1‑150 (MQ=255)
  agcagcTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTtcg                                                                                                                                                  >  1:217645/1‑146 (MQ=255)
    cagcTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGc                                                                                                                                           <  1:305563/151‑1 (MQ=255)
    cagcTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGc                                                                                                                                           <  2:647517/151‑1 (MQ=255)
      gctggctgACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTTCGCTGGCGCt                                                                                                                                          >  1:412160/1‑150 (MQ=255)
      gctggctgACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCCCCGGGTTTGCGCTggc                                                                                                                                             >  2:1196529/1‑145 (MQ=255)
        tggctgACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCCGGGTTTTCGCTGCCGcccc                                                                                                                                        >  2:81185/1‑150 (MQ=255)
         ggctgACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCtt                                                                                                                                      <  1:830988/151‑1 (MQ=255)
         ggctgACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCt                                                                                                                                       <  1:1152271/150‑1 (MQ=255)
          gctgACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTa                                                                                                                                     <  1:929125/151‑1 (MQ=255)
          gctgACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTa                                                                                                                                     <  2:511697/151‑1 (MQ=255)
           ctgACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGCCGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGTTTACGCTCATTa                                                                                                                                     <  2:337294/150‑1 (MQ=255)
           ctgACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGGCGCCCCCGGGTTTTCGCTGACGCCCCtt                                                                                                                                      >  2:570982/1‑149 (MQ=255)
            tgACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGGCGCCGCCGGGTTTTCGCTGGCGCTcc                                                                                                                                        >  2:526677/1‑146 (MQ=255)
            tgACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCCGGGTTTGCGCTGGCGCTCCTTACg                                                                                                                                   >  2:674497/1‑151 (MQ=255)
             gACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTAc                                                                                                                                    <  2:860812/149‑1 (MQ=255)
             gACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGc                                                                                                                                  <  2:805425/151‑1 (MQ=255)
             gACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCCCAGGGTTTTTGCTGACGCTCCCTACGc                                                                                                                                  >  2:1162234/1‑151 (MQ=255)
               ccGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTg                                                                                                                                <  1:1230315/151‑1 (MQ=255)
               ccGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTGCGCTGACGCTCCTTACg                                                                                                                                   >  1:682384/1‑148 (MQ=255)
                cGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGAGGCCGCCGCAGGGTTTGCGCCGACGCTcc                                                                                                                                        >  1:1025302/1‑142 (MQ=255)
                 gCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGcc                                                                                                                              <  2:243718/151‑1 (MQ=255)
                    atcttatctGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGGCGCCGCCGGGTTTTCGCTGACGCCCCTTACGCTGCCgcg                                                                                                                           >  1:843469/1‑149 (MQ=255)
                    atcttatctGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCg                                                                                                                             <  2:1013752/149‑1 (MQ=255)
                        tatctGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCCGGGTTTTTGCTGACGCTCCTTAAGCTGGCgggg                                                                                                                          >  2:440871/1‑148 (MQ=255)
                           ctGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGGCGCCGCCGGGTTTTTGCTGACGCTCCCTACGCTGGCCTGGCTggg                                                                                                                     >  1:18096/1‑150 (MQ=255)
                            tGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCCGGGTTTTCGCTGACGCTCCCTTCGCTGCCGTggggggg                                                                                                                     >  2:906825/1‑144 (MQ=255)
                             ggCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTggg                                                                                                                     <  1:1287755/148‑1 (MQ=255)
                             ggCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTaa                                                                                                                  <  1:987641/151‑1 (MQ=255)
                             ggCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTa                                                                                                                   <  1:1072105/150‑1 (MQ=255)
                             ggCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTggg                                                                                                                     >  1:107507/1‑148 (MQ=255)
                               cAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTGCGCTGACGGCCCTTATCCTGCCGGGGCTGGGGAAt                                                                                                                 >  1:895819/1‑150 (MQ=255)
                                aGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGGCGCCGCAGGGTTTGGGGGGAGGCTCCCTTCGCTGGCGGGGTTGGGGAAta                                                                                                                >  2:1337258/1‑149 (MQ=255)
                                aGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCa                                                                                                               <  1:104023/151‑1 (MQ=255)
                                 ggTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCCGGGTTTTCGCTGACGCTCCTTACGCCGCCGGGGCTGGGTAATCa                                                                                                               >  1:88561/1‑150 (MQ=255)
                                    gATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGc                                                                                                            <  1:556918/150‑1 (MQ=255)
                                    gATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCt                                                                                                           <  2:649708/151‑1 (MQ=255)
                                    gATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCt                                                                                                           <  2:1333391/151‑1 (MQ=255)
                                        ttGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTtg                                                                                                        <  2:1179208/150‑1 (MQ=255)
                                        ttGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTGCCCTGACGCCCCTTAAGCTGCCGCGGCTGGGGAATAATGGTTtgt                                                                                                       >  1:374547/1‑151 (MQ=255)
                                            cGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCCGGGTTTTTGCTGACGCCCCTTTCGCTGCCGTGGGTGGGGAAACATGCTTTTtgg                                                                                                     >  2:1438046/1‑149 (MQ=255)
                                             ggtggtTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGtggctg                                                                                                  <  1:967681/151‑1 (MQ=255)
                                              gtggtTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGtggctgg                                                                                                 <  2:405351/151‑1 (MQ=255)
                                                    ggCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGGCGCCGCAGGGTTTGCGCCGACGCCCCTTACGCCGCCGGGGGTGGGGAATCAAGCTTTGTTGGTGGCAtt                                                                                             >  2:888957/1‑149 (MQ=255)
                                                      cGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTTCGCTGACGCCCCTTACGCTGCCGTGGCTGGGTAATTATGCTTTTTGGCTGGCATTTAcc                                                                                         >  1:246875/1‑151 (MQ=255)
                                                      cGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTGCGCTGACGCTCCCTACGCTGCCGTGGGTGGGGAATCATGCTTTGTGGCCGGCAtt                                                                                             >  1:1166128/1‑147 (MQ=255)
                                                       gTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGGCGCCGCAGGGTTTTCGCTGACGCTCCTTAAGCTGCCGGGGCTGGGGAAATATGCTTTTTGGGGGGCATTaa                                                                                           >  2:907012/1‑148 (MQ=255)
                                                        tCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGGCGCCGCCGGGTTTTTGCTGACGCCCCTTTCGCCGGCGGGGGTGGGGGATTATTCTTTTTTGGTGGGAtttaac                                                                                         >  2:358031/1‑145 (MQ=255)
                                                           ggTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTTCGCTGACGCTCCTTTAGCTGCCGGGGGTGGGGAATTATGCTTTGTGGGTGGGATTTACCGTCtt                                                                                    >  1:1238993/1‑151 (MQ=255)
                                                             tgtTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGGCCCCGCAGGGTTTGCGCTGACGCGCCCTATGCTGCCGCGGGTGGGGAATCATTCTTTTTGGCTGGGCTTAACCggcctt                                                                                   >  2:181113/1‑145 (MQ=255)
                                                              gtTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCt                                                                                 <  1:729231/151‑1 (MQ=255)
                                                               ttATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCt                                                                                 <  2:615644/150‑1 (MQ=255)
                                                                tATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTg                                                                                <  1:1009720/150‑1 (MQ=255)
                                                                 aTCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTgg                                                                               <  1:195876/150‑1 (MQ=255)
                                                                    tgctgGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCCGGGTTTTGGCTGGCGCTCCTTTCCCTGCCGGGGCTGGGGGATCATGCTTTTTGGCGGGGATTAACCggcctttt                                                                                 >  2:1047125/1‑138 (MQ=255)
                                                                     gctgGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGtt                                                                           <  1:1221912/150‑1 (MQ=255)
                                                                      ctgGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCCCAGGGTTTGCGCTGACCCCCCTTACGCTGCCGGGGCTGGGTAAATATACTTTTTGGCTGGGATTTACCTTCTTTTTTGCGtt                                                                           >  2:672919/1‑149 (MQ=255)
                                                                        ggACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTgcg                                                                        <  1:1124319/150‑1 (MQ=255)
                                                                        ggACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTgc                                                                         <  2:1278973/149‑1 (MQ=255)
                                                                        ggACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTgc                                                                         <  2:783928/149‑1 (MQ=255)
                                                                              cATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCtt                                                                 <  1:889680/151‑1 (MQ=255)
                                                                               aTGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCCGGGTTTGCGCTGACGCTCCTTACGCTGCCGCGGCTGGGGAATCATGCTTTTTGGCTGGCATTAAACCTCTTTTTTGCGGTGCGCGGGCgtt                                                                >  2:17502/1‑151 (MQ=255)
                                                                               aTGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCttt                                                                <  2:1025277/151‑1 (MQ=255)
                                                                                tGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCCGGGTTTTCGCTGACGCCCCTTACACTGCCGCGGCCGGGGAATCCTGCTTTTTGGCTGGGCTTAACCCTTCTTCTTGCGTTGCGCggggctt                                                                >  2:395782/1‑146 (MQ=255)
                                                                                 gTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCCGGGTTTTCGCTGACGCTCCTTACGCTGCCGTGGCTGGGGAATCATGCTTTTTGGCTGGGCTTTACCCTCTTTCTTGCGTTGGGCGGGGttt                                                                >  2:1309009/1‑149 (MQ=255)
                                                                                 gTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTTCGCCGACGCTCCTTACGCTGCCGGGGCTGGGTAATCATTCTTTTTGGCTGGCATTTACCCTCTTTCTGGGGTTGGGCGGGGttttt                                                              >  1:642900/1‑149 (MQ=255)
                                                                                     ataGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTgg                                                          <  1:906847/151‑1 (MQ=255)
                                                                                     ataGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTg                                                           <  1:20328/150‑1 (MQ=255)
                                                                                     ataGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTTCGCTGACGCCCCTTACGCTGCCGTGGCTGGGTAATTATGCTTTTTGGCTGGCATTAAACGTCTTTCTTGGGGTGCGCGGGCttt                                                                >  1:1268541/1‑145 (MQ=255)
                                                                                                gcgTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCCCCGCAGGGTTTTCGCTGACGCTCCCTAAGCCGCCGGGGGTGGGGAATCATGCTTTTTGGGTGGGATTAAACCTCTTTTTGGGGGTGCGCGGGGGTTCTCTGGCGGGTAt                                                   >  2:1100882/1‑147 (MQ=255)
                                                                                                 cgTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCCGGGTTTTCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTTACCGTCTTTCTTGCGTTGCGCGGGCTTTTTCTTGCGGCTATCttg                                               >  1:249053/1‑148 (MQ=255)
                                                                                                  gTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCCGGGTTTTTGCCGACGCCCCTTAAGCCGCCCTGGGTGGGTAATCATGCTTTTTGGCTGGCATTAACCGTCTTTTTGGCGGTGGGCGGGGTTTTTCTCGCGGGTATCt                                                 >  2:738445/1‑147 (MQ=255)
                                                                                                    gCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTggcggc                                           <  1:910924/151‑1 (MQ=255)
                                                                                                      cACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGGTTTTCGCTGACGCTCCTTACGCTGCCGTGGGTGGGGAATCATTCCTTTTGGGCGGGCTTAAACCTCTTTCTTGCGGTGTGCGGGCCTTTTTTTGCGGCCATCTggggg                                            >  2:1320969/1‑148 (MQ=255)
                                                                                                      cACCGAAATGCGTAACAGTATGGCGGTGGCCCCCGCCGGTTTTGCGCTGACGCTCCTTACGCTGCCGCGGCTGGGGAATCATGCTTTTTGGGTGGGATTAACCGTCTTTCTGGCGGTGCGCGGGCTTTCTCTGGCGGGTATCTGgcggcg                                          >  2:74417/1‑150 (MQ=255)
                                                                                                        ccGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCa                                       <  1:860811/151‑1 (MQ=255)
                                                                                                        ccGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGGGGGTGGGGAAACATGCTTTGTGGCTGGCATTTACCGTCTTTTTGGGGTTGCGCGGGGTTTCTCTGGCGGCTATCTggcgggggcc                                       >  1:1251541/1‑146 (MQ=255)
                                                                                                         cGAAATGTGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTTCGCTGACGCCCCCTTCGCTGCCGTGGCTGGGGAATCATGGTTTTTGGGTGGCATTTAACGGCTTTTTTGCGTTGCGCGGGGTTTTTCTGGCGGCTATCTTGCGGCGTc                                        >  1:644634/1‑149 (MQ=255)
                                                                                                         cGAAATGCGTAACAGTATGGCGGTGGGCGCCGCAGGGTTTTCGCTGACGCTCCTTACGCTGCCGGGGGTGGGGAATCATTCTTTTTGGGTGGCATTTAACGCCTTTTTTGCGGTGCGCGGGGTTTTTCTGGCGGCTATCTTGCGGCGGCAc                                      >  2:1254675/1‑151 (MQ=255)
                                                                                                         cGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTTCGCTGACGCCCCTTACGCTGCCGGGGGTGGGGAATCATGCTTTTTGGGTGGCATTAACCGGCTTTTTTGCGGTGGGCGGGGCTTCTCTTGCGGCTATCTGGCGGCGGCAc                                      >  2:17397/1‑151 (MQ=255)
                                                                                                            aaTGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTgg                                   <  1:1347739/151‑1 (MQ=255)
                                                                                                            aaTGCGTAACAGTATGGCGGTGGCCGCCCCAGGTTTTTCGCTGACGCTCCTTACGCTGCCGTGGCTGGGGAATCATGCTTTTTGGCTGGGATTTAACGTCTTTCTTGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGGGGCGGCACTgg                                   >  1:559207/1‑151 (MQ=255)
                                                                                                                    aCAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATgg                           <  1:358037/151‑1 (MQ=255)
                                                                                                                     cAGTATGGCGGTGGGCCCCGCAGGTTTTTCGGTGACGCTCCTTACGCTGCCGGGGCTGGGCGATCATGCTTTGTGGCTGGCATTAACCCTCCTTCTGGCGTTTCGCCGGCCTTCTCTCGCGGGTATCTGGCGGCGTCACTgg                                   >  2:525696/1‑142 (MQ=255)
                                                                                                                     cAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGt                          <  1:526685/151‑1 (MQ=255)
                                                                                                                         aTGGCGGTGGCCGCCGCCGGGTTTTCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGGTTTCTCTGGCGGCTATCTTGCGGCGTCACTGGCGCAATgggaac                       >  2:297839/1‑146 (MQ=255)
                                                                                                                              ggTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGttt                 <  1:575484/151‑1 (MQ=255)
                                                                                                                                tGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGc               >  1:1117425/1‑151 (MQ=255)
                                                                                                                                 ggCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGcc              >  1:1085595/1‑151 (MQ=255)
                                                                                                                                   ccgccgCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGCCGc            >  1:1103357/1‑151 (MQ=255)
                                                                                                                                           ggTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGCCGCAACGTGAc    >  2:808108/1‑151 (MQ=255)
                                                                                                                                           ggTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGCCGCAACGTGAc    <  1:1196560/151‑1 (MQ=255)
                                                                                                                                            gTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGCCGCAACGTGACg   <  2:107504/151‑1 (MQ=255)
                                                                                                                                            gTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTC‑TTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGCCGCAACGTGACgg  <  1:1162265/151‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                      
TCAGCAGCTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGCCGCAACGTGACGG  >  NC_000913/4258833‑4259124

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑