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.A99.2% 1676.4 / 1.8 402coding (542/1020 nt)PA14_RS16675ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base . (2/1);  new base A (183/216);  total (185/217)
Fisher's exact test for biased strand distribution p-value = 5.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 422 total aligned reads displayed.

GACGTGCTCGCCCGGGTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTTCCTGGAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTCCGCGCCGAGTTCCTCCGCGGG  >  NC_008463/3670391‑3670687
                                                                                                                                                        |                                                                                                                                                 
gACGTGCTCGCCCGGGTGATTTACGGTTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGGTTCCTgg                                                                                                                                                    >  1:10446074/1‑151 (MQ=255)
gACGTGCTCGCCCGGGTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTgg                                                                                                                                                    >  1:2793997/1‑151 (MQ=255)
 aCGTGCTCGCCCGGGTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGa                                                                                                                                                  >  1:11028597/1‑151 (MQ=255)
 aCGTGCTCGACCGGGTGATTTACTGCTTCCGCGTGTCGATCCTCTTCGCCCTGAACCTGACCGAGCTCAGTTATCTGATCTGCGTCACCGTCGGCGCCCTGCAAGGATAATACTGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGa                                                                                                                                                  >  2:6141291/1‑151 (MQ=255)
  cGTGCTCGCCCGGGTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGaa                                                                                                                                                 <  2:3492535/151‑2 (MQ=255)
   gTGCTCGCCCGGGTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCAACTACGGCGGCACGGTCGACCTGCTCGGCCAGCGTTTCCTGGaag                                                                                                                                                >  2:3295847/1‑149 (MQ=255)
    tGCTCGCCCGGGTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGaagt                                                                                                                                               >  1:10505521/1‑148 (MQ=255)
      cTCGCCCGGGTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGaagtct                                                                                                                                             >  2:4949351/1‑146 (MQ=255)
           ccGGGTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCg                                                                                                                                         <  1:3315716/151‑1 (MQ=255)
           ccGGGTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCg                                                                                                                                         >  2:5998675/1‑151 (MQ=255)
             gggTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCggg                                                                                                                                       >  1:1469042/1‑151 (MQ=255)
              ggTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCgggg                                                                                                                                      <  2:3945386/151‑1 (MQ=255)
              ggTGATTTACCGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCgggg                                                                                                                                      >  2:8363566/1‑151 (MQ=255)
               gTGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCCAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGc                                                                                                                                     >  2:1536574/1‑151 (MQ=255)
                tGATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCt                                                                                                                                    <  1:2574645/151‑1 (MQ=255)
                 gATTTACGGCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTg                                                                                                                                   <  1:7676694/151‑1 (MQ=255)
                         gCTTCCGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGct                                                                                                                           >  2:1937323/1‑151 (MQ=255)
                             ccGCGTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCACGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTAc                                                                                                                       >  2:119491/1‑151 (MQ=255)
                               gcgTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACct                                                                                                                     <  1:10015023/151‑1 (MQ=255)
                                cgTGTCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACctg                                                                                                                    <  1:1552546/151‑1 (MQ=255)
                                 gtgtCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACctgc                                                                                                                   >  1:314961/1‑151 (MQ=255)
                                 gtgtCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACctgc                                                                                                                   >  2:1718125/1‑151 (MQ=255)
                                  tgtCGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACctgct                                                                                                                  >  2:3217898/1‑151 (MQ=255)
                                     cGATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGat                                                                                                               <  2:844588/151‑1 (MQ=255)
                                      gctCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGatc                                                                                                              <  1:5296460/149‑1 (MQ=255)
                                      gATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGatc                                                                                                              <  1:10212750/151‑1 (MQ=255)
                                      gATCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGatc                                                                                                              <  1:4235175/151‑1 (MQ=255)
                                       aTCCTCTTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGatca                                                                                                             >  1:8734144/1‑151 (MQ=255)
                                          ctctTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATcc                                                                                                          >  1:7402112/1‑151 (MQ=255)
                                          ctctTCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATcc                                                                                                          >  1:6231625/1‑151 (MQ=255)
                                             ttCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTgg                                                                                                       <  2:9698148/151‑1 (MQ=255)
                                              tCGCCCTGACCCTGACCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCAGCCTGGc                                                                                                      >  2:10507250/1‑151 (MQ=255)
                                                 ccctgaccctgaCCGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAg                                                                                                   >  2:3397599/1‑151 (MQ=255)
                                                 ccctgaccctgaCCGTGCTCAGTTCACTGATCGGCGTCTCCGTCGGAGCCCTGCAAGGATACTACGTCGTCTGGGTCGACCTGATCTGCCAGAGTTTCCTGGAAGTCTGGTCGGGGCTGCCTGTGCTCTACCTGCTGATCATCCTGGCCAg                                                                                                   <  2:6737372/151‑1 (MQ=255)
                                                    tgaccctgaccGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCtt                                                                                                <  1:7070948/151‑1 (MQ=255)
                                                    tgaccctgaccGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCtt                                                                                                <  2:3682810/151‑1 (MQ=255)
                                                     gaccctgaccGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTc                                                                                               <  1:5967070/151‑1 (MQ=255)
                                                     gaccctgaccGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTc                                                                                               >  1:5549299/1‑151 (MQ=255)
                                                       ccctgaccGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGt                                                                                             <  1:8518497/151‑1 (MQ=255)
                                                         ctgaccGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCGCGGCAAGCTTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGTTCATCCTGGCCAGCTTCGTcc                                                                                           >  2:2657224/1‑151 (MQ=255)
                                                          tgaccGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCa                                                                                          >  2:4314840/1‑151 (MQ=255)
                                                          tgaccGTGCTCAGTTCTCTGACCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCa                                                                                          >  1:4705604/1‑151 (MQ=255)
                                                           gaccGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAg                                                                                         >  1:637947/1‑151 (MQ=255)
                                                            accGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGc                                                                                        >  2:2959295/1‑151 (MQ=255)
                                                             ccGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGcc                                                                                       >  1:3441332/1‑151 (MQ=255)
                                                             ccGTGCTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGcc                                                                                       >  1:4908161/1‑151 (MQ=255)
                                                                  cTCAGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACt                                                                                  >  1:3679959/1‑151 (MQ=255)
                                                                     aGTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCGTCGTCCAGCCGAACTTCt                                                                               >  2:3049493/1‑151 (MQ=255)
                                                                     aGTTCTCTGATCGGCGTCACCGTCGGCGCCCAGCAAGGCTACTACGGCGGCTGGGTCGACCCGCTCGGCCAGCGTGTCCTCGAAGTCTGGTCGGGGCTGCCGGCGCACGACCTGCTGATCATCCTGGCCACCATCGTCCAGCCGAACTTCt                                                                               <  2:509368/151‑1 (MQ=255)
                                                                      gTTCTCTGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCtg                                                                              <  1:5668156/151‑1 (MQ=255)
                                                                            tGATCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGgct                                                                        <  2:7299190/151‑1 (MQ=255)
                                                                              aTCGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGgctgc                                                                      >  2:8640886/1‑151 (MQ=255)
                                                                                cGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGctgctg                                                                    >  2:6121702/1‑151 (MQ=255)
                                                                                cGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGctgctg                                                                    <  1:3582442/151‑1 (MQ=255)
                                                                                cGGCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGatgatg                                                                    >  2:1013281/1‑151 (MQ=255)
                                                                                 ggCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTgg                                                                   >  1:9628297/1‑151 (MQ=255)
                                                                                  gCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTggg                                                                  <  2:8226639/151‑1 (MQ=255)
                                                                                  gCGTCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTggg                                                                  <  1:2909445/151‑1 (MQ=255)
                                                                                     tCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGcat                                                               <  2:2781112/151‑1 (MQ=255)
                                                                                     tCACCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGcat                                                               <  2:2186377/151‑1 (MQ=255)
                                                                                       aCCGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGcatca                                                             >  1:8027178/1‑151 (MQ=255)
                                                                                        ccGTCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGcatcat                                                            <  2:10988070/151‑1 (MQ=255)
                                                                                           tCGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCAtgct                                                         >  1:5835002/1‑151 (MQ=255)
                                                                                            cGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCAtgctg                                                        <  2:1270663/151‑1 (MQ=255)
                                                                                            cGGCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCAtgctg                                                        >  2:958730/1‑151 (MQ=255)
                                                                                             ggCGCCCTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCAtgctgc                                                       >  1:1085809/1‑151 (MQ=255)
                                                                                                  ccTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGttt                                                  <  2:9061984/151‑1 (MQ=255)
                                                                                                  ccTGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGttt                                                  <  1:1987065/151‑1 (MQ=255)
                                                                                                    tGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTc                                                <  2:7959552/151‑1 (MQ=255)
                                                                                                    tGCAAGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTc                                                <  2:6476947/151‑1 (MQ=255)
                                                                                                        aGGCTACTTCGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTgg                                            <  1:10627763/151‑1 (MQ=255)
                                                                                                        aGGCTACTACGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTgg                                            >  2:9501703/1‑151 (MQ=255)
                                                                                                             actacGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAg                                       >  1:6926051/1‑151 (MQ=255)
                                                                                                             actacGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAg                                       >  1:11141528/1‑151 (MQ=255)
                                                                                                               tacGGCGGCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGcc                                     >  1:5010192/1‑151 (MQ=255)
                                                                                                                 cggcggCTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCtg                                   <  1:2319432/151‑1 (MQ=255)
                                                                                                                  ggcggcTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCtgg                                  >  1:8264724/1‑151 (MQ=255)
                                                                                                                   gcggcTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCtggt                                 <  1:9305431/151‑1 (MQ=255)
                                                                                                                   gcggcTGGGTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCtggt                                 <  1:1283030/151‑1 (MQ=255)
                                                                                                                         gggTCGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGt                           <  2:576074/151‑1 (MQ=255)
                                                                                                                             cGACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTc                       >  1:6082863/1‑151 (MQ=255)
                                                                                                                              gACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTcc                      >  2:5273272/1‑151 (MQ=255)
                                                                                                                              gACCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTcc                      >  1:9926205/1‑151 (MQ=255)
                                                                                                                               aCCTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTCcg                     >  1:6542377/1‑151 (MQ=255)
                                                                                                                                ccTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTCcgc                    >  2:4081480/1‑151 (MQ=255)
                                                                                                                                 cTGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTCcgcg                   >  1:8921358/1‑151 (MQ=255)
                                                                                                                                  tGCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTCCgcgc                  >  2:4792456/1‑151 (MQ=255)
                                                                                                                                   gCTCGGCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCCGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTCCGCGcc                 <  1:6008441/151‑1 (MQ=255)
                                                                                                                                     tCGGCCAGCGTGTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTCCGCGCCGa               <  1:10795557/151‑1 (MQ=255)
                                                                                                                                       ggCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTCCGCGCCGAGt             >  1:9038084/1‑151 (MQ=255)
                                                                                                                                       ggCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTCCGCGCCGAGt             <  1:9439948/151‑1 (MQ=255)
                                                                                                                                       ggCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTCCGCGCCGAGt             <  2:10717903/151‑1 (MQ=255)
                                                                                                                                        gCCAGCGTTTCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTCCGCGCCGAGtt            >  1:3928702/1‑151 (MQ=255)
                                                                                                                                             gcgtttCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGACCTTCTGGTGGCTGCTGGGCCTCTGGCTGCTGTGTGCCTGGAGGAGCCTGGTGTACGTGGTCCGCGCCGAGTtcctc        >  2:4219851/4‑151 (MQ=255)
                                                                                                                                              cgtttCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTCCGCGCCGAGTtcctcc       >  1:7813150/3‑151 (MQ=255)
                                                                                                                                               gtttCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTCCGCGCCGAGTTCCTCcg      <  1:2730104/150‑1 (MQ=255)
                                                                                                                                                tttCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCGGTTTTCCTGGATGAGCCTGGTGGACGTGGTCCGCGCCGAGTTCCTCcgc     >  2:6356943/1‑151 (MQ=255)
                                                                                                                                                 ttCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTCCGCGCCGAGTTCCTCcgcg    >  2:6959691/1‑151 (MQ=255)
                                                                                                                                                 ttCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTCCGCGCCGAGTTCCTCcgcg    >  2:4411189/1‑151 (MQ=255)
                                                                                                                                                  tCCTGGAAGTCTGGTCGGGGCTGCCGGTGCTCTACCTGCTGATCATCCTGGCCAGCTTCGTCCAGCCGAACTTCTGGTGGCTGCTGGGCATCATGCTGCTGTTTTCCTGGATGAGCCTGGTGGACGTGGTCCGCGCCGAGTTCCTCCGCgg   <  1:11001173/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: