Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 5,028,197 Δ1 bp intergenic (‑119/+18) 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,1950G.97.2% 2143.4 / 51.0 465intergenic (‑117/+20)PA14_RS22975/PA14_RS22980Nramp family divalent metal transporter/ribbon‑helix‑helix domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base G (7/6);  new base . (230/220);  total (238/227)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 500 total aligned reads displayed.

TGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGGTGGGTTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTCGCCGCGATAGTCC  >  NC_008463/5028045‑5028338
                                                                                                                                                      |                                                                                                                                               
tGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCg                                                                                                                                                 >  2:5513157/1‑150 (MQ=255)
   gcAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCggtc                                                                                                                                              <  2:8579388/150‑3 (MQ=255)
        ggcgAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCggtcgcgtc                                                                                                                                         <  1:7003410/150‑8 (MQ=255)
        ggcgAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCggtcgcgtc                                                                                                                                         >  1:5986084/1‑143 (MQ=255)
                 tCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTtccctccc                                                                                                                               <  2:1531719/150‑1 (MQ=255)
                   ggccggcACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTccctcccct                                                                                                                             >  1:10389264/1‑150 (MQ=255)
                   ggccggcACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTccctcccct                                                                                                                             <  1:5646512/150‑1 (MQ=255)
                    gccggcACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTccctcccctc                                                                                                                            <  2:162850/150‑1 (MQ=255)
                     ccggcACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCa                                                                                                                           <  1:2831449/150‑1 (MQ=255)
                     ccggcACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCa                                                                                                                           <  1:10978921/150‑1 (MQ=255)
                     ccggcACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCa                                                                                                                           <  1:2596554/150‑1 (MQ=255)
                      cggcACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAg                                                                                                                          >  1:6169408/1‑150 (MQ=255)
                        gcACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGGGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGc                                                                                                                        >  2:3070357/1‑150 (MQ=255)
                         cACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGcc                                                                                                                       <  2:2070450/150‑1 (MQ=255)
                           ccGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGCAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTg                                                                                                                     >  2:7509834/1‑150 (MQ=255)
                           ccGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTg                                                                                                                     <  2:2596243/150‑1 (MQ=255)
                           ccGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTg                                                                                                                     <  1:8888320/150‑1 (MQ=255)
                           ccGCCACGCTGTGGCTGACCGCCGAGGGACAGAAAGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTg                                                                                                                     >  2:3687133/1‑150 (MQ=255)
                            cGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGTGGGGCAGGAATTATTGCACGAGGCGACAACTGAGCGGGGAGGAAGAGCTAACAGCCAAAACGTTTCACACAGGGAGCGTCTTCACTCCCCTAAGGCCTgg                                                                                                                    >  2:9843047/1‑150 (MQ=255)
                             gCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGGGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCTGGCCTggg                                                                                                                   >  1:2461173/1‑150 (MQ=255)
                                  gCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGgcgcgc                                                                                                              >  1:4538618/1‑150 (MQ=255)
                                  gCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGgcgcgc                                                                                                              <  1:10857847/150‑1 (MQ=255)
                                  gCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTAATTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGgcgcgc                                                                                                              <  2:1823522/150‑1 (MQ=255)
                                   cTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCt                                                                                                             <  1:10062843/150‑1 (MQ=255)
                                      tgGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGc                                                                                                          >  2:105591/1‑150 (MQ=255)
                                      tgGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGc                                                                                                          >  1:585709/1‑150 (MQ=255)
                                           gACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGAcc                                                                                                     >  1:9285770/1‑150 (MQ=255)
                                           gACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGAcc                                                                                                     <  2:6104463/150‑1 (MQ=255)
                                                ccgAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCt                                                                                                >  1:9043990/1‑150 (MQ=255)
                                                  gAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCa                                                                                              <  2:3495269/150‑1 (MQ=255)
                                                  gAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCa                                                                                              <  1:1803917/150‑1 (MQ=255)
                                                    gggACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGc                                                                                            <  2:2730431/150‑1 (MQ=255)
                                                    gggACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGc                                                                                            >  1:1473659/1‑150 (MQ=255)
                                                       aCAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCg                                                                                         <  2:1195934/150‑1 (MQ=255)
                                                        cAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCgg                                                                                        >  2:10143332/1‑150 (MQ=255)
                                                        cAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCgg                                                                                        >  1:3633936/1‑150 (MQ=255)
                                                        cAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCgg                                                                                        >  1:7993896/1‑150 (MQ=255)
                                                         aGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAATCGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGTCCAGGCTCAGCGCGGc                                                                                       >  2:5815135/1‑150 (MQ=255)
                                                            aCGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGAc                                                                                    <  1:3254645/150‑1 (MQ=255)
                                                               gcgcGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCtt                                                                                 >  1:8386814/1‑150 (MQ=255)
                                                                cgcgGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTtc                                                                                >  2:1375062/1‑150 (MQ=255)
                                                                     gtggccGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTcgc                                                                            >  1:10268234/5‑150 (MQ=255)
                                                                       ggccGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCgcgc                                                                          <  1:5276512/148‑1 (MQ=255)
                                                                        gccGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGcgcc                                                                         <  2:2334454/149‑1 (MQ=255)
                                                                            tggtggGCAGGGAGTATGTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTg                                                                     <  2:8405132/150‑1 (MQ=255)
                                                                             ggtggGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTgc                                                                    >  1:6710500/1‑150 (MQ=255)
                                                                             ggtggGCAGGGAGTATATCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTgc                                                                    >  2:9498994/1‑150 (MQ=255)
                                                                              gtggGCAGGGAGTCTTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTgcg                                                                   <  2:6970120/150‑1 (MQ=255)
                                                                              gtggGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCGCGGCGACCTTCTCGCGCCGCTgcg                                                                   >  2:5310693/1‑150 (MQ=255)
                                                                                gggCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGcc                                                                 <  1:4805532/150‑1 (MQ=255)
                                                                                 ggCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGCGGAAGTGCGACCAGCATAAACTTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCa                                                                <  1:5051424/150‑1 (MQ=255)
                                                                                  gCAGGGAGTATTTCACGAGGCGACGCCGGATCAGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAg                                                               >  2:1922384/1‑150 (MQ=255)
                                                                                   cAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAgg                                                              <  1:2925278/150‑1 (MQ=255)
                                                                                    aGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGt                                                             <  1:787040/150‑1 (MQ=255)
                                                                                    aGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCGGCGCCTGGt                                                             >  1:6836873/1‑150 (MQ=255)
                                                                                     gggAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTa                                                            >  2:8735254/1‑150 (MQ=255)
                                                                                     gggAGTAGTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTa                                                            <  1:6564163/150‑1 (MQ=255)
                                                                                       gcgTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCGTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCTACCTTCTCGCGCCGCTGCGCCAGGTAGc                                                          <  2:505396/148‑1 (MQ=255)
                                                                                            tttCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTg                                                     >  2:4189701/1‑150 (MQ=255)
                                                                                            tttCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTg                                                     <  1:7564572/150‑1 (MQ=255)
                                                                                              tCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCa                                                   <  2:9037720/150‑1 (MQ=255)
                                                                                               cACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTCGCGGGTGCaa                                                  >  2:7985396/1‑150 (MQ=255)
                                                                                               cACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCaa                                                  >  2:5665548/1‑150 (MQ=255)
                                                                                                  gAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTg                                               <  1:8275880/150‑1 (MQ=255)
                                                                                                   aGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGa                                              <  2:6293905/150‑1 (MQ=255)
                                                                                                   aGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCAGCTGCGCCAGGTAGCGGGTGCAACTGa                                              >  1:10567960/1‑150 (MQ=255)
                                                                                                   aGGCGACGCCGGATCGGGGAGGAAGAGCGACCAGCCTACACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGa                                              <  2:3786458/150‑1 (MQ=255)
                                                                                                     gCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGaca                                            <  1:5471692/150‑1 (MQ=255)
                                                                                                      cGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGacac                                           <  1:3758115/150‑1 (MQ=255)
                                                                                                        aCGCCGGATCGGGGAGGACGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCCGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGc                                         >  2:3899471/1‑150 (MQ=255)
                                                                                                        aCGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGc                                         >  1:1962343/1‑150 (MQ=255)
                                                                                                         cGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCa                                        >  1:4081244/1‑150 (MQ=255)
                                                                                                          gCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCTCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCag                                       <  1:7075116/150‑1 (MQ=255)
                                                                                                           ccGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCagg                                      >  2:10790567/1‑150 (MQ=255)
                                                                                                            cGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCagga                                     >  1:2160320/1‑150 (MQ=255)
                                                                                                             ggATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCaggaa                                    >  2:6449978/1‑150 (MQ=255)
                                                                                                              gATCGGGGAGGACGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCaggaag                                   >  2:4907589/1‑150 (MQ=255)
                                                                                                                tCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCaggaagga                                 <  1:9428750/150‑1 (MQ=255)
                                                                                                                tCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCGTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCACCTGACACGCaggaagga                                 >  2:10441735/1‑150 (MQ=255)
                                                                                                                tCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCATCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCaggaagga                                 >  2:8854350/1‑150 (MQ=255)
                                                                                                                  ggggAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCCGGGAGCGGGTGCAACTGACACGCAGGGaggagg                               >  2:6747198/1‑150 (MQ=255)
                                                                                                                  ggggAGGAAGTGCGACCAGCCTAAACGTTACACCCCGGACGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTACTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAaggagg                               <  2:7817694/150‑1 (MQ=255)
                                                                                                                   gggAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGc                              >  1:10328073/1‑150 (MQ=255)
                                                                                                                   gggAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGc                              <  1:3511456/150‑1 (MQ=255)
                                                                                                                    ggaggaAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCg                             >  2:5959997/1‑150 (MQ=255)
                                                                                                                    ggaggaAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCg                             >  2:8070351/1‑150 (MQ=255)
                                                                                                                    ggaggaAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCg                             <  1:7292532/150‑1 (MQ=255)
                                                                                                                    ggaggaAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCg                             <  1:1274648/150‑1 (MQ=255)
                                                                                                                     gaggaAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGa                            >  1:8147726/1‑150 (MQ=255)
                                                                                                                      aggaAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGaa                           <  1:1074229/150‑1 (MQ=255)
                                                                                                                        gaAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCGGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGt                         <  2:7670197/150‑1 (MQ=255)
                                                                                                                        gaAGTGCGACCAGCCTAAACGTGTCACCCCGGTCGCGTCGTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCACGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGt                         <  2:7221263/150‑1 (MQ=255)
                                                                                                                          aGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTg                       >  2:4777841/1‑150 (MQ=255)
                                                                                                                          aGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTg                       >  1:569474/1‑150 (MQ=255)
                                                                                                                           gTGCGACCAGCCAAAACGTTTCACCACGGTCGCGACTTCCCTCCCCTCAGCCCTGTGCGCGCTGGCCTACCAGGAGCAGCACGCCGACAATCTCGCGCCGATCCGCCAGGTAGCGGGTGCAACTGACACGAAGTAAGGAGGCGAAGATGa                      >  2:7056818/1‑150 (MQ=255)
                                                                                                                                    gCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTcgc             <  1:9571841/150‑1 (MQ=255)
                                                                                                                                      cTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTcgccg           >  2:4392553/1‑150 (MQ=255)
                                                                                                                                         aaCGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTCGCCGCGa        <  1:6319778/150‑1 (MQ=255)
                                                                                                                                             tttCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTCGCCGCGATAGt    >  2:3350238/1‑150 (MQ=255)
                                                                                                                                             tttCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTCGCCGCGATAGt    >  1:8738378/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                               
TGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGGTGGGTTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTCGCCGCGATAGTCC  >  NC_008463/5028045‑5028338

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