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.8% 2351.3 / 15.4 485pseudogene (1014/1275 nt)PA14_RS16215fatty acid desaturase
Reads supporting (aligned to +/- strand):  ref base C (0/6);  new base . (266/213);  total (266/219)
Fisher's exact test for biased strand distribution p-value = 8.16e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 510 total aligned reads displayed.

CGCTTCCTCGGCCGCCCCTACCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACAAGGCGCGCCCCTCAAAGTTCGGCCCGATGACCAAAACCTGGGGCATGTGGGAAAACTTCCTGATCCTGCGCAACCACCTCTGCCGCGGCCAG  >  NC_008463/3566514‑3566808
                                                                                                                                                     |                                                                                                                                                 
      cTCGGCCGCCCCTACCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAg                                                                                                                                            <  1:10096268/150‑1 (MQ=255)
        cGGCCGCCCCTACCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGcc                                                                                                                                          >  1:4716102/1‑150 (MQ=255)
          gccgccCCTACCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCgg                                                                                                                                        >  2:3642967/1‑150 (MQ=255)
            cgccCCTACCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCgggg                                                                                                                                      <  2:2050811/150‑1 (MQ=255)
            cgccCCTACCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCgggg                                                                                                                                      <  1:8414508/150‑1 (MQ=255)
            cgccCCTACCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCgggg                                                                                                                                      <  2:6507080/150‑1 (MQ=255)
             gccCCTACCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCggggg                                                                                                                                     >  1:9469005/1‑150 (MQ=255)
               cccTACCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTg                                                                                                                                   <  2:10821648/150‑1 (MQ=255)
               cccTACCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTg                                                                                                                                   <  1:2968245/150‑1 (MQ=255)
                  tACCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCACGCCGGGGGTGAAc                                                                                                                                >  1:2242502/1‑150 (MQ=255)
                  tACCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAAc                                                                                                                                >  1:11025585/1‑150 (MQ=255)
                   aCCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACt                                                                                                                               >  1:6598198/1‑150 (MQ=255)
                      cGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCt                                                                                                                            >  1:4234808/1‑150 (MQ=255)
                       ggCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCTTGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTg                                                                                                                           >  2:9463749/1‑150 (MQ=255)
                       ggCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGACCTTCTg                                                                                                                           >  2:2444007/1‑150 (MQ=255)
                       ggCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTg                                                                                                                           >  1:3673116/1‑150 (MQ=255)
                         ccGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGgc                                                                                                                         <  1:313474/150‑1 (MQ=255)
                             ccaGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCAt                                                                                                                     <  2:224745/148‑1 (MQ=255)
                              aaGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATc                                                                                                                    >  1:52984/1‑150 (MQ=255)
                                 cAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGACCTTCTGGCGCATcgcg                                                                                                                 >  1:2756034/1‑150 (MQ=255)
                                   gggCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGcc                                                                                                               >  1:9192725/1‑150 (MQ=255)
                                   gggCCTGGCCTTCGCCCTGGCGTTCGTGCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCCTGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGcc                                                                                                               <  2:9093164/150‑1 (MQ=255)
                                    ggcctggcctTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCg                                                                                                              >  2:1700874/1‑150 (MQ=255)
                                    ggcctggcctTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACGTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCg                                                                                                              >  2:8722565/1‑150 (MQ=255)
                                       ctggcctTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAg                                                                                                           <  1:2853730/150‑1 (MQ=255)
                                       ctggcctTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAg                                                                                                           >  2:10192090/1‑150 (MQ=255)
                                       ctggcctTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAg                                                                                                           >  2:6982511/1‑150 (MQ=255)
                                       ctggcctTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCCCCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAg                                                                                                           <  1:1807111/150‑1 (MQ=255)
                                        tggcctTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAgc                                                                                                          <  1:9757372/150‑1 (MQ=255)
                                             ttCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGcgcgcg                                                                                                     >  2:7226433/1‑150 (MQ=255)
                                                gCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCCCCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGcc                                                                                                  <  2:8161492/150‑1 (MQ=255)
                                                  ccTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCcaa                                                                                                >  2:4719449/1‑150 (MQ=255)
                                                  ccTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCcaa                                                                                                <  2:8974466/150‑1 (MQ=255)
                                                        gttcgttcGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGgc                                                                                          <  1:1077747/150‑1 (MQ=255)
                                                           cgttcgCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGgcgcg                                                                                       <  1:7659945/150‑1 (MQ=255)
                                                             ttcgCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGcc                                                                                     >  1:7244219/1‑150 (MQ=255)
                                                             ttcgCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGcc                                                                                     >  2:5008715/1‑150 (MQ=255)
                                                              tcgCTGGAGCCTGGGCGTGCTGCTGTTCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGccc                                                                                    <  1:3502940/150‑1 (MQ=255)
                                                               cgCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGcccc                                                                                   <  1:2040827/150‑1 (MQ=255)
                                                                gCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCt                                                                                  >  1:5447237/1‑150 (MQ=255)
                                                                 cTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTc                                                                                 <  2:5475398/150‑1 (MQ=255)
                                                                  tGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTcg                                                                                >  2:4616200/1‑149 (MQ=255)
                                                                   ggAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGTGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTcga                                                                               <  2:3027317/150‑3 (MQ=255)
                                                                   ggAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGCTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTcga                                                                               <  1:4412152/150‑3 (MQ=255)
                                                                    gAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCgaa                                                                              >  2:370530/1‑150 (MQ=255)
                                                                        cTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTc                                                                          >  1:1427496/1‑150 (MQ=255)
                                                                        cTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTc                                                                          <  2:5289342/150‑1 (MQ=255)
                                                                        cTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTc                                                                          <  2:3831722/150‑1 (MQ=255)
                                                                        cTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTc                                                                          <  1:4845160/150‑1 (MQ=255)
                                                                         tGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCg                                                                         >  2:10693575/1‑150 (MQ=255)
                                                                         tGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCg                                                                         >  1:8137966/1‑150 (MQ=255)
                                                                         tGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCg                                                                         >  2:2491830/1‑150 (MQ=255)
                                                                             cGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGccc                                                                     <  2:2639120/150‑1 (MQ=255)
                                                                              gTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCg                                                                    >  1:9909226/1‑150 (MQ=255)
                                                                              gTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCg                                                                    >  1:6358257/1‑150 (MQ=255)
                                                                               tgctgcTGATCGACATCCCGCTGCGGGTGTTCTCCTTCCTGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGa                                                                   <  1:5819917/150‑1 (MQ=255)
                                                                               tgctgcTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGa                                                                   <  1:1554860/150‑1 (MQ=255)
                                                                                gctgctGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGAt                                                                  >  1:4157279/1‑150 (MQ=255)
                                                                                 ctgctgATCGACATCCCGCTGCGGGTGTTCTCCTTCCTGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATg                                                                 <  2:7559414/150‑1 (MQ=255)
                                                                                 ctgctgATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATg                                                                 >  2:1765610/1‑150 (MQ=255)
                                                                                  tgctgATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGa                                                                <  1:9361902/150‑1 (MQ=255)
                                                                                    ctgATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGTCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGAcc                                                              >  2:1514218/1‑150 (MQ=255)
                                                                                    ctgATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGAcc                                                              <  1:2536174/150‑1 (MQ=255)
                                                                                      gATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGAccga                                                            <  1:3840120/150‑3 (MQ=255)
                                                                                           aCATCCCGCTGCGGGTGTTCTCCTTCCTGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCt                                                       <  2:5965835/150‑1 (MQ=255)
                                                                                           aCATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCt                                                       >  2:978069/1‑150 (MQ=255)
                                                                                             aTCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTgg                                                     >  1:9020897/1‑150 (MQ=255)
                                                                                               cccGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTgggg                                                   >  1:3411892/1‑150 (MQ=255)
                                                                                                 cGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCa                                                 >  2:923987/1‑150 (MQ=255)
                                                                                                  gCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCAt                                                >  1:10321909/1‑150 (MQ=255)
                                                                                                  gCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCAt                                                >  2:1874659/1‑150 (MQ=255)
                                                                                                   cTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCAtg                                               >  2:7830422/1‑150 (MQ=255)
                                                                                                   cTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCAtg                                               >  1:6286696/1‑150 (MQ=255)
                                                                                                    tGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCAtgt                                              <  1:693932/150‑1 (MQ=255)
                                                                                                       gggTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTggg                                           >  1:7538918/1‑150 (MQ=255)
                                                                                                        ggTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGa                                          >  2:4211629/1‑150 (MQ=255)
                                                                                                        ggTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGa                                          >  2:4882957/1‑150 (MQ=255)
                                                                                                          tgtTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGaga                                        >  1:1055996/1‑148 (MQ=255)
                                                                                                            ttCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAAc                                      <  2:5775630/150‑1 (MQ=255)
                                                                                                               tcCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTc                                   >  2:10954306/1‑150 (MQ=255)
                                                                                                               tcCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTc                                   >  2:1670901/1‑150 (MQ=255)
                                                                                                                ccgtCAAGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGCGACCTGGGGCATGTGGGAGAACTGcc                                  <  2:1055735/147‑1 (MQ=255)
                                                                                                                ccTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTcc                                  >  1:684848/1‑150 (MQ=255)
                                                                                                                   tCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGa                               <  2:4016131/150‑1 (MQ=255)
                                                                                                                    cATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGAt                              <  2:7300536/150‑1 (MQ=255)
                                                                                                                     aTGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATc                             >  2:3277497/1‑150 (MQ=255)
                                                                                                                     aTGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATc                             <  2:6874492/150‑1 (MQ=255)
                                                                                                                         aGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTgc                         >  1:872956/1‑150 (MQ=255)
                                                                                                                            aCCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTgcgcg                      >  1:4030445/1‑149 (MQ=255)
                                                                                                                            aCCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTgcgcg                      >  2:6746214/1‑149 (MQ=255)
                                                                                                                             ccTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTgcgcga                     <  2:6277755/150‑3 (MQ=255)
                                                                                                                                gCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGACCTTCCTGATCCTGCGCGacca                  >  1:6749179/1‑150 (MQ=255)
                                                                                                                                gCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGacca                  >  1:1220565/1‑150 (MQ=255)
                                                                                                                                  ccAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCAGGCGCATCGCGCCGGAGCGCGCTGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTACCTGATCCTGCGCGaccacc                >  2:1400842/1‑150 (MQ=255)
                                                                                                                                         cGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCcg         >  2:395743/1‑150 (MQ=255)
                                                                                                                                          gACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCcgc        <  1:10682764/150‑1 (MQ=255)
                                                                                                                                            cTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCgg      >  1:5997785/1‑150 (MQ=255)
                                                                                                                                             ttCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGc     <  1:9703010/150‑1 (MQ=255)
                                                                                                                                              tCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGcc    >  2:8555383/1‑150 (MQ=255)
                                                                                                                                               ccaccaCCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCa   <  2:8362040/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
CGCTTCCTCGGCCGCCCCTACCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACAAGGCGCGCCCCTCAAAGTTCGGCCCGATGACCAAAACCTGGGGCATGTGGGAAAACTTCCTGATCCTGCGCAACCACCTCTGCCGCGGCCAG  >  NC_008463/3566514‑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: