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.9% 1707.6 / 3.5 448intergenic (‑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/3);  new base C (198/245);  total (200/248)
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
*NC_0084635,028,2020TG98.9% 1716.6 / 6.5 448intergenic (‑124/+13)PA14_RS22975/PA14_RS22980Nramp family divalent metal transporter/ribbon‑helix‑helix domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  major base G (197/245);  minor base . (2/3);  total (200/248)
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 448 total aligned reads displayed.

CGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGGTGGGTTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTCGCCGCGATAGTCC  >  NC_008463/5028060‑5028338
                                                                                                                                             ||                                                                                                                                        
  tCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTtccctccc                                                                                                                               <  2:5420347/150‑1 (MQ=255)
    ggccggcACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTccctcccct                                                                                                                             >  2:1168429/1‑150 (MQ=255)
    ggccggcACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTccctcccct                                                                                                                             <  2:3775952/150‑1 (MQ=255)
     gccggcACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTccctcccctc                                                                                                                            >  2:530224/1‑150 (MQ=255)
      ccggcACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCa                                                                                                                           <  2:11020956/150‑1 (MQ=255)
        ggcACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAgg                                                                                                                         >  2:4930278/1‑150 (MQ=255)
           aCCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCt                                                                                                                      <  2:3455156/150‑1 (MQ=255)
            ccGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTg                                                                                                                     >  2:1773594/1‑150 (MQ=255)
              gCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTggg                                                                                                                   >  1:3406623/1‑150 (MQ=255)
               ccACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGgc                                                                                                                  >  2:7163642/1‑150 (MQ=255)
                 ccGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGgcgc                                                                                                                >  2:8027077/2‑150 (MQ=255)
                 aCGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGG‑CCGTGGTGGGCCGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTGAGGCCTGGgcgc                                                                                                                >  2:2530927/1‑150 (MQ=255)
                 aCGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGgcgc                                                                                                                >  1:1974412/1‑150 (MQ=255)
                  cGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGgcgcg                                                                                                               <  1:2342893/150‑1 (MQ=255)
                   gCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGgcgcgc                                                                                                              <  1:3233855/150‑1 (MQ=255)
                   gCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGgcgcgc                                                                                                              <  1:8190410/150‑1 (MQ=255)
                    cTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCt                                                                                                             <  1:5175456/150‑1 (MQ=255)
                      gtgGCTGACCGCCGAGGGACAGAACGGCGCGGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTgg                                                                                                           <  2:4600675/150‑1 (MQ=255)
                            gACCGCCGAGGGACAGAACGGCGCGGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGAcc                                                                                                     <  2:9693748/150‑1 (MQ=255)
                             aCCGCCGAGGGACAGAACGGCGCGGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCa                                                                                                    <  2:10321545/150‑1 (MQ=255)
                               cgccgAGGGACAGAACGGCGCGGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGATTCACCCC‑GGTCGCGTCTTCCCACCCCTCAGGCCTGGGCGCGCTGCCCGACCAgg                                                                                                  <  1:206235/150‑1 (MQ=255)
                                   gAGGGACAGAACGGCGCGGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCa                                                                                              <  2:6363227/150‑1 (MQ=255)
                                    aGGGACAGAACGGCGCGGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGCTCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAg                                                                                             <  1:8622692/150‑1 (MQ=255)
                                    aGGGACAGAACGGCGCGGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAg                                                                                             >  2:2475168/1‑150 (MQ=255)
                                     gggACAGAACGGCGCGGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGc                                                                                            <  2:3583385/150‑1 (MQ=255)
                                     gggACAGAACGGCGCGGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGc                                                                                            >  1:1703881/1‑150 (MQ=255)
                                       gACAGAACGGCGCGGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTc                                                                                          >  1:5655720/1‑150 (MQ=255)
                                           gAACGGCGCGGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCg                                                                                      <  1:4593231/150‑1 (MQ=255)
                                           gAACGGCGCGGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCg                                                                                      <  1:1676685/150‑1 (MQ=255)
                                            aaCGGCGCGGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGa                                                                                     <  1:7024768/150‑1 (MQ=255)
                                            aaCGGCGCGGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGa                                                                                     <  1:9219096/150‑1 (MQ=255)
                                             aCGGCGCGGTGG‑CCGTGGTGGGCAGGGAGTATTTCCCGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGGGTCACCCC‑GGTCGCGTCTTCCCTCCCCTCCGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGAc                                                                                    <  2:7613523/150‑1 (MQ=255)
                                              cGGCGCGGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGAcc                                                                                   <  1:6308599/150‑1 (MQ=255)
                                               ggCGCGGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCt                                                                                  <  1:10987501/150‑1 (MQ=255)
                                               cgcgcGGTGG‑CCGGGCTGGGCAGGGAGGATTTCACGCGGCGACGCCGGATCGGGGAGGACGTGCGACCATCCTAAACGTGTCACCCC‑GGTCGCGACTGCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCt                                                                                  <  2:6490089/149‑1 (MQ=255)
                                                 cgcgGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTtc                                                                                <  2:4424043/150‑1 (MQ=255)
                                                   cgGTGG‑CCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTtctc                                                                              <  1:10239967/150‑1 (MQ=255)
                                                    ggTGG‑CCGTGGTGGGCAGTGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTcg                                                                             >  2:336964/1‑150 (MQ=255)
                                                      gtggccGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTcgc                                                                            >  1:10548545/5‑150 (MQ=255)
                                                        ggccGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCgcgc                                                                          <  2:8715611/148‑1 (MQ=255)
                                                        ggccGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCgcgc                                                                          <  2:9100222/148‑1 (MQ=255)
                                                         gccGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGcgcc                                                                         <  1:10009547/149‑1 (MQ=255)
                                                                 gggCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCGCGGCGACCTTCTCGCGCCGCTGCGcc                                                                 >  1:351652/1‑150 (MQ=255)
                                                                  ggCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCa                                                                >  2:1455666/1‑150 (MQ=255)
                                                                   gCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGCCCAGCCTAACCGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAg                                                               <  2:4274643/150‑1 (MQ=255)
                                                                     aGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGt                                                             >  2:9374488/1‑150 (MQ=255)
                                                                      gggAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGGCAGGtc                                                            >  1:2418590/1‑149 (MQ=255)
                                                                      gggAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTa                                                            <  2:7302055/150‑1 (MQ=255)
                                                                       ggAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAg                                                           <  1:4401491/150‑1 (MQ=255)
                                                                       ggAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAg                                                           >  1:7312062/1‑150 (MQ=255)
                                                                       ggAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAg                                                           >  1:8949267/1‑150 (MQ=255)
                                                                        gAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGc                                                          <  2:7370557/150‑1 (MQ=255)
                                                                          gTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCgg                                                        <  1:8360444/150‑1 (MQ=255)
                                                                          gTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCgg                                                        >  1:3667998/1‑150 (MQ=255)
                                                                            aTTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGt                                                      <  1:9608743/150‑1 (MQ=255)
                                                                                cACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCaa                                                  >  2:10134501/1‑150 (MQ=255)
                                                                                 aCGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAAc                                                 <  2:7737120/150‑1 (MQ=255)
                                                                                 aCGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAAc                                                 <  2:621946/150‑1 (MQ=255)
                                                                                   gAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTg                                               >  1:8988173/1‑150 (MQ=255)
                                                                                   gAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCGTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTg                                               >  2:6053650/1‑150 (MQ=255)
                                                                                     ggCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGac                                             <  1:1389573/150‑1 (MQ=255)
                                                                                     ggCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGac                                             >  1:7974454/1‑150 (MQ=255)
                                                                                     ggCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGac                                             <  1:1137534/150‑1 (MQ=255)
                                                                                         aCGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGc                                         >  2:5711984/1‑150 (MQ=255)
                                                                                         aCGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCGTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCTACTGACACGc                                         >  1:10387249/1‑150 (MQ=255)
                                                                                          cGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCa                                        <  1:2091651/150‑1 (MQ=255)
                                                                                           gCCGGATCGGGGAGCCAGTTCGACCTGCCTAAACGTGTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCTGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCag                                       <  2:9253157/150‑1 (MQ=255)
                                                                                            ccGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCGTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCagg                                      >  1:7559990/1‑150 (MQ=255)
                                                                                             cGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCagga                                     <  1:9114822/150‑1 (MQ=255)
                                                                                             cGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCagga                                     <  2:8148970/150‑1 (MQ=255)
                                                                                               gATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTCGCGGGTGCAACTGACACGCaggaag                                   >  2:2034367/1‑150 (MQ=255)
                                                                                                aTCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCaggaagg                                  <  1:1603629/150‑1 (MQ=255)
                                                                                                aTCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCGTCTCGCGCCGCGGCGCCAGGGAGCGGGTGCAACTGACACGCgggaagg                                  >  1:279590/1‑150 (MQ=255)
                                                                                                 tCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTCGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCaggaagga                                 <  1:5107445/150‑1 (MQ=255)
                                                                                                    gggAGGACGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGc                              >  2:1592319/1‑150 (MQ=255)
                                                                                                    gggAGGAAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGc                              <  2:3936548/150‑1 (MQ=255)
                                                                                                     ggaggaAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCg                             <  1:8741297/150‑1 (MQ=255)
                                                                                                         gaAGTGCGACCAGCCTACACGTTTCACCCC‑GGTCGCGTCGTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGt                         <  2:10794/150‑1 (MQ=255)
                                                                                                         gaAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGt                         >  2:8582292/1‑150 (MQ=255)
                                                                                                         gaAGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCGTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGt                         >  1:5959896/1‑150 (MQ=255)
                                                                                                          aaGTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGtt                        <  2:6805559/150‑1 (MQ=255)
                                                                                                            gTGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGa                      <  1:4012477/150‑1 (MQ=255)
                                                                                                             tGCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAg                     <  2:8988173/150‑1 (MQ=255)
                                                                                                              gCGACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGc                    >  2:8344640/1‑150 (MQ=255)
                                                                                                               cGACCAGCCTACACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCa                   >  1:9708017/1‑150 (MQ=255)
                                                                                                                gACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCAc                  <  2:10669776/150‑1 (MQ=255)
                                                                                                                gACCAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCGTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCAc                  >  2:9806464/1‑150 (MQ=255)
                                                                                                                   cAGCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTc               <  2:4098404/150‑1 (MQ=255)
                                                                                                                     gCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTcgc             >  2:7938684/1‑150 (MQ=255)
                                                                                                                     gCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTcgc             <  1:3053528/150‑1 (MQ=255)
                                                                                                                     gCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTcgc             >  2:10885591/1‑150 (MQ=255)
                                                                                                                     gCCTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTcgc             <  1:7776196/150‑1 (MQ=255)
                                                                                                                      ccTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTcgcc            <  1:6528974/150‑1 (MQ=255)
                                                                                                                      ccTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTcgcc            <  1:6933979/150‑1 (MQ=255)
                                                                                                                       cTAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTcgccg           >  1:9893226/1‑150 (MQ=255)
                                                                                                                        tAAACGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTcgccgc          <  1:4750382/150‑1 (MQ=255)
                                                                                                                         aaaCGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGGGCAACTGACACGCAGGAAGGAGGCGAAGTTGCGCACCTCGCcgcg         >  1:10759457/1‑150 (MQ=255)
                                                                                                                          aaCGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTCGCCGCGa        <  2:10467351/150‑1 (MQ=255)
                                                                                                                           aCGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCCGGAAGGAGGCGAAGTTGAGCACCTCGCCGCGAt       >  2:288830/1‑150 (MQ=255)
                                                                                                                           aCGTTTCACCCC‑GGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTCGCCGCGAt       >  2:2375099/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: ‑