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.2% 2336.4 / 27.4 492pseudogene (1014/1275 nt)PA14_RS16215fatty acid desaturase
Reads supporting (aligned to +/- strand):  ref base C (1/8);  new base . (225/256);  total (227/265)
Fisher's exact test for biased strand distribution p-value = 4.24e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 519 total aligned reads displayed.

CGCTTCCTCGGCCGCCCCTACCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACAAGGCGCGCCCCTCAAAGTTCGGCCCGATGACCAAAACCTGGGGCATGTGGGAAAACTTCCTGATCCTGCGCAACCACCTCTGCCGCGGCCAG  >  NC_008463/3566514‑3566808
                                                                                                                                                     |                                                                                                                                                 
  cTTCCTCGGCCGCCCCTACCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCaccaccg                                                                                                                                                 >  2:9016518/1‑149 (MQ=255)
  cTTCCTCGGCCGCCCCTACCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCaccaccg                                                                                                                                                 >  2:10115885/1‑149 (MQ=255)
     ccTCGGCCGCCCCTACCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCaccaccgcaa                                                                                                                                              <  1:2913434/150‑5 (MQ=255)
      cTCGGCCGCCCCTACCCGGCCGACAATCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAg                                                                                                                                            >  2:775568/1‑150 (MQ=255)
       tCGGCCGCCCCTACCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGc                                                                                                                                           <  1:211167/150‑1 (MQ=255)
       tCGGCCGCCCCTACCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGc                                                                                                                                           <  1:10019816/150‑1 (MQ=255)
        cGGCCGCCCCTACCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGcc                                                                                                                                          >  2:9712745/1‑150 (MQ=255)
          gccgccCCTACCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCgg                                                                                                                                        >  2:5116640/1‑150 (MQ=255)
          gccgccCCTACCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCgg                                                                                                                                        >  1:7969502/1‑150 (MQ=255)
          gccgccCCTACCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGAAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCgg                                                                                                                                        >  2:5465076/1‑150 (MQ=255)
           ccgccCCTACCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCggg                                                                                                                                       <  2:6391383/150‑1 (MQ=255)
            cgccCCTACCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCgggg                                                                                                                                      >  2:4294025/1‑150 (MQ=255)
              ccccTACCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGt                                                                                                                                    >  2:9290217/1‑150 (MQ=255)
              ccccTACCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGt                                                                                                                                    <  2:8857739/150‑1 (MQ=255)
               cccTACCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTg                                                                                                                                   >  2:10563442/1‑150 (MQ=255)
                ccTACCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGa                                                                                                                                  >  2:699226/1‑150 (MQ=255)
                   aCCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACt                                                                                                                               <  2:712749/150‑1 (MQ=255)
                   aCCCGGCCGACAAGCAGGGCCTGGCCGTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACt                                                                                                                               <  1:3153870/150‑1 (MQ=255)
                     ccGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTc                                                                                                                             <  1:3923520/150‑1 (MQ=255)
                           gACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGgcgc                                                                                                                       >  1:5946631/1‑150 (MQ=255)
                               aGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATcg                                                                                                                   <  1:9065893/150‑1 (MQ=255)
                                 cAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATcgcg                                                                                                                 >  2:2369069/1‑150 (MQ=255)
                                  aGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCgcgc                                                                                                                <  1:9920044/150‑1 (MQ=255)
                                   gggCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGcc                                                                                                               >  1:2565285/1‑150 (MQ=255)
                                   gggCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGAGCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGcc                                                                                                               >  2:3818734/1‑150 (MQ=255)
                                     gcctggcctTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCAACGCGCCgg                                                                                                             <  1:9560935/150‑1 (MQ=255)
                                       ctggcctTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAg                                                                                                           <  2:1866927/150‑1 (MQ=255)
                                         ggcctTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCACGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAgcg                                                                                                         >  2:7527113/1‑150 (MQ=255)
                                                gCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGcc                                                                                                  >  1:7753696/1‑150 (MQ=255)
                                                gCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGACGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGcc                                                                                                  >  1:10879811/1‑150 (MQ=255)
                                                  ccTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCcaa                                                                                                <  1:4624314/150‑1 (MQ=255)
                                                     ggCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCcaacga                                                                                             >  2:1850333/1‑148 (MQ=255)
                                                      gCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAg                                                                                            >  2:3302477/1‑150 (MQ=255)
                                                       cgttcgttCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAgg                                                                                           <  1:6859912/150‑1 (MQ=255)
                                                       cgttcgttCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAgg                                                                                           >  1:5143222/1‑150 (MQ=255)
                                                        gttcgttcGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGgc                                                                                          >  2:1450329/1‑150 (MQ=255)
                                                         ttcgttcgCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGgcg                                                                                         <  1:1982449/150‑1 (MQ=255)
                                                         ttcgttcgCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGgcg                                                                                         >  2:8754976/1‑150 (MQ=255)
                                                          tcgttcgCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGgcgc                                                                                        <  1:7195066/150‑1 (MQ=255)
                                                          tcgttcgCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCGTCATGCAGGACCTGCCCAGCCACGACTTCCACCCCCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGgcgc                                                                                        <  1:3785954/150‑1 (MQ=255)
                                                            ggtcgCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAACCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGgcgcgc                                                                                      <  2:8568595/148‑1 (MQ=255)
                                                             ttcgCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGcc                                                                                     >  2:4609435/1‑150 (MQ=255)
                                                               cgCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGcccc                                                                                   <  2:2118341/150‑1 (MQ=255)
                                                                gCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCt                                                                                  >  2:2630786/1‑150 (MQ=255)
                                                                gCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCt                                                                                  >  2:2839822/1‑150 (MQ=255)
                                                                 cTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTc                                                                                 >  1:7258735/1‑150 (MQ=255)
                                                                        cTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCACGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTc                                                                          <  2:9613160/150‑1 (MQ=255)
                                                                        cTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTc                                                                          <  1:9213754/150‑1 (MQ=255)
                                                                         tGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCg                                                                         >  2:10897480/1‑150 (MQ=255)
                                                                          gggCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCgg                                                                        >  2:8440919/1‑150 (MQ=255)
                                                                           ggCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGc                                                                       >  2:5350979/1‑150 (MQ=255)
                                                                           ggCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGc                                                                       >  1:9768784/1‑150 (MQ=255)
                                                                           ggCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGc                                                                       <  1:6559229/150‑1 (MQ=255)
                                                                            gCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGcc                                                                      <  2:5592233/150‑1 (MQ=255)
                                                                                gctgctGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGAt                                                                  >  2:8845620/1‑150 (MQ=255)
                                                                                gctgctGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGAt                                                                  <  1:2798186/150‑1 (MQ=255)
                                                                                  tgctgATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGa                                                                <  2:3616127/150‑1 (MQ=255)
                                                                                    ctgATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGAcc                                                              >  1:5737161/1‑150 (MQ=255)
                                                                                     tgATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGAccg                                                             <  1:6619160/150‑2 (MQ=255)
                                                                                      gATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGAccga                                                            >  2:4077682/1‑148 (MQ=255)
                                                                                          gACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGAcc                                                        <  2:1486665/150‑1 (MQ=255)
                                                                                          gACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGAcc                                                        <  1:3729114/150‑1 (MQ=255)
                                                                                           aCATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCt                                                       <  1:5317565/150‑1 (MQ=255)
                                                                                             aTCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTgg                                                     >  1:6394229/1‑150 (MQ=255)
                                                                                                 cGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCa                                                 <  2:6732008/150‑1 (MQ=255)
                                                                                                 cGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCa                                                 >  2:1650239/1‑150 (MQ=255)
                                                                                                  gCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCAt                                                >  1:1289990/1‑150 (MQ=255)
                                                                                                   cTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCCCGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCAtg                                               >  2:5755804/1‑150 (MQ=255)
                                                                                                    tGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCAtgt                                              <  1:4053233/150‑1 (MQ=255)
                                                                                                    tGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCAtgt                                              <  2:9449415/150‑1 (MQ=255)
                                                                                                       gggTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTggg                                           >  1:2283979/1‑150 (MQ=255)
                                                                                                         gtgtTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGag                                         >  2:6164712/1‑149 (MQ=255)
                                                                                                          tgtTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGaga                                        <  2:5908895/150‑3 (MQ=255)
                                                                                                                ccgtCATGCAGGACCTGCCCAGCCACGACTTCCACCACCCCCAGCCGGGGGTGAACTTCAGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGAGGACCGAGACCTGGGGCATGTGGGAGAACTTcc                                  <  2:8105890/147‑1 (MQ=255)
                                                                                                                ccTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCACGCCGGGGGAGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTcc                                  <  2:7310043/150‑1 (MQ=255)
                                                                                                                ccTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTcc                                  >  1:11032270/1‑150 (MQ=255)
                                                                                                                   tCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGa                               <  2:1148563/150‑1 (MQ=255)
                                                                                                                    cctGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGAt                              <  2:1263034/148‑1 (MQ=255)
                                                                                                                    cATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGAt                              <  2:10307421/150‑1 (MQ=255)
                                                                                                                    cATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGAt                              <  1:8325789/150‑1 (MQ=255)
                                                                                                                    cATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGAt                              >  1:4376114/1‑150 (MQ=255)
                                                                                                                     aTGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATc                             >  1:8858924/1‑150 (MQ=255)
                                                                                                                      tGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATcc                            <  1:4892048/150‑1 (MQ=255)
                                                                                                                        cAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTg                          <  1:1705929/150‑1 (MQ=255)
                                                                                                                         aGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTgc                         <  1:1548056/150‑1 (MQ=255)
                                                                                                                         aGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTgc                         <  2:4778187/150‑1 (MQ=255)
                                                                                                                           gACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTgcgc                       >  1:3595026/1‑150 (MQ=255)
                                                                                                                            aCCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTgcgcg                      >  1:2031078/1‑149 (MQ=255)
                                                                                                                                gCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGacca                  <  1:694817/150‑1 (MQ=255)
                                                                                                                                gCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGacca                  <  1:260253/150‑1 (MQ=255)
                                                                                                                                gCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGacca                  >  2:104005/1‑150 (MQ=255)
                                                                                                                                   cAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACct               <  2:610893/150‑1 (MQ=255)
                                                                                                                                    aGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACctc              <  2:4521194/150‑1 (MQ=255)
                                                                                                                                      ccACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTg            >  1:8982286/1‑150 (MQ=255)
                                                                                                                                       cACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGc           <  1:6195559/150‑1 (MQ=255)
                                                                                                                                        aCGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGcc          <  2:4676970/150‑1 (MQ=255)
                                                                                                                                          gACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCcgc        <  1:10545012/150‑1 (MQ=255)
                                                                                                                                           aCTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCcgcg       <  1:7581403/150‑1 (MQ=255)
                                                                                                                                              tCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGcc    <  2:764209/150‑1 (MQ=255)
                                                                                                                                                caccacCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCAg  >  2:1729458/1‑150 (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: