Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 5,028,119 Δ1 bp intergenic (‑41/+96) PA14_RS22975 ← / ← PA14_RS22980 Nramp family divalent metal transporter/ribbon‑helix‑helix domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084635,028,1170C.99.8% 2643.5 / ‑3.5 533intergenic (‑39/+98)PA14_RS22975/PA14_RS22980Nramp family divalent metal transporter/ribbon‑helix‑helix domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base . (215/317);  total (216/317)
Fisher's exact test for biased strand distribution p-value = 4.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 560 total aligned reads displayed.

CGATGGCGGTGGCCCAGTTGCCGGGGTCCATGTAGCCGACCGAGACCAGCAGGCCGGGGCCGACGAACAGCGCCAGCTTGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGGTGGGTTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGC  >  NC_008463/5027967‑5028264
                                                                                                                                                      |                                                                                                                                                   
cGATGGCGGTGGCCCAGTTGCCGGGGTCCATGTAGCCGACCGAGACCAGCAGGCCGGGGCCGACGAACAGCGCCAGCTTGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGc                                                                                                                                                     >  1:10066168/1‑151 (MQ=255)
  aTGGCGGTGGCCCAGTTGCCGGGGTCCATGTAGCCGACCGAGACCAGCAGGCCGGGGCCGACGAACAGCGCCAGCTTGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGccg                                                                                                                                                   <  1:7022316/151‑2 (MQ=255)
   tGGCGGTGGCCCAGTTGCCGGGGTCCATGTAGCCGACCGAGACCAGCAGGCCGGGGCCGACGAACAGCGCCAGCTTGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGccgt                                                                                                                                                  >  2:1395264/1‑149 (MQ=255)
   tGGCGGTGGCCCAGTTGCCGGGGTCCATGTAGCCGACCGAGACCAGCAGGCCGGGGCCGACGAACAGCGCAAGCTTGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGccgt                                                                                                                                                  <  2:5391617/151‑3 (MQ=255)
     gCGGTGGCCCAGTTGCCGGGGTCCATGTAGCCGACCGAGACCAGCAGGCCGGGGCCGACGAACAGCGCCAGCTTGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGccgtgg                                                                                                                                                <  2:552471/151‑5 (MQ=255)
       ggTGGCCCAGTTGCCGGGGTCCATGTAGCCGACCGAGACCAGCAGGCCGGGGCCGACGAACAGCGCCAGCTTGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCgtggtg                                                                                                                                             >  1:4540514/1‑151 (MQ=255)
        gTGGCCCAGTTGCCGGGGTCCATGTAGCCGACCGAGACCAGCAGGCCGGGGCCGACGCACAGCGCCAGCGTGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGACCGGCGCGGTGGCCGtggtgg                                                                                                                                            <  2:5203750/151‑1 (MQ=255)
           gCCCAGTTGCCGGGGTCCATGTAGCCGACCGAGACCAGCAGGCCGGGGCCGACGAACAGCGCCAGCTTGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCa                                                                                                                                         <  1:998796/151‑1 (MQ=255)
            cccAGTTGCCGGGGTCCATGTAGCCGACCGAGACCAGCAGGCCGGGGCCGACGAACAGCGCCAGCTTGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAg                                                                                                                                        <  2:5637970/151‑1 (MQ=255)
             ccAGTTGCCGGGGTCCATGTAGCCGACCGAGACCAGCAGGCCGGGGCCGACGAACAGCGCCAGCTTGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAgg                                                                                                                                       <  2:11111923/151‑1 (MQ=255)
              cAGTTGCCGGGGTCCATGTAGCCGACCGAGACCAGCAGGCCGGGGCCGACGAACAGCGCCAGCTTGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAggg                                                                                                                                      >  2:4658690/1‑151 (MQ=255)
               aGTTGCCGGGGTCCATGTAGCCGACCGAGACCAGCAGGCCGGGGCCGACGAACAGCGCCAGCTTGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGa                                                                                                                                     <  1:11104029/151‑1 (MQ=255)
                  tGCCGGGGTCCATGTAGCCGACCGAGACTAGCAGGCCGGGGCCGACGAACAGCGCCAGCTTGCGCAGCGGCGAGGCGTCGGCCGGCTCCGCCACGCTGTGGCTGACCGCCGAGGGAAAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTa                                                                                                                                  <  2:4525575/151‑1 (MQ=255)
                   gCCGGGGTCCATGTAGCCGACCGAGACCAGCAGGCCGGGGCCGACGAACAGCGCCAGCTTGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTAt                                                                                                                                 <  1:3919599/151‑1 (MQ=255)
                          tCCATGTAGCCGACCGAGACCAGCAGGCCGGGGCCGACGAACAGCGCCAGCTTGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGa                                                                                                                          <  1:4768648/151‑1 (MQ=255)
                           ccATGTAGCCGACCGAGACCAGCAGGCCGGGGCCGACGAACAGCGCCAGCTTGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAg                                                                                                                         <  2:5081730/151‑1 (MQ=255)
                           ccATGTAGCCGACCGAGACCAGCAGGCCGGGGCCGACGAACAGCGCCAGCTTGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAg                                                                                                                         <  2:2311497/151‑1 (MQ=255)
                           ccATGTAGCCGACCGAGACCAGCAGGCCGGGGCCGACGAACAGCGCCAGCTTGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAg                                                                                                                         <  1:6257498/151‑1 (MQ=255)
                            cATGTAGCCGACCGAGACCAGCAGGCCGGGGCCGACGAACAGCGCCAGCTTGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAgg                                                                                                                        <  1:416312/151‑1 (MQ=255)
                             aTGTAGCCGACCGAGACCAGCAGGCCGGGGCCGACGAACAGCGCCAGCTTGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGc                                                                                                                       <  2:2455700/151‑1 (MQ=255)
                               gTAGCCGACCGAGACCAGCAGGCCGGGGCCGACGAACAGCGCCAGCTTGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGa                                                                                                                     <  1:5115333/151‑1 (MQ=255)
                                tAGCCGACCGAGACCAGCAGGCCGGGGCCGACGAACAGCGCCAGCTTGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGAc                                                                                                                    <  1:7857133/151‑1 (MQ=255)
                                 aGCCGACCGAGACCAGCAGGCCGGGGCCGACGAACAGCGCCAGCTTGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACg                                                                                                                   >  2:4839821/1‑151 (MQ=255)
                                 aGCCGACCGAGACCAGCAGGCCGGGGCCGACGAACAGCGCCAGCTTGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACg                                                                                                                   >  1:6845846/1‑151 (MQ=255)
                                   ccgaccgaGACCAGCAGGCCGGGGCCGACGAACAGCGCCAGCTTGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGcc                                                                                                                 >  2:7835788/1‑151 (MQ=255)
                                     gaccgaGACCAGCAGGCCGGGGCCGACGAACAGCGCCAGCTTGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCTGGGAGTATTTCACGAGGCGACGCCgg                                                                                                               >  2:9917831/1‑151 (MQ=255)
                                     gaccgaGACCAGCAGGCCGGGGCCGACGAACAGCGCCAGCTTGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCgg                                                                                                               >  2:8618492/1‑151 (MQ=255)
                                      accgaGACCAGCAGGCCGGGGCCGACGAACAGCGCCAGCTTGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGa                                                                                                              <  2:703682/151‑1 (MQ=255)
                                      accgaGACCAGCAGGCCGGGGCCGACGAACAGCGCCAGCTTGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGa                                                                                                              <  1:7122185/151‑1 (MQ=255)
                                       ccgaGACCAGCAGGCCGGGGCCGACGAACAGCGCCAGCTTGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGAt                                                                                                             <  1:8971002/151‑1 (MQ=255)
                                        cgaGACCAGCAGGCCGGGGCCGACGAACAGCGCCCCCTTGCGCAGCGGCGAGGCGTCGCCCGGCACCGCCACCCAGTGCCTGACCGCCCGCGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGAAAGACACGCGGCGACGCCGGaac                                                                                                            >  2:10612510/1‑149 (MQ=255)
                                        cgaGACCAGCAGGCCGGGGCCGACGAACAGCGCCAGCTTGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATc                                                                                                            <  2:9686396/151‑1 (MQ=255)
                                         gagaCCAGCAGGCCGGGGCCGACGAACAGCGCCAGCTTGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCg                                                                                                           <  2:2985931/151‑1 (MQ=255)
                                         gagaCCAGCAGGCCGGGGCCGACGAACAGCGCCAGCTTGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCg                                                                                                           <  1:10561318/151‑1 (MQ=255)
                                         gagaCCAGCAGGCCGGGGCCGACGAACAGCGCCAGCTTGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCg                                                                                                           <  1:6646868/151‑1 (MQ=255)
                                          agaCCAGCAGGCCGGGGCCGACGAACAGCGCCAGCTTGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCgg                                                                                                          <  2:2784550/151‑1 (MQ=255)
                                             ccAGCAGGCCGGGGCCGACGAACAGCGCCAGCTTGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGgga                                                                                                       <  1:7277810/151‑1 (MQ=255)
                                             ccAGCAGGCCGGGGCCGACGAACAGCGCCAGCTTGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGgga                                                                                                       <  2:10230340/151‑1 (MQ=255)
                                              cagcagGCCGGGGCCGACGAACAGCGCCAGCTTGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGggag                                                                                                      <  2:9048307/151‑1 (MQ=255)
                                              cagcagGCCGGGGCCGACGAACAGCGCCAGCTTGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGggag                                                                                                      <  2:3468099/151‑1 (MQ=255)
                                               agcagGCCGGGGCCGACGAACAGCGCCAGCTTGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGggagg                                                                                                     >  2:8009139/1‑151 (MQ=255)
                                               agcagGCCGGGGCCGACGAACAGCGCCAGCTTGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGggagg                                                                                                     <  2:10310278/151‑1 (MQ=255)
                                                 cagGCCGGGGCCGACGAACAGCGCCAGCTTGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGaa                                                                                                   >  1:10961620/1‑151 (MQ=255)
                                                 cagGCCGGGGCCGACGAACAGCGCCAGCTTGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGaa                                                                                                   <  1:9735179/151‑1 (MQ=255)
                                                   ggCCGGGGCCGACGAACAGCGCCAGCTTGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGt                                                                                                 >  2:6169093/1‑151 (MQ=255)
                                                       ggggCCGACGAACAGCGCCAGCTTGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCCGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGa                                                                                             >  1:3713399/1‑151 (MQ=255)
                                                       ggggCCGACGAACAGCGCCAGCTTGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGa                                                                                             >  2:1630061/1‑151 (MQ=255)
                                                       ggggCCGACGAACAGCGCCAGCTTGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGa                                                                                             >  1:7495826/1‑151 (MQ=255)
                                                           ccGACGAACAGCGCCAGCTTGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAg                                                                                         <  2:5861007/151‑1 (MQ=255)
                                                           ccGACGAACAGCGCCAGCTTGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAg                                                                                         >  2:193321/1‑151 (MQ=255)
                                                           ccGACGAACAGCGCCAGCTTGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAg                                                                                         <  1:5363438/151‑1 (MQ=255)
                                                           ccGACGAACAGCGCCAGCTTGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAg                                                                                         <  2:814481/151‑1 (MQ=255)
                                                            cgacgaACAGCGCCAGCTTGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGc                                                                                        >  1:6513879/1‑151 (MQ=255)
                                                              acgaACAGCGCCAGCTTGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTAGTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCt                                                                                      >  1:3002058/1‑151 (MQ=255)
                                                              aagaACAGCGACAGCTTGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGCGCAGGGTGTAATTCACGTTGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCt                                                                                      >  2:6324658/3‑151 (MQ=255)
                                                               cgaACAGCGCCAGCTTGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTa                                                                                     >  2:4057502/1‑151 (MQ=255)
                                                                     gcgcCAGCTTGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGtt                                                                               <  1:2365960/151‑1 (MQ=255)
                                                                      cgcCAGCTTGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGttt                                                                              <  1:193349/151‑1 (MQ=255)
                                                                       gcCAGCTTGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTc                                                                             >  1:3895702/1‑151 (MQ=255)
                                                                       gcCAGCTTGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCGGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGCCCTGCCTATTCGtttg                                                                             >  1:10502307/1‑150 (MQ=255)
                                                                          aGCTTGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCAcc                                                                          >  2:10488624/1‑151 (MQ=255)
                                                                             gtGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCTCCCCg                                                                       >  1:408181/2‑151 (MQ=255)
                                                                              tGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCgg                                                                      >  1:3226147/1‑151 (MQ=255)
                                                                               gcgcAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCggt                                                                     <  2:2598800/151‑2 (MQ=255)
                                                                                cgcAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCggtc                                                                    >  1:10295939/1‑149 (MQ=255)
                                                                                cgcAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCggtc                                                                    >  1:8271940/1‑149 (MQ=255)
                                                                                 gcAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCggtcg                                                                   >  2:9482458/1‑148 (MQ=255)
                                                                                 gcAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCggtcg                                                                   <  2:3735362/151‑4 (MQ=255)
                                                                                 gcAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCggtcg                                                                   >  1:1741971/1‑148 (MQ=255)
                                                                                   aGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCCGGGAGTATTTCTCGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCGAAACGTTTCACCCCggtcgcg                                                                 >  2:3261943/1‑146 (MQ=255)
                                                                                                cggccggcACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTccctcccct                                                   >  1:2599606/1‑151 (MQ=255)
                                                                                                 ggccggcACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTccctcccctc                                                  <  1:7196008/151‑1 (MQ=255)
                                                                                                  gccggcACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCa                                                 <  1:5701178/151‑1 (MQ=255)
                                                                                                   ccggcACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAg                                                <  2:10955494/151‑1 (MQ=255)
                                                                                                   ccggcACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAg                                                >  1:147527/1‑151 (MQ=255)
                                                                                                    cggcACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAgg                                               <  1:3527097/151‑1 (MQ=255)
                                                                                                    cggcACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAgg                                               <  2:882710/151‑1 (MQ=255)
                                                                                                           gCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGgc                                        <  1:9465184/151‑1 (MQ=255)
                                                                                                           gCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGgc                                        <  2:6482318/151‑1 (MQ=255)
                                                                                                            ccACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGgcg                                       >  2:9844810/1‑151 (MQ=255)
                                                                                                              aCGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGgcgcg                                     >  1:10633958/1‑151 (MQ=255)
                                                                                                              aCGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGgcgcg                                     <  2:6878834/151‑1 (MQ=255)
                                                                                                               cGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGgcgcgc                                    <  1:8174711/151‑1 (MQ=255)
                                                                                                                gCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCt                                   <  1:555964/151‑1 (MQ=255)
                                                                                                                 cTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTg                                  <  2:2851773/151‑1 (MQ=255)
                                                                                                                  tgtgGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTgg                                 >  2:8492732/1‑151 (MQ=255)
                                                                                                                  tgtgGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCAACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTgg                                 >  1:10818726/1‑151 (MQ=255)
                                                                                                                   gtgGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGc                                >  2:8407962/1‑151 (MQ=255)
                                                                                                                    tgGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGcc                               <  1:9162544/151‑1 (MQ=255)
                                                                                                                    tgGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGcc                               >  1:7925885/1‑151 (MQ=255)
                                                                                                                    tgGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGcc                               >  1:2198780/1‑151 (MQ=255)
                                                                                                                    tgGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGGCTTCCCTCCCCGCAGGCCTGGGCGCGCGGGcc                               >  1:6076900/1‑151 (MQ=255)
                                                                                                                                 aGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGc                  >  2:10861749/1‑151 (MQ=255)
                                                                                                                                  gggACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCt                 <  2:8936505/151‑1 (MQ=255)
                                                                                                                                  gggACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCt                 <  2:734929/151‑1 (MQ=255)
                                                                                                                                  gggACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCt                 <  1:4870/151‑1 (MQ=255)
                                                                                                                                     aCAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCgg              >  2:7676801/1‑151 (MQ=255)
                                                                                                                                          aCGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGAcc         <  2:9221569/151‑1 (MQ=255)
                                                                                                                                               gcgGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTtctc    <  1:3069962/151‑1 (MQ=255)
                                                                                                                                               gcgGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTtctc    >  1:4306964/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
CGATGGCGGTGGCCCAGTTGCCGGGGTCCATGTAGCCGACCGAGACCAGCAGGCCGGGGCCGACGAACAGCGCCAGCTTGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGGTGGGTTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGC  >  NC_008463/5027967‑5028264

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: