Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 2,040,296 (G)10→11 coding (148/1017 nt) PA14_RS09495 → glycosyltransferase family 4 protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084632,040,2861.G94.9% 1106.0 / 43.4 295coding (138/1017 nt)PA14_RS09495glycosyltransferase family 4 protein
Reads supporting (aligned to +/- strand):  ref base . (11/4);  new base G (135/145);  total (146/149)
Fisher's exact test for biased strand distribution p-value = 6.72e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 339 total aligned reads displayed.

AGCATCAACAATAATGAACTTCTGGTTTTGGCTTCCCGTAGTTATGGGGCTTTCTTTAGTCTTGACCGGAGGTCTGCGTCGATACGCTTTGGCGCGTAGCTTGATGGACGTTCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGT‑‑GGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTGGTTTTGCTGATTCTGGCGGCCCTTGGTCAAGTGAGCTGGGACTATACGTGGGCGTTGCTAGGTGCTGGTATTGGCATTGCGGTTCTTGGCTTTCTAGACGATCATGGCC  >  NC_008463/2040136‑2040431
                                                                                                                                                        |                                                                                                                                                   
      aacaaTAATGAACTTCTGGTTTTGGCTTCCCGTAGTTATGGGGCTTTCTTTAGTCTTGACCGGAGGTCTGCGTCGATACGCTTTGGCGCGTAGCTTGATGGACGTTCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGT‑‑gggggg                                                                                                                                              <  2:9542281/151‑1 (MQ=255)
       acaaTAATGAACTTCTGGTTTTGGCTTCCCGTAGTTATGGGGCTTTCTTTAGTCTTGACCGGAGGTCTGCGTCGATACGCTTTGGCGCGTAGCTTGATGGACGTTCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGT‑‑ggggggg                                                                                                                                             <  2:5592500/151‑1 (MQ=255)
       acaaTAATGAACTTCTGGTTTTGGCTTCCCGTAGTTATGGGGCTTTCTTTAGTCTTGACCGGAGGTCTGCGTCGATACGCTTTGGCGCGTAGCTTGATGGACGTTCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGT‑‑ggggggg                                                                                                                                             >  1:3515755/1‑151 (MQ=255)
        caaTAATGAAATTCTGGTTTTGGCGTCCAGTAGTTATGGGGCTTTCTTTAGTCTTGACAGGAGTTCTGCGAAGATACGCTTATGCTCGAAGATTGATGGACGTTCCCAATGCCCGTATCTCTCACTATGTACCTACACCCCGT‑‑gggtgggg                                                                                                                                            >  2:5193048/1‑151 (MQ=255)
         aataatGAACTTCTGGTTTTGGCTTCCCGTAGTTATGGGGCTTTCTTTAGTCTTGACCGGAGGTCTGCGTCGATACGCTTTGGCGCGTAGCTTGATGGACGTTCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGT‑‑ggggggggg                                                                                                                                           <  2:552407/151‑1 (MQ=255)
         aataatGAACTTCTGGTTTTGGCTTCCCGTAGTTATGGGGCTTTCTTTAGTCTTGACCGGAGGTCTGCGTCGATACGCTTTGGCGCGTAGCTTGATGGACGTTCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGT‑‑ggggggggg                                                                                                                                           >  2:8593442/1‑151 (MQ=255)
          ataatGAACTTCTGGTTTTGGCTTCCCGTAGTTATGGGGCTTTCTTTAGTCTTGACCGGAGGTCTGCGTCGATACGCTTTGGCGCGTAGCTTGATGGACGTTCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGT‑‑gggggggggg                                                                                                                                          <  2:935671/151‑1 (MQ=255)
             atGAACTTCTGGTTTTGGCTTCCCGTAGTTATGGGGCTTTCTTTAGTCTTGACCGGAGGTCTGCGTCGATACGCTTTGGCGCGTAGCTTGATGGACGTTCCCAATGACCGTAGCTCTCACTCGGTACCTAAACCACGT‑‑gggggggggggtg                                                                                                                                       >  2:7057587/1‑148 (MQ=255)
              tGAACTTCTGGTTTTGGCTTCCCGTAGTTATGGGGCTTTCTTTAGTCTTGACCGGAGGTCTGCGTCGATACGCTTTGGCGCGTAGCTTGATGGACGTTCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGT‑‑gggggggggggtgg                                                                                                                                      >  1:6116805/1‑147 (MQ=255)
                 aCTTCTGGTTTTGGCTTCCCGTAGTTATGGGGCTTTCTTTAGTCTTGACCGGAGGTCTGCGTCGATACGCTTTGGCGCGTAGCTGGATGGACGTTCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTa                                                                                                                                   >  1:3012509/1‑151 (MQ=255)
                  cTTCTGGTTTTGGCTTCCCGTAGTTATGGGGCTTTCTTTAGTCTTGACCGGAGGTCTGCGTCGATACGCTTTGGCGCGTAGCTTGATGGACGTTCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTAt                                                                                                                                  >  1:7252516/1‑151 (MQ=255)
                    tCTGGTTTTGGCTTCCCGTAGTTATGGGGCTTTCTTTAGTCTTGACCGGAGGTCTGCGTCGCTACGCTTTGGCGCGTAGCTTGATGGACGTTACCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCg                                                                                                                                >  2:3847166/1‑151 (MQ=255)
                    tCTGGTTTTGGCTTCCCGTAGTTATGGGGCTTTCTTTAGTCTTGACCGGAGGTCTGCGTCGATACGCTTTGGCGCGTAGCTTGATGGACGTTCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCg                                                                                                                                >  1:9815202/1‑151 (MQ=255)
                    tCTGGTTTTGGCTTCCCGTAGTTATGGGGCTTTCTTTAGTCTTGACCGGAGGTCTGCGTCGATACGCTTTGGCGCGTAGCTTGATGGACGTTCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCg                                                                                                                                >  2:8208436/1‑151 (MQ=255)
                    tCTGGTTTTGGCTTCCCGTAGTTATGGGGCTTTCTTTAGTCTTGACCGGAGGTCTGCGTCGATACGCTTTGGCGCGTAGCTTGATGGACGTTCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCg                                                                                                                                >  1:629834/1‑151 (MQ=255)
                        gTTTTGGCTTCCCGTAGTTATGGGGCTTTCTTTAGTCTTGACCGGAGGTCTGCGTCGATACGCTTTGGCGCGTAGCTTGATGGACGTTCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTAt                                                                                                                            >  2:9088408/1‑151 (MQ=255)
                        gTTTTGGCTTCCCGTAGTTATGGGGCTTTCTTTAGTCTTGACCGGAGGTCTGCGTCGATACGCTTTGGCGCGTAGCTTGATGGACGTTCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTAt                                                                                                                            >  1:4060153/1‑151 (MQ=255)
                         ttttGGCTTCCCGTAGTTATGGGGCTTTCTTTAGTCTTGACCGGAGGTCTGCGTCGATACGCTTTGGCGCGTAGCTTGATGGACGTTCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATg                                                                                                                           >  2:7478504/1‑151 (MQ=255)
                           ttGGCTTCCCGTAGTTATGGGGCTTTCTTTAGTCTTGACCGGAGGTCTGCGTCGATACGCTTTGGCGCGTAGCTTGATGGACGTTCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAg                                                                                                                         >  1:1526252/1‑151 (MQ=255)
                             ggCTTCCCGTAGTTATGGGGCTTTCATTAGTCTTGACCGGAGGTCTGCGTCGATACGCTTTGGCGCGTAGCTTGATGGACGTTCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGtt                                                                                                                       <  2:10095599/151‑1 (MQ=255)
                               cTTCCCGTAGTTATGGGGCTTTCTTTAGTCTTGACCGGAGGTCTGCGTCGATACGCTTTGGCGCGTAGCTTGATGGACGTTCCCAATGACCGTAGCTCTCACTCGGTACCTACACACCGTGGGGGGGGGGGTGGGCAATCGTTATGAGtta                                                                                                                      >  2:7217932/1‑150 (MQ=255)
                                 tCCCGTAGTTATGGGGCTTTCTTTAGTCTTGACCGGAGGTCTGCGTCGATACGCTTTGGCGCGTAGCTTGATGGACGTTCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTAATGAGTTTCCt                                                                                                                   >  2:5732763/1‑151 (MQ=255)
                                 tCCCGGAGTTATGGGGCTTTCTTTAGTCTTGACCGGAGGTCTGCGTCGATACGCTGTGGCGCGTAGCTTGATGGACGGTCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTTTACt                                                                                                                   >  2:10385196/1‑151 (MQ=255)
                                   ccGTATTTATGGGGCTTTCATTAGTCTTGACCGGCGGTCTGCGTCGTTACGCGTTGGCGCGTTGCTTGATGGACGTTCCCAATGCCCGTAGCTCTCACACGGTACCTACACCCCGGG‑‑GGGGGGGGGGTGGCTATCGTTATGAGTTTCCTAc                                                                                                                 <  2:2539368/151‑1 (MQ=255)
                                    cGTAGTTATGGGGCTTTCTTTAGTCTTGACCGGAGGTCTGCGTCGATACGCTTTGGCGCGTAGCTTGATGGACGTTCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACt                                                                                                                >  1:3348741/1‑151 (MQ=255)
                                    cGTAGTTATGGGGCTTTCTTTAGTCTTGACCGGAGGTCTGCGTCGATACGCTTTGGCGCGTAGCTTGATGGACGTTCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACt                                                                                                                >  1:6977885/1‑151 (MQ=255)
                                     gTAGTTATGGGGCTTTCTTTAGTCTTGACCGGAGGTCTGCGTCGATACGCTTTGGCGCGTAGCTTGATGGACGTTCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTg                                                                                                               >  2:4299220/1‑151 (MQ=255)
                                     gTAGTTATGGGGCTTTCTTTAGTCTTGACCGGAGGTCTGCGTCGATACGCTTTGGCGCGTAGCTTGATGGACGTTCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTg                                                                                                               >  1:494625/1‑151 (MQ=255)
                                     gTAGTTATGGGGCTTTCTTTAGTCTTGACCGGAGGTCTGCGTCGATACGCTTTGGCGCGTAGCTTGATGGACGTTCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTg                                                                                                               >  2:9032328/1‑151 (MQ=255)
                                      tAGTTATGGGGCTTTCTTTAGTCTTGACCGGAGGTCTGCGTCGATACGCTTTGGCGCGTAGCTTGATGGACGTTCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTgg                                                                                                              <  2:5048237/151‑1 (MQ=255)
                                       aGTTATGGGGCTTTCTTTAGTCTTGACCGGAGGTCTGCGTCGATACGCTTTGGCGCGTAGCTTGATGGACGTTCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTGGt                                                                                                             <  2:7638143/151‑1 (MQ=255)
                                       aGTTATGGGGCTTTCTTTAGTCTTGACCGGAGGTCTGCGTCGATACGCTTTGGCGCGTAGCTTGATGGACGTTCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTGGt                                                                                                             >  1:5972048/1‑151 (MQ=255)
                                         ttATGGGGCTTTCTTTAGTCTTGACCGGAGGTCTGCGTCGATACGCTTTGGCGCGTAGCTTGATGGACGTTCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTGGttt                                                                                                           >  1:9550824/1‑151 (MQ=255)
                                            tGGGGCTTTCTTTAGTCTTGACCGGAGGTCTGCGTCGATACGCTTTGGCGCGTAGCTTGCTGGACGTTCCCAATGCCCGTAGCTCTCCCTCGGTACCTCCACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTGGTTTTGc                                                                                                        <  2:2080347/151‑1 (MQ=255)
                                            tGGGGCTTTCTTTAGTCTTGACCGGAGGTCTGCGTCGATACGCTTTGGCGCGTAGCTTGCAGGACGTTCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTGGTTTTGc                                                                                                        <  2:9113032/151‑1 (MQ=255)
                                                gCTTTCTTTAGTCTTGACCGGAGGTCTGCGTCGATACGCTTTGGCGCGTAGCTTGATGGACGTTCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTGGTTTTGCTGAt                                                                                                    <  2:605139/151‑1 (MQ=255)
                                                 ctttctttAGTCTTTACCGGAGGTCTGCGTCGATCCGCTTGGGCGCGTAGCTTGATGGACGTTCCCATGGCCCTTAGCTCCCACCCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTGGTTTTGCTGAtt                                                                                                   <  2:4537329/151‑1 (MQ=255)
                                                 ctttctttAGTCTTGACCGGAGGTCTGCGTCGATACGCTTTGGCGCGTAGCTTGATGGACGTTCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTGGTTTTGCTGAtt                                                                                                   <  1:2225784/151‑1 (MQ=255)
                                                  tttctttAGTCTTGACCGGAGGTCTGCGTCGATACGCTTTGGCGCGTAGCTTGATGGACGTTCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTGGTTTTGCTGATTc                                                                                                  >  1:1920543/1‑151 (MQ=255)
                                                    tctttAGTCTTGACCGGAGGTCTGCGTCGATACGCTTTGGCGCGTAGCTTGATGGACGTTCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTGGTTTTGCTGATTCTg                                                                                                >  1:7421548/1‑151 (MQ=255)
                                                        tAGTCTTGACCGGAGGTCTGCGTCGATACGCTTTGGCGCGTAGCTTGATGGACGTTCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTGGTTTTGCTGATTCTggcgg                                                                                            <  1:476347/151‑1 (MQ=255)
                                                         aGTCTTGACCGGAGGTCTGCGTCGATACGCTTTGGCGCGTAGCTTGATGGACGTTCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTGGTTTTGCTGATTCTgggggc                                                                                           >  2:4998445/1‑151 (MQ=255)
                                                             ttGACCGGAGGTCTGCGTCGATACGCTTTGGCGCGTAGCTTGATGGACGTTCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTGGTTTTGCTGATTCTGGCGGCCCtt                                                                                       <  1:2614055/151‑1 (MQ=255)
                                                              tGACCGGAGGTCTGCGTCGATACGCTTTGGCGCGTAGCTTGATGGACGTTCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGGG‑‑GGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTGGTTTTGCTGATTCTGGCGGCCCTTg                                                                                      >  1:948824/1‑151 (MQ=255)
                                                              tGACCGGAGGTCTGCGTCGATACGCTTTGGCGCGTAGCTTGATGGACGTTCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTGGTTTTGCTGATTCTGGCGGCCCTTg                                                                                      <  2:638396/151‑1 (MQ=255)
                                                                aCCGGAGGTCTGCGTCGATACGCTTTGGCGCGTAGCTTGATGGACGTTCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTGGTTTTGCTGATTCTGGCGGCCCTTGGt                                                                                    >  1:5148743/1‑151 (MQ=255)
                                                                 ccGGAGGTCTGCGTCGATACGCTTTGGCGCGTAGCTTGATGGACGTTCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTTTCCTAATGGTTTTGCTGATTATGGCGGCCCTTGGTc                                                                                   >  2:9464565/1‑151 (MQ=255)
                                                                   ggAGGTCTGCGTCGATACGCTTTGGCGCGTAGCTTGATGGACGTTCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTGGTTTTGCTGATTCTGGCGGCCCTTGGTCaa                                                                                 >  1:1467477/1‑151 (MQ=255)
                                                                     aGGTCTGCGTCGATACGCTTTGGCGCGTAGCTTGATGGACGTTCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTGGTTTTGCTGATTCTGGCGGCCCTTGGTCAAGt                                                                               >  2:3285073/1‑151 (MQ=255)
                                                                     aGGTCTGCGTCGATACGCTTTGGCGCGTAGCTTGATGGACGTTCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTGGTTTTGCTGATTCTGGCGGCCCTTGGTCAAGt                                                                               <  2:9768061/151‑1 (MQ=255)
                                                                      ggTCTGCGTCGATACGCTTTGGCGCGTAGCTTGATGGACGTTCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTGGTTTTGCTGATTCTGGCGGCCCTTGGTCAAGTg                                                                              <  1:9101608/151‑1 (MQ=255)
                                                                       gTCTGCGTCGATACGCTTTGGCGCGTAGCTTGATGGACGTTCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTGGTTTTGCTGATTCTGGCGGCCCTTGGTCAAGTGa                                                                             >  1:883592/1‑151 (MQ=255)
                                                                        tCTGCGTCGATACGCTTTGGCGCGTAGCTTGATGGACGTTCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTGGTTTTGCTGATTCTGGCGGCCCTTGGTCAAGTGAg                                                                            >  1:995543/1‑151 (MQ=255)
                                                                         cTGCGTCGATACGCTTTGGCGCGTAGCTTGATGGACGTTCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTGGTTTTGCTGATTCTGGCGGCCCTTGGTCAAGTGAGc                                                                           >  1:7922651/1‑151 (MQ=255)
                                                                          tGCGTCTATACGCTTTGGCGCGTATCTTGATGGACGTTCCCATTGCCCGTAGCTCTACATCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTGGTTTTGCTGATTCTGGCGGCCCTTGGTCAAGTGAGCt                                                                          <  2:10543538/151‑1 (MQ=255)
                                                                          tGCGTCGATACGCTTTGGCGCGTAGCTTGATGGACGTTCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTGGTTTTGCTGATTCTGGCGGCCCTTGGTCAAGTGAGCt                                                                          >  1:1078289/1‑151 (MQ=255)
                                                                          tGCGTCGATACGCTTTGGCGCGTAGCTTGATGGACGTTCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTGGTTTTGCTGATTCTGGCGGCCCTTGGTCAAGTGAGCt                                                                          <  2:1341981/151‑1 (MQ=255)
                                                                          tGCGTCGATACGCTTTGGCGCGTAGCTTGATGGACGTTCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTGGTTTTGCTGATTCTGGCGGCCCTTGGTCAAGTGAGCt                                                                          >  1:10212379/1‑151 (MQ=255)
                                                                            cGTCGATACGCTTTGGCGCGTAGCTTGATGGACGTTCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTGGTTTTGCTGATTCTGGCGGCCCTTGGTCAAGTGAGCTgg                                                                        >  1:4663427/1‑151 (MQ=255)
                                                                             gTCGATACGCTTTGGCGCGTAGCTTGATGGACGTTCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATTAGTTTACTACTGGTTTTGCTGATTATTGCGGCCCATGGTCAAGTGAGCTggg                                                                       >  1:8764217/1‑151 (MQ=255)
                                                                             gTCGATACGCTTTGGCGCGTAGCTTGATGGACGTTCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTGGTTTTGCTGATTCTGGCGGCCCTTGGTCAAGTGAGCTggg                                                                       >  1:8165992/1‑151 (MQ=255)
                                                                                gATACGCTTTGGCGCGTAGCTTGATGGACGTTCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGGG‑‑GGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTGGTTTTGCTGATTCTGGCGGCCCTTGGTCAAGTGAGCTGGGACt                                                                    >  1:1650610/1‑151 (MQ=255)
                                                                                gATACGCTTTGGCGCGTAGCTTGATGGACGTTCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTGGTTTTGCTGATTCTGGCGGCCCTTGGTCAAGTGAGCTGGGACt                                                                    >  1:5037050/1‑151 (MQ=255)
                                                                                 aTACGCTTTGGCGCGTAGCTTGATGGACGTTCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTTTCATACTGGTTTTGCTGATTCTGGCGGCCCTTGGTCAAGTGAGCTGGGACta                                                                   >  2:9256557/1‑151 (MQ=255)
                                                                                  tccGCTTTGGCGCGTAGCTTGATGGACGTTCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTGGTTTTGCTGATTCTGGCGGCCCTTGGTCAAGTGAGCTGGGACtat                                                                  <  2:6883577/149‑1 (MQ=255)
                                                                                  tACGCTTTGGCGCGTAGCTTGATGGACGTTCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTGGTTTTGCTGATTCTGGCGGCCCTTGGTCAAGTGAGCTGGGACtat                                                                  >  1:4376138/1‑151 (MQ=255)
                                                                                   aCGCTTTGGCGCGTCGCTTGAGGGACGGTCCCACTGCCCGTAGCTCACAACCGGTAACTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTGGTTTTGCTGATTCTGGCGGCCCTTGGTCAAGTGAGCTGGGACtata                                                                 <  2:6570771/151‑1 (MQ=255)
                                                                                    cGCTTTGGCGCGTAGCTTGATGGACGTTCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTGGTTTTGCTGATTCTGGCGGCCCTTGGTCAAGTGAGCTGGGACTATAc                                                                <  1:3826394/151‑1 (MQ=255)
                                                                                      cTTTGTCTCGTAGCTTGATGGCCGGTCCAATTTCCCGTAGCTCTCTCTCGTTAACTAAACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTGTCCTACTGGTTTGGCTGATTCTGGCGGCCCTTGGTCAAGTGAGCTGGGACTATACGt                                                              <  2:473831/151‑1 (MQ=255)
                                                                                      cTTTGGCGCGTAGCTTGATGGACGTTCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTGGTTTTGCTGATTCTGGCGGCCCTTGGTCAAGTGAGCTGGGACTATACGt                                                              >  1:3211886/1‑151 (MQ=255)
                                                                                       tttGGCGCGTAGCTTGATGTAGGTTCCCAAGGCCCGTTGCTCACACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTGGTTTTGCTGATTCTGGCGGCCCTTGGTCAAGTGAGCTGGGACTATACGTg                                                             <  2:7966345/151‑1 (MQ=255)
                                                                                       tttGGCGCGTAGCTTGATGGACGTTCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTGGTTTTGCTGATTCTGGCGGCCCTTGGTCAAGTGAGCTGGGACTATACGTg                                                             <  1:9052392/151‑1 (MQ=255)
                                                                                            cgcgTAGCTTGATGGACGTTCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGGTATGAGTTTCCTACTGGGTTTGCTGATTCTGGCGGCCCTTGGTCAAGGGAGCTGGGACTATAGGTGGGCGt                                                        >  2:4356837/1‑151 (MQ=255)
                                                                                              cgTAGCTTGATGGACGTTCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTGGTTTTGCTGATTCTGGCGGCCCTTGGTCAAGTGAGCTGGGACTATACGTGGGCGTTg                                                      >  2:3702068/1‑151 (MQ=255)
                                                                                                tAGCTTGATGGCCGTTCCCAATGCCCGTTGCTCCCCCCCGGTACCCACACCCCGTGGGGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTGGTTTTGCTGATTCTGGCGGCCCTTGGTCAAGTGAGCTGGGACTATACGTGGGCGTTGc                                                     <  1:690026/151‑1 (MQ=255)
                                                                                                 aGCTTGATGGCCGTTCCCAATGACCGTAGCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTGGTTTTGCTGATTCTGGCGGCCCTTGGTCAAGTGAGCTGGGACTATACGTGGGCGTTGCTa                                                   <  2:3201855/151‑1 (MQ=255)
                                                                                                     tGATGGACGTTCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTGGTTTTGCTGATTCTGGCGGCCCTTGGTCAAGTGAGCTGGGACTATACGTGGGCGTTGATAGGTg                                               <  1:8246037/151‑1 (MQ=255)
                                                                                                         ggACGTTCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTGGTTTTGCTGATTCTGGCGGCCCTTGGTCAAGTGAGCTGGGACTATACGTGGGCGTTGCTAGGTGCTgg                                           >  1:6033546/1‑151 (MQ=255)
                                                                                                          gACGTTCCCAATGCCCGTAGCTCTCACGCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTGGTTTTGCTGATTCTGGCGGCCCTTGGTCAAGTGAGCTGGGACTATACGTGGGCGTTGCTAGGTGCTGGt                                          <  2:2300202/151‑1 (MQ=255)
                                                                                                            cGTTCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTGGTTTTGCTGATTCTGGCGGCCCTTGGTCAAGTGAGCTGGGACTATACGTGGGCGTTGCTAGGTGCTGGTAt                                        >  1:4245000/1‑151 (MQ=255)
                                                                                                               tCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTGGTTTTGCTGATTCTGGCGGCCCTTGGTCAAGTGAGCTGGGACTATACGTGGGCGTTGCTAGGTGCTGGTATTgg                                     <  1:8541756/151‑1 (MQ=255)
                                                                                                                cccAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTGGTTTTGCTGATTCTGGCGGCCCTTGGTCAAGTGAGCTGGGACTATACGTGGGCGTTGCTAGGTGCTGGTATTGGc                                    >  1:6639716/1‑151 (MQ=255)
                                                                                                                cccAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTGGTTTTGCTGATTCTGGCGGCCCTTGGTCAAGTGAGCTGGGACTATACGTGGGCGTTGCTAGGTGCTGGTATTGGc                                    <  1:7587551/151‑1 (MQ=255)
                                                                                                                  cAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCGATAGTTATGAGTTTCCTACTGGTTTTGCTGATTCTGGCGGCCCTTGGTCAAGTGAGCTGGGACTATACGTGGGCGTTGCTAGGTGCTGGTATTGGCAt                                  >  2:1172390/1‑151 (MQ=255)
                                                                                                                   aaTGCCCGTAGCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTGGTTTTGCTGATTCTGGCGGCCCTTGGTCAAGTGAGCTGGGACTATACGTGGGCGTTGCTAGGTGCTGGTATTGGCAtt                                 <  2:1521725/151‑1 (MQ=255)
                                                                                                                   aaTGCCCGTAGCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTGGTTTTGCTGATTCTGGCGGCCCTTGGTCAAGTGAGCTGGGACTATACGTGGGCGTTGCTAGGTGCTGGTATTGGCAtt                                 <  2:1671564/151‑1 (MQ=255)
                                                                                                                    aTGCCCGTAGCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTGGTTTTGCTGATTCTGGCGGCCCTTGGTCAAGTGAGCTGGGACTATACGTGGGCGTTGCTAGGTGCTGGTATTGGCATTg                                >  1:6484123/1‑151 (MQ=255)
                                                                                                                    aTGCCCGTAGCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTGGTTTTGCTGATTCTGGCGGCCCTTGGTCAAGTGAGCTGGGACTATACGTGGGCGTTGCTAGGTGCTGGTATTGGCATTg                                <  2:8781236/151‑1 (MQ=255)
                                                                                                                     tGCCCGTAGCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTGGTTTTGCTGATTCTGGCGGCCCTTGGTCAAGTGAGCTGGGACTATACGTGGGCGTTGCTAGGTGCTGGTATTGGCATTGc                               <  2:2881584/151‑1 (MQ=255)
                                                                                                                             gCTCTCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTGGTTTTGCTGTTTCTGGCGGCCCTTGGTCAAGTGAGCTGGGACTATACGTGGGCGTTGCTAGGTGCTGGTATTGGCATTGCGGTTCTTg                       >  2:10660772/1‑151 (MQ=255)
                                                                                                                              ctctCACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTGGTTTTGCTGATTCTGGCGGCCCTTGGTCAAGTGAGCTGGGACTATACGTGGGCGTTGCTAGGTGCTGGTATTGGCATTGCGGTTCTTgg                      >  1:10359153/1‑151 (MQ=255)
                                                                                                                                ctcACTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTAGGAGTTTGCTACTGGTTTTGCTGATTCTGGCGGCCCTTGGTCAAGGGAGCTGGGACTATACGTGGGCGTTGCTAGGTGGTGGTATTGGCATTGCGGTTCTTGGCt                    >  2:2178306/1‑151 (MQ=255)
                                                                                                                                 tcACGCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTGGTTTTGCTGATTCTGGCGGCCCTTGGTCAAGTGAGCTGGGACTATACGTGGGCGTTGCTAGGTGCTGGTATTGGCATTGCGGTTCTTGGCtt                   <  2:6636102/151‑1 (MQ=255)
                                                                                                                                   aCTCGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTGGTTTTGCTGATTCTGGCGGCCCTTGGTCAAGTGAGCTGGGACTATACGTGGGCGTTGCTAGGTGCTGGTATTGGCATTGCGGTTCTTGGCTTTc                 <  2:7737259/151‑1 (MQ=255)
                                                                                                                                      cgttACCTACACCCCGT‑‑GGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTGGTTTTGCTGATTCTGGCGGCCCTTGGTCAAGTGAGCTGGGACTATACGTGGGCGTTGCTAGGTGCTGGTATTGGCATTGCGGTTCTTGGCTTTCTAGa             <  2:178071/148‑1 (MQ=255)
                                                                                                                                      cGGTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTGGTTTTGCTGATTCTGGCGGCCCTTGGTCAAGTGAGCTGGGACTATACGTGGGCGTTGCTAGGTGCTGGTATTGGCATTGCGGTTCTTGGCTTTCTAg              >  1:5768549/1‑151 (MQ=255)
                                                                                                                                       ggTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTGGTTTTGCTGATTCTGGCGGCCCTTGGTCAAGTGAGCTGGGACTATACGTGGGCGTTGCTAGGTGCTGGTATTGGCATTGCGGTTCTTGGCTTTCTAGa             <  1:3926416/151‑1 (MQ=255)
                                                                                                                                        gTACCTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTGGTTTTGCTGATTCTGGCGGCCCTTGGTCAAGTGAGCTGGGACTATACGTGGGCGTTGCTAGGTGCTGGTATTGGCATTGCGGTTCTTGGCTTTCTAGAc            >  1:5328702/1‑151 (MQ=255)
                                                                                                                                          aacTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTGGTTTTGCTGATTCTGGCGGCCCTTGGTCAAGTGAGCTGGGACTATACGTGGGCGTTGCTAGGTGCTGGTATTGGCATTGCGGTTCTTGGCTTTCTAGACGa          <  2:399475/149‑1 (MQ=255)
                                                                                                                                            cTACACCCCGTGGGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTGGTTTTGCTGATTCTGGCGGCCCTTGGTCAAGTGAGCTGGGACTATACGTGGGCGTTGCTAGGTGCTGGTATTGGCATTGCGGTTCTTGGCTTTCTAGACGATc        >  2:1948776/1‑151 (MQ=255)
                                                                                                                                                        |                                                                                                                                                   
AGCATCAACAATAATGAACTTCTGGTTTTGGCTTCCCGTAGTTATGGGGCTTTCTTTAGTCTTGACCGGAGGTCTGCGTCGATACGCTTTGGCGCGTAGCTTGATGGACGTTCCCAATGCCCGTAGCTCTCACTCGGTACCTACACCCCGT‑‑GGGGGGGGGGTGGCTATCGTTATGAGTTTCCTACTGGTTTTGCTGATTCTGGCGGCCCTTGGTCAAGTGAGCTGGGACTATACGTGGGCGTTGCTAGGTGCTGGTATTGGCATTGCGGTTCTTGGCTTTCTAGACGATCATGGCC  >  NC_008463/2040136‑2040431

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: