Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 3,670,543 +A coding (543/1020 nt) PA14_RS16675 → ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084633,670,5421.A98.4% 1836.3 / 15.6 446coding (542/1020 nt)PA14_RS16675ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base . (5/2);  new base A (232/204);  total (239/207)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 472 total aligned reads displayed.

GACGTGCTCGCCCGGGTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTTCCTGGAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTCCGCGCCGAGTTCCTCCGCGGG  >  NC_008463/3670391‑3670687
                                                                                                                                                        |                                                                                                                                                  
gACGTGCTCGCCCGGGTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTgg                                                                                                                                                     >  1:4861200/1‑151 (MQ=255)
gACGTGCTCGCCCGGGTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTgg                                                                                                                                                     <  1:3918714/151‑1 (MQ=255)
   gTGCTCGCCCGGGTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGaag                                                                                                                                                 <  1:8148250/151‑3 (MQ=255)
    tGCTCGCCCGGGTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTATCCTGGaagt                                                                                                                                                >  1:9604760/1‑148 (MQ=255)
      cTCGCCCGGGTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGaagtct                                                                                                                                              >  1:2365228/1‑146 (MQ=255)
       tCGCCCGGGTGATTTACGGCTACCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTg                                                                                                                                             >  2:6165078/1‑151 (MQ=255)
        cGCCCGGGTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTgg                                                                                                                                            >  2:178228/1‑151 (MQ=255)
           ccGGGTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCg                                                                                                                                         >  2:4646812/1‑151 (MQ=255)
            cGGGTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCgg                                                                                                                                        >  1:10321421/1‑151 (MQ=255)
             gggTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCggg                                                                                                                                       >  1:919715/1‑151 (MQ=255)
             gggTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCggg                                                                                                                                       <  2:6475815/151‑1 (MQ=255)
              ggagATTTACGGCGTCCGCGAGTCGCTCCTCTTCGCCCTGACCCGGACCGTGCTCAGTCCTCTGTTCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCgggg                                                                                                                                      <  2:10804011/148‑1 (MQ=255)
              ggTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCgggg                                                                                                                                      <  1:8286243/151‑1 (MQ=255)
              ggTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTGTCCTGGAAGTCTGGTCgggg                                                                                                                                      >  2:2461039/1‑151 (MQ=255)
               gTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGc                                                                                                                                     >  1:11291563/1‑151 (MQ=255)
                tGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCt                                                                                                                                    >  2:4004867/1‑151 (MQ=255)
                   tttACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGcc                                                                                                                                 <  1:1078820/151‑1 (MQ=255)
                    ttACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTACCTGGAAGTCTGGTCGGGGCTGCCg                                                                                                                                >  1:4654558/1‑151 (MQ=255)
                     tACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCgg                                                                                                                               <  1:2511753/151‑1 (MQ=255)
                         gCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGct                                                                                                                           >  2:2779819/1‑151 (MQ=255)
                           ttCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGctct                                                                                                                         <  2:4505815/151‑1 (MQ=255)
                           ttCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGctct                                                                                                                         >  2:4861135/1‑151 (MQ=255)
                               gcgTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACct                                                                                                                     <  1:355052/151‑1 (MQ=255)
                               gcgTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACct                                                                                                                     <  1:8959706/151‑1 (MQ=255)
                                      gATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGatc                                                                                                              >  2:4213992/1‑151 (MQ=255)
                                      gATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCCGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGatc                                                                                                              <  2:3146074/151‑1 (MQ=255)
                                       aTCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGatca                                                                                                             <  2:9238182/151‑1 (MQ=255)
                                         ccTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGatcatc                                                                                                           >  2:5737568/1‑151 (MQ=255)
                                         ccTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGatcatc                                                                                                           >  2:5290766/1‑151 (MQ=255)
                                           tctTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCt                                                                                                         >  2:7112082/1‑151 (MQ=255)
                                              tCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGc                                                                                                      <  1:2044826/151‑1 (MQ=255)
                                               cGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGcc                                                                                                     <  1:10903418/151‑1 (MQ=255)
                                                     gaccctgaccGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTc                                                                                               <  2:2897237/151‑1 (MQ=255)
                                                     gaccctgaccGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTc                                                                                               >  2:7039090/1‑151 (MQ=255)
                                                       ccctgaccGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGt                                                                                             <  2:7472632/151‑1 (MQ=255)
                                                        cctgaccGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTc                                                                                            >  1:9759213/1‑151 (MQ=255)
                                                           gaccGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAg                                                                                         >  1:7040942/1‑151 (MQ=255)
                                                             ccGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGcc                                                                                       <  1:1587438/151‑1 (MQ=255)
                                                             ccGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGcc                                                                                       >  2:1264039/1‑151 (MQ=255)
                                                             ccGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGcc                                                                                       <  1:7839848/151‑1 (MQ=255)
                                                              cGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCg                                                                                      <  2:3397437/151‑1 (MQ=255)
                                                              cGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCg                                                                                      <  2:2200456/151‑1 (MQ=255)
                                                               gTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGa                                                                                     <  2:7713232/151‑1 (MQ=255)
                                                               gTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGa                                                                                     >  2:7850265/1‑151 (MQ=255)
                                                                tGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGaa                                                                                    <  1:5479536/151‑1 (MQ=255)
                                                                   tCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACtt                                                                                 >  1:15337/1‑151 (MQ=255)
                                                                   tCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACtt                                                                                 >  1:10586667/1‑151 (MQ=255)
                                                                      gTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCtg                                                                              >  2:6770717/1‑151 (MQ=255)
                                                                         ctcgGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCtggtg                                                                           <  2:5641916/151‑1 (MQ=255)
                                                                          tctGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCtggtgg                                                                          <  2:9619114/151‑1 (MQ=255)
                                                                          tctGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCtggtgg                                                                          <  1:9323654/151‑1 (MQ=255)
                                                                           ctGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGGCGAACTGCTCGGCCAGCGTAGCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCAAGCTTCGTCCAGCCGAACTTCTGGTGgc                                                                         >  2:8462300/1‑151 (MQ=255)
                                                                           ctGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACGACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGgc                                                                         <  2:9336177/151‑1 (MQ=255)
                                                                            tGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGgct                                                                        <  2:5879443/151‑1 (MQ=255)
                                                                            tGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGgct                                                                        <  1:8618096/151‑1 (MQ=255)
                                                                             gATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGgctg                                                                       <  1:421232/151‑1 (MQ=255)
                                                                             gATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGgctg                                                                       <  1:586282/151‑1 (MQ=255)
                                                                             gATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGgctg                                                                       <  1:6380383/151‑1 (MQ=255)
                                                                              aTCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGgctgc                                                                      <  1:4394219/151‑1 (MQ=255)
                                                                                cGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGctgctg                                                                    <  2:7325745/151‑1 (MQ=255)
                                                                                 ggCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTgg                                                                   <  2:10508602/151‑1 (MQ=255)
                                                                                 ggCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTgg                                                                   >  1:1321728/1‑151 (MQ=255)
                                                                                 ggCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTgg                                                                   <  1:3064925/151‑1 (MQ=255)
                                                                                 ggCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTgg                                                                   <  2:2099721/151‑1 (MQ=255)
                                                                                  gCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTggg                                                                  >  1:91079/1‑151 (MQ=255)
                                                                                  gCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTggg                                                                  <  1:3676932/151‑1 (MQ=255)
                                                                                     tCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGcat                                                               >  2:1133884/1‑151 (MQ=255)
                                                                                      cACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGcatc                                                              <  2:5823638/151‑1 (MQ=255)
                                                                                         cGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCAtg                                                           <  2:8783157/151‑1 (MQ=255)
                                                                                         cGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCAtg                                                           >  1:3340834/1‑151 (MQ=255)
                                                                                          gTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCAtgc                                                          >  2:11078474/1‑151 (MQ=255)
                                                                                            cGGCGCCCTGCAAGGCTACGACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCAtgctg                                                        <  2:11567731/151‑1 (MQ=255)
                                                                                              gcgcCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATgctgct                                                      <  1:8059368/151‑1 (MQ=255)
                                                                                              gcgcCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATgctgct                                                      <  1:8792589/151‑1 (MQ=255)
                                                                                               cgcCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGctgctg                                                     <  2:8122179/151‑1 (MQ=255)
                                                                                                gcCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGt                                                    <  1:9081430/151‑1 (MQ=255)
                                                                                                gcCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGt                                                    >  1:11541146/1‑151 (MQ=255)
                                                                                                ccccAGCACGGCTACTACGGCGGCTGGGTCGACCCGCCCGGCCAGCGTTTCCTGGACGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCAGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGt                                                    <  2:2975921/150‑1 (MQ=255)
                                                                                                 cccTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGtt                                                   >  1:501398/1‑151 (MQ=255)
                                                                                                   cTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGACGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGtttt                                                 >  2:10118172/1‑151 (MQ=255)
                                                                                                      cAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCt                                              <  1:1752006/151‑1 (MQ=255)
                                                                                                      cAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCt                                              <  2:136779/151‑1 (MQ=255)
                                                                                                       aaGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTg                                             >  1:9940573/1‑151 (MQ=255)
                                                                                                           ctactaCGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATg                                         >  1:6687479/1‑151 (MQ=255)
                                                                                                            tactacGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGa                                        >  2:1617253/1‑151 (MQ=255)
                                                                                                            tactacGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGa                                        <  1:1238209/151‑1 (MQ=255)
                                                                                                               tacGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGcc                                     <  1:6164614/151‑1 (MQ=255)
                                                                                                                      gcTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGa                              <  1:3254891/151‑1 (MQ=255)
                                                                                                                               aCCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTCcg                     >  1:2611448/1‑151 (MQ=255)
                                                                                                                                  tGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTCCgcgc                  >  2:5112350/1‑151 (MQ=255)
                                                                                                                                   gCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTCCGCGcc                 >  1:6690701/1‑151 (MQ=255)
                                                                                                                                         ccAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTCCGCGCCGAGTtc           >  1:5743591/1‑151 (MQ=255)
                                                                                                                                            agcgtttCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTCCGCGCCGAGTtcct         >  1:10164357/5‑151 (MQ=255)
                                                                                                                                              cgtttCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTCCGCGCCGAGTtcctcc       <  2:5484365/149‑1 (MQ=255)
                                                                                                                                              cgtttCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTCCGCGCCGAGTtcctcc       <  1:5157126/149‑1 (MQ=255)
                                                                                                                                                tttCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTCCGCGCCGAGTTCCTCcgc     <  1:2886402/151‑1 (MQ=255)
                                                                                                                                                 ttCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTCCGCGCCGAGTTCCTCcgcg    <  2:381456/151‑1 (MQ=255)
                                                                                                                                                 ttCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTCCGCGCCGAGTTCCTCcgcg    >  1:7527066/1‑151 (MQ=255)
                                                                                                                                                  tCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTCCGCGCCGAGTTCCTCCGCgg   <  1:5253855/151‑1 (MQ=255)
                                                                                                                                                   ccTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTCCGCGCCGAGTTCCTCCGCggg  <  2:10241403/151‑1 (MQ=255)
                                                                                                                                                        |                                                                                                                                                  
GACGTGCTCGCCCGGGTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTTCCTGGAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTCCGCGCCGAGTTCCTCCGCGGG  >  NC_008463/3670391‑3670687

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