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% 2695.5 / 16.4 556pseudogene (1014/1275 nt)PA14_RS16215fatty acid desaturase
Reads supporting (aligned to +/- strand):  ref base C (0/7);  new base . (262/287);  total (262/294)
Fisher's exact test for biased strand distribution p-value = 1.61e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 580 total aligned reads displayed.

CGCTTCCTCGGCCGCCCCTACCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACAAGGCGCGCCCCTCAAAGTTCGGCCCGATGACCAAAACCTGGGGCATGTGGGAAAACTTCCTGATCCTGCGCAACCACCTCTGCCGCGGCCAG  >  NC_008463/3566514‑3566808
                                                                                                                                                     |                                                                                                                                                 
  cTTCCTCGGCCGCCCCTACCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCAGGCGGTCGTTCGCAGGAGCCTGGGCGTCCAGCGGAACGACATCCCCCTCCGCGTCTTGTCCTTCATGCAGGACCTCCCCCGCCACGACTACCACCAccc                                                                                                                                                 >  2:5015522/1‑150 (MQ=255)
   ttCCTCGGCCGCCCCTACCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCaccaccgc                                                                                                                                                >  2:9258331/1‑148 (MQ=255)
    tCCTCGGCCGCCCCTACCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCaccaccgca                                                                                                                                               >  2:1974374/1‑147 (MQ=255)
    tCCTCGGCCGCCCCTACCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCaccaccgca                                                                                                                                               >  1:7277132/1‑147 (MQ=255)
    tCCTCGGCCGCCCCTACCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTGCCaccaccgca                                                                                                                                               >  2:679315/1‑147 (MQ=255)
         ggCCGCCCCTACCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCg                                                                                                                                         <  1:10382883/150‑1 (MQ=255)
          gccgccCCTACCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCgg                                                                                                                                        <  1:7610522/150‑1 (MQ=255)
                ccTACCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGa                                                                                                                                  <  1:248793/150‑1 (MQ=255)
                 cTACCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGaa                                                                                                                                 <  2:8056297/150‑1 (MQ=255)
                   aCCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACt                                                                                                                               >  1:8296734/1‑150 (MQ=255)
                   aCCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACt                                                                                                                               <  1:9733196/150‑1 (MQ=255)
                   aCCCGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACt                                                                                                                               >  1:2691396/1‑150 (MQ=255)
                      cGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCt                                                                                                                            >  1:545217/1‑150 (MQ=255)
                      cGGCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCt                                                                                                                            <  2:5551832/150‑1 (MQ=255)
                        gCCGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTgg                                                                                                                          <  1:8470851/150‑1 (MQ=255)
                          cGACCAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGgcg                                                                                                                        <  1:509302/150‑1 (MQ=255)
                          cGACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGgcg                                                                                                                        <  2:9191427/150‑1 (MQ=255)
                           gACAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGgcgc                                                                                                                       >  2:4846801/1‑150 (MQ=255)
                            aCAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCa                                                                                                                      >  1:10098930/1‑150 (MQ=255)
                             cAAGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCAt                                                                                                                     >  2:1631033/1‑150 (MQ=255)
                              aaGCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATc                                                                                                                    >  2:7435909/1‑150 (MQ=255)
                                gCAGGGCCTGGCCTTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATcgc                                                                                                                  >  2:7665348/1‑150 (MQ=255)
                                         ggcctTCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAgcg                                                                                                         <  2:2685952/150‑1 (MQ=255)
                                             ttCGCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGcgcgcg                                                                                                     <  2:7948642/150‑1 (MQ=255)
                                                gCCCTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGcc                                                                                                  <  1:198020/150‑1 (MQ=255)
                                                  ccTGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCcaa                                                                                                <  2:7235160/150‑1 (MQ=255)
                                                    tGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCcaacg                                                                                              >  1:132257/1‑149 (MQ=255)
                                                    tGGCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCcaacg                                                                                              >  2:6567979/1‑149 (MQ=255)
                                                      gCGTTCGTTCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAg                                                                                            >  2:33613/1‑150 (MQ=255)
                                                       cgttcgttCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAgg                                                                                           <  1:372973/150‑1 (MQ=255)
                                                       cgttcgttCGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAgg                                                                                           <  1:5916035/150‑1 (MQ=255)
                                                        gttcgttcGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGgc                                                                                          >  1:9462123/1‑150 (MQ=255)
                                                        gttcgttcGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGgc                                                                                          <  2:8790921/150‑1 (MQ=255)
                                                        gttcgttcGCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCCGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGgc                                                                                          <  1:10248366/150‑1 (MQ=255)
                                                             ttcgCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGcc                                                                                     <  1:3936327/150‑1 (MQ=255)
                                                              tcgCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGccc                                                                                    >  2:3763549/1‑150 (MQ=255)
                                                              tcgCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGccc                                                                                    <  1:2398326/150‑1 (MQ=255)
                                                               cgCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGcccc                                                                                   >  2:2942993/1‑150 (MQ=255)
                                                                gCTGGAGCCTGGGCGTGCTGCTGCTCGACATCCCGCTGCGGGTGTTCTCCCTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGACCTTCTGGCGCATCGCGCCGGAGCGCGCGGCCACCGAGGCGCGCCCCt                                                                                  <  2:5425377/150‑1 (MQ=255)
                                                                gCTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCt                                                                                  <  2:5968580/150‑1 (MQ=255)
                                                                 cTGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTc                                                                                 >  2:5211395/1‑150 (MQ=255)
                                                                  tGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTcg                                                                                >  2:10702131/1‑149 (MQ=255)
                                                                  tGGAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTcg                                                                                >  2:5666455/1‑149 (MQ=255)
                                                                    gAGCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCgaa                                                                              >  2:10046296/1‑150 (MQ=255)
                                                                      gCCTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGt                                                                            <  2:10310109/150‑1 (MQ=255)
                                                                       ccTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGtt                                                                           <  2:7773156/150‑1 (MQ=255)
                                                                       ccTGGGCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGtt                                                                           <  1:4176995/150‑1 (MQ=255)
                                                                           ggCGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGc                                                                       <  1:3343027/150‑1 (MQ=255)
                                                                             cGTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGccc                                                                     <  2:3232488/150‑1 (MQ=255)
                                                                              gTGCTGCTGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCg                                                                    <  2:2319781/150‑1 (MQ=255)
                                                                                gctgctGATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGAt                                                                  <  1:7839710/150‑1 (MQ=255)
                                                                                gctgctGATCGACATCCCGCGGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGAt                                                                  >  2:698275/1‑150 (MQ=255)
                                                                                  tgctgATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGa                                                                <  1:4683962/150‑1 (MQ=255)
                                                                                   gctgATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGAc                                                               <  2:9013814/150‑1 (MQ=255)
                                                                                      gATCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGAccga                                                            <  2:1240381/150‑3 (MQ=255)
                                                                                        tCGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGAccgaga                                                          >  1:924008/1‑146 (MQ=255)
                                                                                         cGACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGAccgagac                                                         >  2:8518966/1‑145 (MQ=255)
                                                                                          gACATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGAcc                                                        >  1:2920603/1‑150 (MQ=255)
                                                                                           aCATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCt                                                       <  2:10917694/150‑1 (MQ=255)
                                                                                            cATCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTg                                                      <  1:10535208/150‑1 (MQ=255)
                                                                                              tCCCGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTggg                                                    >  1:8801015/1‑150 (MQ=255)
                                                                                               cccGCTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTgggg                                                   >  1:1492808/1‑150 (MQ=255)
                                                                                                   cTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCAtg                                               >  1:4406489/1‑150 (MQ=255)
                                                                                                   cTGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCAtg                                               <  1:504930/150‑1 (MQ=255)
                                                                                                    tGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCAtgt                                              >  2:7526263/1‑150 (MQ=255)
                                                                                                    tGCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCAtgt                                              <  2:1038842/150‑1 (MQ=255)
                                                                                                     gCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCAtgtg                                             <  2:8318682/150‑1 (MQ=255)
                                                                                                     gCGGGTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCAtgtg                                             >  2:6329313/1‑150 (MQ=255)
                                                                                                        ggTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCACGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGa                                          <  1:6885384/150‑1 (MQ=255)
                                                                                                        ggTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGa                                          <  2:4072752/150‑1 (MQ=255)
                                                                                                        ggTGTTCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGa                                          >  1:1081741/1‑150 (MQ=255)
                                                                                                            ttCTCCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAAc                                      <  1:4766721/150‑1 (MQ=255)
                                                                                                              ctcCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACtt                                    >  1:1828550/1‑150 (MQ=255)
                                                                                                               tcCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTc                                   <  1:11057020/150‑1 (MQ=255)
                                                                                                               tcCTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTc                                   >  1:6136954/1‑150 (MQ=255)
                                                                                                                ccTTCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCCAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTcc                                  <  2:3578706/150‑1 (MQ=255)
                                                                                                                  ttCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTg                                >  1:10874521/1‑150 (MQ=255)
                                                                                                                  ttCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTg                                <  1:5514503/150‑1 (MQ=255)
                                                                                                                   tCATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGa                               >  1:6391955/1‑150 (MQ=255)
                                                                                                                    cATGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGAt                              <  1:6577282/150‑1 (MQ=255)
                                                                                                                     aTGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGACCTTCCTGATc                             >  2:9468962/1‑150 (MQ=255)
                                                                                                                     aTGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATc                             <  1:5684457/150‑1 (MQ=255)
                                                                                                                      tGCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATcc                            <  2:9752319/150‑1 (MQ=255)
                                                                                                                       gCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCt                           <  2:9902099/150‑1 (MQ=255)
                                                                                                                       gCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCt                           <  1:3452145/150‑1 (MQ=255)
                                                                                                                       gCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCt                           >  1:491169/1‑150 (MQ=255)
                                                                                                                       gCAGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCt                           >  1:6741610/1‑150 (MQ=255)
                                                                                                                         aGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTgc                         >  1:3199028/1‑150 (MQ=255)
                                                                                                                         aGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTgc                         <  2:4667642/150‑1 (MQ=255)
                                                                                                                         aGGACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTgc                         >  2:2465371/1‑150 (MQ=255)
                                                                                                                           gACCTGCCCAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTgcgc                       >  1:9119079/1‑150 (MQ=255)
                                                                                                                                 cccAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAAGGTGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGCTCCTGCGCGaccac                 >  2:8415418/1‑150 (MQ=255)
                                                                                                                                  ccAGCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGaccacc                >  2:4377643/1‑150 (MQ=255)
                                                                                                                                     gCCACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACctct             <  2:1119412/150‑1 (MQ=255)
                                                                                                                                      ccACGACTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGAGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTg            <  2:4209417/150‑1 (MQ=255)
                                                                                                                                            cTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCgg      >  1:2072214/1‑150 (MQ=255)
                                                                                                                                            cTTCCACCACCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCgg      >  1:535566/1‑150 (MQ=255)
                                                                                                                                               ccaccaCCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCa   <  1:2986027/150‑1 (MQ=255)
                                                                                                                                               ccaccaCCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCa   <  2:1352814/150‑1 (MQ=255)
                                                                                                                                               ccaccaCCGCAAGCCGGGGGTGAACTTCTGGCGCATCGCGCCGGAGCGCGCGGCCAACGAGGCGCGCCCCTCGAAGTTCGGCCCGATGACCGAGACCTGGGGCATGTGGGAGAACTTCCTGATCCTGCGCGACCACCTCTGCCGCGGCCa   <  2:10570776/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: