Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 4,707,787 G→C A208A (GCG→GCC)  PA14_RS21625 → zinc‑binding alcohol dehydrogenase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084634,707,7870GC99.7% 2465.8 / NA 622A208A (GCG→GCC) PA14_RS21625zinc‑binding alcohol dehydrogenase family protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base C (284/336);  minor base A (2/0);  total (286/336)
Fisher's exact test for biased strand distribution p-value = 2.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 648 total aligned reads displayed.

TGGCCAAGGCCCTGGGGGCGGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCGGAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAGGCGGCGGTCCCGGCGCTGGCCGCCTTCGGTCGCCTGGCGATCATCGCCGCACCAGCCGACGGCCATGTGCGCCTGCCGGCGTTGAACCTGTACCGTCGCGGTGGCTCGGTG  >  NC_008463/4707639‑4707934
                                                                                                                                                    |                                                                                                                                                   
tggcCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGa                                                                                                                                                   <  1:9792284/151‑1 (MQ=255)
 ggcCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAg                                                                                                                                                  <  2:1310567/151‑1 (MQ=255)
 ggcCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAg                                                                                                                                                  <  2:1874733/151‑1 (MQ=255)
  gcCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAgg                                                                                                                                                 <  1:8103354/151‑1 (MQ=255)
   ccAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCACCGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGt                                                                                                                                                >  2:710978/1‑151 (MQ=255)
   ccAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCNAGGt                                                                                                                                                <  2:10164982/151‑1 (MQ=255)
      aGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGAt                                                                                                                                             <  2:6152600/151‑1 (MQ=255)
       ggCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATc                                                                                                                                            <  2:6456638/151‑1 (MQ=255)
        gCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCt                                                                                                                                           >  1:3726094/1‑151 (MQ=255)
          ccTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTc                                                                                                                                         >  1:1141134/1‑151 (MQ=255)
            tGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGa                                                                                                                                       <  2:5324854/151‑1 (MQ=255)
             gggggCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGAt                                                                                                                                      <  1:6894745/151‑1 (MQ=255)
             gggggCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGAt                                                                                                                                      <  1:9962417/151‑1 (MQ=255)
              ggggCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGATa                                                                                                                                     >  1:7431749/1‑151 (MQ=255)
              ggggCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGATa                                                                                                                                     <  2:6286612/151‑1 (MQ=255)
                ggCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACg                                                                                                                                   <  1:8577636/151‑1 (MQ=255)
                    gAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACGACCg                                                                                                                               <  1:7698561/151‑1 (MQ=255)
                     aaGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACGACCgg                                                                                                                              >  1:472226/1‑151 (MQ=255)
                        cGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACGACCGGttt                                                                                                                           >  2:1150460/1‑151 (MQ=255)
                        cGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACGACCGGttt                                                                                                                           >  1:7956614/1‑151 (MQ=255)
                          aTCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACGACCGGTTTCt                                                                                                                         <  2:6681276/151‑1 (MQ=255)
                          aTCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACGACCGGTTTCt                                                                                                                         >  2:5119967/1‑151 (MQ=255)
                                gcgGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACGACCGGTTTCTGGCTCg                                                                                                                   <  1:323741/151‑1 (MQ=255)
                                 cgGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACGACCGGTTTCTGGCTCGa                                                                                                                  <  1:9657274/151‑1 (MQ=255)
                                   gTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAgg                                                                                                                <  1:822044/151‑1 (MQ=255)
                                     gcggcgCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAggcg                                                                                                              <  1:9194233/151‑1 (MQ=255)
                                       ggcgCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAggcggc                                                                                                            <  1:9319586/151‑1 (MQ=255)
                                         cgcgAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAGGcggcgg                                                                                                          >  2:3691764/1‑151 (MQ=255)
                                          gcgAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAGGCGGCGGt                                                                                                         >  2:6835249/1‑151 (MQ=255)
                                            gaggagCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAGGCGGCGGTcc                                                                                                       <  1:2960316/151‑1 (MQ=255)
                                             aggagCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAGGCGGCGGTccc                                                                                                      <  2:2448275/151‑1 (MQ=255)
                                              ggagCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAGGCGGCGGTCCCg                                                                                                     <  2:9745877/151‑1 (MQ=255)
                                               gagCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAGGCGGCGGTCCCgg                                                                                                    >  2:7639076/1‑151 (MQ=255)
                                                 gCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAGGCGGCGGTCCCGgcg                                                                                                  <  1:1513257/151‑1 (MQ=255)
                                                  cAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAGGCGGCGGTCCCGgcgc                                                                                                 >  1:10009304/1‑151 (MQ=255)
                                                  cAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAGGCGGCGGTCCCGgcgc                                                                                                 >  1:5644330/1‑151 (MQ=255)
                                                   aaCGCCAGGCGCTGGTAGCGCAGGGCATCGCCGCGCGGCTGCTGGGCGAGCCGGCGGGACTGGCGGCACAGGTGGAGGGGATCTTCCCCGGCGGGGCCGAGGTGTTCTTCGTTGCGACCGGGTTCTGGCTCGAGGCGGCGGTCCCGgggcg                                                                                                >  2:8152791/1‑150 (MQ=255)
                                                   aaCGCCAGGCGCTGGAAGCGCAGGGCCTCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAGGGGCTCTTCCCCGGCGGGGCAGAGGTGATCTTCGAGACGACCGGTTTCTGGCTCGTGGCGGCGGTCCCGgcgcg                                                                                                >  2:10276721/1‑150 (MQ=255)
                                                   aaCGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAGGCGGCGGTCCCGGCGCt                                                                                                >  1:10441730/1‑151 (MQ=255)
                                                      gCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAGGCGGCGGTCCCGGCGCTGgc                                                                                             >  1:5826430/1‑151 (MQ=255)
                                                       ccAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAGGCGGCGGTCCCGGCGCTGgcc                                                                                            >  2:6952402/1‑151 (MQ=255)
                                                       ccAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAGGCGGCGGTCCCGGCGCTGgcc                                                                                            <  1:8801868/151‑1 (MQ=255)
                                                       ccAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAGGCGGCGGTCCCGGCGCTGgcc                                                                                            >  1:8933396/1‑151 (MQ=255)
                                                          ggCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAGGCGGCGGTCCCGGCGCTGgccgcc                                                                                         >  2:10451366/1‑151 (MQ=255)
                                                           gcgcTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAGGCGGCGGTCCCGGCGCTGGCCGCCt                                                                                        >  1:7228755/1‑151 (MQ=255)
                                                            cgcTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAGGCGGCGGTCCCGGCGCTGGCCGCCtt                                                                                       >  2:7819388/1‑151 (MQ=255)
                                                               tGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAGGCGGCGGTCCCGGCGCTGGCCGCCTGCgg                                                                                    >  2:2864684/1‑151 (MQ=255)
                                                                     cgcAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAGGCGGCGGTCCCGGCGCTGGCCGCCTTCGGTCGCCt                                                                              <  1:5180098/151‑1 (MQ=255)
                                                                       cAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAGGCGGCGGTCCCGGCGCTGGCCGCCTTCGGTCGCCTgg                                                                            <  2:4318152/151‑1 (MQ=255)
                                                                         gggCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAGGCGGCGGTCCCGGCGCTGGCCGCCTTCGGTCGCCTGGCg                                                                          <  2:7456749/151‑1 (MQ=255)
                                                                            cATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAGGCGGCGGTCCCGGCGCTGGCCGCCTTCGGTCGCCTGGCGatc                                                                       >  1:4474510/1‑151 (MQ=255)
                                                                             aTCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAGGCGGCGGTCCCGGCGCTGGCCGCCTTCGGTCGCCTGGCGatca                                                                      <  1:5484035/151‑1 (MQ=255)
                                                                               cgccgcGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAGGCGGCGGTCCCGGCGCTGGCCGCCTTCGGTCGCCTGGCGatcatc                                                                    <  1:836663/151‑1 (MQ=255)
                                                                               cgccgcGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAGGCGGCGGTCCCGGCGCTGGCCGCCTTCGGTCGCCTGGCGatcatc                                                                    <  1:6729278/151‑1 (MQ=255)
                                                                               cgccgcGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAGGCGGCGGTCCCGGCGCTGGCCGCCTTCGGTCGCCTGGCGatcatc                                                                    <  1:4982506/151‑1 (MQ=255)
                                                                                   gcgcTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAGGCGGCGGTCCCGGCGCTGGCCGCCTTCGGTCGCCTGGCGATCATcgccg                                                                >  2:3251256/1‑151 (MQ=255)
                                                                                    cgcTGCTGCTGTGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCAGAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAGGCGGCGGTCCCGGCGCTGGCCGCCTTCGGTCGCCTGTCGATCATcgccgc                                                               >  2:471817/1‑151 (MQ=255)
                                                                                     gctgctgctGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAGGCGGCGGTCCCGGCGCTGGCCGCCTTCGGTCGCCTGGCGATCATCGCCGCa                                                              >  2:221782/1‑151 (MQ=255)
                                                                                      ctgctgctgGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAGGCGGCGGTCCCGGCGCTGGCCGCCTTCGGTCGCCTGGCGATCATCGCCGCAc                                                             <  1:201495/151‑1 (MQ=255)
                                                                                         ctgctgGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAGGCGGCGGTCCCGGCGCTGGCCGCCTTCGGTCGCCTGGCGATCATCGCCGCACCAg                                                          <  2:357904/151‑1 (MQ=255)
                                                                                         ctgctgGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAGGCGGCGGTCCCGGCGCTGGCCGCCTTCGGTCGCCTGGCGATCATCGCCGCACCAg                                                          <  2:11065688/151‑1 (MQ=255)
                                                                                         ctgctgGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAGGCGGCGGTCCCGGCGCTGGCCGCCTTCGGTCGCCTGGCGATCATCGCCGCACCAg                                                          <  2:2630238/151‑1 (MQ=255)
                                                                                              gggCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAGGCGGCGGTCCCGGCGCTGGCCGCCTTCGGTCGCCTGGCGATCATCGCCGCACCAGCCGAc                                                     <  2:5918374/151‑1 (MQ=255)
                                                                                              gggCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACGACCGGTGTCTGGCTCGAGGCGGCGGTCCCGGCGCTGGCCGCCTTCGGTCGCCTGGCGATCATCGCCGCACCAGCCGAc                                                     >  1:1232511/1‑151 (MQ=255)
                                                                                                        gccgAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAGGCGGCGGTCCCGGCGCTGGCCGCCTTCGGTCGCCTGGCGATCATCGCCGCACCAGCCGACGGCCATGTgc                                           <  1:3091260/151‑1 (MQ=255)
                                                                                                        gccgAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAGGCGGCGGTCCCGGCGCTGGCCGCCTTCGGTCGCCTGGCGATCATCGCCGCACCAGCCGACGGCCATGTgc                                           >  2:4583262/1‑151 (MQ=255)
                                                                                                         ccgAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAGGCGGCGGTCCCGGCGCTGGCCGCCTTCGGTCGCCTGGCGATCATCGCCGCACCAGCCGACGGCCATGTgcg                                          >  1:615599/1‑151 (MQ=255)
                                                                                                          cgAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAGGCGGCGGTCCCGGCGCTGGCCGCCTTCGGTCGCCTGGCGATCATCGCCGCACCAGCCGACGGCCATGTgcgc                                         >  2:9666129/1‑151 (MQ=255)
                                                                                                            aaCTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAGGCGGCGGTCCCGGCGCTGGCCGCCTTCGGTCGCCTGGCGATCATCGCCGCACCAGCCGACGGCCATGTGCGCCt                                       <  2:8456398/151‑1 (MQ=255)
                                                                                                            aaCTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAGGCGGCGGTCCCGGCGCTGGCCGCCTTCGGTCGCCTGGCGATCATCGCCGCACCAGCCGACGGCCATGTGCGCCt                                       >  1:2749198/1‑151 (MQ=255)
                                                                                                            aaCTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAGGCGGCGGTCCCGGCGCTGGCCGCCTTCGGTCGCCTGGCGATCATCGCCGCACCAGCCGACGGCCATGTGCGCCt                                       <  1:4272484/151‑1 (MQ=255)
                                                                                                               tGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAGGCGGCGGTCCCGGCGCTGGCCGCCTTCGGTCGCCTGGCGATCATCGCCGCACCAGCCGACGGCCATGTGCGCCTGcc                                    >  1:6338523/1‑151 (MQ=255)
                                                                                                                 gcggcACAGGTGGAAGGGATCTGCCCCGGCGGGGCCGAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAGGCGGCGGTCCCGGCGCTGGCCGCCTTCGGTCGCCTGGCGATCATCGCCGCACCAGCCGACGGCCATGTGCGCCTGCCgg                                  <  2:6766573/151‑1 (MQ=255)
                                                                                                                  cggcACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAGGCGGCGGTCCCGGCGCTGGCCGCCTTCGGTCGCCTGGCGATCATCGCCGCACCAGCCGACGGCCATGTGCGCCTGCCGGc                                 <  2:855032/151‑1 (MQ=255)
                                                                                                                     cacaGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAGGCGGCGGTCCCGGCGCTGGCCGCCTTCGGTCGCCTGGCGATCATCGCCGCACCAGCCGACGGCCATGTGCGCCTGCCGGCGtt                              >  1:2218752/1‑151 (MQ=255)
                                                                                                                      acaGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAGGCGGCGGTCCCGGCGCTGGCCGCCTTCGGTCGCCTGGCGATCATCGCCGCACCAGCCGACGGCCATGTGCGCCTGCCGGCGTTg                             <  1:10268526/151‑1 (MQ=255)
                                                                                                                      acaGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAGGCGGCGGTCCCGGCGCTGGCCGCCTTCGGTCGCCTGGCGATCATCGCCGCACCAGCCGACGGCCATGTGCGCCTGCCGGCGTTg                             >  2:7552146/1‑151 (MQ=255)
                                                                                                                       caGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAGGCGGCGGTCCCGGCGCTGGCCGCCTTCGGTCGCCTGGCGATCATCGCCGCACCAGCCGACGGCCATGTGCGCCTGCCGGCGTTGa                            <  2:8798887/151‑1 (MQ=255)
                                                                                                                        aGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAGGCGGCGGTCCCGGCGCTGGCCGCCTTCGGTCGCCTGGCGATCATCGCCGCACCAGCCGACGGCCATGTGCGCCTGCCGGCGTTGaa                           >  2:9226442/1‑151 (MQ=255)
                                                                                                                        aGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAGGCGGCGGTCCCGGCGCTGGCCGCCTTCGGTCGCCTGGCGATCATCGCCGCACCAGCCGAAGGCCATGTGCGCCTGCCGGCGTTGaa                           >  1:1006145/1‑151 (MQ=255)
                                                                                                                           tGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAGGCGGCGGTCCCGGCGCTGGCCGCCTTCGGTCGCCTGGCGATCATCGCCGCACCAGCCGACGGCCATGTGCGCCTGCCGGCGTTGAACCt                        >  2:1698804/1‑151 (MQ=255)
                                                                                                                            ggAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAGGCGGCGGTCCCGGCGCTGGCCGCCTTCGGTCGCCTGGCGATCATCGCCGCACCAGCCGACGGCCATGTGCGCCTGCCGGCGTTGAACCTg                       >  2:4885753/1‑151 (MQ=255)
                                                                                                                            ggAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAGGCGGCGGTCCCGGCGCTGGCCGCCTTCGGTCGCCTGGCGATCATCGCCGCACCAGCCGACGGCCATGTGCGCCTGCCGGCGTTGAACCTg                       <  1:2589128/151‑1 (MQ=255)
                                                                                                                               aGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAGGCGGCGGTCCCGGCGCTGGCCGCCTTCGGTCGCCTGGCGATCATCGCCGCACCAGCCGACGGCCATGTGCGCCTGCCGGCGTTGAACCTGTAc                    <  2:5751115/151‑1 (MQ=255)
                                                                                                                                 ggATCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACGCCCGGTTTCTGGCTCGAGGCGGCGGTCCCGGCGCTGGCCGCCTTCGGTCGCCTGGCGATCATCGCCGCACCAGCCGACGGCCATGTGCGCCTGCCGGCGTTGAACCTGTACCg                  <  2:10727163/151‑1 (MQ=255)
                                                                                                                                   aTCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAGGCGGCGGTCCCGGCGCTGGCCGCCTTCGGTCGCCTGGCGATCATCGCCGCACCAGCCGACGGCCATGTGCGCCTGCCGGCGTTGAACCTGTACCGTc                >  2:3966073/1‑151 (MQ=255)
                                                                                                                                   aTCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAGGCGGCGGTCCCGGCGCTGGCCGCCTTCGGTCGCCTGGCGATCATCGCCGCACCAGCCGACGGCCATGTGCGCCTGCCGGCGTTGAACCTGTACCGTc                <  2:3377486/151‑1 (MQ=255)
                                                                                                                                   aTCTTCCCCGGCGGGGCCGAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAGGCGGCGGTCCCGGCGCTGGCCGCCTTCGGTCGCCTGGCGATCATCGCCGCACCAGCCGACGGCCATGTGCGCCTGCCGGCGTTGAACCTGTACCGTc                <  1:10773076/151‑1 (MQ=255)
                                                                                                                                     cTTCCCCGGCGGGGCCGAGGTGATCTTCGATCCGACCGGTTTCTGGCTCGAGGCGGCGGTCCCGGCGCTGGCCGCCTTCGGTCGCCTGGCGAGCATCGCCGCACCAGCCGACGGCCATGTGCGCCTGCCGGCGTTGAACCTGTACCGTcgc              >  2:9557253/1‑151 (MQ=255)
                                                                                                                                         cccGGCGGGGCCGAGGTGATCTTCGATTCGTCCGGGTTCTGGCTCGAGGCGGCGGTCCCGGCGCGGTCCGCCTTCGGTCGACTGGCGATCATCGCCGCACCAGCCGACGTCCATGTGCGCCTGCCGGCGTTGAACCTGTACCGTCGCGGTg          <  2:7996092/151‑1 (MQ=255)
                                                                                                                                         cccGGCGGGGCCGAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAGGCGGCGGTCCCGGCGCTGGCCGCCTTCGGTCGCCTGGCGATCATCGCCGCACCAGCCGACGGCCATGTGCGCCTGCCGGCGTTGAACCTGTACCGTCGCGGTg          >  2:5603661/1‑151 (MQ=255)
                                                                                                                                          ccGGCGGGGCCGAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAGGCGGCGGTCCCGGCGCTGGCCGCCTTCGGTCGCCTGGCGATCATCGCCGCACCAGCCGACGGCCATGTGCGCCTGCCGGCGTTGAACCTGTACCGTCGCGGTgg         <  1:6991722/151‑1 (MQ=255)
                                                                                                                                          ccGGCGGGGCCGAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAGGCGGCGGTCCCGGCGCTGGCCGCCTTCGGTCGCCTGGCGATCATCGCCGCACCAGCCGACGGCCATGTGCGCCTGCCGGCGTTGAACCTGTACCGTCGCGGTgg         <  1:6497431/151‑1 (MQ=255)
                                                                                                                                          ccGGCGGGGCCGAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAGGCGGCGGTCCCGGCGCTGGCCGCCTTCGGTCGCCTGGCGATCATCGCCGCACCAGCCGACGGCCATGTGCGCCTGCCGGCGTTGAACCTGTACCGTCGCGGTgg         <  2:2096351/151‑1 (MQ=255)
                                                                                                                                           cggcggGGCCGAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAGGCGGCGGTCCCGGCGCTGGCCGCCTTCGGTCGCCTGGCGATCATCGCCGCACCAGCCGACGGCCATGTGCGCCTGCCGGCGTTGAACCTGTACCGTCGCGGTGGc        <  2:3119525/151‑1 (MQ=255)
                                                                                                                                             gcggggccgAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAGGCGGCGGTCCCGGCGCTGGCCGCCTTCGGTCGCCTGGCGATCATCGCCGCACCAGCCGACGGCCATGTGCGCCTGCCGGCGTTGAACCTGTACCGTCGCGGTGGCTc      <  1:413920/151‑1 (MQ=255)
                                                                                                                                              cggggccgAGGTGCTCTTCGATACGACCGGTTTCTGGCTCGAGGCGGCGGTCCCGGCGCTGGCCGCCTTCGGTCGCCTGGCGATCATCGCCGCACCAGCCGACGGCCATGTGCGCCTGCCGGCGTTGAACCTGTACCGTCGCGGTGGCTCg     <  1:10656217/151‑1 (MQ=255)
                                                                                                                                              cggggccgAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAGGCGGCGGTCCCGGCGCTGGCCGCCTTCGGTCGCCTGGCGATCATCGCCGCACCAGCCGACGGCCATGTGCGCCTGCCGGCGTTGAACCTGTACCGTCGCGGTGGCTCg     >  2:10821069/1‑151 (MQ=255)
                                                                                                                                              cggggccgAGGTGATCTTCGATACGACCGGTTGCTGGCTCGAGGCGGCGGTCCCGGCGCTGGCCGCCTTCGGTCGCCTGGCGATCATCGCCGCACCAGCCGACGGCCATGTGCGCCTGCCGGCGATGAACCTGTACCGTCGCGGTGGCTCg     <  2:7246806/151‑1 (MQ=255)
                                                                                                                                                 ggccgAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAGGCGGCGGTCCCGGCGCTGGCCGCCTTCGGTCGCCTGGCGATCATCGCCGCACCAGCCGACGGCCATGTGCGCCTGCCGGCGTTGAACCTGTACCGTCGCGGTGGCTCggtg  >  2:9052115/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TGGCCAAGGCCCTGGGGGCGGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCGGAGGTGATCTTCGATACGACCGGTTTCTGGCTCGAGGCGGCGGTCCCGGCGCTGGCCGCCTTCGGTCGCCTGGCGATCATCGCCGCACCAGCCGACGGCCATGTGCGCCTGCCGGCGTTGAACCTGTACCGTCGCGGTGGCTCGGTG  >  NC_008463/4707639‑4707934

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: ‑