Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,995,572:1 +G 6.2% 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 6.2% 454.9 / 22.2 112coding (1188/1272 nt)bkdBbranched‑chain alpha‑keto acid dehydrogenase complex lipoamide acyltransferase
Reads supporting (aligned to +/- strand):  ref base . (51/54);  new base G (3/4);  total (54/58)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.67e-01

Only 100 of 114 total aligned reads displayed.

GCCCTGGGTGGCATTGTCAGCACGCCGGTGGTCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGATCATCGGCGGCGGCCCTGGCGGCTATGTAGCA  >  NC_002947/4995423‑4995721
                                                                                                                                                      |                                                                                                                                                     
gcCCTGGGTGGCATTGTCAGCACGCCGGTGGTCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACcgc                                                                                                                                                        >  2:927293/1‑150 (MQ=255)
  ccTGGGTGGCATTGTCAGCACGCCGGTGGTCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGt                                                                                                                                                     <  2:1431619/150‑1 (MQ=255)
  ccTGGGTGGCATTGTCAGACCGCCGGTGGTCAACACACCGGAATTGACAATGGTCGGGGTAAACCGCATGGTCGAACGGCGAGCAGCAATCGACGACCAGATCGTCGTGCGCAAGATGATGAACCTGTACAACTCGTTCGACCACCGCGt                                                                                                                                                     >  2:863530/1‑150 (MQ=255)
        tGGCATTGTCAGCACGCCGGTGGTCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACGCGGTGGTCGa                                                                                                                                               <  1:1903329/150‑1 (MQ=255)
          gCATTGTCAGCATGACGGTGGTCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGCACGCCCAGTCGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACGCGGTGGTCGaac                                                                                                                                             <  1:495692/150‑3 (MQ=255)
          gCATTGTCAGCACGCCGGTGGTCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGGCCACGCGGTGGTCGaac                                                                                                                                             <  1:1858170/150‑3 (MQ=255)
          gCACTGCCAGCACGCCGGTGGTCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTTGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTTGTCGaac                                                                                                                                             <  2:477121/150‑3 (MQ=255)
             ttGTCAGCACGCCGGTGGTCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCa                                                                                                                                          <  2:1537398/150‑1 (MQ=255)
             ttGTCAGCACGCCGGTGGTCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCa                                                                                                                                          >  2:2098556/1‑150 (MQ=255)
                  aGCACGCCGGTGGTCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGAt                                                                                                                                     >  1:1672318/1‑150 (MQ=255)
                    cACGCCGGTGGTCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATgc                                                                                                                                   >  2:915870/1‑150 (MQ=255)
                     aCGCCGGTGGTCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATgcc                                                                                                                                  >  1:2319537/1‑150 (MQ=255)
                     aCGCCGGTGGTCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATgcc                                                                                                                                  <  2:2296188/150‑1 (MQ=255)
                       gCCGGTGGTCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATgccgc                                                                                                                                <  2:2058191/150‑1 (MQ=255)
                         cGGTGGTCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGccc                                                                                                                              >  2:1976054/1‑150 (MQ=255)
                          ggtggtCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCt                                                                                                                             <  1:1528119/150‑1 (MQ=255)
                              gtCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTc                                                                                                                         <  1:1350892/150‑1 (MQ=255)
                              gtCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTc                                                                                                                         >  1:1905730/1‑150 (MQ=255)
                              gtCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCCGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTc                                                                                                                         <  1:2031256/150‑1 (MQ=255)
                               tCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCa                                                                                                                        >  1:1369323/1‑150 (MQ=255)
                                  acacCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATcc                                                                                                                     >  1:938204/1‑150 (MQ=255)
                                   cacCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTggc                                                                                                                                                  >  1:1521557/1‑119 (MQ=255)
                                   cacCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTggc                                                                                                                                                  <  2:1521557/120‑2 (MQ=255)
                                     ccccGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAgg                                                                                                                  <  2:2124323/150‑1 (MQ=255)
                                      cccGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGtt                                                                                                                          >  2:1179677/1‑141 (MQ=255)
                                      cccGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGtt                                                                                                                          <  1:1179677/141‑1 (MQ=255)
                                      cccGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGc                                                                                                                 <  1:2048172/150‑1 (MQ=255)
                                      cccGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGc                                                                                                                 >  2:601910/1‑150 (MQ=255)
                                       ccGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCt                                                                                                                             <  1:1065066/137‑1 (MQ=255)
                                       ccGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCt                                                                                                                             >  2:1065066/1‑137 (MQ=255)
                                       ccGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGcc                                                                                                                >  2:1177941/1‑150 (MQ=255)
                                              tGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGTCGCCCTGTTCATCCAGGCCGTACGTg                                                                                                         <  2:1781518/150‑1 (MQ=255)
                                              tGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTg                                                                                                         <  1:826890/150‑1 (MQ=255)
                                                gCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGc                                                                                                       <  1:1456708/150‑1 (MQ=255)
                                                gCAATCGTCGGGGTCAACCGCATGGTCGAACGCAAAGTGGTGTTCGACGTTCATATCGTCGTGCGCAAGCTGCTGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGc                                                                                                       <  1:540575/150‑1 (MQ=255)
                                                   aTCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTg                                                                                                    >  1:933788/1‑150 (MQ=255)
                                                    tcgtcgGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGc                                                                                                   >  1:904123/1‑150 (MQ=255)
                                                      gtcgGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTc                                                                                                 <  2:2347544/150‑1 (MQ=255)
                                                          gggTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGaac                                                                                             >  2:353215/1‑150 (MQ=255)
                                                            gTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGaacaa                                                                                           <  2:896389/150‑1 (MQ=255)
                                                               aaCCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAAccc                                                                                        <  2:933788/150‑1 (MQ=255)
                                                                  cGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCgcc                                                                                     <  2:1703467/150‑1 (MQ=255)
                                                                        gTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAgg                                                                                                                  >  1:196618/1‑115 (MQ=255)
                                                                        gTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAgg                                                                                                                  <  2:196618/115‑1 (MQ=255)
                                                                        gTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGcctgcctg                                                                               >  1:1692824/1‑150 (MQ=255)
                                                                           gAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTc                                                                            <  2:1265569/150‑1 (MQ=255)
                                                                           gAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTc                                                                            >  1:108650/1‑150 (MQ=255)
                                                                            aaCGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGt                                                                                                                           >  1:788004/1‑102 (MQ=255)
                                                                            aaCGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGt                                                                                                                           <  2:788004/102‑1 (MQ=255)
                                                                               ggCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACTACCGCGTGGTCGATGTCATGGATACCGCCCTGTTGCTCCAGGCCGTACGTGGGCTGCTCGAACAACCGGCCTGCCTGTTCGTgg                                                                        >  2:461751/1‑150 (MQ=255)
                                                                               ggCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTgg                                                                        <  2:1905730/150‑1 (MQ=255)
                                                                                 ccAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTg                                                                         >  2:1928954/1‑147 (MQ=255)
                                                                                 ccAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTg                                                                         <  1:1928954/147‑1 (MQ=255)
                                                                                    gtggtgATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGAt                                                                                                                                              <  2:696279/75‑1 (MQ=255)
                                                                                    ggggtgATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGAt                                                                                                                                              >  1:696279/3‑75 (MQ=255)
                                                                                                ggCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGcc                                                                                                                >  1:1786121/1‑93 (MQ=255)
                                                                                                ggCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGcc                                                                                                                <  2:1786121/93‑1 (MQ=255)
                                                                                                ggCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGa                                                       >  2:1282390/1‑150 (MQ=255)
                                                                                                ggCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGa                                                       <  1:1282390/150‑1 (MQ=255)
                                                                                                ggCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGa                                                       <  2:1692824/150‑1 (MQ=255)
                                                                                                  ccAGATCGTCGTGCGCAACATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGAtt                                                     <  1:461751/150‑1 (MQ=255)
                                                                                                   cAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACATGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTa                                                    >  2:161955/1‑150 (MQ=255)
                                                                                                     gATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATc                                                  <  1:353215/150‑1 (MQ=255)
                                                                                                     gATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATc                                                  >  1:1546539/1‑150 (MQ=255)
                                                                                                       tcgtcgTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCa                                                <  1:2229729/150‑1 (MQ=255)
                                                                                                         gtcgtGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGt                                                                                                              >  2:265032/1‑86 (MQ=255)
                                                                                                         gtcgtGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGt                                                                                                              <  1:265032/86‑1 (MQ=255)
                                                                                                          tcgtGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGAc                                             <  1:1177941/150‑1 (MQ=255)
                                                                                                          tcgtGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACATGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGAc                                             <  1:161955/150‑1 (MQ=255)
                                                                                                             tGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTAc                                          >  2:1582460/1‑150 (MQ=255)
                                                                                                                gcAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGt                                                                              >  2:654507/1‑111 (MQ=255)
                                                                                                                gcAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACa                                                         <  2:1031571/132‑1 (MQ=255)
                                                                                                                gcAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACa                                                         >  1:1031571/1‑132 (MQ=255)
                                                                                                                gcAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCt                                       >  2:539425/1‑150 (MQ=255)
                                                                                                                gcAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCt                                       >  1:624244/1‑150 (MQ=255)
                                                                                                                gcAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCt                                       <  2:904123/150‑1 (MQ=255)
                                                                                                                  aaGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGt                                     <  2:1546539/150‑1 (MQ=255)
                                                                                                                   aGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGtt                                    <  1:2318030/150‑1 (MQ=255)
                                                                                                                   aGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGtt                                    >  2:1067851/1‑150 (MQ=255)
                                                                                                                    gatgatGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGc                                                                                                       >  2:1934889/1‑82 (MQ=255)
                                                                                                                    gatgatGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGc                                                                                                       <  1:1934889/82‑1 (MQ=255)
                                                                                                                     atgatgAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGa                                                       <  2:1120103/129‑1 (MQ=255)
                                                                                                                     atgatgAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGa                                                       >  1:1120103/1‑129 (MQ=255)
                                                                                                                      tgatgaACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAAATGATCCAGACTACCCTGTTGat                                 >  2:1947201/1‑150 (MQ=255)
                                                                                                                           aaCCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGATCATcg                            <  1:1776222/150‑1 (MQ=255)
                                                                                                                             ccTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGt                                                                                                                           >  2:1617895/1‑53 (MQ=255)
                                                                                                                             ccTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGt                                                                                                                           <  1:1617895/53‑1 (MQ=255)
                                                                                                                              cTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGATCATcggcg                         >  1:1609217/1‑150 (MQ=255)
                                                                                                                               tGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGATCATcggcgg                        <  2:1813437/150‑1 (MQ=255)
                                                                                                                                       tCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCGGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGATCATCGGCGGCGGCCCTg                >  2:575627/1‑150 (MQ=255)
                                                                                                                                         catCGACCACGCGGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGa                                                                       <  2:1807207/91‑1 (MQ=255)
                                                                                                                                         catCGACCACGCGGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGa                                                                       >  1:1807207/3‑93 (MQ=255)
                                                                                                                                          ttCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGATCATCGGCGGCGGCCCTggcg             >  1:1177731/1‑150 (MQ=255)
                                                                                                                                          atCGACCACGCGGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGATCATCGGCGGCGGCCCTggcg             >  1:2361536/2‑150 (MQ=255)
                                                                                                                                          atCGACCACGCGGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGATCATCGGCGGCGGCCCTggcg             >  1:2124061/2‑150 (MQ=255)
                                                                                                                                             gACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGATCATCGGCGGCGGCCCTGGCGGCt          >  1:207396/1‑150 (MQ=255)
                                                                                                                                               ccaccGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACa                                                         >  2:1944591/1‑101 (MQ=255)
                                                                                                                                               ccaccGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACa                                                         <  1:1944591/101‑1 (MQ=255)
                                                                                                                                               ccaccGCCTGCTCGATAGCATGGATGCCGCCTTGTTCATCCAGGCCGTACGTGGCCTGCGCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGTATGCAACAGATTATCCAGACTACCCTGTTGATCATCGGCGGCGGCCCTGGCGGCTAt        <  1:863530/150‑1 (MQ=255)
                                                                                                                                                  ccGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGATCATCGGCGGCGGCCCTGGCGGCTATGTa     <  2:1672318/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
GCCCTGGGTGGCATTGTCAGCACGCCGGTGGTCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGATCATCGGCGGCGGCCCTGGCGGCTATGTAGCA  >  NC_002947/4995423‑4995721

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