Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 3,390,474 Δ1 bp pseudogene (440/1997 nt) PA14_RS15455 → alkaline phosphatase D family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084633,390,4720G.99.6% 2319.0 / ‑0.4 469pseudogene (438/1997 nt)PA14_RS15455alkaline phosphatase D family protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base . (212/255);  minor base C (1/1);  total (213/256)
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 493 total aligned reads displayed.

CCTGCGCCTCGCCGCCGACGGCGGGGAGCCCGCGGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGGCGGCCGATGGGCCTGCTGCAGGTAGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCTGGTGGCCATCCATGGG  >  NC_008463/3390322‑3390618
                                                                                                                                                      |                                                                                                                                                  
ccTGCGCCTCGCCGCCGACGGCGGGGAGCCCGCGGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCg                                                                                                                                                    <  1:5553165/151‑1 (MQ=255)
ccTGCGCCTCGCCGCCGACGGCGGGGAGCCCGCGGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCg                                                                                                                                                    <  1:1654252/151‑1 (MQ=255)
 cTGCGCCTCGCCGCCGACGGCGGGGAGCCCGCGGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCgg                                                                                                                                                   <  2:10327170/151‑1 (MQ=255)
 cTGCGCCTCGCCGCCGACGGCGGGGAGCCCGCGGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCgg                                                                                                                                                   >  1:8014589/1‑151 (MQ=255)
     gcCTCGCCGCCGACGGCGGGGAGCCCGCGGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCggcggc                                                                                                                                               >  2:583217/1‑147 (MQ=255)
      ccTCGCCGCCGACGGCGGGGAGCCCGCGGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGcc                                                                                                                                             <  1:7030613/151‑1 (MQ=255)
       cTCGCCGCCGACGGCGGGGAGCCCGCGGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCCGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCTAGCCGCATCACCCGGCGGCCg                                                                                                                                            >  2:2567569/1‑151 (MQ=255)
       cTCGCCGCCGACGGCGGGGAGCCCGCGGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCg                                                                                                                                            >  1:2767327/1‑151 (MQ=255)
        tCGCCGCCGACGGCGGGGAGCCCGCGGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGa                                                                                                                                           >  1:9866798/1‑151 (MQ=255)
               cgACGGCGGGGAGCCCGCGGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTCGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATggcctg                                                                                                                                    >  1:8573043/1‑147 (MQ=255)
                  cggcggGGAGCCCGCGGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCctgctg                                                                                                                                <  2:2924027/151‑1 (MQ=255)
                   ggcggGGAGCCCGCGGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCtgctgc                                                                                                                               >  2:4378434/1‑151 (MQ=255)
                    gcggGGAGCCCGCGGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCa                                                                                                                              <  1:10728592/151‑1 (MQ=255)
                      ggggAGCCCGCGGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAgg                                                                                                                            >  2:6576505/1‑151 (MQ=255)
                       gggAGCCCGCGGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGt                                                                                                                           >  1:9629624/1‑151 (MQ=255)
                        ggAGCCCGCGGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGtc                                                                                                                          <  1:10978002/151‑2 (MQ=255)
                        ggAGCCCGCGGGAATCGCCGAAGGGGCGACGCCTCTGCTCGTCCCCTGTCGCGAGCGCGCCAGCCTGGTGCTCTAGCCTCTTCTCGACCAGTTGCTCCACGGCTCCTGCCGCTAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGtc                                                                                                                          <  1:7846381/151‑2 (MQ=255)
                         gAGCCCGCGGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGtcg                                                                                                                         <  1:2025622/151‑3 (MQ=255)
                          aGCCCGCGGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGa                                                                                                                        >  2:3546603/1‑151 (MQ=255)
                           gCCCGCGGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGAc                                                                                                                       <  1:8444697/151‑1 (MQ=255)
                            cccGCGGGACTCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGAcc                                                                                                                      >  2:7909399/1‑151 (MQ=255)
                            cccGCGGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGAcc                                                                                                                      <  2:3346436/151‑1 (MQ=255)
                            cccGCGGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGAcc                                                                                                                      >  1:4966586/1‑151 (MQ=255)
                            cccGCGGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGAcc                                                                                                                      <  2:8807502/151‑1 (MQ=255)
                             ccGCGGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCg                                                                                                                     <  2:7182653/151‑1 (MQ=255)
                             ccGCGGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCg                                                                                                                     <  1:617039/151‑1 (MQ=255)
                              cgcgGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCgg                                                                                                                    >  1:9448426/1‑151 (MQ=255)
                              cgcgGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCgg                                                                                                                    <  2:1925466/151‑1 (MQ=255)
                                cgGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGgct                                                                                                                  <  2:7563948/151‑1 (MQ=255)
                                cgGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGgct                                                                                                                  <  1:6264095/151‑1 (MQ=255)
                                   gAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGgctgct                                                                                                               <  1:3367685/151‑1 (MQ=255)
                                    aaTCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGctgctg                                                                                                              >  1:2545786/1‑151 (MQ=255)
                                       cGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTggag                                                                                                           <  1:5703515/151‑1 (MQ=255)
                                       cGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTggag                                                                                                           >  2:9099132/1‑151 (MQ=255)
                                        gCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCAGCAGGTCGAGCGGCTGCaggagg                                                                                                          >  1:1080191/1‑151 (MQ=255)
                                         ccGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTggagga                                                                                                         <  1:755354/151‑1 (MQ=255)
                                         ccGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTggagga                                                                                                         >  1:9552736/1‑151 (MQ=255)
                                          cGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGaa                                                                                                        >  1:10570454/1‑151 (MQ=255)
                                            aaTGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGaaaa                                                                                                      >  2:6253082/1‑151 (MQ=255)
                                            aaTGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGaaaa                                                                                                      >  2:8043732/1‑151 (MQ=255)
                                               gggCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACct                                                                                                   <  2:3800646/151‑1 (MQ=255)
                                               gggCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACct                                                                                                   >  2:1342831/1‑151 (MQ=255)
                                                ggCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACctg                                                                                                  <  2:10174423/151‑1 (MQ=255)
                                                ggCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACctg                                                                                                  <  2:5031929/151‑1 (MQ=255)
                                                  cgccgcATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACctgct                                                                                                <  1:7276744/151‑1 (MQ=255)
                                                  cgccgcATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACctgct                                                                                                >  1:9743439/1‑151 (MQ=255)
                                                    ccgcATCTGCTCTCCCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCtgctgc                                                                                              <  2:9487785/151‑1 (MQ=255)
                                                    ccgcATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCACGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCCTCTCCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGACAACCtgctgc                                                                                              >  2:3997305/1‑151 (MQ=255)
                                                       cATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAgg                                                                                           >  1:2986969/1‑151 (MQ=255)
                                                        aTCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTTCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGc                                                                                          <  1:3763428/151‑1 (MQ=255)
                                                        aTCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGc                                                                                          >  2:1684519/1‑151 (MQ=255)
                                                           tGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGa                                                                                       <  1:5312187/151‑1 (MQ=255)
                                                                  ccccGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAgccgc                                                                                <  1:5995460/151‑1 (MQ=255)
                                                                   cccGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAgccgcc                                                                               >  1:4510868/1‑151 (MQ=255)
                                                                      ggTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCgg                                                                            >  2:4206074/1‑151 (MQ=255)
                                                                       gTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGc                                                                           <  2:8164919/151‑1 (MQ=255)
                                                                          gCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGtt                                                                        >  1:11059369/1‑151 (MQ=255)
                                                                          gCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGtt                                                                        >  2:1819542/1‑151 (MQ=255)
                                                                              gcgcgcCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCtg                                                                    >  2:3155249/1‑151 (MQ=255)
                                                                                    cAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCtgatgatg                                                              <  2:8687310/151‑1 (MQ=255)
                                                                                       ccTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAgc                                                           <  1:8770365/151‑1 (MQ=255)
                                                                                         tGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAgcgg                                                         >  2:10959835/1‑151 (MQ=255)
                                                                                           gTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAgcggcg                                                       <  1:210442/151‑1 (MQ=255)
                                                                                            tGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGa                                                      <  1:3307342/151‑1 (MQ=255)
                                                                                            tGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGa                                                      <  1:11250904/151‑1 (MQ=255)
                                                                                              cTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGAcc                                                    <  1:5632966/151‑1 (MQ=255)
                                                                                               tCAAGCCACGTCTCGACCTGTTGCTCCACGGCTCCTGCCGCTAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGTCCGGCGGCTGGAGGAACTCCTGCTGGAGGCCGAAGGCCGCCCGGCGGTGCTGATGAGGGgcggcgggcg                                                   >  2:8396573/1‑147 (MQ=255)
                                                                                                cAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCCGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAg                                                  <  2:5668690/151‑1 (MQ=255)
                                                                                                    ccACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCt                                              <  1:4170200/151‑1 (MQ=255)
                                                                                                    ccACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCt                                              >  1:910286/1‑151 (MQ=255)
                                                                                                      aCGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTAc                                            <  1:8956076/151‑1 (MQ=255)
                                                                                                           tcGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGccta                                       <  1:2691734/151‑3 (MQ=255)
                                                                                                            cGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCcgac                                      <  2:1189044/151‑1 (MQ=255)
                                                                                                              aCCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCcgacga                                    >  2:10764959/1‑151 (MQ=255)
                                                                                                                 aGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGt                                 <  2:1456897/151‑1 (MQ=255)
                                                                                                                     gCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCGGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGCTGAGCGGCGCCCAGGTCTACGCCGACGACGTCGcc                             <  2:5521354/151‑1 (MQ=255)
                                                                                                                      cTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCg                            >  1:5859600/1‑151 (MQ=255)
                                                                                                                       tCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCgg                           <  1:7773047/151‑1 (MQ=255)
                                                                                                                          aCGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGCTGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGcc                        <  1:6868057/151‑1 (MQ=255)
                                                                                                                             gCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGAt                     <  1:3606232/151‑1 (MQ=255)
                                                                                                                              cTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATg                    >  2:8557840/1‑151 (MQ=255)
                                                                                                                               tCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGc                   >  1:10196115/1‑151 (MQ=255)
                                                                                                                                ccTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCt                  >  2:6763832/1‑151 (MQ=255)
                                                                                                                                ccTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCt                  >  2:3926177/1‑151 (MQ=255)
                                                                                                                                 cTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCtg                 <  1:1767017/151‑1 (MQ=255)
                                                                                                                                     cGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCtggtgg             >  2:8862681/1‑151 (MQ=255)
                                                                                                                                     cGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGGTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCtggtgg             >  2:7311515/1‑151 (MQ=255)
                                                                                                                                       cAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCTGGTGGcc           <  1:1418557/151‑1 (MQ=255)
                                                                                                                                        aaGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGTTGAGCGGCGACCCGGTCTACGCCGACGCCGTCGCCGGGCCGATGCTGGTGGcca          >  2:6032586/1‑151 (MQ=255)
                                                                                                                                        aaGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCTGGTGGcca          <  2:4731918/151‑1 (MQ=255)
                                                                                                                                        aaGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCAGGTGGcca          <  1:10079899/151‑1 (MQ=255)
                                                                                                                                         aGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCTGGTGGccat         <  2:7828197/151‑1 (MQ=255)
                                                                                                                                          gCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCTGGTGGccatc        <  2:10566485/151‑1 (MQ=255)
                                                                                                                                           ccGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCTGGTGGccatcc       >  2:924084/1‑151 (MQ=255)
                                                                                                                                            cGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCTGGTGGccatcca      <  1:2298132/151‑1 (MQ=255)
                                                                                                                                             gCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGTTGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCTGGTGGccatccat     >  1:6531456/1‑151 (MQ=255)
                                                                                                                                             gCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGCCGTCGCCGGGCCGATGCTGGTGGccatccat     >  2:7071072/1‑151 (MQ=255)
                                                                                                                                              cATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCTGGTGGCCATCCATg    <  1:1335397/151‑1 (MQ=255)
                                                                                                                                              cATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCTGGTGGCCATCCATg    <  2:2787344/151‑1 (MQ=255)
                                                                                                                                               aTCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCTGGTGGCCATCCATgg   <  2:2394307/151‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
CCTGCGCCTCGCCGCCGACGGCGGGGAGCCCGCGGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGGCGGCCGATGGGCCTGCTGCAGGTAGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCTGGTGGCCATCCATGGG  >  NC_008463/3390322‑3390618

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