Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,995,572:1 +G 9.0% coding (1188/1272 nt) 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_0029474,995,5721.G 9.0% 506.8 / 45.2 133coding (1188/1272 nt)bkdBbranched‑chain alpha‑keto acid dehydrogenase complex lipoamide acyltransferase
Reads supporting (aligned to +/- strand):  ref base . (56/65);  new base G (9/3);  total (65/68)
Fisher's exact test for biased strand distribution p-value = 7.26e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.18e-01

Only 100 of 136 total aligned reads displayed.

CCTGGGTGGCATTGTCAGCACGCCGGTGGTCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGATCATCGGCGGCGGCCCTGGCGGCTATGTAGCA  >  NC_002947/4995425‑4995721
                                                                                                                                                    |                                                                                                                                                      
ccTGGGTGGCATTGTCAGCACGCCGGTGGTCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGt                                                                                                                                                      <  2:2138034/150‑1 (MQ=255)
 cTGGGTGGCATTGTCAGCACGCCGGTGGTCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCCCGTg                                                                                                                                                     >  2:2133244/1‑150 (MQ=255)
  tGGGTGGCATTGTCAGCACGCCGGTGGTCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTgg                                                                                                                                                    <  2:1187746/150‑1 (MQ=255)
       ggCATTGTCAGCACGCCGGTGGTCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTTGGACTGCCAGTGGTTATTGATGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACGCGGTGGTCGaa                                                                                                                                               <  2:675125/150‑2 (MQ=255)
       ggCATTGTCAGCACGCCGGTGGTCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGAt                                                                                                                                               <  2:174739/150‑1 (MQ=255)
       ggCATTGTCAGCACGCCGGTGGTCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGAt                                                                                                                                               <  1:253320/150‑1 (MQ=255)
       ggCATTGTCAGCACGCCGGTGGTCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACGCGGTGGTCGaa                                                                                                                                               <  1:2394909/150‑2 (MQ=255)
       ggCATTGTCAGCACGCCGGGGGTCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGAt                                                                                                                                               <  2:914048/150‑1 (MQ=255)
        gCATTGTCAGCACGCCGGTGGTCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATg                                                                                                                                              <  2:2384407/150‑1 (MQ=255)
         cATTGTCCGAACGCTGGTGGCCAGGTCTCCGGTGTTGTCAAGCGGCGGGGTCAACCGCATGGTCAGCCAGCGCTGGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATgg                                                                                                                                             <  1:2253772/150‑1 (MQ=255)
          aTTGTCAGCACGCCGGTGGTCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGc                                                                                                                                            >  2:1002412/1‑150 (MQ=255)
           ttGTCAGCACGCCGGTGGTCAACACCCCGGGAGTGGCAATCGTCGGGGGCAACCGCATGGTCGAACGGCCAGGGGTGATCGACGGACAGATCGTCGTGCGCAAGATGATGAACCTGGCCAGCTCGTTGGACCACCGCGTGGTGGATGGCa                                                                                                                                           >  2:79325/1‑150 (MQ=255)
           ttGTCAGCACGCCGGTGGTCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCa                                                                                                                                           <  1:676010/150‑1 (MQ=255)
           ttGTCAGCACGCCGGTGGTCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCa                                                                                                                                           <  2:1795920/150‑1 (MQ=255)
             gTCAGCACGCCGGTGGTCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATg                                                                                                                                         <  2:2192424/150‑1 (MQ=255)
                     gCCGGTGGTCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGCTCGAACGGCCAGTGGTGATCGACGGCCAAAGGATCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATgccgc                                                                                                                                 >  1:411905/1‑150 (MQ=255)
                       cGGTGGTCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGccc                                                                                                                               <  1:988602/150‑1 (MQ=255)
                            gtCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCTCAAGATGATGAACGTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCACCCTGTTc                                                                                                                          >  2:1307059/1‑150 (MQ=255)
                             tCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGt                                                                                                                            >  2:1497454/1‑147 (MQ=255)
                             tCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGt                                                                                                                            <  1:1497454/147‑1 (MQ=255)
                               aaCACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATc                                                                                                                       <  1:1781890/150‑1 (MQ=255)
                                acacCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATcc                                                                                                                      <  1:1162175/150‑1 (MQ=255)
                                acacCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATcc                                                                                                                      <  1:1630874/150‑1 (MQ=255)
                                    cccGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAGCCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGc                                                                                                                  <  1:1215224/150‑1 (MQ=255)
                                    cccGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGc                                                                                                                  >  1:2204916/1‑150 (MQ=255)
                                       ggAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGt                                                                                                               >  1:1838572/1‑150 (MQ=255)
                                       ggAAGTGGCAATCGTCGGGGTCAACCGCATGGTCAAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGt                                                                                                               >  1:1678667/1‑150 (MQ=255)
                                            tGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTActga                                                                                                          <  2:1838572/150‑4 (MQ=255)
                                             ggCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTgg                                                                                                         >  1:91958/1‑150 (MQ=255)
                                                aaTCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCt                                                                                                      >  1:1550525/1‑150 (MQ=255)
                                                aaTCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGCGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCt                                                                                                      >  2:857748/1‑150 (MQ=255)
                                                    gtcgGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGCGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTc                                                                                                  <  1:857748/150‑1 (MQ=255)
                                                     tcgGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCg                                                                                                 >  1:143039/1‑150 (MQ=255)
                                                     tcgGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCACTGTTCATCCAGGCCGTACGTGGCCTGCTCg                                                                                                 >  2:1248310/1‑150 (MQ=255)
                                                        gggTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGaac                                                                                              <  1:1950499/150‑1 (MQ=255)
                                                          gTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGaacaa                                                                                            <  1:1406495/150‑1 (MQ=255)
                                                           ccAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGaacaac                                                                                           <  2:670990/149‑1 (MQ=255)
                                                            cAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAAcc                                                                                          <  1:1390812/150‑1 (MQ=255)
                                                             aaCCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCa                                                                                                                                           >  2:780361/1‑100 (MQ=255)
                                                             aaCCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCa                                                                                                                                           <  1:780361/100‑1 (MQ=255)
                                                               ccGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGa                                                                                                                                                >  2:2293639/1‑93 (MQ=255)
                                                               ccGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGa                                                                                                                                                <  1:2293639/93‑1 (MQ=255)
                                                               ccGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGc                                                                                                        >  2:1338948/1‑133 (MQ=255)
                                                               ccGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGc                                                                                                        <  1:1338948/133‑1 (MQ=255)
                                                               ccGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCgc                                                                                       <  2:1567816/150‑1 (MQ=255)
                                                               ccGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCgc                                                                                       <  1:2110230/150‑1 (MQ=255)
                                                                      gTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGcctgcctg                                                                                >  2:1195043/1‑150 (MQ=255)
                                                                         gAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTCTCGAACAACCCGCCTGCCTGTTCg                                                                            >  1:1869803/1‑150 (MQ=255)
                                                                           aCGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGt                                                                           <  2:1550525/150‑1 (MQ=255)
                                                                               ccAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAg                                                                       >  1:1710341/1‑150 (MQ=255)
                                                                                cAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGc                                                                                                                                            <  1:503560/80‑1 (MQ=255)
                                                                                cAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGc                                                                                                                                            >  2:503560/1‑80 (MQ=255)
                                                                                cAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGt                                                                      <  2:1678667/150‑1 (MQ=255)
                                                                                  gtggtgATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGa                                                                    <  2:143039/150‑1 (MQ=255)
                                                                                   tggtgATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAg                                                                   >  2:38815/1‑150 (MQ=255)
                                                                                                 cAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTg                                                                          <  1:733024/129‑1 (MQ=255)
                                                                                                   gATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATc                                                  <  1:38815/150‑1 (MQ=255)
                                                                                                   gATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATc                                                  <  2:326382/150‑1 (MQ=255)
                                                                                                        tcgtGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGAc                                             >  2:1662364/1‑150 (MQ=255)
                                                                                                              gcAAGCTGATGAACCTGTCCAGCTCGTATGACCACGGCTTGGGCGATGGCATGGATGCCGCATTGTAACGCCGATCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCt                                       <  1:1002412/150‑1 (MQ=255)
                                                                                                                aaGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACa                                                         >  1:388563/1‑130 (MQ=255)
                                                                                                                 aGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGtt                                    >  1:2193792/1‑150 (MQ=255)
                                                                                                                 aGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGACCTGCTCGAACAACCCGCCTGCCTGTTCGTGGa                                                                        <  2:2292262/115‑1 (MQ=255)
                                                                                                                 aGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGACCTGCTCGAACAACCCGCCTGCCTGTTCGTGGa                                                                        >  1:2292262/1‑115 (MQ=255)
                                                                                                                    tgatgaACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGat                                 <  2:2149283/150‑1 (MQ=255)
                                                                                                                     gatgaACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGatc                                <  2:860358/150‑1 (MQ=255)
                                                                                                                        gaACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGatcatc                             <  1:2415830/150‑1 (MQ=255)
                                                                                                                        gaACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGatcatc                             <  2:2372064/150‑1 (MQ=255)
                                                                                                                            cTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGATCATcggcg                         >  2:344859/1‑150 (MQ=255)
                                                                                                                            cTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGATCATcggcg                         >  2:956588/1‑150 (MQ=255)
                                                                                                                                 cAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGATCATCggcggcggc                    >  1:2014413/1‑150 (MQ=255)
                                                                                                                                 cAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGATCATCggcggcggc                    >  1:2006213/1‑150 (MQ=255)
                                                                                                                                 cAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGATCATCggcggcggc                    >  1:2176747/1‑150 (MQ=255)
                                                                                                                                 cAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGCTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCt                                       >  1:531233/1‑131 (MQ=255)
                                                                                                                                  aGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTAcc                                         <  2:1284325/128‑1 (MQ=255)
                                                                                                                                  aGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTAcc                                         >  1:1284325/1‑128 (MQ=255)
                                                                                                                                    tgccatCGACCACGCGGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGATCATCGGCGGCGGCCCt                 >  1:2106349/6‑150 (MQ=255)
                                                                                                                                    cTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGATCATCGGCGGCGGcc                   >  1:1050349/1‑148 (MQ=255)
                                                                                                                                     gccatCGACCACGCGGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGATCATCGGCGGCGGCCCTg                <  2:2106349/146‑1 (MQ=255)
                                                                                                                                       gTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGc                                                                  <  1:795269/99‑1 (MQ=255)
                                                                                                                                       gTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGa                                              <  2:2073537/118‑1 (MQ=255)
                                                                                                                                       catCGACCACGCGGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGATCATCGGCGGCGGCCCTggc              >  2:2337864/3‑150 (MQ=255)
                                                                                                                                        ttCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGaacaac                                                                                           >  2:2223634/1‑73 (MQ=255)
                                                                                                                                        atCGACCACGCGGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCTTGGCTGTGAGGGCGGTAGCAAGACAAGTGGTTATCCAGACTACCCTGTTGATCCTCGGCGGCGGCCCTggcg             >  1:876807/2‑150 (MQ=255)
                                                                                                                                        atCGACCACGCGGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGATCATCGGCGGCGGCCCTggcg             >  2:1436848/2‑150 (MQ=255)
                                                                                                                                         tCGACCACGCGGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGt                                                                           >  2:1181482/1‑88 (MQ=255)
                                                                                                                                            accaccGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGATCATCGGCGGCGGCCCTGGCGGCTa         <  2:1710341/150‑1 (MQ=255)
                                                                                                                                             ccaccGCGTGGTCGCTGGCACACGTTCCGCCCTGTTCATCCAGGCCGTACGTGGCCTCCTCGCCCAACCCGCCTGCCTGTTCGTCGAGTGTGCATGCAACAGATTATCCAGACTACCCTGTTGATCATCGGCGGCGGCCCTGGCGGCTAt        <  1:1307059/150‑1 (MQ=255)
                                                                                                                                              caccGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGATCATCGGCGGCGGCCCTGGCGGCTATg       >  1:1326954/1‑150 (MQ=255)
                                                                                                                                               accGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGtt                                    >  2:841130/1‑120 (MQ=255)
                                                                                                                                               accGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGtt                                    <  1:841130/120‑1 (MQ=255)
                                                                                                                                               accGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGATCATCGGCGGCGGCCCTGGCGGCTATGt      >  2:1746524/1‑150 (MQ=255)
                                                                                                                                               accGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGATCATCGGCGGCGGCCCTGGCGGCTATGt      >  2:2108077/1‑150 (MQ=255)
                                                                                                                                               accGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGATCATCGGCGGCGGCCCTGGCGGCTATGt      <  1:1195043/150‑1 (MQ=255)
                                                                                                                                                ccGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGATCATCGGCGGCGGCCCTggcg             >  1:644199/1‑142 (MQ=255)
                                                                                                                                                ccGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGATCATCGGCGGCGGCCCTggcg             <  2:644199/142‑1 (MQ=255)
                                                                                                                                                ccGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGATCATCGGCGGCGGCCCTGGCGGCTATGTag    <  2:1869803/150‑1 (MQ=255)
                                                                                                                                                   cgTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGATCATCGGCGGCGGCCCTGGCGGCTATGTagca  <  2:750869/150‑1 (MQ=255)
                                                                                                                                                   cgTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGATCATCGGCGGCGGCCCTGGCGGCTATGTagca  <  2:2193792/150‑1 (MQ=255)
                                                                                                                                                   cgTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGGTGATCATCGGGGGCGGCCCTGGCGGCTATGTagca  <  2:411905/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
CCTGGGTGGCATTGTCAGCACGCCGGTGGTCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGATCATCGGCGGCGGCCCTGGCGGCTATGTAGCA  >  NC_002947/4995425‑4995721

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