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% 2224.1 / 16.0 462pseudogene (1014/1275 nt)PA14_RS16215fatty acid desaturase
Reads supporting (aligned to +/- strand):  ref base C (0/6);  new base . (229/226);  total (230/232)
Fisher's exact test for biased strand distribution p-value = 3.03e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 492 total aligned reads displayed.

GCTTCCTCGGCCGCCCCTACCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACAAGGCGCGCCCCTCAAAGTTCGGCCCGATGACCAAAACCTGGGGCATGTGGGAAAACTTCCTGATCCTGCGCAACCACCTCTGCCGCGGCCAG  >  NC_008463/3566515‑3566808
                                                                                                                                                    |                                                                                                                                                 
gCTTCCTCGGCCGCCCCTACCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCaccacc                                                                                                                                                  >  1:7967435/1‑150 (MQ=255)
gCTTCCTCGGCCGCCCCTACCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCaccacc                                                                                                                                                  >  1:318322/1‑150 (MQ=255)
gCTTCCTCGGCCGCCCCTACCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCaccacc                                                                                                                                                  >  1:3417202/1‑150 (MQ=255)
 cTTCCTCGGCCGCCCCTACCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCaccaccg                                                                                                                                                 >  1:3682482/1‑149 (MQ=255)
  ttCCTCGGCCGCCCCTACCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCaccaccgc                                                                                                                                                <  2:902399/150‑3 (MQ=255)
   tCCTCGGCCGCCCCGACCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTACGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCaccaccgca                                                                                                                                               <  2:5354718/150‑4 (MQ=255)
    ccTCGGCCGCCCCTACCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCaccaccgcaa                                                                                                                                              >  2:1774428/1‑146 (MQ=255)
      tCGGCCGCCCCTACCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGc                                                                                                                                           >  1:8152531/1‑150 (MQ=255)
      tCGGCCGCCCCTACCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGc                                                                                                                                           <  1:4031571/150‑1 (MQ=255)
            gccCCTACCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCggggg                                                                                                                                     <  2:1659569/150‑1 (MQ=255)
            gccCCTACCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCggggg                                                                                                                                     <  2:2063498/150‑1 (MQ=255)
              cccTACCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTg                                                                                                                                   >  2:9806024/1‑150 (MQ=255)
                  aCCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACt                                                                                                                               >  2:10329280/1‑150 (MQ=255)
                  aCCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACt                                                                                                                               >  2:4958759/1‑150 (MQ=255)
                    ccGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTc                                                                                                                             <  1:10351970/150‑1 (MQ=255)
                      ggCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTg                                                                                                                           >  2:4076145/1‑150 (MQ=255)
                      ggCCGACAAGCAGGGCCTGGCCTTCGCCCCGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTg                                                                                                                           <  1:2395952/150‑1 (MQ=255)
                       gCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTgg                                                                                                                          >  1:8397747/1‑150 (MQ=255)
                        ccGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGgc                                                                                                                         >  1:2730072/1‑150 (MQ=255)
                          gACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGgcgc                                                                                                                       >  2:7360126/1‑150 (MQ=255)
                            cAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCAt                                                                                                                     <  2:331719/150‑1 (MQ=255)
                            cAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCAt                                                                                                                     >  2:2609875/1‑150 (MQ=255)
                             aaGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATc                                                                                                                    <  1:6160117/150‑1 (MQ=255)
                               gCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATcgc                                                                                                                  <  1:1415574/150‑1 (MQ=255)
                                 aGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCgcgc                                                                                                                <  1:8814445/150‑1 (MQ=255)
                                  gggCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGcc                                                                                                               >  2:8010801/1‑150 (MQ=255)
                                   ggcctggcctTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCg                                                                                                              <  2:8757471/150‑1 (MQ=255)
                                   ggcctggcctTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCg                                                                                                              >  2:8339531/1‑150 (MQ=255)
                                    gcctggcctTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCgg                                                                                                             >  2:2514273/1‑150 (MQ=255)
                                     cctggcctTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGa                                                                                                            >  1:5009136/1‑150 (MQ=255)
                                       tggcctTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAgc                                                                                                          >  2:10547206/1‑150 (MQ=255)
                                        ggcctTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACGTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAgcg                                                                                                         >  2:5775865/1‑150 (MQ=255)
                                           ctTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAgcgcgc                                                                                                      >  1:6860969/1‑150 (MQ=255)
                                           ctTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAgcgcgc                                                                                                      >  2:5195410/1‑150 (MQ=255)
                                            ttCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGcgcgcg                                                                                                     >  1:4734359/1‑150 (MQ=255)
                                             tCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGATCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCAGCGCGCCGGAGCGCGCgg                                                                                                    >  2:6827943/1‑150 (MQ=255)
                                                cccTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCca                                                                                                 >  2:6490324/1‑150 (MQ=255)
                                                    ggCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCcaacga                                                                                             <  2:4143462/150‑3 (MQ=255)
                                                      cgttcgttCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAgg                                                                                           <  1:7719640/150‑1 (MQ=255)
                                                      cgttcgttCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAgg                                                                                           <  2:4508055/150‑1 (MQ=255)
                                                           gttcgCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGgcgcgc                                                                                      >  1:6645401/1‑150 (MQ=255)
                                                            ttcgCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGcc                                                                                     >  2:3530059/1‑150 (MQ=255)
                                                              cgCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGcccc                                                                                   <  2:1078525/150‑1 (MQ=255)
                                                               gCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCt                                                                                  <  2:7900285/150‑1 (MQ=255)
                                                               gCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCt                                                                                  >  1:6210371/1‑150 (MQ=255)
                                                                caggAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTc                                                                                 >  1:10429649/3‑150 (MQ=255)
                                                                cTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTc                                                                                 <  2:1309350/150‑1 (MQ=255)
                                                                    aGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAg                                                                             >  1:4479927/1‑150 (MQ=255)
                                                                    aGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAg                                                                             >  1:8852949/1‑150 (MQ=255)
                                                                     gCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGt                                                                            >  2:6086352/1‑150 (MQ=255)
                                                                     gCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGt                                                                            <  1:5604935/150‑1 (MQ=255)
                                                                     gCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGt                                                                            <  1:8181358/150‑1 (MQ=255)
                                                                     gCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGt                                                                            >  2:8509153/1‑150 (MQ=255)
                                                                      ccTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGGTCTCCTTCATGCAGGACCTGCCCCGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGtt                                                                           >  2:2369744/1‑150 (MQ=255)
                                                                        tGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCAACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCg                                                                         >  2:965568/1‑150 (MQ=255)
                                                                        tGGGCGTGCTGCTGATCGACATCCCGCGGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCTCCGCCAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGTGGCGCGCCCCGCGAAGTGCg                                                                         >  2:7099107/1‑150 (MQ=255)
                                                                           gCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGcc                                                                      <  2:2807173/150‑1 (MQ=255)
                                                                             gTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCg                                                                    >  2:3789545/1‑150 (MQ=255)
                                                                              tgctgcTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGa                                                                   >  1:9535096/1‑150 (MQ=255)
                                                                               gctgctGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGAt                                                                  <  1:5217406/150‑1 (MQ=255)
                                                                               gctgctGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGAt                                                                  <  2:1049373/150‑1 (MQ=255)
                                                                               gctgctGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGAt                                                                  <  1:10059928/150‑1 (MQ=255)
                                                                                 tgctgATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGa                                                                >  2:3636660/1‑150 (MQ=255)
                                                                                  gctgATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGAc                                                               >  2:7589315/1‑150 (MQ=255)
                                                                                      aTCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGAccgag                                                           >  2:5548168/1‑147 (MQ=255)
                                                                                       tCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGAccgaga                                                          <  1:4338845/150‑5 (MQ=255)
                                                                                        cGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGAccgagac                                                         <  2:2141550/150‑6 (MQ=255)
                                                                                         gACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGAcc                                                        <  1:1751176/150‑1 (MQ=255)
                                                                                         gACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGAcc                                                        >  2:7001430/1‑150 (MQ=255)
                                                                                          aCATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCt                                                       >  1:5721881/1‑150 (MQ=255)
                                                                                             tCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTggg                                                    <  1:848700/150‑1 (MQ=255)
                                                                                                   tGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCAtgt                                              >  2:4007446/1‑150 (MQ=255)
                                                                                                       ggTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGa                                          >  2:5940294/1‑150 (MQ=255)
                                                                                                       ggTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGa                                          >  1:9209156/1‑150 (MQ=255)
                                                                                                            tctcCTTCATGCAGGACCTGCCCAGCCTCGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACt                                     <  2:5494874/150‑1 (MQ=255)
                                                                                                              tcCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTc                                   >  2:8233718/1‑150 (MQ=255)
                                                                                                               cctacctGCAGGACCTGCACAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTcc                                  <  2:1400355/144‑1 (MQ=255)
                                                                                                               ccTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTcc                                  <  2:9313215/150‑1 (MQ=255)
                                                                                                                   cATGCAGGACCAGCCCAGCCACGACTTCCACCACCCCACGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGGTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGAt                              <  1:6479815/150‑1 (MQ=255)
                                                                                                                      gCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCt                           >  1:9856844/1‑150 (MQ=255)
                                                                                                                      gCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCt                           >  1:2224071/1‑150 (MQ=255)
                                                                                                                          gACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTgcgc                       <  1:10093976/150‑1 (MQ=255)
                                                                                                                              tGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGacc                   <  1:9931203/150‑1 (MQ=255)
                                                                                                                              tGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGacc                   >  2:4435550/1‑150 (MQ=255)
                                                                                                                               gCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGacca                  >  1:606347/1‑150 (MQ=255)
                                                                                                                                cccAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGaccac                 <  2:4750379/150‑1 (MQ=255)
                                                                                                                                cccAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGaccac                 >  1:8684312/1‑150 (MQ=255)
                                                                                                                                 ccAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGaccacc                <  1:194633/150‑1 (MQ=255)
                                                                                                                                  cAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACct               >  2:861502/1‑150 (MQ=255)
                                                                                                                                    gCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCAGGTGGGAGAACTTCCTGATCCTGCGCGACCACctct             >  2:2737091/1‑150 (MQ=255)
                                                                                                                                      cACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGc           <  1:8989454/150‑1 (MQ=255)
                                                                                                                                      cACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGc           <  1:7581878/150‑1 (MQ=255)
                                                                                                                                       aCGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGcc          <  1:1125935/150‑1 (MQ=255)
                                                                                                                                          aCTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCcgcg       >  2:3077460/1‑150 (MQ=255)
                                                                                                                                          aCTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCcgcg       <  2:2312353/150‑1 (MQ=255)
                                                                                                                                          aCTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCcgcg       >  1:6952781/1‑150 (MQ=255)
                                                                                                                                           cTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCgg      >  1:7151442/1‑150 (MQ=255)
                                                                                                                                            ttCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGc     >  2:7199129/1‑150 (MQ=255)
                                                                                                                                             tCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGcc    <  1:2538842/150‑1 (MQ=255)
                                                                                                                                               caccacCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAg  <  2:2869716/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GCTTCCTCGGCCGCCCCTACCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACAAGGCGCGCCCCTCAAAGTTCGGCCCGATGACCAAAACCTGGGGCATGTGGGAAAACTTCCTGATCCTGCGCAACCACCTCTGCCGCGGCCAG  >  NC_008463/3566515‑3566808

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: