Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 5,028,201 2 bp→CG intergenic (‑123/+13) 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,2010GC99.1% 1669.1 / 3.1 430intergenic (‑123/+14)PA14_RS22975/PA14_RS22980Nramp family divalent metal transporter/ribbon‑helix‑helix domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base G (4/0);  new base C (212/214);  total (216/214)
Fisher's exact test for biased strand distribution p-value = 1.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01
*NC_0084635,028,2020TG99.1% 1692.3 / 4.4 430intergenic (‑124/+13)PA14_RS22975/PA14_RS22980Nramp family divalent metal transporter/ribbon‑helix‑helix domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base G (211/214);  minor base . (4/0);  total (216/214)
Fisher's exact test for biased strand distribution p-value = 1.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 430 total aligned reads displayed.

CGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGGTGGGTTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTCGCCGCGATAGTCC  >  NC_008463/5028060‑5028338
                                                                                                                                             ||                                                                                                                                        
cGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTtccctc                                                                                                                                 <  1:7210854/150‑1 (MQ=255)
cGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTtccctc                                                                                                                                 >  1:2157418/1‑150 (MQ=255)
 gTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGTCAGAACGGCGCGGTGGCCGGGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCCGGCTAAACGTGGCACCCCGGTCGCGTCTtccctcc                                                                                                                                >  2:771590/1‑150 (MQ=255)
 gTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTtccctcc                                                                                                                                >  2:7565687/1‑150 (MQ=255)
  tCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTtccctccc                                                                                                                               >  2:6317998/1‑150 (MQ=255)
    ggccggcACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTccctcccct                                                                                                                             >  1:7419990/1‑150 (MQ=255)
      ccggcACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCa                                                                                                                           <  2:1539798/150‑1 (MQ=255)
        ggcACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAgg                                                                                                                         <  1:3415275/150‑1 (MQ=255)
         gcACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGACGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGc                                                                                                                        >  2:2263398/1‑150 (MQ=255)
          cACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGcc                                                                                                                       <  2:8188325/150‑1 (MQ=255)
             cGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTgg                                                                                                                    <  2:9643884/150‑1 (MQ=255)
             cGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTgg                                                                                                                    <  1:10039318/150‑1 (MQ=255)
              gCCACGCTGTGGCTGACCGCCGTGGGACAGAACGGCGCGGTGGCCGTGGTGGGCTGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGAGACCAGCCTTAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTggg                                                                                                                   <  2:3111334/150‑1 (MQ=255)
              gCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTggg                                                                                                                   >  2:5912944/1‑150 (MQ=255)
                cACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGgcg                                                                                                                 <  2:4795322/150‑1 (MQ=255)
                   gCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGgcgcgc                                                                                                              <  2:678055/150‑1 (MQ=255)
                   gCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGgcgcgc                                                                                                              <  1:10631870/150‑1 (MQ=255)
                    cTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCATGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCt                                                                                                             <  2:2744289/150‑1 (MQ=255)
                     tgtgGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTg                                                                                                            <  2:4771319/150‑1 (MQ=255)
                      gtgGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGACGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTgg                                                                                                           >  1:2747894/1‑150 (MQ=255)
                      gtgGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTgg                                                                                                           >  2:8245316/1‑150 (MQ=255)
                        ggCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGcc                                                                                                         <  2:7235833/150‑1 (MQ=255)
                         gCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCg                                                                                                        >  1:8346536/1‑150 (MQ=255)
                              ccgccgAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAg                                                                                                   <  1:5458848/150‑1 (MQ=255)
                               cgccgAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAgg                                                                                                  <  1:6877527/150‑1 (MQ=255)
                                gccgAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGc                                                                                                 <  1:10500573/150‑1 (MQ=255)
                                 ccgAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTACACGTATCACCCCGGTCGCGTCTTCCCTCCGCTCAGGCCTGGGCGCGCTGGCCGACCAGGCt                                                                                                >  2:6857975/1‑150 (MQ=255)
                                 ccgAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCt                                                                                                <  1:16592/150‑1 (MQ=255)
                                  cgAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGGTGCACCCCGGTCGCGTCTTCCCTCCCCTCTGGCCTGGGCGCTGTCGCCGTGCAGGCTc                                                                                               >  2:4085228/1‑150 (MQ=255)
                                    aGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAg                                                                                             >  1:1501619/1‑150 (MQ=255)
                                     gggACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGc                                                                                            <  2:1187432/150‑1 (MQ=255)
                                      ggACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCCGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCGCAGcg                                                                                           >  2:9142613/1‑149 (MQ=255)
                                      ggACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCt                                                                                           <  2:4210796/150‑1 (MQ=255)
                                       gACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTc                                                                                          >  2:17929/1‑150 (MQ=255)
                                         cAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCgg                                                                                        <  1:8712123/150‑1 (MQ=255)
                                          aGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGc                                                                                       >  2:10282313/1‑150 (MQ=255)
                                           gAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTGGGCg                                                                                      >  2:7969979/1‑150 (MQ=255)
                                            aaCGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGa                                                                                     <  2:453441/150‑1 (MQ=255)
                                              cGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGAcc                                                                                   >  1:6965828/1‑150 (MQ=255)
                                              cGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGAcc                                                                                   <  1:9210216/150‑1 (MQ=255)
                                              cGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCCTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGAcc                                                                                   >  1:9615335/1‑150 (MQ=255)
                                               ggCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCt                                                                                  <  1:7096451/150‑1 (MQ=255)
                                                   cgGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTtctc                                                                              >  1:9444058/1‑150 (MQ=255)
                                                    ggTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTcg                                                                             >  2:5106596/1‑150 (MQ=255)
                                                      gtggccGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTcgc                                                                            >  1:7932311/5‑150 (MQ=255)
                                                      gtggccGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTcgc                                                                            >  2:9034708/5‑150 (MQ=255)
                                                      gtggccGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTcgc                                                                            <  2:8564489/146‑1 (MQ=255)
                                                         gccGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGcgcc                                                                         <  1:6317998/149‑1 (MQ=255)
                                                          ccGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGcgccg                                                                        >  1:5905254/1‑150 (MQ=255)
                                                                 gggCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGcc                                                                 <  1:4499551/150‑1 (MQ=255)
                                                                   gCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAg                                                               >  2:8713759/1‑150 (MQ=255)
                                                                   gCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAg                                                               >  2:4401126/1‑150 (MQ=255)
                                                                     aGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGt                                                             <  1:3982040/150‑1 (MQ=255)
                                                                     aGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGt                                                             <  1:3066080/150‑1 (MQ=255)
                                                                     aGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGt                                                             <  1:1966017/150‑1 (MQ=255)
                                                                      gggAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTa                                                            >  1:9932422/1‑150 (MQ=255)
                                                                      gggAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTa                                                            >  1:175535/1‑150 (MQ=255)
                                                                       ggAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAg                                                           <  1:9777323/150‑1 (MQ=255)
                                                                       ggAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAg                                                           <  2:4943671/150‑1 (MQ=255)
                                                                         aGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCg                                                         >  2:3599478/1‑150 (MQ=255)
                                                                              ttCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGc                                                    >  1:7847729/1‑150 (MQ=255)
                                                                               tCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCGCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCa                                                   <  2:8474188/150‑1 (MQ=255)
                                                                                   gAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTg                                               <  1:4565219/150‑1 (MQ=255)
                                                                                     ggCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGac                                             >  1:1292193/1‑150 (MQ=255)
                                                                                      gCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGaca                                            <  2:7006437/150‑1 (MQ=255)
                                                                                        gACGCCGGATCGGGGAGGACGTGCGTCCAGCCGAATCGTTTCACCCCGGTCGCGTCGTCCCTCCCCTCAGGCGTGGGCGCGCTGGCCGACCAGGCTCAGCGCGGCGGGCGTCTCGCGCCGCTGCGCCGGGTAGCGGGTGCGAGTGACCCg                                          >  2:10663392/1‑150 (MQ=255)
                                                                                         aCGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACGGACACGc                                         >  1:4802255/1‑150 (MQ=255)
                                                                                          cGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCa                                        <  1:3545642/150‑1 (MQ=255)
                                                                                           gCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCag                                       <  2:3292857/150‑1 (MQ=255)
                                                                                                 tCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCaggaaggt                                 <  1:251215/150‑2 (MQ=255)
                                                                                                    gggAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGc                              >  1:5381656/1‑150 (MQ=255)
                                                                                                    gggAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGc                              >  2:6039785/1‑150 (MQ=255)
                                                                                                     ggaggaAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCCGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGGGCAACTGACACGCAGGCAGGAGGCg                             >  2:7377489/1‑150 (MQ=255)
                                                                                                     ggaggaAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCg                             >  2:3771559/1‑150 (MQ=255)
                                                                                                       aggaAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGaa                           >  1:5291278/1‑150 (MQ=255)
                                                                                                          aaGTGCGACCAGCCAAAACGTTTCACCCCGGTCGCGTCTTCCCTCCGCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGtt                        <  2:2901941/150‑1 (MQ=255)
                                                                                                           aGTGCGACCAGCCTAACCGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTg                       <  2:715233/150‑1 (MQ=255)
                                                                                                           aGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTg                       <  2:3387034/150‑1 (MQ=255)
                                                                                                           aGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTg                       <  1:1258600/150‑1 (MQ=255)
                                                                                                           aGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTg                       >  1:7599964/1‑150 (MQ=255)
                                                                                                           aGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTg                       <  1:8869787/150‑1 (MQ=255)
                                                                                                              gCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGc                    >  1:5592937/1‑150 (MQ=255)
                                                                                                                  ccAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCt                <  1:5942680/150‑1 (MQ=255)
                                                                                                                  ccAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCt                >  2:9399379/1‑150 (MQ=255)
                                                                                                                  ccAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAAGAGGCGAAGTTGAGCACCt                <  1:4710153/150‑1 (MQ=255)
                                                                                                                    aGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTcg              >  2:5411159/1‑150 (MQ=255)
                                                                                                                      ccTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTcgcc            >  1:6633275/1‑150 (MQ=255)
                                                                                                                       cTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCTGGAAGGAGGCGAAGTTGAGCACCTcgccg           >  2:1272756/1‑150 (MQ=255)
                                                                                                                        tAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTcgccgc          >  1:10327868/1‑150 (MQ=255)
                                                                                                                        tAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTcgccgc          <  1:2390687/150‑1 (MQ=255)
                                                                                                                        tAAACGTGTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCACGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTcgccgc          <  2:6555626/150‑1 (MQ=255)
                                                                                                                           aCGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTCGCCGCGAt       >  2:2100974/1‑150 (MQ=255)
                                                                                                                           aCGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTCGCCGCGAt       >  2:6442905/1‑150 (MQ=255)
                                                                                                                            cGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTCGCCGCGATa      >  2:5376964/1‑150 (MQ=255)
                                                                                                                             gTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTCGCCGCGATAg     <  2:8759564/150‑1 (MQ=255)
                                                                                                                             gTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTCGCCGCGATAg     >  1:8251624/1‑150 (MQ=255)
                                                                                                                              tttCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTCGCCGCGATAGt    <  2:952582/150‑1 (MQ=255)
                                                                                                                               ttCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTCGCCGCGATAGTc   >  1:4200155/1‑150 (MQ=255)
                                                                                                                               ttCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTCGCCGCGATAGTc   <  2:5834281/150‑1 (MQ=255)
                                                                                                                                tCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTCGCCGCGATAGTcc  >  1:380949/1‑150 (MQ=255)
                                                                                                                                             ||                                                                                                                                        
CGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGGTGGGTTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTCGCCGCGATAGTCC  >  NC_008463/5028060‑5028338

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