Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 4,707,658 G→C A165A (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,6580GC99.9% 2612.0 / NA 686A165A (GCG→GCCPA14_RS21625zinc‑binding alcohol dehydrogenase family protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base C (305/380);  total (305/381)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 686 total aligned reads displayed.

CCCAGGCCGCCAGTTGCGGGGTGCCCTACAGCACCGCCTGGGACGCCCTGCAACGCAGCCAGGTCAAGGCCGGCACCCGGCTCCTGGTAATCGGCGCCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCGGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCGGAGGTGATCTTCGATACG  >  NC_008463/4707510‑4707805
                                                                                                                                                    |                                                                                                                                                   
cccAGGCCGCCAGTTGCGGGGTGCCCTACAGCACCGCCTGGGACGCCCTGCAACGCAGCCAGGTCAAGGCCGGCACCCGGCTCCTGGTAATCGGCGCCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGa                                                                                                                                                   >  1:4247272/1‑151 (MQ=255)
 ccAGGCCGCCAGTTGCGGGGTGCCCTACAGCACCGCCTGGGACGCCCTGCAACGCAGCCAGGTCAAGGCCGGCACCCGGCTCCTGGTAATCGGCGCCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGaa                                                                                                                                                  >  2:766316/1‑151 (MQ=255)
  cAGGCCGCCAGTTGCGGGGTGCCCTACAGCACCGCCTGGGACGCCCTGCAACGCAGCCAGGTCAAGGCCGGCACCCGGCTCCTGGTAATCGGCGCCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAg                                                                                                                                                 <  1:8888417/151‑1 (MQ=255)
  cAGGCCGCCAGTTGCGGGGTGCCCTACAGCACCGCCTGGGACGCCCTGCAACGCAGCCAGGTCAAGGCCGGCACCCGGCTCCTGGTAATCGGCGCCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAg                                                                                                                                                 <  2:1766238/151‑1 (MQ=255)
   aGGCCGCCAGTTGCGGGGTGCCCTACAGCACCGCCTGGGACGCCCTGCAACGCAGCCAGGTCAAGGCCGGCACCCGGCTCCTGGTAATCGGCGCCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGc                                                                                                                                                >  1:3507569/1‑151 (MQ=255)
   aGGCCGCCACTTGCGGGGTGCCCTACAGCACCGCCAGGGACGCCCTGCAACGCAGCCAGGACAAGGCCGGCACCCGGCTCCTGGTAATCGGCGCCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGc                                                                                                                                                <  2:5957232/151‑1 (MQ=255)
       cgccAGTTGCGGGGTGCCCTACAGCACCGCCTGGGACGCCCTGCAACGCAGCCAGGTCAAGGCCGGCACCCGGCTCCTGGTAATCGGCGCCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATc                                                                                                                                            >  2:7730421/1‑151 (MQ=255)
       cgccAGTTGCGGGGTGCCCTACAGCACCGCCTGGGACGCCCTGCAACGCAGCCAGGTCAAGGCCGGCACCCGGCTCCTGGTAATCGGCGCCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATc                                                                                                                                            <  2:36328/151‑1 (MQ=255)
         ccAGTTGCGGGGTGCCCTACAGCACCGCCTGGGACGCCCTGCAACGCAGCCAGGTCAAGGCCGGCACCCGGCTCCTGGTAATCGGCGCCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCgg                                                                                                                                          >  2:8436024/1‑151 (MQ=255)
         ccAGTTGCGGGGTGCCCTACAGCACCGCCTGGGACGCCCTGCAAAGCAGCCAGGTCAAGGCCGGCACCCGGCTCCTGGTAATCGGCGCCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCgg                                                                                                                                          >  1:8605016/1‑151 (MQ=255)
          cAGTTGCGGGGTGCCCTACAGCACCGCCTGGGACGCCCTGCAACGCAGCCAGGTCAAGGCCGGCACCCGGCTCCTGGTAATCGGCGCCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGgc                                                                                                                                         >  1:6710639/1‑151 (MQ=255)
          cAGTTGCGGGGTGCCCTACAGCACCGCCTGGGACGCCCTGCAACGCAGCCAGGTCAAGGCCGGCACCCGGCTCCTGGTAATCGGCGCCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGgc                                                                                                                                         <  1:2870829/151‑1 (MQ=255)
            gTTGCGGGGTGCCCTACAGCACCGCCTGGGACGCCCTGCAACGCAGCCAGGTCAAGGCCGGCACCCGGCTCCTGGTAATCGGCGCCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGgcgc                                                                                                                                       <  2:11154053/151‑1 (MQ=255)
               gCGGGGTGCCCTACAGCACCGCCTGGGACGCCCTGCAACGCAGCCAGGTCAAGGCCGGCACCCGGCTCCTGGTAATCGGCGCCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGt                                                                                                                                    <  1:3192814/151‑1 (MQ=255)
                cGGGGTGCCCTACAGCACCGCCTGGGACGCCCTGCAACGCAGCCAGGTCAAGGCCGGCACCCGGCTCCTGGTAATCGGCGCCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTg                                                                                                                                   <  2:5015812/151‑1 (MQ=255)
                 ggggTGCCCTACAGCACCGCCTGGGACGCCCTGCAACGCAGCCAGGTCAAGGCCGGCACCCGGCTCCTGGTAATCGGCGCCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTgc                                                                                                                                  >  1:3975629/1‑151 (MQ=255)
                   ggTGCCCTACAGCACCGCCTGGGACGCCCTGCAACGCAGCCAGGTCAAGGCCGGCACCCGGCTCCTGGTAATCGGCGCCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTgcgg                                                                                                                                <  1:8515417/151‑1 (MQ=255)
                   ggTGCCCTACAGCACCGCCTGGGACGCCCTGCAACGCAGCCAGGTCAAGGCCGGCACCCGGCTCCTGGTAATCGGCGCCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTgcgg                                                                                                                                >  2:3232405/1‑151 (MQ=255)
                   ggTGCCCTACAGCACCGCCTGGGACGCCCTGCAACGCAGCCAGGTCAAGGCCGGCACCCGGCTCCTGGTAATCGGCGCCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTgcgg                                                                                                                                <  2:110085/151‑1 (MQ=255)
                     tGCCCTACAGCACCGCCTGGGACGCCCTGCAACGCAGCCAGGTCAAGGCCGGCACCCGGCTCCTGGTAATCGGCGCCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTgcggcg                                                                                                                              >  2:8350802/1‑151 (MQ=255)
                      gCCCTACAGCACCGCCTGGGACGCCCTGCAACGCAGCCAGGTCAAGGCCGGCACCCGGCTCCTGGTAATCGGCGCCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGgcgc                                                                                                                             <  2:3430480/151‑1 (MQ=255)
                         cTACAGCACCGCCTGGGACGCCCTGCAACGCAGCCAGGTCAAGGCCGGCACCCGGCTCCTGGTAATCGGCGCCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCgag                                                                                                                          <  1:10073591/151‑1 (MQ=255)
                         cTACAGCACCGCCTGGGACGCCCTGCAACGCAGCCAGGTCAAGGCCGGCACCCGGCTCCTGGTAATCGGCGCCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCgag                                                                                                                          >  1:1819905/1‑151 (MQ=255)
                            cAGCACCGCCTGGGACGCCCTGCAACGCAGCCAGGTCAAGGCCGGCACCCGGCTCCTGGTAATCGGCGCCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCgaggag                                                                                                                       <  1:809390/151‑1 (MQ=255)
                            cAGCACCGCCTGGGACGCCCTGCAACGCAGCCAGGTCAAGGCCGGCACCCGGCTCCTGGTAATCGGCGCCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCgaggag                                                                                                                       >  2:614814/1‑151 (MQ=255)
                             aGCACCGCCTGGGACGCCCTGCAACGCAGCCAGGTCAAGGCCGGCACCCGGCTCCTGGTAATCGGCGCCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGc                                                                                                                      >  1:8406886/1‑151 (MQ=255)
                                aCCGCCTGGGACGCCCTGCAACGCAGCCAGGTCAAGGCCGGCACCCGGCTCCTGGTAATCGGCGCCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAAc                                                                                                                   <  2:6353532/151‑1 (MQ=255)
                                  cGCCTGGGACGCCCTGCAACGCAGCCAGGTCAAGGCCGGCACCCGGCTCCTGGTAATCGGCGCCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGc                                                                                                                 <  2:3113127/151‑1 (MQ=255)
                                   gCCTGGGACGCCCTGCAACGCAGCCAGGTCAAGGCCGGCACCCGGCTCCTGGTAATCGGCGCCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGcc                                                                                                                <  2:559691/151‑1 (MQ=255)
                                   gCCTGGGACGCCCTGCAACGCAGCCAGGTCAAGGCCGGCACCCGGCTCCTGGTAATCGGCGCCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGcc                                                                                                                >  1:200436/1‑151 (MQ=255)
                                       gggACGCCCTGCAACGCAGCCAGGTCAAGGCCGGCACCCGGCTCCTGGTAATCGGCGCCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGgc                                                                                                            >  1:6040923/1‑151 (MQ=255)
                                          aCGCCCTGCAACGCAGCCAGGTCAAGGCCGGCACCCGGCTCCTGGTAATCGGCGCCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCt                                                                                                         <  2:9111357/151‑1 (MQ=255)
                                           cGCCCTGCAACGCAGCCAGGTCAAGGCCGGCACCCGGCTCCTGGTAATCGGCGCCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTg                                                                                                        <  2:10740874/151‑1 (MQ=255)
                                            gCCCTGCAACGCAGCCAGGTCAAGGCCGGCACCCGGCTCCTGGTAATCGGCGCCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTgg                                                                                                       <  1:10308877/151‑1 (MQ=255)
                                            gCCCTGCAACGCAGCCAGGTCAAGGCCGGCACCCGGCTCCTGGTAATCGGCGCCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCGGGCCAAGGCCCGGGGGGCCGACGCGATCGGCGCGGTGCGGCGCGAGGCGCAACGCCAGGCGCTgg                                                                                                       >  2:660055/1‑151 (MQ=255)
                                              ccTGCAACGCAGCCAGGTCAAGGCCGGCACCCGGCTCCTGGTAATCGGCGCCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGaa                                                                                                     >  1:384410/1‑151 (MQ=255)
                                               cTGCAACGCAGCCAGGTCAAGGCCGGCACCCGGCTCCTGGTAATCGGCGCCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAg                                                                                                    >  1:6352161/1‑151 (MQ=255)
                                                 gCAACGCAGCCAGGTCAAGGCCGGCACCCGGCTCCTGGTAATCGGCGCCACCGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAgcg                                                                                                  <  2:1578490/151‑1 (MQ=255)
                                                   aaCGCAGCCAGGTCAAGGCCGGCACCCGGCTCCTGGTAATCGGCGCCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCa                                                                                                <  2:8908405/151‑1 (MQ=255)
                                                      gCAGCCAGGTCAAGGCCGGCACCCGGCTCCTGGTAATCGGCGCCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAggg                                                                                             <  1:4931713/151‑1 (MQ=255)
                                                       cAGCCAGGTCAAGGCCGGCACCCGGCTCCTGGTAATCGGCGCCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGc                                                                                            >  2:3326894/1‑151 (MQ=255)
                                                       cAGCCAGGTCAAGGCCGGCACCCGGCTCCTGGTAATCGGCGCCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGc                                                                                            <  1:215074/151‑1 (MQ=255)
                                                         gCCAGGTCAAGGCCGGCACCCGGCTCCTGGTAATCGGCGCCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCAt                                                                                          <  2:9272218/151‑1 (MQ=255)
                                                           cAGGTCAAGGCCGGCACCCGGCTCCTGGTAATCGGCGCCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATcg                                                                                        >  1:9657451/1‑151 (MQ=255)
                                                           cAGGTCAAGGCCGGCACCCGGCTCCTGGTAATCGGCGCCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATcg                                                                                        <  2:9585223/151‑1 (MQ=255)
                                                              gTCAAGGCCGGCACCCGGCTCCTGGTAATCGGCGCCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATcgccg                                                                                     >  2:10972727/1‑151 (MQ=255)
                                                                cAAGGCCGGCACCCGGCTCCTGGTAATCGGCGCCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCcgcg                                                                                   <  2:10048918/151‑1 (MQ=255)
                                                                  cggCCGGCACCCGGCTCCTGGTAATCGGCGCCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCgct                                                                                 <  2:1138980/150‑1 (MQ=255)
                                                                  aGGCCGGCACCCGGCTCCTGGTAATCGGCGCCACCGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCgct                                                                                 <  1:7296429/151‑1 (MQ=255)
                                                                  aGGCCGGCACCCGGCTCCTGGTAATCGGCGCCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCgct                                                                                 >  1:5349944/1‑151 (MQ=255)
                                                                        gCACCCGGCTCCTGGTAATCGGCGCCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCgctgctgct                                                                           >  1:9251986/1‑151 (MQ=255)
                                                                        gCACCCGGCTCCTGGTAATCGGCGCCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCgctgctgct                                                                           <  1:1001669/151‑1 (MQ=255)
                                                                        gCACCCGGCTCCTGGTAATCGGCGCCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCgctgctgct                                                                           <  2:5501045/151‑1 (MQ=255)
                                                                        gCACCCGGCTCCTGGTAATCGGCGCCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCgctgctgct                                                                           >  1:8711959/1‑151 (MQ=255)
                                                                         cACCCGGCTCCTGGTAATCGGCGCCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGctgctgctg                                                                          >  2:4721221/1‑151 (MQ=255)
                                                                             cGGCTCCTGGTAATCGGCGCCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCg                                                                      >  2:6450840/1‑151 (MQ=255)
                                                                             cGGCTCCTGGTAAGCGGCGCCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCg                                                                      <  1:3038117/151‑1 (MQ=255)
                                                                              ggCTCCTGGTAATCGGCGCCACCGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGa                                                                     <  1:4649357/151‑1 (MQ=255)
                                                                               gCTCCTGGTAATCGGCGCCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAg                                                                    >  1:10838090/1‑151 (MQ=255)
                                                                               gCTCCTGGTAATCGGCGCCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAg                                                                    <  2:8719099/151‑1 (MQ=255)
                                                                                 tCCTGGTAATCGGCGCCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAgcg                                                                  >  2:7221582/1‑150 (MQ=255)
                                                                                  ccTGGTAATCGGCGCCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAgccg                                                                 <  2:136630/151‑1 (MQ=255)
                                                                                   cTGGTAATCGGCGCCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCACGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAgccgg                                                                <  2:8515441/151‑1 (MQ=255)
                                                                                   cTGGTAATCGGCGCCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAgccgg                                                                >  1:11174738/1‑151 (MQ=255)
                                                                                    tGGTAATCGGCGCCAACGGCGCGGTCGGCTAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAgccggc                                                               >  1:1300717/1‑151 (MQ=255)
                                                                                    tGGTAATCGGCGCCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAgccggc                                                               <  2:4346694/151‑1 (MQ=255)
                                                                                      gTAATCGGCGCCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAgccggccg                                                             >  1:6907106/1‑151 (MQ=255)
                                                                                       tAATCGGCGCCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGa                                                            >  1:5819899/1‑151 (MQ=255)
                                                                                            ggCGCCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTgg                                                       <  1:567412/151‑1 (MQ=255)
                                                                                            ggCGCCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCAGCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTgg                                                       >  2:2421263/1‑151 (MQ=255)
                                                                                              cgcCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTggcg                                                     <  2:1668372/151‑1 (MQ=255)
                                                                                              cgcCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTggcg                                                     <  1:10615464/151‑1 (MQ=255)
                                                                                               gcCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTggcgg                                                    >  2:4033595/1‑151 (MQ=255)
                                                                                               gcCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTggcgg                                                    >  1:2238116/1‑151 (MQ=255)
                                                                                                 cAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGACCTGGCGGca                                                  >  1:1432439/1‑151 (MQ=255)
                                                                                                 cAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGca                                                  >  2:9466505/1‑151 (MQ=255)
                                                                                                  aaCGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGcac                                                 <  2:6909854/151‑1 (MQ=255)
                                                                                                      gcgcGGTCGGGAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCGGAAGCGACCGGCGCGGTGCGGCGCGAGGAGCTACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAggc                                             <  1:6465634/151‑2 (MQ=255)
                                                                                                      gcgcGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGt                                             <  1:1613867/151‑1 (MQ=255)
                                                                                                        gcgGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTgg                                           >  2:2783420/1‑151 (MQ=255)
                                                                                                         cgGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGa                                          >  2:7971314/1‑151 (MQ=255)
                                                                                                          ggTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGaa                                         >  1:2381530/1‑151 (MQ=255)
                                                                                                          ggTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGaa                                         >  1:7849722/1‑151 (MQ=255)
                                                                                                              ggCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGa                                     <  1:5074949/151‑1 (MQ=255)
                                                                                                                 aaGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCt                                  >  2:5391386/1‑151 (MQ=255)
                                                                                                                 aaGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCt                                  >  2:4417133/1‑151 (MQ=255)
                                                                                                                 aaGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCt                                  <  2:1913333/151‑1 (MQ=255)
                                                                                                                  aGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCtt                                 <  2:4179610/151‑1 (MQ=255)
                                                                                                                         gctggcgctggcCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCcggc                          >  1:6985396/1‑151 (MQ=255)
                                                                                                                              cgctggcCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCggcggggc                     >  2:10535643/1‑151 (MQ=255)
                                                                                                                                 tggcCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGa                  <  1:7506281/151‑1 (MQ=255)
                                                                                                                                 tggcCAAGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGa                  <  2:2617349/151‑1 (MQ=255)
                                                                                                                                      aaGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGa             >  2:5295537/1‑151 (MQ=255)
                                                                                                                                      aaGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGa             <  1:7251016/151‑1 (MQ=255)
                                                                                                                                      aaGGCCCGGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGACGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACAGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGa             >  1:2770473/1‑151 (MQ=255)
                                                                                                                                       aGGCCCTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGAt            <  1:9055022/151‑1 (MQ=255)
                                                                                                                                           ccTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTc        >  2:10333949/1‑151 (MQ=255)
                                                                                                                                           ccTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTc        <  2:2029339/151‑1 (MQ=255)
                                                                                                                                           ccTGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTc        <  2:1474400/151‑1 (MQ=255)
                                                                                                                                             tGGGGGCCGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCCGAGGTGATCTTCGa      >  1:2537788/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CCCAGGCCGCCAGTTGCGGGGTGCCCTACAGCACCGCCTGGGACGCCCTGCAACGCAGCCAGGTCAAGGCCGGCACCCGGCTCCTGGTAATCGGCGCCAACGGCGCGGTCGGCAAGGCCGCGCTGGCGCTGGCCAAGGCCCTGGGGGCGGAAGCGATCGGCGCGGTGCGGCGCGAGGAGCAACGCCAGGCGCTGGAAGCGCAGGGCATCGCCGCGCTGCTGCTGGGCGAGCCGGCCGAACTGGCGGCACAGGTGGAAGGGATCTTCCCCGGCGGGGCGGAGGTGATCTTCGATACG  >  NC_008463/4707510‑4707805

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: