Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 3,566,665 Δ1 bp pseudogene (1016/1275 nt) PA14_RS16215 → fatty acid desaturase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084633,566,6630C.98.7% 2305.8 / 14.9 478pseudogene (1014/1275 nt)PA14_RS16215fatty acid desaturase
Reads supporting (aligned to +/- strand):  ref base C (1/5);  new base . (231/240);  total (232/246)
Fisher's exact test for biased strand distribution p-value = 2.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 500 total aligned reads displayed.

ACGCTTCCTCGGCCGCCCCTACCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACAAGGCGCGCCCCTCAAAGTTCGGCCCGATGACCAAAACCTGGGGCATGTGGGAAAACTTCCTGATCCTGCGCAACCACCTCTGCCGCGGCCAGT  >  NC_008463/3566513‑3566809
                                                                                                                                                      |                                                                                                                                                  
  gCTTCCTCGGCCGCCCCTACCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCaccaccg                                                                                                                                                  >  2:4112528/1‑150 (MQ=255)
  gCTTCCTCGGCCGCCCCTACCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCaccaccg                                                                                                                                                  >  2:5772372/1‑150 (MQ=255)
    ttCCTCGGCCGCCCCTACCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCaccaccgca                                                                                                                                                >  1:4003431/1‑148 (MQ=255)
    ttCCTCGGCCGCCCCTACCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCaccaccgca                                                                                                                                                >  2:7997032/1‑148 (MQ=255)
         cGGCCGCCCCTACCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCg                                                                                                                                          <  2:417988/151‑1 (MQ=255)
         cGGCCGCCCCTACCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCg                                                                                                                                          <  2:8815251/151‑1 (MQ=255)
          ggCCGCCCCTACCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCgg                                                                                                                                         <  1:6906805/151‑1 (MQ=255)
           gccgccCCTACCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCggg                                                                                                                                        <  1:4354286/151‑1 (MQ=255)
              gccCCTACCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGt                                                                                                                                     >  2:2418087/1‑151 (MQ=255)
               ccccTACCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTg                                                                                                                                    >  1:5771417/1‑151 (MQ=255)
                 ccTACCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGaa                                                                                                                                  >  2:7492901/1‑151 (MQ=255)
                   tACCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACt                                                                                                                                >  2:4939590/1‑151 (MQ=255)
                     cccGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTc                                                                                                                              >  1:3055147/1‑151 (MQ=255)
                      ccGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCt                                                                                                                             <  1:2130636/151‑1 (MQ=255)
                      ccGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCt                                                                                                                             <  1:4141056/151‑1 (MQ=255)
                       cGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTg                                                                                                                            >  2:1983627/1‑151 (MQ=255)
                       cGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTg                                                                                                                            >  1:2300620/1‑151 (MQ=255)
                         gCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGgc                                                                                                                          <  1:1071717/151‑1 (MQ=255)
                         gCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGgc                                                                                                                          <  2:1036330/151‑1 (MQ=255)
                         gCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGgc                                                                                                                          >  2:6938359/1‑151 (MQ=255)
                          ccGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGgcg                                                                                                                         >  2:6413336/1‑151 (MQ=255)
                          ccGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGgcg                                                                                                                         <  1:505334/151‑1 (MQ=255)
                           cGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGgcgc                                                                                                                        <  1:8087294/151‑1 (MQ=255)
                           cGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGgcgc                                                                                                                        >  1:8257228/1‑151 (MQ=255)
                            gACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCa                                                                                                                       >  2:9190866/1‑151 (MQ=255)
                            gACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCa                                                                                                                       <  1:2547833/151‑1 (MQ=255)
                            gACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGGTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCa                                                                                                                       <  2:319537/151‑1 (MQ=255)
                             aCAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCAt                                                                                                                      >  2:8369740/1‑151 (MQ=255)
                                 gCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATcgcg                                                                                                                  >  1:279040/1‑151 (MQ=255)
                                   aGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGcc                                                                                                                >  1:8179798/1‑151 (MQ=255)
                                      gcctggcctTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGa                                                                                                             <  1:8557352/151‑1 (MQ=255)
                                       cctggcctTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAg                                                                                                            <  1:457215/151‑1 (MQ=255)
                                           gcctTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAgcgcg                                                                                                        >  2:6001454/1‑151 (MQ=255)
                                             ctTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCCTGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGcgcgcg                                                                                                      <  2:906085/151‑1 (MQ=255)
                                             ctTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGcgcgcg                                                                                                      <  2:5099701/151‑1 (MQ=255)
                                              ttCGCCCGGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCgg                                                                                                     <  2:7594909/151‑1 (MQ=255)
                                                   ccTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCcaac                                                                                                >  1:7172786/1‑151 (MQ=255)
                                                   ccTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTACCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCcaac                                                                                                >  1:4836353/1‑151 (MQ=255)
                                                    cTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCcaacg                                                                                               >  1:6226035/1‑150 (MQ=255)
                                                            cgttcgCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGgcgcgc                                                                                       >  2:232475/1‑151 (MQ=255)
                                                               tcgCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGcccc                                                                                    <  1:7806329/151‑1 (MQ=255)
                                                                 gCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTc                                                                                  <  2:1326069/151‑1 (MQ=255)
                                                                  cTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTcg                                                                                 <  1:8777131/151‑2 (MQ=255)
                                                                    ggAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCgaa                                                                               >  1:3687674/1‑151 (MQ=255)
                                                                        ccTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTc                                                                           <  2:3505988/151‑1 (MQ=255)
                                                                                tgctgcTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGAt                                                                   <  1:6762276/151‑1 (MQ=255)
                                                                                  ctgctgATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGa                                                                 >  2:5170374/1‑151 (MQ=255)
                                                                                   tgctgATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGAc                                                                >  2:9529302/1‑151 (MQ=255)
                                                                                   tgctgATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGAc                                                                >  1:7089158/1‑151 (MQ=255)
                                                                                    gctgATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGAcc                                                               >  2:657667/1‑151 (MQ=255)
                                                                                      tgATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGAccga                                                             <  1:1448192/151‑3 (MQ=255)
                                                                                      ggATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTGCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGAccga                                                             <  2:7276362/150‑3 (MQ=255)
                                                                                       gATCGACATCCCGCTGCGGGTGGTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGAccgag                                                            <  2:1019651/151‑4 (MQ=255)
                                                                                        aTCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGAccgaga                                                           <  1:5465726/151‑5 (MQ=255)
                                                                                         tCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGAccgagac                                                          <  2:8659560/151‑6 (MQ=255)
                                                                                          cGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGAcc                                                         >  2:3533887/1‑151 (MQ=255)
                                                                                             cATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTgg                                                      >  2:982918/1‑151 (MQ=255)
                                                                                              aTCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTggg                                                     >  1:8978596/1‑151 (MQ=255)
                                                                                                cccGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGc                                                   <  2:2945170/151‑1 (MQ=255)
                                                                                                 ccGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCTACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGct                                                  >  1:3462801/1‑150 (MQ=255)
                                                                                                   gCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCAtg                                                >  2:4586865/1‑151 (MQ=255)
                                                                                                     tGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCAtgtg                                              <  1:9711593/151‑1 (MQ=255)
                                                                                                      gCGGGTGTTCTCCTTCATGCAGGACCTGCCCCGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGCGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTgg                                             <  2:10602285/151‑1 (MQ=255)
                                                                                                       cGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTggg                                            <  2:5086128/151‑1 (MQ=255)
                                                                                                         ggTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGag                                          <  1:5306538/151‑2 (MQ=255)
                                                                                                            gtTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAAc                                       <  1:10251519/151‑1 (MQ=255)
                                                                                                             ttCTCCTTCCTGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGTGGTGATCTTCTGGCGCATCGCGCCGGTGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGTCCTGGGGCATGTGGGAGAACt                                      <  1:10448229/151‑1 (MQ=255)
                                                                                                             ttCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACt                                      >  1:10007718/1‑151 (MQ=255)
                                                                                                              tctcCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACtt                                     >  1:6503530/1‑151 (MQ=255)
                                                                                                               ctcCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTc                                    >  2:1476160/1‑151 (MQ=255)
                                                                                                                  cTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTg                                 <  2:2733741/151‑1 (MQ=255)
                                                                                                                   ttCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGa                                <  1:9363383/151‑1 (MQ=255)
                                                                                                                   ttCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGa                                >  2:2605548/1‑151 (MQ=255)
                                                                                                                    tcctGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGAt                               <  2:794051/148‑1 (MQ=255)
                                                                                                                    tCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCACCGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGAt                               <  2:3922210/151‑1 (MQ=255)
                                                                                                                    tCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGAt                               >  1:329786/1‑151 (MQ=255)
                                                                                                                     cATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATc                              <  1:5850947/151‑1 (MQ=255)
                                                                                                                      aTGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATcc                             <  2:5325107/151‑1 (MQ=255)
                                                                                                                       tGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGGACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCt                            >  2:6171267/1‑151 (MQ=255)
                                                                                                                        gCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTg                           <  2:6282110/151‑1 (MQ=255)
                                                                                                                          aGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTgcg                         <  1:10513503/151‑1 (MQ=255)
                                                                                                                          aGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTgcg                         >  2:3112234/1‑151 (MQ=255)
                                                                                                                               cTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGacc                    <  1:9585632/151‑1 (MQ=255)
                                                                                                                               cTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGacc                    >  1:6631899/1‑151 (MQ=255)
                                                                                                                               cTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGacc                    >  1:8444870/1‑151 (MQ=255)
                                                                                                                                 gCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGaccac                  >  1:566871/1‑151 (MQ=255)
                                                                                                                                 gCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGaccac                  <  1:2666904/151‑1 (MQ=255)
                                                                                                                                  cccAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGaccacc                 >  1:7401496/1‑151 (MQ=255)
                                                                                                                                  cccAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGaccacc                 <  2:8551966/151‑1 (MQ=255)
                                                                                                                                  cccAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGaccacc                 <  1:2944104/151‑1 (MQ=255)
                                                                                                                                   ccAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACct                <  2:1611402/151‑1 (MQ=255)
                                                                                                                                   ccAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACct                <  2:3766225/151‑1 (MQ=255)
                                                                                                                                      gCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTg             <  1:1244714/151‑1 (MQ=255)
                                                                                                                                       ccACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGc            >  2:8222795/1‑151 (MQ=255)
                                                                                                                                         aCGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCcg          >  2:4802559/1‑151 (MQ=255)
                                                                                                                                              ttCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGTGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGcc     <  2:282806/151‑1 (MQ=255)
                                                                                                                                               tCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCa    <  1:1818207/151‑1 (MQ=255)
                                                                                                                                               tCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCa    <  2:7723678/151‑1 (MQ=255)
                                                                                                                                               tCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCa    <  2:10056384/151‑1 (MQ=255)
                                                                                                                                                ccaccaCCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAg   <  1:7898171/151‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
ACGCTTCCTCGGCCGCCCCTACCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACAAGGCGCGCCCCTCAAAGTTCGGCCCGATGACCAAAACCTGGGGCATGTGGGAAAACTTCCTGATCCTGCGCAACCACCTCTGCCGCGGCCAGT  >  NC_008463/3566513‑3566809

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: