Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 3,390,485 Δ1 bp pseudogene (451/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,4830G.99.8% 2344.2 / ‑3.6 473pseudogene (449/1997 nt)PA14_RS15455alkaline phosphatase D family protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base . (245/227);  total (245/228)
Fisher's exact test for biased strand distribution p-value = 4.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 491 total aligned reads displayed.

GCCGCCGACGGCGGGGAGCCCGCGGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGGCGGCCGATGGGCCTGCTGCAGGTAGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCTGGTGGCCATCCATGGGCTGATCCGCCG  >  NC_008463/3390332‑3390629
                                                                                                                                                       |                                                                                                                                                  
gccgccGACGGCGGGGAGCCCGCGGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATg                                                                                                                                                    <  1:5077861/151‑1 (MQ=255)
gccgccGACGGCGGGGAGCCCGCGGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATg                                                                                                                                                    >  1:11236165/1‑151 (MQ=255)
  cgccgACGGCGGGGAGCCCGCGGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATggc                                                                                                                                                  <  1:1253827/151‑2 (MQ=255)
    ccgACGGCGGGGAGCCCGCGGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATggcct                                                                                                                                                >  2:1035085/1‑148 (MQ=255)
       aCGGCGGGGAGCCCGCGGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCctgct                                                                                                                                            >  1:1999421/1‑151 (MQ=255)
          gcggGGAGCCCGCGGGCATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCa                                                                                                                                         <  2:2876429/151‑1 (MQ=255)
           cggGGAGCCCGCGGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAg                                                                                                                                        <  1:385236/151‑1 (MQ=255)
            ggggAGCCCGCGGGAATCGCCGAATGGGCGCCGCAGCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAgg                                                                                                                                       >  2:5851904/1‑151 (MQ=255)
              ggAGCCCGCGGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGtc                                                                                                                                     <  1:11007134/151‑2 (MQ=255)
              ggAGCCCGCGGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGtc                                                                                                                                     <  1:10181314/151‑2 (MQ=255)
              ggAGCCCGCGGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGtc                                                                                                                                     <  1:8147993/151‑2 (MQ=255)
              ggAGCCCGCGGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGtc                                                                                                                                     <  1:2260884/151‑2 (MQ=255)
              ggAGCCCGCGGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGtc                                                                                                                                     <  2:6594826/151‑2 (MQ=255)
               gAGCCCGCGGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGtcg                                                                                                                                    >  2:11331785/1‑149 (MQ=255)
                aGCCCGCGGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGa                                                                                                                                   <  1:9969341/151‑1 (MQ=255)
                aGCCCGCGGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGa                                                                                                                                   <  2:8486186/151‑1 (MQ=255)
                aGCCCGCGGGAATCGCCGAATGGGCGCCGCATCTGCTAGACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGGTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGACCTGCAGCAggacga                                                                                                                                   >  2:4234918/1‑147 (MQ=255)
                  cccGCGGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGAcc                                                                                                                                 <  1:11421680/151‑1 (MQ=255)
                   ccGCGGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCg                                                                                                                                <  1:1433794/151‑1 (MQ=255)
                     gcgGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGgc                                                                                                                              >  2:11080442/1‑151 (MQ=255)
                      cgGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGgct                                                                                                                             >  2:2296902/1‑151 (MQ=255)
                       gggAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGgctg                                                                                                                            <  1:7979755/151‑1 (MQ=255)
                         gAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGgctgct                                                                                                                          >  1:9758704/1‑151 (MQ=255)
                         gAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGgctgct                                                                                                                          <  2:2742627/151‑1 (MQ=255)
                         gAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGgctgct                                                                                                                          >  2:9020369/1‑151 (MQ=255)
                              gCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTggagg                                                                                                                     <  2:5164452/151‑1 (MQ=255)
                              gCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTggagg                                                                                                                     >  2:5207314/1‑151 (MQ=255)
                              gCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCAGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTggagg                                                                                                                     <  1:1069520/151‑1 (MQ=255)
                               ccGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTggagga                                                                                                                    >  2:3585551/1‑151 (MQ=255)
                                cGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGaa                                                                                                                   >  1:8856222/1‑151 (MQ=255)
                                cGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGaa                                                                                                                   >  2:9889525/1‑151 (MQ=255)
                                  aaTGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGaaaa                                                                                                                 >  2:6800305/1‑151 (MQ=255)
                                  aaTGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGaaaa                                                                                                                 <  1:3222564/151‑1 (MQ=255)
                                   aTGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAAc                                                                                                                >  1:9435475/1‑151 (MQ=255)
                                   aTGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAAc                                                                                                                >  2:11767146/1‑151 (MQ=255)
                                    tGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAAcc                                                                                                               <  1:8580287/151‑1 (MQ=255)
                                     gggCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACct                                                                                                              <  1:5297653/151‑1 (MQ=255)
                                      ggCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACctg                                                                                                             <  1:3746761/151‑1 (MQ=255)
                                       gcgcCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACctgc                                                                                                            <  1:565129/151‑1 (MQ=255)
                                         gccgcATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACctgctg                                                                                                          <  2:11630338/151‑1 (MQ=255)
                                               tCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGcc                                                                                                    >  1:6142818/1‑151 (MQ=255)
                                               tCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGcc                                                                                                    <  2:8637017/151‑1 (MQ=255)
                                                cTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCg                                                                                                   <  1:5481651/151‑1 (MQ=255)
                                                     ctACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAgc                                                                                              >  2:9310448/1‑151 (MQ=255)
                                                       aCCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAgccg                                                                                            >  1:6499191/1‑151 (MQ=255)
                                                        ccccGGTCGCCAGCGCGCCAGCCTGGTGCTCCAGCCACGTCTCGACCAGTTGCTCCACGGCTCCATCCGCAAGCCGCATCACCCGCCGGCCCATGGCCTGCAGCAGGTCGACCGGCTGCTGGAGGAACACCTGCTGCAGGCCGAAAgccgc                                                                                           <  2:4715376/151‑1 (MQ=255)
                                                        ccccGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAgccgc                                                                                           <  2:1945535/151‑1 (MQ=255)
                                                          ccGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGccc                                                                                         <  1:4237238/151‑1 (MQ=255)
                                                            ggTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCgg                                                                                       <  2:5629162/151‑1 (MQ=255)
                                                            ggTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCGGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCgg                                                                                       >  1:3078984/1‑151 (MQ=255)
                                                              tCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCg                                                                                     >  1:7058767/1‑151 (MQ=255)
                                                                gCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGtt                                                                                   <  1:10496817/151‑1 (MQ=255)
                                                                 ccAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTg                                                                                  <  1:5745805/151‑1 (MQ=255)
                                                                 ccAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTg                                                                                  >  1:4049437/1‑151 (MQ=255)
                                                                 ccAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTg                                                                                  >  2:10500921/1‑151 (MQ=255)
                                                                  cAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGc                                                                                 >  2:961173/1‑151 (MQ=255)
                                                                       cgcCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCtgatg                                                                            <  1:10260970/151‑1 (MQ=255)
                                                                       cgcCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGGTGCtgatg                                                                            >  1:10019433/1‑151 (MQ=255)
                                                                        gcCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCtgatga                                                                           <  2:10190982/151‑1 (MQ=255)
                                                                        gcCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCtgatga                                                                           >  1:6914604/1‑151 (MQ=255)
                                                                        gcCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCtgatga                                                                           >  2:171821/1‑151 (MQ=255)
                                                                          cAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCAGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCtgatgatg                                                                         >  2:8778578/1‑151 (MQ=255)
                                                                           aGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCtgatgatga                                                                        <  1:3638170/151‑1 (MQ=255)
                                                                           aGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCtgatgatga                                                                        <  1:5914651/151‑1 (MQ=255)
                                                                              cTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAgcg                                                                     <  1:5854122/151‑1 (MQ=255)
                                                                                ggTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAgcggc                                                                   <  2:8893821/151‑1 (MQ=255)
                                                                                    cTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGAcc                                                               >  2:3957839/1‑151 (MQ=255)
                                                                                     tCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGGTGCTGATGATGAGCGGCGACCa                                                              >  1:3901584/1‑151 (MQ=255)
                                                                                       aaGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAgg                                                            >  2:7222180/1‑151 (MQ=255)
                                                                                        aGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAATCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGt                                                           <  1:4698764/151‑1 (MQ=255)
                                                                                        aGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGGGGCTGATGATGAGCGGCGACCAGGt                                                           >  2:7851505/1‑151 (MQ=255)
                                                                                          ccACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCt                                                         >  2:2067514/1‑151 (MQ=255)
                                                                                          ccACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCt                                                         >  1:3493539/1‑151 (MQ=255)
                                                                                           cACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTa                                                        >  2:4989675/1‑151 (MQ=255)
                                                                                            aCGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTAc                                                       <  2:5404945/151‑1 (MQ=255)
                                                                                                ctcGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCcg                                                   >  2:10974759/1‑151 (MQ=255)
                                                                                                 tcGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCcga                                                  >  1:4368126/1‑151 (MQ=255)
                                                                                                     ccAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCgacgac                                              >  1:10819715/1‑151 (MQ=255)
                                                                                                           gCTCCACGGCTCCTGCCGCAAGCCGCCTCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGcc                                        <  2:1444890/151‑1 (MQ=255)
                                                                                                             tCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCgg                                      >  1:9571153/1‑151 (MQ=255)
                                                                                                             tCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCgg                                      >  2:2361951/1‑151 (MQ=255)
                                                                                                             tCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCgg                                      <  2:7452064/151‑1 (MQ=255)
                                                                                                                aCGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGcc                                   >  1:7677632/1‑151 (MQ=255)
                                                                                                                aCGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGcc                                   >  2:7662916/1‑151 (MQ=255)
                                                                                                                    cTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGTTg                               >  2:4549022/1‑151 (MQ=255)
                                                                                                                    cTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATg                               <  1:6675515/151‑1 (MQ=255)
                                                                                                                        tGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCtgg                           <  1:4555028/151‑1 (MQ=255)
                                                                                                                         gCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCtggt                          >  2:6200554/1‑151 (MQ=255)
                                                                                                                           cGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCtggtgg                        <  1:8315409/151‑1 (MQ=255)
                                                                                                                            gCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCTGGTGGc                       >  1:720694/1‑151 (MQ=255)
                                                                                                                              aaGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCGGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCTGGTGGcca                     >  2:11564147/1‑151 (MQ=255)
                                                                                                                               aGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCTGGTGGccat                    >  2:3435666/1‑151 (MQ=255)
                                                                                                                                   gCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCTGGTGGccatccat                >  2:4343201/1‑151 (MQ=255)
                                                                                                                                    cATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCTGGTGGCCATCCATg               <  1:2608492/151‑1 (MQ=255)
                                                                                                                                          cccggcggcCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCTGGTGGCCATCCATGGGCTg          >  2:6929335/4‑151 (MQ=255)
                                                                                                                                          cccggcggcCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCTGGTGGCCATCCATGGGCTg          >  2:8314214/4‑151 (MQ=255)
                                                                                                                                            cggcggcCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCTGGTGGCCATCCATGGGCTGAt        <  1:2803007/150‑1 (MQ=255)
                                                                                                                                                 gcCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCTGGTGGCCATCCATGGGCTGATccgcc   <  1:2504196/151‑1 (MQ=255)
                                                                                                                                                  ccGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCTGGTGGCCATCCATGGGCTGATccgccg  <  2:10750903/151‑1 (MQ=255)
                                                                                                                                                  ccGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCTGGTGGCCATCCATGGGCTGATccgccg  <  1:1912434/151‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                  
GCCGCCGACGGCGGGGAGCCCGCGGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGGCGGCCGATGGGCCTGCTGCAGGTAGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCTGGTGGCCATCCATGGGCTGATCCGCCG  >  NC_008463/3390332‑3390629

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: