Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,918,299 T→C 6.1% intergenic (‑33/+85) PP_1718 ← / ← prc EAL domain‑containing diguanylate phosphodiesterase/tail‑specific protease Prc

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,918,2990TC 6.1% 410.5 / 10.5 131intergenic (‑33/+85)PP_1718/prcEAL domain‑containing diguanylate phosphodiesterase/tail‑specific protease Prc
Reads supporting (aligned to +/- strand):  ref base T (58/65);  new base C (4/4);  total (62/69)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.23e-01

Only 100 of 131 total aligned reads displayed.

CCTCATGGCCGAATACCCTATGTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATCAGGATGCGACCGGTTTCCGCG  >  NC_002947/1918152‑1918447
                                                                                                                                                   |                                                                                                                                                    
ccTCATGGCCGAATACCCTATGTTCACTCAGGCACCTGATTGGCTGGATCAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAATGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGc                                                                                                                                                    <  2:398814/150‑1 (MQ=255)
ccTCATGGCCGAATACCCTATGTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGc                                                                                                                                                    <  2:592201/150‑1 (MQ=255)
  tCATGGCCGAATACCCTATGTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCAt                                                                                                                                                  >  2:1467293/1‑150 (MQ=255)
   cATGGCCGAATACCCTATGTTCACTCAGGCACATGATCGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATg                                                                                                                                                 >  1:3090600/1‑150 (MQ=255)
      ggCCGAATACCCTATGTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGccc                                                                                                                                              <  2:2677343/150‑1 (MQ=255)
        ccGAATACCCTATGTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATGCAGATGCCCgg                                                                                                                                            <  2:328539/150‑1 (MQ=255)
        ccGAATACCCTATGTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCgg                                                                                                                                            <  2:3087678/150‑1 (MQ=255)
          gAATACCCTATGTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATGCAGATGCCCGGcc                                                                                                                                          <  1:3352201/150‑1 (MQ=255)
          gAATACCCTATGTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGcc                                                                                                                                          <  1:406673/150‑1 (MQ=255)
            aTACCCTATGTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCt                                                                                                                                        <  2:1336445/150‑1 (MQ=255)
             tACCCTATGTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATGCAGATGCCCGGCCCtt                                                                                                                                       <  2:4205178/150‑2 (MQ=255)
                ccTATGTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCata                                                                                                                                    >  1:543816/1‑150 (MQ=255)
                ccTATGTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCata                                                                                                                                    >  1:2361530/1‑150 (MQ=255)
                  tATGTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATAtt                                                                                                                                  >  2:2185309/1‑150 (MQ=255)
                    tGTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTc                                                                                                                                <  2:4749555/150‑1 (MQ=255)
                     gTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCa                                                                                                                               >  2:482402/1‑150 (MQ=255)
                     gTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCa                                                                                                                               <  1:4750311/150‑1 (MQ=255)
                      ttCACTCAGGCTCATGTTTGTCTGGAACAGACTTTTTATACTGCGCTGAGCCAGGATGGAACTAAAGGAAGTTAGGTGTTCGGATACGGTCTTGTTGGATTCATTAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCATCATATTTCAc                                                                                                                              <  1:2050538/150‑1 (MQ=255)
                      ttCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCAc                                                                                                                              <  1:3060733/150‑1 (MQ=255)
                         aCTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGAc                                                                                                                           >  2:3449932/1‑150 (MQ=255)
                         aCTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGAc                                                                                                                           >  1:3409119/1‑150 (MQ=255)
                          cTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACa                                                                                                                          >  2:687847/1‑150 (MQ=255)
                          cTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACa                                                                                                                          <  1:5387362/150‑1 (MQ=255)
                              ggCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTTACTTTTCAGGaatccat                                                                                                                      <  1:451656/150‑7 (MQ=255)
                               gCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACg                                                                                                                     >  1:159956/1‑150 (MQ=255)
                                cacaTGATGGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGa                                                                                                                    <  1:604021/150‑1 (MQ=255)
                                 acaTGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGAt                                                                                                                   <  2:1603972/150‑1 (MQ=255)
                                 acaTGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGAt                                                                                                                   <  2:1415765/150‑1 (MQ=255)
                                 acaTGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGAt                                                                                                                   <  1:3847107/150‑1 (MQ=255)
                                 acaTGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGAt                                                                                                                   <  1:3846996/150‑1 (MQ=255)
                                    tGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGAc                                                                                                                >  2:5111442/1‑150 (MQ=255)
                                     gATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACt                                                                                                               >  2:72942/1‑150 (MQ=255)
                                       ttGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGt                                                                                                             >  1:417766/1‑150 (MQ=255)
                                       ttGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGt                                                                                                             >  2:1271919/1‑150 (MQ=255)
                                              gAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTtgatga                                                                                                      >  2:3808381/1‑150 (MQ=255)
                                              gAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCAGGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTtgatga                                                                                                      <  1:1965527/150‑1 (MQ=255)
                                                aCAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACa                                                                                                    <  1:5109799/150‑1 (MQ=255)
                                                  aGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGa                                                                                                  <  2:164862/150‑1 (MQ=255)
                                                      tgtgAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATcac                                                                                              >  1:2107324/1‑150 (MQ=255)
                                                       gtgAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACt                                                                                                               >  2:4532797/1‑132 (MQ=255)
                                                       gtgAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATcaca                                                                                             <  2:4099503/150‑1 (MQ=255)
                                                       gtgAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATcaca                                                                                             >  1:2289943/1‑150 (MQ=255)
                                                       gcgAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATcaca                                                                                             >  1:4787171/3‑150 (MQ=255)
                                                        tgAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACat                                                                                            >  1:3710014/1‑150 (MQ=255)
                                                            tACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACattatt                                                                                        <  2:2861368/150‑1 (MQ=255)
                                                              cTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGCCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCt                                                                                      <  2:2211984/150‑1 (MQ=255)
                                                               tGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTg                                                                                     >  1:1977808/1‑150 (MQ=255)
                                                                    tGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATGATTCTGCcatc                                                                                >  2:2387340/1‑150 (MQ=255)
                                                                      aGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCcatcat                                                                              >  1:365044/1‑150 (MQ=255)
                                                                        ccAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTg                                                                            >  1:1096582/1‑150 (MQ=255)
                                                                          aGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCt                                                                          >  2:1663921/1‑150 (MQ=255)
                                                                             aTGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAg                                                                       <  1:3527227/150‑1 (MQ=255)
                                                                                  aCTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCata                                                                                                                                    <  1:5360527/84‑1 (MQ=255)
                                                                                  aCTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCata                                                                                                                                    >  2:5360527/1‑84 (MQ=255)
                                                                                  aCTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGcc                                                                  >  2:153819/1‑150 (MQ=255)
                                                                                  aCTCAAGGCACTGAGCTGTTCGGCAACGGCCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGcc                                                                  <  2:3378099/150‑1 (MQ=255)
                                                                                    tCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCtt                                                                >  1:5299504/1‑150 (MQ=255)
                                                                                       aGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGa                                                                                                                 <  2:1972941/98‑1 (MQ=255)
                                                                                       aGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGa                                                                                                                 >  1:1972941/1‑98 (MQ=255)
                                                                                       aGGCACTGAGCTGTTCGGCAACGGTCACGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCt                                                                 >  2:743135/1‑146 (MQ=255)
                                                                                            cTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTg                                                        <  2:3937146/150‑1 (MQ=255)
                                                                                            cTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTg                                                        <  2:5128269/150‑1 (MQ=255)
                                                                                            cTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTg                                                        <  1:842089/150‑1 (MQ=255)
                                                                                            cTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTg                                                        <  1:123476/150‑1 (MQ=255)
                                                                                            cTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTg                                                        <  1:1582890/150‑1 (MQ=255)
                                                                                             tGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGttt                                                                                                           <  1:4009414/98‑1 (MQ=255)
                                                                                                  tGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCAcc                                                  >  1:2770994/1‑150 (MQ=255)
                                                                                                    ttCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCtt                                                >  1:1104773/1‑150 (MQ=255)
                                                                                                    ttCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGGTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCtt                                                >  2:5255904/1‑150 (MQ=255)
                                                                                                            cGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTgg                                        <  2:1100016/150‑1 (MQ=255)
                                                                                                             ggTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGc                                       >  2:2437362/1‑150 (MQ=255)
                                                                                                               tCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGttt                                                                                                           >  2:2104654/1‑80 (MQ=255)
                                                                                                               tCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGttt                                                                                                           <  1:2104654/80‑1 (MQ=255)
                                                                                                               tCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCAc                                     >  2:5404851/1‑150 (MQ=255)
                                                                                                                  tgatgGATTCCTTAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGa                                  <  1:4852152/150‑1 (MQ=255)
                                                                                                                  tgatgGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGa                                  >  2:3970136/1‑150 (MQ=255)
                                                                                                                  tgatgGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGa                                  >  1:4470106/1‑150 (MQ=255)
                                                                                                                   gatgGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAg                                 <  2:1941117/150‑1 (MQ=255)
                                                                                                                   gatgGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAg                                 >  2:456858/1‑150 (MQ=255)
                                                                                                                   gatgGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAg                                 <  1:5046197/150‑1 (MQ=255)
                                                                                                                     tgGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGt                               >  2:2870751/1‑150 (MQ=255)
                                                                                                                        aTTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGt                            <  2:1741939/150‑1 (MQ=255)
                                                                                                                             tGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATc                       <  1:1467293/150‑1 (MQ=255)
                                                                                                                             tGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATc                       >  1:778442/1‑150 (MQ=255)
                                                                                                                             tGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATc                       >  1:656167/1‑150 (MQ=255)
                                                                                                                             tGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCAAATTTCACGACAGCACGATGACTGTTTGATGACAGATCAAATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATc                       >  1:1142785/1‑150 (MQ=255)
                                                                                                                              gaaatatAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATCa                      >  1:2951464/7‑150 (MQ=255)
                                                                                                                               aaaaGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATCAg                     >  2:3627650/1‑150 (MQ=255)
                                                                                                                               aaaaGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATCAg                     >  2:1409952/1‑150 (MQ=255)
                                                                                                                                     gaGGGCCGGGCATGCAGATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGt                                                                     >  1:2339279/2‑96 (MQ=38)
                                                                                                                                       gggCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATCAGGATGCGAc             <  2:2361530/150‑1 (MQ=255)
                                                                                                                                        ggCCGGGCATGCAGATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATCAGGATGCGAcc            >  1:2258140/1‑150 (MQ=255)
                                                                                                                                        ggCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATCAGGATGCGAcc            >  1:658464/1‑150 (MQ=255)
                                                                                                                                        ggCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATCAGGATGCGAcc            >  1:3833091/1‑150 (MQ=255)
                                                                                                                                        ggCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATCAGGATGCGAcc            >  1:4958801/1‑150 (MQ=255)
                                                                                                                                         gCCGGGCATGCAGATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATCAGGATGCGACCg           >  1:3926172/1‑150 (MQ=255)
                                                                                                                                            gggCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATCAGGATGCGACCGGtt        >  2:1313537/1‑150 (MQ=255)
                                                                                                                                              gCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATcaca                                                                                             <  1:3758465/63‑1 (MQ=255)
                                                                                                                                              gCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATcaca                                                                                             >  2:3758465/1‑63 (MQ=255)
                                                                                                                                               cATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATCAGGATGCGACCGGTTTcc     <  2:4809843/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CCTCATGGCCGAATACCCTATGTTCACTCAGGCACATGATTGGCTGGAACAGACTGTGAATACTGCGCTGAGCCAGGATGGCACTCAAGGCACTGAGCTGTTCGGCAACGGTCATGATGGATTCCTGAAAAGTAAGGGCCGGGCATCTGCATGCCCGGCCCTCATATTTCACGACAGCACGATGACTGTTTGATGACAGATCACATTATTCTGCCATCATTGCTGAGGTGCCTTACTTCTTGGCCACCTTGGTGCTGGCACTGAGGTAGTCGATCAGGATGCGACCGGTTTCCGCG  >  NC_002947/1918152‑1918447

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