Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,000,521 C→G 6.0% R396R (CGC→CGG) ‡ bkdB → branched‑chain alpha‑keto acid dehydrogenase complex lipoamide acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,000,5210CG 6.0% 727.7 / 27.4 215R396R (CGC→CGG) ‡bkdBbranched‑chain alpha‑keto acid dehydrogenase complex lipoamide acyltransferase
Reads supporting (aligned to +/- strand):  ref base C (81/121);  new base G (8/5);  total (89/126)
Fisher's exact test for biased strand distribution p-value = 1.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

Only 100 of 221 total aligned reads displayed.

GCCCTGGGTGGCATTGTCAGCACGCCGGTGGTCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGATCATCGGCGGCGGCCCTGGCGGCTATGTAGCA  >  NC_002947/5000372‑5000670
                                                                                                                                                     |                                                                                                                                                     
   cTGGGTGGCATTGTCAGCACGCCGGTGGTCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTg                                                                                                                                                    <  1:4536039/150‑1 (MQ=255)
      ggTGGCATTGTCAGCACGCCGGTGGTCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTc                                                                                                                                                 <  2:2790621/150‑1 (MQ=255)
        tGGCATTGTCAGCACGCCGGTGGTCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGa                                                                                                                                               <  1:4011683/150‑1 (MQ=255)
         ggCATTGTCAGCACGCCGGTGGTCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGAt                                                                                                                                              <  1:2627860/150‑1 (MQ=255)
          gCATTGTCATCTCGCCGGTGGTCTACACCCCGGAAGTGGCAATCGTCGTGGTCTAGCGCATGGTCGAACGGGCAGTGGTGATCGACGGCGAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCTCGTGGTCGATg                                                                                                                                             <  2:3014427/150‑1 (MQ=255)
          gCATTGTCAGCACGGCGGTGGTCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGATCGACCACGCGGTGGTCGaac                                                                                                                                             <  2:4711309/150‑3 (MQ=255)
          gCATTGTCAGCACGCCGGTGGTCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACGCGGTGGTCGaac                                                                                                                                             <  1:1085615/150‑3 (MQ=255)
            aTTGTCAGCACGCCGGTGGTCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGc                                                                                                                                           >  1:3796511/1‑150 (MQ=255)
            aTTGACAGCACGCCGGTGGTCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGATCGACCACCGCGTGCTCGATGGc                                                                                                                                           >  2:4971043/1‑150 (MQ=255)
             ttGTCAGCACGCCGGTGGTCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCa                                                                                                                                          >  1:1580009/1‑150 (MQ=255)
             ttGTCAGCACGCCGGTGGTCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCa                                                                                                                                          <  2:5300704/150‑1 (MQ=255)
               gTCAGCACGCCGGTGGTCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATg                                                                                                                                        <  1:3680667/150‑1 (MQ=255)
                 cAGCACGCCGGTGGTCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGa                                                                                                                                      >  1:2700004/1‑150 (MQ=255)
                   gCACGCCGGTGGTCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCAtccatg                                                                                                                                    <  1:435038/150‑6 (MQ=255)
                     aCGCCGGTGGTCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATgcc                                                                                                                                  <  1:3308100/150‑1 (MQ=255)
                      cGCCGGTGGTCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATgccg                                                                                                                                 <  1:430239/150‑1 (MQ=255)
                      cGCCGGTGGTCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATgccg                                                                                                                                 <  2:140833/150‑1 (MQ=255)
                          ggtggtCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCCTGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCt                                                                                                                             <  2:4137191/150‑1 (MQ=255)
                          ggtggtCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCt                                                                                                                             >  1:2982756/1‑150 (MQ=255)
                            tggtCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGt                                                                                                                           <  2:1200023/150‑1 (MQ=255)
                              gtCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTc                                                                                                                         <  2:457264/150‑1 (MQ=255)
                               tCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCa                                                                                                                        <  2:5607025/150‑1 (MQ=255)
                                   cacCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCa                                                                                                                    >  1:1057682/1‑150 (MQ=255)
                                   cacCCCGGAAGTGGCAATCGTCGGGGTCAACCGCAGGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCa                                                                                                                    <  2:3404386/150‑1 (MQ=255)
                                    acCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAg                                                                                                                   <  2:4890698/150‑1 (MQ=255)
                                    acCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAg                                                                                                                   <  2:3972938/150‑1 (MQ=255)
                                      cccGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGc                                                                                                                 >  1:640270/1‑150 (MQ=255)
                                      cccGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGc                                                                                                                 <  1:5050593/150‑1 (MQ=255)
                                      cccGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGc                                                                                                                 >  2:2151295/1‑150 (MQ=255)
                                      cccGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGc                                                                                                                 >  1:4083215/1‑150 (MQ=255)
                                       ccGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGcc                                                                                                                >  2:2497704/1‑150 (MQ=255)
                                       ccGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGcc                                                                                                                <  2:3138031/150‑1 (MQ=255)
                                        cGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCg                                                                                                               >  2:4277234/1‑150 (MQ=255)
                                             gTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGt                                                                                                          <  1:202399/150‑1 (MQ=255)
                                              tGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTg                                                                                                         <  1:3182151/150‑1 (MQ=255)
                                                 cAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGcc                                                                                                      <  1:3377515/150‑1 (MQ=255)
                                                  aaTCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCt                                                                                                     >  2:4521192/1‑150 (MQ=255)
                                                  aaTCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCt                                                                                                     <  1:3158820/150‑1 (MQ=255)
                                                  aaTCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCt                                                                                                     <  1:3811197/150‑1 (MQ=255)
                                                    tcgtcgGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGc                                                                                                   <  1:5214407/150‑1 (MQ=255)
                                                         ggggTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGaa                                                                                              <  1:5276528/150‑1 (MQ=255)
                                                         ggggTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGaa                                                                                              <  2:2645979/150‑1 (MQ=255)
                                                               aaCCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAAccc                                                                                        >  1:2901746/1‑150 (MQ=255)
                                                                  cGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCgcc                                                                                     >  1:4894030/1‑150 (MQ=255)
                                                                     aTGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCgcctgc                                                                                  <  1:1096248/150‑1 (MQ=255)
                                                                     aTGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCgcctgc                                                                                  <  2:5059271/150‑1 (MQ=255)
                                                                        gTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTTGCCTGCTCGAACAACCCGcctgcctg                                                                               <  1:131293/150‑1 (MQ=255)
                                                                           gAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTc                                                                            >  2:1921033/1‑150 (MQ=255)
                                                                              cGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGTATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTg                                                                         <  2:1077636/150‑1 (MQ=255)
                                                                              cGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTg                                                                         <  2:4905406/150‑1 (MQ=255)
                                                                              cGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTg                                                                         <  2:2404046/150‑1 (MQ=255)
                                                                                  cAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGt                                                                     >  2:3937351/1‑150 (MQ=255)
                                                                                    gtggtgATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGa                                                                   <  2:2700004/150‑1 (MQ=255)
                                                                                     tggtgATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAg                                                                  <  2:1660592/150‑1 (MQ=255)
                                                                                           tCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCa                                                            >  2:4606034/1‑150 (MQ=255)
                                                                                             gACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAAc                                                          >  1:3248854/1‑150 (MQ=255)
                                                                                                ggCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGa                                                       <  1:3834245/150‑1 (MQ=255)
                                                                                                ggCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGa                                                       <  1:4421379/150‑1 (MQ=255)
                                                                                                   cAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTa                                                    >  1:5061047/1‑150 (MQ=255)
                                                                                                      aTCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATcc                                                 >  2:4369498/1‑150 (MQ=255)
                                                                                                       tcgtcgTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCa                                                <  1:710926/150‑1 (MQ=255)
                                                                                                       tcgtcgTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCa                                                <  2:400278/150‑1 (MQ=255)
                                                                                                          tcgtGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGAc                                             <  2:1580009/150‑1 (MQ=255)
                                                                                                          tcgtGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGAc                                             >  2:2779933/1‑150 (MQ=255)
                                                                                                              gcgcAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTAcc                                         <  1:283696/150‑1 (MQ=255)
                                                                                                                gcAAGATGATGAACCTGTCCATCTCGTTCGACCACCGCGTTGTGGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCt                                       <  2:718537/150‑1 (MQ=255)
                                                                                                                   aGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGtt                                    <  1:200804/150‑1 (MQ=255)
                                                                                                                   aGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGtt                                    <  1:4583627/150‑1 (MQ=255)
                                                                                                                   aGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGtt                                    <  1:2297503/150‑1 (MQ=255)
                                                                                                                    gatgagGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTg                                   >  2:1025326/1‑150 (MQ=255)
                                                                                                                     atgatgAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGa                                  <  1:2790386/150‑1 (MQ=255)
                                                                                                                       gatgaACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGatc                                >  2:653404/1‑150 (MQ=255)
                                                                                                                         tgaACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGatcat                              >  1:2389596/1‑150 (MQ=255)
                                                                                                                          gaACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGatcatc                             <  1:2218505/150‑1 (MQ=255)
                                                                                                                             ccTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGATCATcggc                          >  2:589267/1‑150 (MQ=255)
                                                                                                                              cTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGa                                                       >  2:5170090/1‑120 (MQ=255)
                                                                                                                              cTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGa                                                       <  1:5170090/120‑1 (MQ=255)
                                                                                                                              cTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGATCATcggcg                         <  2:3344001/150‑1 (MQ=255)
                                                                                                                                   cAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCAGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGATCATCggcggcggc                    >  2:5501020/1‑150 (MQ=255)
                                                                                                                                      cTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGATCATCGGCGGCGGCCCt                 >  2:4089819/1‑150 (MQ=255)
                                                                                                                                      cTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGATCATCGGCGGCGGCCCt                 <  2:1324968/150‑1 (MQ=255)
                                                                                                                                      cTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGATCATCGGCGGCGGCCCt                 >  1:2452754/1‑150 (MQ=255)
                                                                                                                                         gTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGATCATCGGCGGCGGCCCTggc              <  2:917487/150‑1 (MQ=255)
                                                                                                                                         catCGACCACGCGGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGATCATCGGCGGCGGCCCTggc              >  1:4783693/3‑150 (MQ=255)
                                                                                                                                          ttCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGATCATCGGCGGCGGCCCTggcg             <  2:3866759/150‑1 (MQ=255)
                                                                                                                                          atCGACCACGCGGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTgg                                                                                                        <  1:3905028/58‑1 (MQ=37)
                                                                                                                                          atCGACCACGCGGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGATCATCGGCGGCGGCCCTggcg             >  1:572030/2‑150 (MQ=255)
                                                                                                                                           tCGACCACGCGGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCACCAGATTATCCAGACTACCCTGTTGATCATCGGCGGCGGCCCTggcgg            >  2:3621221/1‑150 (MQ=255)
                                                                                                                                             gACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGATCATCGGCGGCGGCCCTGGCGGCt          >  1:1936570/1‑150 (MQ=255)
                                                                                                                                             gACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGATCATCGGCGGCGGCCCTGGCGGCt          >  2:3282556/1‑150 (MQ=255)
                                                                                                                                              accaccGCGTGGTCGCTGGCATGGATGCCGCCCTGGTCATCCAGGCCGTACGTGGCCTGCTCGAACACCCCGCCTGCCTGTTCGTGGAGTGAGCCTGCACCAGCGTATCCAGACTCCCCTGTTGATCCTCGGCGGCGGCCCTGGCGGCTa         <  2:633045/150‑1 (MQ=255)
                                                                                                                                              accaccGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTCTCCAGACTACCCTGTTGATCATCGGCGGCGGCCCTGGCGGCTa         >  2:5371590/1‑150 (MQ=255)
                                                                                                                                              accaccGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGATCATCGGCGGCGGCCCTGGCGGCTa         <  1:882232/150‑1 (MQ=255)
                                                                                                                                               ccaccGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGATCATcgg                           <  2:2236369/131‑1 (MQ=255)
                                                                                                                                                caccGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGATCATCGGCGGCGGCCCTGGCGGCTATg       >  1:5157657/1‑150 (MQ=255)
                                                                                                                                                caccGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGATCATCGGCGGCGGCCCTGGCGGCTATg       >  1:4235975/1‑150 (MQ=255)
                                                                                                                                                caccGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGATCATCGGCGGCGGCCCTGGCGGCTATg       <  1:3937351/150‑1 (MQ=255)
                                                                                                                                                caccGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGATCATCGGCGGCGGCCCTGGCGGCTATg       >  1:3543815/1‑150 (MQ=255)
                                                                                                                                                 accGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGATCATCGGCGGCGGCCCTGGCGGCTATGt      >  2:3809601/1‑150 (MQ=255)
                                                                                                                                                     cgTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGATCATCGGCGGCGGCCCTGGCGGCTATGTagca  <  1:2151295/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCCCTGGGTGGCATTGTCAGCACGCCGGTGGTCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGATCATCGGCGGCGGCCCTGGCGGCTATGTAGCA  >  NC_002947/5000372‑5000670

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: