Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,649,537 T→C 9.2% H522R (CAT→CGT)  PP_1446 ← TonB‑dependent receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,649,5370TC 9.2% 186.7 / 7.5 65H522R (CAT→CGT) PP_1446TonB‑dependent receptor
Reads supporting (aligned to +/- strand):  ref base T (30/29);  new base C (3/3);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.31e-01

TTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATGTAGTTGCGAAAGTGGCTGTAGAACAACCCAACACTGCCCTTGTGGGTGCCGTTGTCGAAGCGCAGCGCCAGG  >  NC_002947/1649390‑1649683
                                                                                                                                                   |                                                                                                                                                  
ttCTTGGCTCGGGTATAGTCACCGTACAACTCCCGCGCGAAGCCAACGCAACGGTTCTCGGCCAACTGCCAGCGGTGCTGTGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTCGTATTCCGGGAAGCTCtcatgatgatga                                                                                                                                                  <  2:2077123/150‑1 (MQ=255)
        tCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCatgatcatgatcatgatc                                                                                                                                          >  1:2174311/1‑150 (MQ=255)
        tCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCatgatcatgatcatgatc                                                                                                                                          >  2:1356182/1‑150 (MQ=255)
        tCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCatgatcatgatcatgatc                                                                                                                                          <  2:1433365/150‑1 (MQ=255)
            gTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCtcatgatcatgatcatgatcatga                                                                                                                                      <  1:1308794/150‑1 (MQ=255)
            gTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGa                                                                                                                                      >  1:984097/1‑150 (MQ=255)
                    aCCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGAtcctcctc                                                                                                                              >  1:623396/1‑150 (MQ=255)
                    aCCGGAAAACTCCAGCGCGAAGCTGCCGTAACGGTTGGCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCTCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGAtcctcctc                                                                                                                              <  1:991464/150‑1 (MQ=255)
                      cGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCACGATCatgatcatgatcatgagg                                                                                                                              <  1:1609135/148‑7 (MQ=255)
                      cGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCACGATCatgatcatgatcatgagg                                                                                                                              >  2:1609135/1‑142 (MQ=255)
                       ggACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCctcctcctc                                                                                                                           >  2:1378018/1‑150 (MQ=255)
                       ggACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCCCACGCCTTGGTACTGGTATTCCGGGAAGCCCtcatgatcatgatc                                                                                                                                                <  1:1238131/129‑1 (MQ=255)
                       ggACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCCCACGCCTTGGTACTGGTATTCCGGGAAGCCCtcatgatcatgatc                                                                                                                                                >  2:1238131/1‑129 (MQ=255)
                          atcttCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTTGGCCTGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATg                                                                                                                        <  1:1870374/146‑1 (MQ=255)
                              tcgaGCGCGAAGCGGCCGTAACGGGTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCa                                                                                                                    >  2:1241564/4‑150 (MQ=255)
                              tCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCa                                                                                                                    >  1:1189445/1‑150 (MQ=255)
                                    gcgAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCg                                                                                                              >  1:1080388/1‑150 (MQ=255)
                                         gCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTc                                                                                                         >  2:2153714/1‑150 (MQ=255)
                                            gCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTcgcg                                                                                                      <  1:637671/150‑1 (MQ=255)
                                                tAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCACGATCatgatcatgatcatga                                                                                                                                >  2:2110754/1‑116 (MQ=255)
                                                tAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCACGATCatgatcatgatcatga                                                                                                                                <  1:2110754/120‑5 (MQ=255)
                                                  aCGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTAATCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGtt                                                                                                >  1:2098672/1‑150 (MQ=255)
                                                          cGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGc                                                                                        <  2:2098672/150‑1 (MQ=255)
                                                          cGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCACGATCatgatcatgatcatga                                                                                                                                >  1:2354638/1‑106 (MQ=255)
                                                          cGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCACGATCatgatcatgatcatga                                                                                                                                <  2:2354638/110‑5 (MQ=255)
                                                             ccAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCaa                                                                                     >  1:1285397/1‑150 (MQ=255)
                                                               aaCTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCtcatgatcat                                                                                                                                                    >  2:1004408/1‑85 (MQ=255)
                                                               aaCTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCtcatgatcat                                                                                                                                                    <  1:1004408/85‑1 (MQ=255)
                                                                aCTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCa                                                                                  <  1:363526/150‑1 (MQ=255)
                                                                aCTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCa                                                                                  >  2:1515159/1‑150 (MQ=255)
                                                                  tGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAgg                                                                                >  2:2103756/1‑150 (MQ=255)
                                                                    ccAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGcc                                                                              <  2:1285397/150‑1 (MQ=255)
                                                                    ccAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGcc                                                                              >  1:97144/1‑150 (MQ=255)
                                                                           tCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCtcatgatcatga                                                                                                                                                  <  2:1349292/75‑1 (MQ=255)
                                                                           tCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCtcatgatcatga                                                                                                                                                  >  1:1349292/1‑75 (MQ=255)
                                                                             cTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGACCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATGTAGtt                                                                     >  2:823851/1‑150 (MQ=255)
                                                                                   cTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGAtcc                                                                                                                                   >  2:1934228/1‑82 (MQ=255)
                                                                                   cTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGAtcc                                                                                                                                   <  1:1934228/82‑1 (MQ=255)
                                                                                     cGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGc                                                                                        <  1:1528/123‑1 (MQ=255)
                                                                                     cGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGc                                                                                        >  2:1528/1‑123 (MQ=255)
                                                                                     cGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATGTAGTTGCGAAAGt                                                             <  2:2174311/150‑1 (MQ=255)
                                                                                       aTGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATGTAGTTGCGAAAGTgg                                                           <  1:1515159/150‑1 (MQ=255)
                                                                                        tGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATGTAGTTGCGAAAGTGGc                                                          <  2:2106354/150‑1 (MQ=255)
                                                                                         gCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATGTAGTTGCGAAAGTGGCt                                                         <  1:1241564/150‑1 (MQ=255)
                                                                                           cGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATGTAGTTGc                                                                   >  1:566780/1‑138 (MQ=255)
                                                                                           cGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATGTAGTTGc                                                                   <  2:566780/138‑1 (MQ=255)
                                                                                            gTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGACCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATGTAGTTGCGAAAGTTGCTGTa                                                      <  1:823851/150‑1 (MQ=255)
                                                                                            gTAGAAGCGCGCCCGCACGACTTGGTAATGGTATTCAGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGATGCCAAACAGGACGATGTAGTTGCGAAAGTGGCTGTa                                                      >  2:142744/1‑150 (MQ=255)
                                                                                             tAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGAtcctcctcc                                                                                                                             <  2:160430/78‑1 (MQ=255)
                                                                                             tAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGAtcctcctcc                                                                                                                             >  1:160430/1‑78 (MQ=255)
                                                                                             tAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATGTAGTTGCGAAAGTGGCTGTAg                                                     <  2:793003/150‑1 (MQ=255)
                                                                                                    gcgcCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATGTAGTTGCGAAAGTGGCTGTAGAACAAcc                                              <  1:1356182/150‑1 (MQ=255)
                                                                                                       cccGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATGTAGTTGCGAAAGTGGCTGTAGAACAACCCaa                                           >  2:180959/1‑150 (MQ=255)
                                                                                                        ccGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATGTAGTTGCGAAAGTGGCTGTAGAACAACCCAac                                          >  2:400722/1‑150 (MQ=255)
                                                                                                        ccGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATGTAGTTGCGAAAGTGGCTGTAGAACAACCCAac                                          >  1:1683301/1‑150 (MQ=255)
                                                                                                         cGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATGTAGTTGCGAAAGTGGCTGTAGAACAACCCAaca                                         <  2:623396/150‑1 (MQ=255)
                                                                                                          gCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCa                                                                                  >  2:489890/1‑108 (MQ=255)
                                                                                                          gCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCa                                                                                  <  1:489890/108‑1 (MQ=255)
                                                                                                            aCGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATGTAGTTGCGAAAGTGGCTGTAGAACAAACCAACACTg                                      <  2:1189445/150‑1 (MQ=255)
                                                                                                             cGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATGTAGTTGCGAAAGTGGCTGTAGAACAACCCAACACTGc                                     >  2:146712/1‑150 (MQ=255)
                                                                                                              gCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATGTAGTTGCGAAAGTGGCTGTAGAACAACCCAACACTGcc                                    >  2:1974310/1‑150 (MQ=255)
                                                                                                                     tACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATGTAGTTGCGAAAGTGGCTGTAGAACAACCCAACACTGCCCTTGTgg                             <  1:180959/150‑1 (MQ=255)
                                                                                                                         ggTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATAAGGCCGATGTAGTTGCGAAAGTGGCTGTAGAACAACCCAACACTGCCCTTGTGGGTGc                         <  1:2103756/150‑1 (MQ=255)
                                                                                                                             ttCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCaa                                                                                     <  2:1957132/86‑1 (MQ=255)
                                                                                                                             ttCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATGTAGTTGCGAAAGTGGCTGTAGAACAACCCAACACTGCCCTTGTGGGTGCCGtt                     >  1:577152/1‑150 (MQ=255)
                                                                                                                             ttCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCAGGTGCCaa                                                                                     >  1:1957132/1‑86 (MQ=255)
                                                                                                                               ccGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATGTAGTTGCGAAAGTGGCTGTAGAACAACCCAACACTGCCCTTGTGGGTGCCGTTGt                   <  2:413193/150‑1 (MQ=255)
                                                                                                                               ccGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATGTAGCTGCGAAAGTGGCTGTAGAACAACCCAACACTGCCCTTGTGGGTGCCGTTGt                   <  1:1378018/150‑1 (MQ=255)
                                                                                                                                cGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATGTAGTTGCGAAAGTGGCTGTAGAACAACCCAACACTGCCCTTGTGGGTGCCGTTGTc                  <  1:146712/150‑1 (MQ=255)
                                                                                                                                cGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATGTAGTTGCGAAAGTGGCTGTAGAACAACCCAACACTGCCCTTGTGGGTGCCGTTGTc                  >  2:730414/1‑150 (MQ=255)
                                                                                                                                     aGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATGTAGTTGCGAAAGTGGCTGTAGAACAACCCAACACTGCCCTTGTGGGTGCCGTTGTCGAagc             <  2:1080388/150‑1 (MQ=255)
                                                                                                                                       cccTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATGTAGTTGCGAAAGTGGCTGTAGAACAACCCAACACTGCCCTTGTGGGTGCCGTTGTCGAagcgc           <  1:730414/150‑1 (MQ=255)
                                                                                                                                       cccTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATGTAGTTGCGAAAGTGGCTGTAGAACAACCCAAAACTCGCCTTGTGGGTGCCGTTGTCGAagcgc           >  1:403/1‑150 (MQ=255)
                                                                                                                                        gatcatgatcatgatcatgaTCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATGTAGTTGCGaaa                                                               <  2:259881/95‑1 (MQ=255)
                                                                                                                                        gatcatgatcatgatcatgaTCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATGTAGTTGCGaaa                                                               >  1:259881/3‑97 (MQ=255)
                                                                                                                                        gatcatgatcatgatcatgaTCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATGTAGTTGCGAAAGTGGCTGTAGAACAACCCAACACTGCCCTtgt                               <  2:991725/127‑1 (MQ=255)
                                                                                                                                        gatcatgatcatgatcatgaTCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATGTAGTTGCGAAAGTGGCTGTAGAACAACCCAACACTGCCCTtgt                               >  1:991725/3‑129 (MQ=255)
                                                                                                                                        ccTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATGTAGTTGCGAAAGTGGCTGTAGAACAACCCAACACTGCCCTTGTGGGTGCCGTTGTCGAagcgca          >  2:2233308/1‑150 (MQ=255)
                                                                                                                                         cTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATGTAGtt                                                                     >  2:2175171/1‑90 (MQ=255)
                                                                                                                                         cTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATGTAGtt                                                                     <  1:2175171/90‑1 (MQ=255)
                                                                                                                                          tcatgatcatgatcatgaTCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAAt                                                                                    <  1:1061904/74‑1 (MQ=255)
                                                                                                                                          tcatgatcatgatcatgaTCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAAt                                                                                    >  2:1061904/1‑74 (MQ=255)
                                                                                                                                          tcatgatcatgatcatgaTCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGAt                                                                           >  2:1529590/1‑83 (MQ=255)
                                                                                                                                          tcatgatcatgatcatgaTCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGAt                                                                           <  1:1529590/83‑1 (MQ=255)
                                                                                                                                          tcatgatcatgatcatgaTCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATGTAGTTGCGAAAGTGGCTGTAGAACAACCCAACACTGCCCTTGTGGGTGCCGTTGTCGAagcgcagc        >  2:1037181/1‑150 (MQ=255)
                                                                                                                                               atcatgatcatgatcGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATGTAGTTGCGAAAGTGGCTGTAGAACAACCCAACACTGCCCTTGTGGGTGCCGTTGTCGAAGCGCAGCGCCAg   >  2:934782/1‑150 (MQ=255)
                                                                                                                                                tcatgatcatgatcGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATGTAGTTGCGAAAGTGGCTGTAGa                                                    >  2:2026681/1‑100 (MQ=255)
                                                                                                                                                tcatgatcatgatcGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATGTAGTTGCGAAAGTGGCTGTAGa                                                    <  1:2026681/100‑1 (MQ=255)
                                                                                                                                                tcatgatcatgatcGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATGTAGTTGCGAAAGTGGCTGTAGAACAACCCAACACTGCCCTTGTGGGTGCCGTTGTCGAAGCGCAGCGCCAgg  >  2:2084305/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TTCTTGGCTCGGGTATAGTCACCGGACAACTCCAGCGCGAAGCTGCCGTAACGGTTCTCGGCCAACTGCCAGCGGTCCTGGGCCTCGATGCCGTAGAAGCGCGCCCGCACGCCTTGGTACTGGTATTCCGGGAAGCCCTCATGATCATGATCATGATCGTGATCCTCCTCCTCATGGTCATGCTCGCCTTCGCGGGTGTTGCCGGTGCCAATCAGGCCGATGTAGTTGCGAAAGTGGCTGTAGAACAACCCAACACTGCCCTTGTGGGTGCCGTTGTCGAAGCGCAGCGCCAGG  >  NC_002947/1649390‑1649683

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: