Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,992,977 Δ1 bp 5.6% coding (886/1233 nt) bkdAA → branched‑chain alpha‑keto acid dehydrogenase complex subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,992,9770G. 5.6% 226.8 / 7.9 71coding (886/1233 nt)bkdAAbranched‑chain alpha‑keto acid dehydrogenase complex subunit alpha
Reads supporting (aligned to +/- strand):  ref base G (27/40);  new base . (2/2);  total (29/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.80e-01

TTTCCACCTTCCAGGCCATCGCCGGTGGCGAGTCGACCACCTTTGCCGGCCGTGGCGTGGGTTGCGGTATTGCTTCGCTGCGGGTTGACGGCAACGACTTCGTCGCCGTGTACGCTGCCTCGCGCTGGGCGGCCGAGCGCGCCCGCCGCGGCCTGGGCCCAAGCCTGATCGAGTGGGTCACCTACCGTGCCGGCCCGCACTCGACGTCGGACGACCCCTCCAAGTACCGCCCTGCCGACGACTGGAGCCACTTCCCGCTGGGTGACCCGATCGCCCGCCTGAAGCAGCACCTGATCAA  >  NC_002947/4992828‑4993125
                                                                                                                                                     |                                                                                                                                                    
tttCCACCTTCCAGGCCATCGCCGGTGGCGAGTCGACCACCTTTGCCGGCCGTGGCGTGGGTTGCGGTATTGCTTCGCTGCGGGTTGACGGCAACGACTTCGTCGCCGTGTACGCTGCCTCGCGCTGGGCGGCCGAGCGCGCCCGCcgcg                                                                                                                                                      >  2:1417890/1‑150 (MQ=255)
    cACCTTCCAGGCCATCGCCGGTGGCGAGTCGACCACCTTTGCCGGCCGTGGCGTGGGTTGCGGTATTGCTTCGCTGCGGGTTGACGGCAACGACTTCGTCGCCGTGTACGCTGCCTCGCGCTGGGCGGCCGAGCGCGCCCGCCGCGGCCt                                                                                                                                                  >  2:660515/1‑150 (MQ=255)
       cTTCCAGGCCATCGCCGGTGGCGAGTCGACCACCTTTGCCGGCCGTGGCGTGGGTTGCGGTATTGCTTCGCTGCGGGTTGACGGCAACGACTTCGTCGCCGTGTACGCTGCCTCGCTCTGGGCGGCCGAGCGCGCCCGCCGCGGCCTggg                                                                                                                                               <  1:1620267/150‑1 (MQ=255)
        ttCCAGGCCATCGCCGGTGGCGAGTCGACCACCTTTGCCGGCCGTGGCGTGGGTTGCGGTATTGCTTCGCTGCGGGTTGACGGCAACGACTTCGTCGCCGTGTACGCTGCCTCGCGCTGGGCGGCCGAGCGCGCCCGCCGCGGCCTGGGc                                                                                                                                              <  1:1587471/150‑1 (MQ=255)
          ccAGGCCTTCGCCGGTGGAGAGTCGATCACCTTTCCCGGCCGTGGCGTGGTTTGCGGTATTGCTTCGCTGCGGGTTGACGGCAACGACTTCGTCGCCGTGTACGCTGCCTCGCGCTGGGCGGCCGAGCGCGCCCGCCGCGGCCTGGGccc                                                                                                                                            <  2:1250339/150‑1 (MQ=255)
            aGGCCATCGCCGGTGGCGAGTCGACCACCTTTGCCGGCCGTGGCGTGGGTTGCGGTATTGCTTCGCTGCGGGTTGACGGCAACGACTTCGTCGCCGTGTACGCTGCCTCGCGCTGGGCGGCCGAGCGCGCCCGCCGCGGCCTGGGCCCaa                                                                                                                                          >  2:1268281/1‑150 (MQ=255)
              gCCATCGCCGGTGGCGAGTCGACCACCTTTGCCGGCCGTGGCGTGGGTTGCGGTATTGCTTCGCTGCGGGTTGACGGCAACGACTTCGTCGCCGTGTACGCTGCCTCGCGCTGGGCGGCCGAGCGCGCCCGCCGCGGCCTGGGCCCAAGc                                                                                                                                        <  2:1323023/150‑1 (MQ=255)
              gCCATCGCCGGTGGCGAGTCGACCACCTTTGCCGGCCGTGGCGTGGGTTGCGGTATTGCTTCGCTGCGGGTTGACGGCAACGACTTCGTCGCCGTGTACGCTGCCTCGCGCTGGGCGGCCGAGCGCGCCCGCCGCGGCCTGGGCCCAAGc                                                                                                                                        <  1:1268281/150‑1 (MQ=255)
              gCCATCGCCGGTGGCGAGTCGACCACCTTTGCCGGCCGTGGCGTGGGTTGCGGTATTGCTTCGCTGCGGGTTGACGGCAACGACTTCGTCGCCGTGTACGCTGCCTCGCGCTGGGCGGCCGAGCGCGCCCGCCGCGGCCTGGGCCCAAGc                                                                                                                                        >  2:724259/1‑150 (MQ=255)
              gCCATCGCCGGTGGCGAGTCGACCACCTTTGCCGGCCGTGGCGTGGGTTGCGGTATTGCTTCGCTGCGGGTTGACGGCAACGACTTCGTCGCCGTGTACGCTGCCTCGCGCTGGGCGGCCGAGCGCGCCCGCCGCGGCCTGGGCCCAAGc                                                                                                                                        <  2:1535779/150‑1 (MQ=255)
                  tCGCCGGTGGCGAGTCGACCACCTTTGCCGGCCGTGGCGTGGGTTGCGGTATTGCTTCGCTGCGGGTTGACGGCAACGACTTCGTCGCCGTGTACGCTGCCTCGCGCTGGGCGGCCGAGCGCGCCCGCCGCGGCCTGGGCTCAAGCCTGa                                                                                                                                    <  2:297384/150‑1 (MQ=255)
                   cGCCGGTGGCGAGTCGACCACCTTTGCCGGCCGTGGCGTGGGTTGCGGTATTGCTTCGCTGCGGGTTGACGGCAACGACTTCGTCGCCGTGTACGCTGCCTCGCGCTGGGCGGCCGAGCGCGCCCGCCGCGGCCTGGGCCCAAGCCTGAt                                                                                                                                   <  1:273240/150‑1 (MQ=255)
                    gCCGGTGGCGAGTCGACCACCTTTGCCGGCCGTGGCGTGGGTTGCGGTATTGCTTCGCTGCGGGTTGACGGCAACGACTTCGTCGCCGTGTACGCTGCCTCGCGCTGGGCGGCCGAGCGCGCCCGCCGCGGCCTGGGCCCAAGCCTGATc                                                                                                                                  <  2:1279237/150‑1 (MQ=255)
                    gCCGGTGGCGAGTCGACCACCTTTGCCGGCCGTGGCGTGGGTTGCGGTATTGCTTCGCTGCGGGTTGACGGCAACGACTTCGTCGCCGTGTACGCTGCCTCGCGCTGGGCGGCCGAGCGCGCCCGCCGCGGCCTGGGCCCAAGCCTGATc                                                                                                                                  <  2:551326/150‑1 (MQ=255)
                      cGGTGGCGAGTCGACCACCTTTGCCGGCCGTGGCGTGGGTTGCGGTATTGCTTCGCTGCGGGTTGACGGCAACGACTTCGTCGCCGTGTACGCTGCCTCGCGCTGGGCGGCCGAGCGCGCCCGCCGCGGCCTGGGCCCAAGCCTGATCGa                                                                                                                                >  2:1482997/1‑150 (MQ=255)
                        gTGGCGAGTCGACCACCTTTGCCGGCCGTGGCGTGGGTTGCGGTATTGCTTCGCTGCGGGTTGACGGCAACGACTTCGTCGCCGTGTACGCTGCCTCGCGCTGGGCGGCCGAGCGCGCCCGCCGCGGCCTGGGCCCAAGCCTGATCGAGt                                                                                                                              <  1:660515/150‑1 (MQ=255)
                        gTGGCGAGTCGACCACCTTTGCCGGCCGTGGCGTGGGTTGCGGTATTGCTTCGCTGCGGGTTGACGGCAACGACTTCGTCGCCGTGTACGCTGCCTCGCGCTGGGCGGCCGAGCGCGCCCGCCGCGGCCTGGGCCCAAGCCTGATCGAGt                                                                                                                              <  2:963812/150‑1 (MQ=255)
                         tGGCGAGTCGACCACCTTTGCCGGCCGTGGCGTGGGTTGCGGTATTGCTTCGGTGCGGGTTGACGGCAACGACTTCGGCGCCGTGTACGCTCCCACGCGCGGGGCGGCCGAGCGCGCCCGCCGCGGCGTGGGCCCAAGCCTGATCGGGTg                                                                                                                             <  2:1187155/150‑1 (MQ=255)
                          ggCGAGTCGACCACCTTTGCCGGCCGTGGCGTGGGTTGCGGTATTGCTTCGCTGCGGGTTGACGGCAACGACTTCGTCGCCGTGTACGCTGCCTCGCGCTGGGCGGCCGAGCGCGCCCGCCGCGGCCTGGGCCCAAGCCTGATCGAGTgg                                                                                                                            <  2:972856/150‑1 (MQ=255)
                              aGTCGACCACCTTTGCCGGCCGTGGCGTGGGTTGCGGTATTGCTTCGCTGCGGGTTGACGGCAACGACTTCGTCGCCGTGTACGCTGCCTCGCACTGGGCGGCCGAGCGCGCCCGCCGCGGCCTGGGCCCAAGCCTGATCGAGTGGGTCa                                                                                                                        <  2:1059926/150‑1 (MQ=255)
                                tCGACCACCTTTGCCGGCCGTGGCGTGGGTTGCGGTATTGCTTCGCTGCGGGTTGACTGCAACGACTTCGTCGCCGTTTACGCGGCCTCGCGCTGGGCGGCCGAGCGCGCCCGCCGCGGCCTGGGCCCAAGCCTGATCGAGTGGGTCAcc                                                                                                                      <  2:220928/150‑1 (MQ=255)
                                   accaccTTTGCCGGCCGTGGCGTGGGTTGCGGTATTGCTTCGCTGCGGGTTGACGGCAACGACTTCGTCGCCGTGTACGCTGCCTCGCGCTGGGCGGCCGAGCGACAATGCCGCGGACTGGGCCGAAGCCTGATCGAGTGGGTCACCTAc                                                                                                                   >  2:999127/1‑150 (MQ=255)
                                    ccaccTTTGCCGGCCGTGGCGTGGGTTGCGGTATTGCTTCGCTGCGGGTTGACGGCAACGACTTCGTCGCCGTGTACGCTGCCTCGCGCTGGGCGGCCGAGCGCGCCCGCCGCGGCCTGGGCCCAAGCCTGATCGAGTGGGTCACCTAcc                                                                                                                  <  2:938635/150‑1 (MQ=255)
                                    ccaccTTTGCCGGCCGTGGCGTGGGTTGCGGTATTGCTTCGCTGCGGGTTGACGGCAACGACTTCGTCGCCGTGTACGCTGCCTCGCGCTGGGCGGCCGAGCGCGCCCGCCGCGGCCTGGGCCCAAGCCTGATCGAGTGGGTCACCTAcc                                                                                                                  <  1:999127/150‑1 (MQ=255)
                                       ccTTTGCCGGCCGTGGCGTGGGTTGCGGTATTGCTTCGCTGCGGGTTGACGGCAACGACTTCGTCGCCGTGTACGCTGCCTCGCGCTGGGCGGCCGAGCGCGCCCGCCGCGGCCTGGGCCCAAGCCTGATCGAGTGGGTCACCTACCGTg                                                                                                               >  1:264524/1‑150 (MQ=255)
                                        cTTTGCCGGCCGTGGCGTGGGTTGCGGTATTGCTTCGCTGCGGGTTGACAGCAACGACTTCGTCGCCGTGTACGCTGCCTCGCGCTGGGCGGCCGAGCGCGCCCGCCGCGGCCTGGGCCCAAGCCTGATCGAGTGGGTCACCTACCGTGc                                                                                                              >  2:1423867/1‑150 (MQ=255)
                                        cTTTGCCGGCCGTGGCGTGGGGGGCGGTATTGCTTCGCTGCGGGTTGACGGCAACGACTTCGTCGCCGTGTACGCTGCCTCGCGCTGGGCGGCCGAGCGCGCCCGCCGCGGCCTGGGCCCAAGCCTGATCGAGTGGGTCACCTACCGTGc                                                                                                              >  2:68304/1‑150 (MQ=255)
                                              cggccgTGGCGTGGGTTGCGGTATTGCTTCGCTGCGGGTTGACGGCAACGACTTCGTCGCCGTGTACGCTGCCTCGCGCTGGGCGGCCGAGCGCGCCCGCCGCGGCCTGGGCCCAAGCCTGATCGAGTGGGTCACCTACCGTGCCGGccc                                                                                                        <  1:1482997/150‑1 (MQ=255)
                                              cggccgTGGCGTGGGTTGCGGTATTGCTTCGCTGCGGGTTGACGGCAACGACTTCGTCGCCGTGTACGCTGCCTCGCGCTGGGCGGCCGAGCGCGCCCGCCGCGGCCTGGGCCCAAGCCTGATCGAGTGGGTCACCTACCGTGCCGGccc                                                                                                        <  2:493825/150‑1 (MQ=255)
                                                gccgTGGCGTGGGTTGCGGTATTGCTTCGCTGCGGGTTGACGGCAACGACTTCGTCGCCGTGTACGCTGCCTCGCGCTGGGCGGCCGAGCGCGCCCGCCGCGGCCAGGGCCCAAGCCTGATCGAGTGGGTCACCTACCGTGCCGGCCCGc                                                                                                      >  2:1184763/1‑150 (MQ=255)
                                                      gcgtggGTTGCGGTATTGCTTCGCTGCGGGTTGACGGCAACGACTTCGTCGCCGTGTACGCTGCCTCGCGCTGGGCGGCCGAGCGCGCCCGCCGCGGCCTGGGCCCAAGCCTGATCGAGTGGGTCACCTACCGTGCCGGCCCGCACTCGa                                                                                                >  2:865934/1‑150 (MQ=255)
                                                                cGGTATTGCTTCGCTGCGGGTTGACGGCAACGACTTCGTCGCCGTGTACGCTGCCTCGCGCTGGGCGGCCGAGCGCGCCCGCCGCGGCCTGGGCCCAAGCCTGATCGAGTGGGTCACCTACCGTGCCGGCCCGCACTCGACGTCGgacga                                                                                      <  1:865934/150‑1 (MQ=255)
                                                                    aTTGCTTCGCTGCGGGTTGACGGCAACGACTTCGTCGCCGTGTACGCTGCCTCGCGCTGGGCGGCCGAGCGCGCCCGCCGCGGCCTGGGCCCAAGCCTGATCGAGTGGGTCACCTACCGTGCCGGCCCGCACTCGACGTCGGACGAcccc                                                                                  <  2:529231/150‑1 (MQ=255)
                                                                      tGCTTCGCTGCGGGTTGACGGCAACGACTTCGTCGCCGTGTACGCTGCCTCGCGCTGGGCGGCCGAGCGCGCCCGCCGCGGCCTGGGCCCAAGCCTGATCGAGTGGGTCACCTACCGTGCCGGCCCGCACTCGACGTCGGACGACCCCTc                                                                                <  1:1417890/150‑1 (MQ=255)
                                                                       gCTTCGCTGCGGGTTGACGGCAACGACTTCGTCGCCGTGTACGCTGCCTCGCGCTGGGCGGCCGAGCGCGCCCGCCGCGGCCTGGGCCCAAGCCTGATCGAGTGGGTCACCTACCGTGCCGGCCCGCACTCGACGTCGGACGACCCCTcc                                                                               <  1:1254867/150‑1 (MQ=255)
                                                                        cTTCGCTGCGGGTTGACGGCAACGACTTCGTCGCCGTGTACGCTGCCTCGCGCTGGGCGGCCGAGCGCGCCCGCCGCGGCCTGGGCCCAAGCCTGATCGAGTGGGTCACCTACCGTGCCGGccc                                                                                                        >  1:86250/1‑124 (MQ=255)
                                                                        cTTCGCTGCGGGTTGACGGCAACGACTTCGTCGCCGTGTACGCTGCCTCGCGCTGGGCGGCCGAGCGCGCCCGCCGCGGCCTGGGCCCAAGCCTGATCGAGTGGGTCACCTACCGTGCCGGccc                                                                                                        <  2:86250/124‑1 (MQ=255)
                                                                         ttCGCTGCGGGTTGACGGCAACGACTTCGTCTCCGTGTACGCTGCCTCGCGCTGGGCGGCCGAGCGCGCCCGCCGCGGCCTGGGCCCAAGCCTGATCGAGTGGGTCACCTACCGTGCCGGCCCGCACTCGACGTCGGACGACCCCTCCaa                                                                             <  2:746028/150‑1 (MQ=255)
                                                                         ttCGCTGCGGGTTGACGGCAACGACTTCGTCGCCGTGTACGCTGCCTCGCGCTGGGCGGCCGAGCGCGCCCGCCGC‑GCCTGGGCCCAAGCCTGATCGAGTGGGTCACCTACCGTGCCGGCCCGCACTCGACGTCGGACGACCCCTCCAAg                                                                            >  1:425521/1‑150 (MQ=255)
                                                                         ttCGCTGCGGGTTGACGGCAACGACTTCGTCGCCGTGTACGCTGCCTCGCGCTGGGCGGCCGAGCGCGCCCGCCGC‑GCCTGGGCCCAAGCCTGATCGAGTGGGTCACCTACCGTGCCGGCCCGCACTCGACGTCGGACGACCCCTCCAAg                                                                            >  1:425326/1‑150 (MQ=255)
                                                                              tGCGGGTTGACGGCAACGACTTCGTCGCCGTGTACGCTGCCTCGCGCTGGGCGGCCGAGCGCGCCCGCCGCGGCCTGGGCCCAAGCCTGAt                                                                                                                                   >  2:1310833/1‑91 (MQ=255)
                                                                              tGCGGGTTGACGGCAACGACTTCGTCGCCGTGTACGCTGCCTCGCGCTGGGCGGCCGAGCGCGCCCGCCGCGGCCTGGGCCCAAGCCTGAt                                                                                                                                   <  1:1310833/91‑1 (MQ=255)
                                                                                 gggTTGACGGCAACGACTTCGTCGCCGTGTACGCTGCCTCGCGCTGGGCGGCCGAGCGCGCCCGCCGCGGCCTGGGCCCAAGCCTGATCGAGTGGGTCACCTACCGTGCCGGCCCGCACTCGACGTCGGACGACCCCTCCAAGTACCGcc                                                                     <  2:264524/150‑1 (MQ=255)
                                                                                    ttGACGGCAACGACTTCGTCGCCGTGTACGCTGCCTCGCGCTGGGCGGCCGAGCGCGCCCGCCGCGGCCTGGGCCCAAGCCTGATCGAGTGGGTCACCTACCGTGCCGGCCCGCACTCGACGTCGGACGACCCCTCCAAGTACCGCCCTg                                                                  >  2:182004/1‑150 (MQ=255)
                                                                                         ggCAACGACTTCGTCGCCGTGTACGCTGCCTCGCGCTGGGCGGCCGAGCGCGCCCGCCGCGGCCTGGGCCCAAGCCTGATCGAGTGGGTCAc                                                                                                                       >  1:1268027/1‑92 (MQ=255)
                                                                                         ggCAACGACTTCGTCGCCGTGTACGCTGCCTCGCGCTGGGCGGCCGAGCGCGCCCGCCGCGGCCTGGGCCCAAGCCTGATCGAGTGGGTCAc                                                                                                                       <  2:1268027/92‑1 (MQ=255)
                                                                                          gCAACGACTTCGTCGCCGTGTACGCTGCCTCGCGCTGGGCGGCCGAGCGCGCCCGCCGCGGCCTGGGCCCAAGCCTGATc                                                                                                                                  <  1:1378964/80‑1 (MQ=255)
                                                                                          gCAACGACTTCGTCGCCGTGTACGCTGCCTCGCGCTGGGCGGCCGAGCGCGCCCGCCGCGGCCTGGGCCCAAGCCTGATc                                                                                                                                  >  2:1378964/1‑80 (MQ=255)
                                                                                           cAACGACTTCGTCGCCGTGTACGCTGCCTCGCGCTGGGCGGCCGAGCGCGCCCGCCGC‑GCCTGGGCCCAAGCCTGATCGAGTGGGTCACCTACCGTGCCGGCCCGCACTCGACGTCGGACGACCCCTCCAAGTACCGCCCTGCCgacgac                                                          <  2:425326/150‑1 (MQ=255)
                                                                                           cAACGACTTCGTCGCCGTGTACGCTGCCTCGCGCTGGGCGGCCGAGCGCGCCCGCCGC‑GCCTGGGCCCAAGCCTGATCGAGTGGGTCACCTACCGTGCCGGCCCGCACTCGACGTCGGACGACCCCTCCAAGTACCGCCCTGCCgacgac                                                          <  2:425521/150‑1 (MQ=255)
                                                                                            aaCGACTTCGTCGCCGTGTACGCTGCCTCGCGCTGGGCGGCCGAGCGCGCCCGCCGCGGCCTGGGCCCAAGCCTGATCGAGTGGGTCACCTACCGTGCCGGCCCGCACTCGACGTCGGACGACCCCTCCAAGTACCGCCCTGCCgacgac                                                          >  1:672576/1‑150 (MQ=255)
                                                                                             aCGACTTCGTCGCCGTGTACGCTGCCTCGCGCTGGGCGGCCGAGCGCGCCCGCCGCGGCCTGGGCCCAAGCCTGATCGAGTGGGTCACCTACCGTGCCGGCCCGCACTCGACGTCGGACGACCCCTCCAAGTACCGCCCTGc                                                                 >  1:660512/1‑142 (MQ=255)
                                                                                             aCGACTTCGTCGCCGTGTACGCTGCCTCGCGCTGGGCGGCCGAGCGCGCCCGCCGCGGCCTGGGCCCAAGCCTGATCGAGTGGGTCACCTACCGTGCCGGCCCGCACTCGACGTCGGACGACCCCTCCAAGTACCGCCCTGc                                                                 <  2:660512/142‑1 (MQ=255)
                                                                                             aCGACTTCGTCGCCGTGTACGCTGCCTCGCGCTGGGCGGCCGAGCGCGCCCGCCGCGGCCTGGGCCCAAGCCTGATCGAGTGGGTCACCTACCGTGCCGGCCCGCACTCGACGTCGGACGACCCCTCCAAGTACCGCCCTGCCGACGACt                                                         <  2:286388/150‑1 (MQ=255)
                                                                                             aCGACTTCGTCGCCGAGTACGCTGCCTCGCGCTGCGCCGCTGAGCGCGCCACCCGCGGCGTGGGCCCAAGCCTGATCGACTGGGTCACCTACCGTGCCGCCCCGCACTCGACGTCCGACGTCCCCTCCAAGTACCGCCCTGCCGACGACt                                                         <  2:107936/150‑1 (MQ=255)
                                                                                                aCTTCGTCGCCGTGTACGCTGCCTCGCGCTGGGCGGCCGAGCGCGCCCGCCGCGGCCTGGGCCCAAGCCTGATCGAGTGGGTCACCTACCGTGCCGGCCCGCACTCGACGTCGGACGACCCCTCCAAGTACCGCCCTGCCGACGACTgtt                                                      <  1:724259/150‑3 (MQ=255)
                                                                                                 cTTCGTCGCCGTGTACGCTGCCTCGCGCTGGGCGGCCGAGCGCGCCCGCCGCGGCCTGGGCCCAAGCCTGATCGAGTGGGTCACCTACCGTGCCGGCCCGCACTCGACGTCGGACGACCCCTCCAAGTACCGCCCTGACGACGACggggg                                                     >  2:1326919/1‑148 (MQ=255)
                                                                                                  ttCGTCGCCGTGTACGCTGCCTCGCGCTGGGCGGCCGAGCGCGCCCGCCGCGGCCTGGGCCCAAGCCTGATCGAGTGGGTCACCTACCGTGCCGGCCCGCACTCGACGTCGGACGACCCCTCCAAGTACCGCCCTGCCGACGACTGGAGc                                                    >  2:319250/1‑150 (MQ=255)
                                                                                                           gtgtACGCTGCCTCGCGCTGGGCGGCCGAGCGCGCCCGCCGCGGCCTGGGCCCAAGCCTGATCGAGTGGGTCACCTACCGTGCCGGCCCGCACTCGACGTCGGACGACCCCTCCAAGTACCGCCCTGCCGACGACTGGAGCCACTTCCCg                                           >  2:82529/1‑150 (MQ=255)
                                                                                                               aCGCTGCCTCGCGCTGGGCGGCCGAGCGCGCCCGCCGCGGCCTGGGCCCAAGCCTGATCGAGTGGGTCACCTACCGTGCCGGCCCGCACTCGACGTCGGACGACCCCTCCAAGTACCGCCCTGCCGACGACTGGAGCCACTTCCCGCTgg                                       <  1:1326919/150‑1 (MQ=255)
                                                                                                                  cTGCCTCGCGCTGGGCGGCCGAGCGCGCCCGCCGCGGCCTGGGCCCAAGCCTGATCGAGTGGGTCACCTACCGTGCCGGCCCGCACTCGACGTCGGCCGACCCCTCCAAGTACCGCCCTGCCGACGACTGGAGCCCCTTCCCGCTGGGTg                                    >  2:585879/1‑150 (MQ=255)
                                                                                                                     ccTCGCGCTGGGCGGCCGAGCGCGCCCGCCGCGGCCTGGGCCCAAGCCTGATCGAGTGGGTCACCTACCGTGCCAGCCCGCACTCGACGTCGGACGACCCCTCCAAGTACCGCCCTGCCGACGACTGGAGCCACTTCCCGCTGGGTGAcc                                 >  1:894265/1‑150 (MQ=255)
                                                                                                                        cgcgCTGGGCGGCCGAGCGCGCCCGCCGCGGCCTGGGCCCAAGCCTGATCGAGTGGGTCACCTACCGTGCCGGCCCGCACTCGACGTCGGACGACCCCTCCAAGTACCGCCCTGCCGACGACTGGAGCCACTTc                                              >  1:1545695/1‑134 (MQ=255)
                                                                                                                        cgcgCTGGGCGGCCGAGCGCGCCCGCCGCGGCCTGGGCCCAAGCCTGATCGAGTGGGTCACCTACCGTGCCGGCCCGCACTCGACGTCGGACGACCCCTCCAAGTACCGCCCTGCCGACGACTGGAGCCACTTc                                              <  2:1545695/134‑1 (MQ=255)
                                                                                                                        cgcgCTGGGCGGCCGAGCGCGCCCGCCGCGGCCTGGGCCCAAGCCTGATCGAGTGGGTCACCTACCGTGCCAGCCCGCACTCGACGTCGGACGACCCCTCCAAGTACCGCCCTGCCGACGACTGGAGCCACTTCCCGCTGGGTGACCCGa                              <  2:894265/150‑1 (MQ=255)
                                                                                                                          cgcTGGGCGGCCGAGCGCGCCCGCCGCGGCCTGGGCCCAAGCCTGATCGAGTGGGTCACCTACCGTGCCGGCCCGCACTCGACGTCGGACGACCCCTCCAAGTACCGCCCTGCCGACGACTGGAGCCACTTCCCGCTGGGTGACCCGATc                            <  1:68304/150‑1 (MQ=255)
                                                                                                                           gcTGGGCGGCCGAGCGCGCCCGCCGCGGCCTGGGCCCAAGCCTGATCGAGTGGGTCACCTACCGTGCCGGCCCGCACTCGACGTCGGACGACCCCTCCAAGTACCGCCCTGCCGACGACTGGAGCCACTTCCCGCTg                                        <  2:576740/137‑1 (MQ=255)
                                                                                                                           gcTGGGCGGCCGAGCGCGCCCGCCGCGGCCTGGGCCCAAGCCTGATCGAGTGGGTCACCTACCGTGCCGGCCCGCACTCGACGTCGGACGACCCCTCCAAGTACCGCCCTGCCGACGACTGGAGCCACTTCCCGCTg                                        >  1:576740/1‑137 (MQ=255)
                                                                                                                           gcTGGGCGGCCGAGCGCGCCCGCCGCGGCCTGGGCCCAAGCCTGATCGAGTGGGTCACCTACCGTGCCGGCCCGCACTCGACGTCGGACGACCCCTCCAAGTACCGCCCTGCCGACGACTGGAGCCACTTCCCGCTGGGTGACCCGATcg                           >  2:462442/1‑150 (MQ=255)
                                                                                                                                     cGAGCGCGCCCGCCGCGGCCTGGGCCCAAGCCTGATCGAGTGGGTCACCTACCGTGCCGGCCCGCACTCGACGTCGGACGACCCCTCCAAGTACCGCCCTGCCGACGACTGGAGCCACTTCCCGCTGGGTGACCCGATCGCCCGCCTGaa                 >  2:1638407/1‑150 (MQ=255)
                                                                                                                                     cGAGCGCGCCCGCCGCGGCCTGGGCCCAAGCCTGATCGAGTGGGTCACCTACCGTGCCGGCCCGCACTCGACGTCGGACGACCCCTCCAAGTACCGCCCTGCCGACGACTGGAGCCACTTCCCGCTGGGTGACCCGATCGCCCGCCTGaa                 <  2:451140/150‑1 (MQ=255)
                                                                                                                                              ccgccgCGGCCTGGGCCCAAGCCTGAGCGAGTGGGTCACCTACCGTGCCGGCCCGCACTCGACGTCGGACGACCCCTCCAAGTACCGCCCTGCCGACGACTGGAGCCACTTCCCGCTGGGTGACCCGATCGCCCGCCTGAAGCAGCACCt        >  1:221051/1‑150 (MQ=255)
                                                                                                                                                    cgGCCTGGGCCCAAGCCTGATCGAGTGGGTCACCTACCGTGCCGGCCCGCACTCGACGTCGGACGACCCCTCCAAGTACCGCCCTGCCGACGACTGGAGCCACTTCCCGCTGGGTGACCCGATCGCCCGCCTGAAGCAGCACCTGATCaa  <  2:221051/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TTTCCACCTTCCAGGCCATCGCCGGTGGCGAGTCGACCACCTTTGCCGGCCGTGGCGTGGGTTGCGGTATTGCTTCGCTGCGGGTTGACGGCAACGACTTCGTCGCCGTGTACGCTGCCTCGCGCTGGGCGGCCGAGCGCGCCCGCCGCGGCCTGGGCCCAAGCCTGATCGAGTGGGTCACCTACCGTGCCGGCCCGCACTCGACGTCGGACGACCCCTCCAAGTACCGCCCTGCCGACGACTGGAGCCACTTCCCGCTGGGTGACCCGATCGCCCGCCTGAAGCAGCACCTGATCAA  >  NC_002947/4992828‑4993125

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: ‑