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,2010GC98.6% 1896.8 / 5.9 485intergenic (‑123/+14)PA14_RS22975/PA14_RS22980Nramp family divalent metal transporter/ribbon‑helix‑helix domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base G (2/5);  new base C (236/241);  total (238/247)
Fisher's exact test for biased strand distribution p-value = 4.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0084635,028,2020TG98.6% 1880.3 / 12.9 485intergenic (‑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 (236/242);  minor base . (2/5);  total (238/247)
Fisher's exact test for biased strand distribution p-value = 4.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 485 total aligned reads displayed.

CGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGGTGGGTTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTCGCCGCGATAGTCCA  >  NC_008463/5028060‑5028339
                                                                                                                                             ||                                                                                                                                         
  tCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTccctcccc                                                                                                                               >  2:6171868/1‑151 (MQ=255)
    ggccggcACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTccctcccctc                                                                                                                             <  2:9677715/151‑1 (MQ=255)
    ggccggcACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTccctcccctc                                                                                                                             >  1:5756274/1‑151 (MQ=255)
       cggcACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAgg                                                                                                                          >  2:10179560/1‑151 (MQ=255)
       cggcACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAgg                                                                                                                          >  1:3211944/1‑151 (MQ=255)
        ggcACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGc                                                                                                                         >  2:4034754/1‑151 (MQ=255)
        ggcACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTC‑CCCCCGGTCGCGTCTTCCCTCCCCTCAGGc                                                                                                                         <  2:7338747/151‑1 (MQ=255)
          cACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCt                                                                                                                       <  2:3528982/151‑1 (MQ=255)
          cACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCt                                                                                                                       >  2:6243313/1‑151 (MQ=255)
           aCCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTg                                                                                                                      >  2:3047562/1‑151 (MQ=255)
           aCCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTg                                                                                                                      >  2:6917818/1‑151 (MQ=255)
             cGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTggg                                                                                                                    >  2:10831347/1‑151 (MQ=255)
                 aCGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGgcgcg                                                                                                                <  1:4108827/151‑1 (MQ=255)
                   gCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCGCCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCt                                                                                                              <  2:1482657/151‑1 (MQ=255)
                    cTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTg                                                                                                             <  1:5501288/151‑1 (MQ=255)
                     tgtgGCTGACCGCCGAGGGACAGAACGGCGCGGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTgg                                                                                                            <  2:10465572/151‑1 (MQ=255)
                          cTGACCGCCGAGGGACAGAACGGCGCGGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGAc                                                                                                       <  2:1892388/151‑1 (MQ=255)
                             aCCGCCGAGGGACAGAACGGCGCGGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAg                                                                                                    <  1:6830939/151‑1 (MQ=255)
                                gccgAGGGACAGAACGGCGCGGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCt                                                                                                 >  1:56259/1‑151 (MQ=255)
                                 ccgAGGGACAGAACGGCGCGGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTc                                                                                                <  1:4651380/151‑1 (MQ=255)
                                  cgAGGGACAGAACGGCGCGGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGACGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCa                                                                                               >  2:7040930/1‑151 (MQ=255)
                                  cgAGGGACAGAACGGCGCGGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCa                                                                                               <  2:8316758/151‑1 (MQ=255)
                                  cgAGGGACAGAACGGCGCGGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCa                                                                                               >  2:4985215/1‑151 (MQ=255)
                                  cgAGGGACAGAACGGCGCGGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCa                                                                                               >  2:9026868/1‑151 (MQ=255)
                                  cgAGGGACAGAACGGCGCGGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCa                                                                                               <  2:7795603/151‑1 (MQ=255)
                                   gAGGGACAGAACGGCGCGGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAg                                                                                              >  1:20801/1‑151 (MQ=255)
                                   gAGGGACAGAACGGCGCGGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAg                                                                                              <  1:706401/151‑1 (MQ=255)
                                     gggACAGAACGGCGCGGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCt                                                                                            <  1:8754243/151‑1 (MQ=255)
                                     gggACAGAACGGCGCGGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCt                                                                                            <  2:6712744/151‑1 (MQ=255)
                                     gggACAGAACGGCGCGGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCt                                                                                            <  1:4913742/151‑1 (MQ=255)
                                       gACAGAACGGCGCGGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCTCCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCg                                                                                          <  2:1395389/151‑1 (MQ=255)
                                         cAGAACGGCGCGGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGc                                                                                        >  2:8491640/1‑151 (MQ=255)
                                          aGAACGGCGCGGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCg                                                                                       <  1:2704897/151‑1 (MQ=255)
                                              cGGCGCGGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCt                                                                                   <  2:3808167/151‑1 (MQ=255)
                                                 cgcgGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTtct                                                                                <  2:4630592/151‑1 (MQ=255)
                                                 cgcgGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCGTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCACGGCGACCTtct                                                                                <  1:6968178/151‑1 (MQ=255)
                                                  gcgGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTtctc                                                                               >  1:11101786/1‑151 (MQ=255)
                                                  gcgGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTtctc                                                                               >  1:11214381/1‑151 (MQ=255)
                                                   cgGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTcg                                                                              <  2:1800109/151‑1 (MQ=255)
                                                      gtggccGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTcgcg                                                                            >  1:3971404/5‑151 (MQ=255)
                                                           cGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTCAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCt                                                                       <  2:7516411/151‑1 (MQ=255)
                                                              ggtggGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTgcg                                                                    <  2:8670513/151‑1 (MQ=255)
                                                               gtggGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTgcgc                                                                   <  1:65563/151‑1 (MQ=255)
                                                               gtggGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTgcgc                                                                   >  1:9996497/1‑151 (MQ=255)
                                                                tggGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGcc                                                                  <  2:5639842/151‑1 (MQ=255)
                                                                  ggCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAg                                                                <  1:1288945/151‑1 (MQ=255)
                                                                         aGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCgg                                                         >  2:10983483/1‑151 (MQ=255)
                                                                          gTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCggg                                                        <  2:534083/151‑1 (MQ=255)
                                                                           tATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGt                                                       >  1:620685/1‑151 (MQ=255)
                                                                              ttCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGCCCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCa                                                    >  2:2351067/1‑151 (MQ=255)
                                                                               tCACGAGGCGACGCGGGATCGGGGAGGACGTGCGACCAGCCTAAACGTTGCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCGGCGCCAGGTCGCGGGGGCaa                                                   >  2:8023236/1‑151 (MQ=255)
                                                                                 aCGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACt                                                 <  2:10634727/151‑1 (MQ=255)
                                                                                  cGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCCGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTg                                                >  1:10345084/1‑151 (MQ=255)
                                                                                  cGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTg                                                <  1:8198954/151‑1 (MQ=255)
                                                                                    aGGCGACGCCGGATCGGGGAGGACGTGCGACCAGCCTAAACGTTTCACCCCGGACTCG‑ACTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGac                                              <  2:2715337/151‑1 (MQ=255)
                                                                                    aGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGGGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGac                                              <  1:8462340/151‑1 (MQ=255)
                                                                                    aGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGac                                              <  2:3216680/151‑1 (MQ=255)
                                                                                       cGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACg                                           >  2:5987068/1‑151 (MQ=255)
                                                                                       cGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCGGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCGTCTCGCGCCGCTGCGCGAGGTAGCGGGTGCAACTGTCACg                                           >  1:9857074/1‑151 (MQ=255)
                                                                                         aCGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCa                                         <  2:6114962/151‑1 (MQ=255)
                                                                                         aCGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCa                                         >  1:8550234/1‑151 (MQ=255)
                                                                                          cGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCag                                        <  2:4797177/151‑1 (MQ=255)
                                                                                          cGCCGGATCGGGGACGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACATTCTCGCGCCGCTGCCCCAGGTAGCGGGTGCAACAGACACGCag                                        >  2:2003951/1‑151 (MQ=255)
                                                                                           gCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCagg                                       >  2:4196560/1‑151 (MQ=255)
                                                                                            ccGGATCGGGGAGGAAGTGCGACCAGCCTAACCGTTTCACCCC‑GGTCGAGTCGTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCagga                                      <  2:1702771/151‑1 (MQ=255)
                                                                                            ccGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCagga                                      <  1:10040407/151‑1 (MQ=255)
                                                                                             cGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCGTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCaggat                                     >  1:1758449/1‑150 (MQ=255)
                                                                                               gATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCaggaagg                                   >  1:665731/1‑151 (MQ=255)
                                                                                               gATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCaggaagg                                   <  1:9083346/151‑1 (MQ=255)
                                                                                               gATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACAGACACGCaggaagg                                   >  1:2515253/1‑151 (MQ=255)
                                                                                                aTCGGGGAGGAAGTGCGACCAGCCTCAACGTTTC‑CCCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCaggaagga                                  <  2:10723414/151‑1 (MQ=255)
                                                                                                  cGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCGGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAaggagg                                <  1:1097031/151‑1 (MQ=255)
                                                                                                  cGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAaggagg                                <  2:5962755/151‑1 (MQ=255)
                                                                                                   ggggAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGc                               <  1:8895762/151‑1 (MQ=255)
                                                                                                    gggAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCg                              <  1:3165310/151‑1 (MQ=255)
                                                                                                         gaAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGtt                         >  1:77186/1‑151 (MQ=255)
                                                                                                           aGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGa                       <  1:4373200/151‑1 (MQ=255)
                                                                                                             tGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGc                     >  1:8652231/1‑151 (MQ=255)
                                                                                                               cGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCAc                   <  1:346024/151‑1 (MQ=255)
                                                                                                                gACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCAcc                  >  1:3240763/1‑151 (MQ=255)
                                                                                                                 aCCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGTAGGAGGCGAAGTTGTGCACCt                 >  1:10642098/1‑151 (MQ=255)
                                                                                                                 aCCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCt                 >  2:7279744/1‑151 (MQ=255)
                                                                                                                 aCCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCt                 <  2:2905202/151‑1 (MQ=255)
                                                                                                                 aCCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCACTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCt                 >  1:3624827/1‑151 (MQ=255)
                                                                                                                  ccAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTc                >  1:2864227/1‑151 (MQ=255)
                                                                                                                  ccAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTc                <  1:9537061/151‑1 (MQ=255)
                                                                                                                    aGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCGCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTcgc              >  1:1934494/1‑151 (MQ=255)
                                                                                                                    aGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCGCAGCTCGGCGACCTTCGCGCGCCGCTGCGCCAGGTAGCGGGTGCCACTGACACGCAGGTAGGAGGCGACGTTGAGCACCTcgc              >  2:11191348/1‑151 (MQ=255)
                                                                                                                     gCCTAAACGTTTCACCCC‑GGTCTCGACTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCTGCTGGGCGACCTTCTCGCGCCGCTGCGCGAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTcgcc             <  2:1192059/151‑1 (MQ=255)
                                                                                                                     gCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTcgcc             <  1:63540/151‑1 (MQ=255)
                                                                                                                        tAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTCGCcgcg          >  1:2993869/1‑151 (MQ=255)
                                                                                                                          aaCGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTCGCCGCGAt        >  1:9241057/1‑151 (MQ=255)
                                                                                                                            cGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTCGCCGCGATAg      >  2:6509703/1‑151 (MQ=255)
                                                                                                                            cGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTCGCCGCGATAg      <  1:691638/151‑1 (MQ=255)
                                                                                                                            cGTTTCACCCC‑GGTCGCGTCGTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTCGCCGCGATAg      <  2:9228065/151‑1 (MQ=255)
                                                                                                                             gTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTCGCCGCGATAGt     <  1:5188736/151‑1 (MQ=255)
                                                                                                                             gTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTCGCCGCGATAGt     >  2:4522005/1‑151 (MQ=255)
                                                                                                                             gTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTCGCCGCGATAGt     >  2:9834779/1‑151 (MQ=255)
                                                                                                                              tttCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTCGCCGCGATAGTc    <  1:7508298/151‑1 (MQ=255)
                                                                                                                               ttCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTCGCCGCGATAGTcc   >  1:1508557/1‑151 (MQ=255)
                                                                                                                                             ||                                                                                                                                         
CGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGGTGGGTTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTCGCCGCGATAGTCCA  >  NC_008463/5028060‑5028339

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: ‑