Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,797,898 C→G 6.1% V47L (GTG→CTG)  vdh ← vanillin dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,797,8980CG 6.1% 311.7 / 9.1 98V47L (GTG→CTG) vdhvanillin dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (48/44);  new base G (4/2);  total (52/46)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.99e-01

AGTTGGCCATGGCCCCGGTTTCACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTGGCGGGGCGCGACTGCCCGCCAATCAGCAAAGGCACCTGCAACATGCTCATTC  >  NC_002947/3797751‑3798044
                                                                                                                                                     |                                                                                                                                                   
aGTTGGCCATGGCCCCGGTTTCACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACcg                                                                                                                                                   >  1:1323572/1‑150 (MQ=255)
 gTTGGCCATGGCCCCGGTTTCACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACcgc                                                                                                                                                  >  2:1548830/1‑150 (MQ=255)
 gTTGGCCATGGCCCCGGTTTCACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCAGCGCTGGAAGCAGACACcgc                                                                                                                                                  >  2:1118287/1‑150 (MQ=255)
     gCCATGGCCCCGGTTTCACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGcggca                                                                                                                                              >  1:1508150/1‑150 (MQ=255)
      ccATGGCCCCGGTTTCACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGcggcat                                                                                                                                             >  2:2171523/1‑150 (MQ=255)
      ccATGGCCCCGGTTTCACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCAGCCGcggcat                                                                                                                                             >  1:1910277/1‑150 (MQ=255)
       cATGGCCCCGGTTTCACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACcg                                                                                                                                                   >  1:2364547/1‑143 (MQ=255)
       cATGGCCCCGGTTTCACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACcg                                                                                                                                                   <  2:2364547/143‑1 (MQ=255)
       cATGGCCCCGGTTTCACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGcggcatc                                                                                                                                            <  2:610130/150‑1 (MQ=255)
       cATGGCCCCGGTTTCACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGcggcatc                                                                                                                                            >  1:1561677/1‑150 (MQ=255)
          ggCCCCGGTTTCACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGcggcatcggc                                                                                                                                         <  2:788953/150‑1 (MQ=255)
            ccccGGTTTCACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGcggcatcggcat                                                                                                                                       >  1:1038003/1‑150 (MQ=255)
             cccGGTTTCACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGcggcatc                                                                                                                                            >  1:1841354/1‑144 (MQ=255)
             cccGGTTTCACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGcggcatc                                                                                                                                            <  2:1841354/144‑1 (MQ=255)
             cccGGTTTCACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCggcatcggcatc                                                                                                                                      >  2:1280257/1‑150 (MQ=255)
             cccGGTTTCACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCggcatcggcatc                                                                                                                                      >  2:1280875/1‑150 (MQ=255)
                  tttCACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCa                                                                                                                                 <  2:2204220/150‑1 (MQ=255)
                     cACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAgtg                                                                                                                              <  2:757229/150‑1 (MQ=255)
                     cACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAgtg                                                                                                                              >  1:1672190/1‑150 (MQ=255)
                            ggcggcgATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCCAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCGACCGCGGCATCGGCATCGGCCAGTGTggcgg‑c                                                                                                                      <  2:847737/150‑1 (MQ=255)
                                 cgATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCNGCCAGCNAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCCGCCGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAg                                                                                                                 >  1:1588701/1‑150 (MQ=255)
                                     gAACTCGGCGGCCGTCGCTTGCAACGGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTg                                                                                                             <  2:1588701/150‑1 (MQ=255)
                                     gAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTg                                                                                                             <  1:689184/150‑1 (MQ=255)
                                      aaCTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGc                                                                                                            <  1:1440466/150‑1 (MQ=255)
                                         tCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGac                                                                                                         >  2:322746/1‑150 (MQ=255)
                                         tCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGac                                                                                                         >  2:1503232/1‑150 (MQ=255)
                                              ggcCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGGTGGCTGCTGCCAGCGcggc                                                                                                                                               >  1:784166/1‑108 (MQ=255)
                                              ggcCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGGTGGCTGCTGCCAGCGcggc                                                                                                                                               <  2:784166/108‑1 (MQ=255)
                                              ggcCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGAcaccac                                                                                                    >  2:609003/1‑150 (MQ=255)
                                              ggcCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGAcaccac                                                                                                    <  1:2171523/150‑1 (MQ=255)
                                                 ccTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGGACAAAAGCCGGGCT‑‑CAGCAGCCACCGCGGCCTCGGAATCGGCCTGCATGTCCG‑CGGCAGCCTGCGACACCACCTc                                                                                                 >  1:656329/1‑150 (MQ=255)
                                                 ccTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACTCCACCTc                                                                                                 <  1:1118287/150‑1 (MQ=255)
                                                 ccTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTc                                                                                                 >  2:673379/1‑150 (MQ=255)
                                                  cTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTTCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCg                                                                                                >  2:155174/1‑150 (MQ=255)
                                                   tCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCGGCCTGCGACACCACCTCGc                                                                                               >  2:2382611/1‑150 (MQ=255)
                                                     gCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCa                                                                                             <  2:1038003/150‑1 (MQ=255)
                                                        tGCAACAGATGAGCGCCGGCCAGCAAGCGGCTGCGCCGCTCGCCCGGGGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTc                                                                                          <  2:1508150/150‑1 (MQ=255)
                                                        tGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTc                                                                                          >  1:1055398/1‑150 (MQ=255)
                                                          cAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGtgtg                                                                                                                            <  1:2355664/115‑1 (MQ=255)
                                                          cAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGtgtg                                                                                                                            >  2:2355664/1‑115 (MQ=255)
                                                          cAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTggcgg‑cggc                                                                                                                   <  1:1777312/123‑1 (MQ=255)
                                                          cAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTggcgg‑cggc                                                                                                                   >  2:1777312/1‑123 (MQ=255)
                                                              aGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCggg                                                                                    <  2:1323572/150‑1 (MQ=255)
                                                              aGATCAGCGCATGCCAGCAACCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCggg                                                                                    <  2:552452/150‑1 (MQ=255)
                                                                 tCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTa                                                                                 >  1:235662/1‑150 (MQ=255)
                                                                   aGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACa                                                                               <  2:1561677/150‑1 (MQ=255)
                                                                   aGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACa                                                                               <  1:997014/150‑1 (MQ=255)
                                                                    gcgcCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAg                                                                              >  1:1910876/1‑150 (MQ=255)
                                                                     cgcCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAgc                                                                             <  1:1280875/150‑1 (MQ=255)
                                                                        cTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCt                                                                          >  1:1367742/1‑150 (MQ=255)
                                                                            cAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGaa                                                                      <  1:343575/150‑1 (MQ=255)
                                                                                 aGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTc                                                                 <  2:1910876/150‑1 (MQ=255)
                                                                                  gCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCg                                                                >  2:2333520/1‑150 (MQ=255)
                                                                                    ggCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACGTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTc                                                              <  1:1952093/150‑1 (MQ=255)
                                                                                    ggCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTc                                                              >  1:1995867/1‑150 (MQ=255)
                                                                                       tGCGCCGCTCGCCCGGTGCCAGGGGGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCAt                                                           <  1:2320528/150‑1 (MQ=255)
                                                                                       tGCGCCGCTCGCCCGGCGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGa                                                                                                          >  2:1210595/1‑103 (MQ=255)
                                                                                       tGCGCCGCTCGCCCGGCGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGa                                                                                                          <  1:1210595/103‑1 (MQ=255)
                                                                                          gccgcTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGc                                                                                               >  2:2075723/1‑111 (MQ=255)
                                                                                          gccgcTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGc                                                                                               <  1:2075723/111‑1 (MQ=255)
                                                                                           ccgcTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCt                                                       <  2:2255141/150‑1 (MQ=255)
                                                                                             gcTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTgg                                                     >  2:1357491/1‑150 (MQ=255)
                                                                                                cGCCCGGTGCCAGGGCGGCCTAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGGATGGTCCATTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTGGCgg                                                  <  2:76912/150‑1 (MQ=255)
                                                                                                cGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAg                                                                     >  1:113521/1‑131 (MQ=255)
                                                                                                cGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAg                                                                     <  2:113521/131‑1 (MQ=255)
                                                                                                   ccGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTGGCGGGgc                                               >  1:2364800/1‑150 (MQ=255)
                                                                                                   ccGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTGGCGGGgc                                               >  1:1637662/1‑150 (MQ=255)
                                                                                                     ggTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACGGCGCTCGAAGGTTCGTCCATCGCTGGCGGGgcgc                                             <  1:1280257/150‑1 (MQ=255)
                                                                                                       tGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTGGCGGGGCGCGa                                           <  1:121453/150‑1 (MQ=255)
                                                                                                        gCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTGGCGGGGCGCGAc                                          <  1:2333520/150‑1 (MQ=255)
                                                                                                          cAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTGGCGGGGCGCGACTg                                        >  1:2142557/1‑150 (MQ=255)
                                                                                                           aGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTGGCGGGGCGCGACTGc                                       >  2:776265/1‑150 (MQ=255)
                                                                                                             ggcggcCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTGGCGGGGCGCGACTGccc                                     <  1:1548830/150‑1 (MQ=255)
                                                                                                                  cccAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCa                                                                                             >  2:1624947/1‑89 (MQ=255)
                                                                                                                  cccAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCa                                                                                             <  1:1624947/89‑1 (MQ=255)
                                                                                                                   ccAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCGATCGCTGTCGGGGCGCGACTGCCCGCCAAt                               >  1:1924584/1‑150 (MQ=255)
                                                                                                                    cAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTGGCGGGGCGCGACTGCCCGCCAATc                              <  1:155174/150‑1 (MQ=255)
                                                                                                                     aGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGacac                                                                                                       <  2:1498946/76‑1 (MQ=255)
                                                                                                                     aGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGacac                                                                                                       >  1:1498946/1‑76 (MQ=255)
                                                                                                                          gg‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTGGCGGGGCGCGACTGCCCGCCAATCAGCaaa                        <  2:1637662/150‑1 (MQ=255)
                                                                                                                             atctagccgaGCTGGCGGCAGCCACCGCGGCATCGGCATCGGCCAGTGTAGAGGACGACGTACGTCCGGTTGAAGT‑GCCAGTCACCGTGTTACGGCGCTCGAAGATACGACAAACGCTGGTTGGGCGCGACTGCCCGCCAATCAGCAAAg                       <  1:2351547/145‑1 (MQ=255)
                                                                                                                              aaaCGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTGGCGGGGCGCGACTGCCCGCCAATCAGCAAAgg                      >  2:1328247/1‑150 (MQ=255)
                                                                                                                                  gccgcGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTGGCGGGGCGCGACTGCCCGCCAATCAGCaaa                        <  2:1156266/144‑1 (MQ=255)
                                                                                                                                  gccgcGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTGGCGGGGCGCGACTGCCCGCCAATCAGCaaa                        >  1:1156266/1‑144 (MQ=255)
                                                                                                                                  gccgcGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTGGCGGGGCGCGACTGCCCGCCAATCAGCAAAggagat                  >  2:2351547/1‑146 (MQ=255)
                                                                                                                                  gccgcGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTGGCGGGGCGCGACTGCCCGCCAATCAGCAAAGGCAcc                  >  2:2216243/1‑150 (MQ=255)
                                                                                                                                  gccgcGCTGGCAGCAGCCACCGCGGCATCGGCATCGCCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTGGCGGGGCGCGACTGCCCGCCAATCAGCAAAGGCAcc                  >  2:439858/1‑150 (MQ=255)
                                                                                                                                   ccgcGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGGGCTCGAAGGTTCGTCCATCGCTGGCGGGGCGCGACTGCCCGCCAATCAGCAAAGGCACCt                 >  2:825429/1‑150 (MQ=255)
                                                                                                                                   ccgcGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTGGCGGGGCGCGACTGCCCGCCAATCAGCAAAGGCACCt                 <  2:1367742/150‑1 (MQ=255)
                                                                                                                                    cgcgGTGGCTGCTGCCAGCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTGGCgg                                                  <  1:1689592/116‑1 (MQ=255)
                                                                                                                                    cgcgGTGGCTGCTGCCAGCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTGGCgg                                                  >  2:1689592/1‑116 (MQ=255)
                                                                                                                                     gcggTGGCTGCTGCCAGCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGacac                                                                                                       >  2:2308188/1‑62 (MQ=25)
                                                                                                                                     gcggTGGCTGCTGCCAGCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTGGCGGGGCGCGACTGCCCGCCAATCAGCAAAGGCACCTGc               >  1:236412/1‑150 (MQ=255)
                                                                                                                                     gcgcTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTGGCGGGGCGCGACTGCCCGCCAATCAGCAAAGGCACCTGc               <  2:2364800/150‑1 (MQ=255)
                                                                                                                                        cTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTGGCGGGGCGCGACTGCCCGCCAATCAGCAAAGGCACCTGCAAc            <  2:656329/150‑1 (MQ=255)
                                                                                                                                           gcagcaGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACAACTTCGCCAGTCAAGGAGGTACAGCGCGCGAAGGTTCGTCCATCGCTGGCGGGGCGCGACTGCCCGCCAATCAGCAAAGGCACCTGCAACATg         >  1:900368/1‑150 (MQ=255)
                                                                                                                                             agcagcCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGACACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTGGCGGGGCGCGACTGCCCGCCAATCAGCAAAGGCACCTGCAACATGCt       >  2:1563842/1‑150 (MQ=255)
                                                                                                                                                  ccACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTGGCGGGGCGCGACTGCCCGCCAATCAGCAAAGGCACCTGCAACATGCTCATTc  <  1:2306328/150‑1 (MQ=255)
                                                                                                                                                     ccGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTGGCGGGGCGCGACTGCCCGCCAAt                               >  1:650926/1‑118 (MQ=255)
                                                                                                                                                     ccGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTGGCGGGGCGCGACTGCCCGCCAAt                               <  2:650926/118‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
AGTTGGCCATGGCCCCGGTTTCACCGGCGGCGGCGATGAACTCGGCGGCCCTCGCCTGCAACAGATCAGCGCCTGCCAGCAAGCGGCTGCGCCGCTCGCCCGGTGCCAGGGCGGCCCAGGCCGG‑‑AAACGCCGCGCTGGCAGCAGCCACCGCGGCATCGGCATCGGCCAGTGTGGCGG‑CGGCAGCCTGCGACACCACCTCGCCAGTCACCGGGTTACAGCGCTCGAAGGTTCGTCCATCGCTGGCGGGGCGCGACTGCCCGCCAATCAGCAAAGGCACCTGCAACATGCTCATTC  >  NC_002947/3797751‑3798044

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑