Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,689,167:1 +C 5.0% coding (421/555 nt) pilI ← twitching motility protein PilI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,689,1671.C 5.0% 416.3 / 12.5 99coding (421/555 nt)pilItwitching motility protein PilI
Reads supporting (aligned to +/- strand):  ref base . (44/50);  new base C (3/2);  total (47/52)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCACGTTCGCAGGTCCTATAACGCCACATCGAGAAAGCCCGGCGCCTGGGCCAGGGCAAAGGGGCTGAAGATCGCCCAGTTACGTTCACGCGGAAAATACCCCTGCACAAAGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGCTGAAGGCCCAGCACCTCGTCGACCAGCAAACCGACGAACAGGTCGTCGTGATCCAGCACCAGCACTCGCCGCTGCTTGCCCGGTGCGGCGTGGCCAAGGCCGAAGAAGCAGCTCAGGT  >  NC_002947/5689018‑5689310
                                                                                                                                                      |                                                                                                                                               
tCACGTTCGCAGGTCCTATAACGCCACATCGAGAAAGCCCGGCGCCTGGGCCAGGGCAAAGGGGCTGAAGATCGCCCAGTTACGTTCACGCGGAAAATACCCCTGCACAAAGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAg                                                                                                                                                  <  1:2183157/150‑1 (MQ=255)
  aCGTTCGCAGGTCCTATAACGCCACATCGAGAAAGCCCGGCGCCTGGGCCAGGGCAAAGGGGCTGAAGATCGCCCAGTTACGTTCACGCGGAAAATACCCCTGCACAAAGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑Ct                                                                                                                                               >  1:563294/1‑150 (MQ=255)
    gTTCGCAGGTCCTATAACGCCACATCGAGAAAGCCCGGCGCCTGGGCCAGGGCAAAGGGGCTGAAGATCGCCCAGTTACGTTCACGCGGAAAATACCCCTGCACAAAGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGc                                                                                                                                             >  1:636862/1‑150 (MQ=255)
         cAGGTCCTATAACGCCACATCGAGAAAGCCCGGCGCCTGGGCCAGGGCAAAGGGGCTGAAGATCGCCCAGTTACGTTCACGCGGAAAATACCCCTGCACAAAGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCt                                                                                                                                        >  2:1602069/1‑150 (MQ=255)
            gTCCTATAACGCCACATCGAGAAAGCCCGGCGCCTGGGCCAGGGCAAAGGGGCTGAAGATCGCCCAGTTACGTTCACGCGGAAAATACCCCTGCACAAAGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCt                                                                                                                                        <  2:300850/147‑1 (MQ=255)
            gTCCTATAACGCCACATCGAGAAAGCCCGGCGCCTGGGCCAGGGCAAAGGGGCTGAAGATCGCCCAGTTACGTTCACGCGGAAAATACCCCTGCACAAAGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCt                                                                                                                                        >  1:300850/1‑147 (MQ=255)
              ccgataACGCCACATCGAGAAAGCCCGGCGCCTGGGCCAGGGCAAAGGGGCTGAAGATCGCCCCGTTACGTTCACGCTGAAAATACCCCTGCACAAAGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGa                                                                                                                                   <  1:1208735/147‑1 (MQ=255)
               cTATAACGCCACATCGAGAAAGCCCGGCGCCTGGGCCAGGGCAAAGGGGCTGAAGATCGCCCAGTTACGTTCACGCGGAAAATACCCCTGCACAAAGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAg                                                                                                                                  <  2:2389140/150‑1 (MQ=255)
               cTATAACGCCACATCGAGAAAGCCCGGCGCCTGGGCCAGGGCAAAGGGGCTGAAGATCGCCCAGTTACGTTCACGCGGAAAATACCCCTGCACAAAGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCGGTCGAg                                                                                                                                  <  2:636862/150‑1 (MQ=255)
                tataACGCCACATCGAGAAAGCCCGGCGCCTGGGCCAGGGCAAAGGGGCTGAAGATCGCCCAGTTACGTTCACGCGGAAAATACCCCTGCACAAAGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAgg                                                                                                                                 <  1:2206759/150‑1 (MQ=255)
                 ataACGCCACATCGAGAAAGCCCGGCGCCTGGGCCAGGGCAAAGGGGCTGAAGATCGCCCAGTTACGTTCACGCGGAAAATACCCCTGCACAAAGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAggg                                                                                                                                <  1:1607488/150‑1 (MQ=255)
                     cGCCACATCGAGAAAGCCCGGCGCCTGGGCCAGGGCAAAGGGGCTGAAGATCGCCCAGTTACGTTCACGCGGAAAATACCCCTGCACAAAGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCaaa                                                                                                                            <  2:822023/150‑1 (MQ=255)
                     cGCCACATCGAGAAAGCCCGGCGCCTGGGCCAGGGCAAAGGGGCTGAAGATCGCCCAGTTACGTTCACGCGGAAAATACCCCTGCACAAAGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCaaa                                                                                                                            <  1:322826/150‑1 (MQ=255)
                         acaTCGAGAAAGCCCGGCGCCTGGGCCAGGGCAAAGGGGCTGAAGATCGACCAGTTACGTTCACGCGGAAAATACCCCTGCACAAAGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGc                                                                                                                        <  2:995515/150‑1 (MQ=255)
                          caTCGAGAAAGCCCGGCGCCTGGGCCAGGGCAAAGGGGCTGAAGATCGCCCAGTTACGTTCACGCGGAAAATACCCCTGCACAAAGCGTGCAGCGGCCAGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGCt                                                                                                                       <  1:1267395/150‑1 (MQ=255)
                           aTCGAGAAAGCCCGGCGCCTGGGCCAGGGCAAAGGGGCTGAAGATCGCCCAGTTACGTTCACGCGGAAAATACCCCTGCACAAAGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGCTg                                                                                                                      >  1:1163010/1‑150 (MQ=255)
                            tCGAGAAAGCCCGGCGCCTGGGCCAGGGCAAAGGGGCTGAAGATCGCCCAGTTACGTTCACGCGGAAAATACCCCTGCACAAAGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGCTGa                                                                                                                     >  2:2280502/1‑150 (MQ=255)
                             cGAGAACGCCCGGCGCCTGGGCCCGGGCAAAGGGGCTGAAGATCGCCCCGTTACGTTCACGCGGAAAATACCCCTGCACAAAGCGTGCAGCGGCAGGGATCCGCGGCTGCGGTGGCGACAG‑CTGCAGGCTGTCCAGCGCAAAGTGCTGaa                                                                                                                    <  2:1005257/150‑1 (MQ=255)
                              gagaAAGCCCGGCGCCTGGGCCAGGGCAAAGGGGCTGAAGATCGCCCAGTTACGTTCACGCGGAAAATACCCCTGCACAAAGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAGCCTGCA‑GCTGTc                                                                                                                                     >  2:2362112/1‑132 (MQ=255)
                              gagaAAGCCCGGCGCCTGGGCCAGGGCAAAGGGGCTGAAGATCGCCCAGTTACGTTCACGCGGAAAATACCCCTGCACAAAGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAGCCTGCA‑GCTGTc                                                                                                                                     <  1:2362112/132‑1 (MQ=255)
                               agaAAGCCCGGCGCCTGGGCCAGGGCAAAGGGGCTGAAGATCGCCCAGTTACGTTCACGCGGAAAATACCCCTGCACAAAGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGCTGAAgg                                                                                                                  <  2:2131154/150‑1 (MQ=255)
                               agaAAGCCCGGCGCCTGGGCCAGGGCAAAGGGGCTGAAGATCGCCCAGTTACGTTCACGCGGAAAATACCCCTGCACAAAGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGCTGAAgg                                                                                                                  <  1:274619/150‑1 (MQ=255)
                                 aaaGCCCGGCGCCTGGGCCAGGGCAAAGGGGCTGAAGATCGCCCAGTTACGTTCACGCGGAAAATACCCCTGCACAAAGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGCTGAAGGcc                                                                                                                <  2:563294/150‑1 (MQ=255)
                                 aaaGCCCGGCGCCTGGGCCAGGGCAAAGGGGCTGAAGATCGCCCAGTTACGTTCACGCGGAAAATACCCCTGCACAAAGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGCTGAAGGcc                                                                                                                >  1:545061/1‑150 (MQ=255)
                                    gCCCGGCGCCTGGGCCAGGGCAAAGGGGCTGAAGATCGCCCAGTTACGTTCACGCGGAAAATACCCCTGCACAAAGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGCTGAAGGCCCAg                                                                                                             >  2:341265/1‑150 (MQ=255)
                                    gCCCGGCGCCTGGGCCAGGGCAAAGGGGCTGAAGATCGCCCAGTTACGTTCACGCGGAAAATACCCCTGCACAAAGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGCTGAAGGCCCAg                                                                                                             <  2:848567/150‑1 (MQ=255)
                                          cgcCTGGGCCAGGGCAAAGGGGCTGAAGATCGCCCAGTTACGTTCACGCGGAAAATACCCCTGCACAAAGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAGAGTGCTGAAGGCCCAGCACCtc                                                                                                       <  2:1756731/150‑1 (MQ=255)
                                             cTGGGCCAGGGCAAAGGGGCTGAAGATCGCCCAGTTACGTTCACGCGGAAAATACCCCTGCACAAAGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAC‑CTGCAGGCTGTGGGGTGCAAAGTGCTGAAGGCCAAGAAAAtcgtc                                                                                                    >  2:2201525/1‑150 (MQ=255)
                                                ggCCAGGGCAAAGGGGCTGAAGATCGCCCAGTTACGTTCACGCGGAAAATACCCCTGCACAAAGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGCTGAAGGCCCAGCACCTCGTCGAc                                                                                                 <  2:1528373/150‑1 (MQ=255)
                                                ggCCAGGGCAAAGGGGCTGAAGATCGCCCAGTTACGTTCACGCGGAAAATACCCCTGCACAAAGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGCTGAAGGCCCAGCACCTCGTCGAc                                                                                                 >  2:1924123/1‑150 (MQ=255)
                                                 gCCAGGGCAAAGGGGCTGAAGATCGCCCAGTTACGTTCACGCGGAAAATACCCCTGCACAAAGCGTGCAGCGGCAGGGATCAGTGGCTGCGGCGGCGACAGCCTGCA‑GCTGTc                                                                                                                                     <  2:611774/113‑1 (MQ=255)
                                                 gCCAGGGCAAAGGGGCTGAAGATCGCCCAGTTACGTTCACGCGGAAAATACCCCTGCACAAAGCGTGCAGCGGCAGGGATCAGTGGCTGCGGCGGCGACAGCCTGCA‑GCTGTc                                                                                                                                     >  1:611774/1‑113 (MQ=255)
                                                   cAGGGCAAAGGGGCTGAAGATCGCCCAGTTACGTTCACGCGGAAAATACCCCTGCACAAAGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGCTGAAGGCCCAGCACCTCGTCGACCAg                                                                                              >  1:2311864/1‑150 (MQ=255)
                                                    aGGGCAAAGGGGCTGAAGATCGCCCAGTTACGTTCACGCGGAAAATACCCCTGCACAAAGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGCTGAAGGCCCAGCACCTCGTCGACCGGc                                                                                             <  1:1924123/150‑1 (MQ=255)
                                                     gggCAAAGGGGCTGAAGATCGCCCAGTTACGTTCACGCGGAAAATACCCCTGCACAAAGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTACAGGCTGTCGAGGGCAAAGTGCTGAAGGCCCAGCACCTCGTCGACCAGCa                                                                                            >  1:801993/1‑150 (MQ=255)
                                                       gCAAAGGGGCTGAAGATCGCCCAGTTACGTTCACGCGGAAAATACCCCTGCACAAAGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGCTGAAGGCCCAGCACCTCGTCGACCAGCaaa                                                                                          <  2:659067/150‑1 (MQ=255)
                                                               gCTGAAGATCGCCCAGTTACGTTCACGCGGAAAATACCCCTGCACAAAGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGc                                                                                                                        <  2:831879/112‑1 (MQ=255)
                                                               gCTGAAGATCGCCCAGTTACGTTCACGCGGAAAATACCCCTGCACAAAGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGc                                                                                                                        >  1:831879/1‑112 (MQ=255)
                                                               gCTGAAGATCGCCCAGTTACGTTCACGCGGAAAATACCCCTGCACAAAGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGCTGAAGGCCCAGCACCTCGTCGACCAGCAAACCGACGaa                                                                                  <  1:1295938/150‑1 (MQ=255)
                                                               gCTGAAGATCGCCCAGTTACGTTCACGCGGAAAATACCCCTGCACAAAGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGCTGAAGGCCCAGCACCTCGTCGACCAGCAAACCGACGaa                                                                                  >  1:2215770/1‑150 (MQ=255)
                                                                cTGAAGATCGCCCAGTTACGTTCACGCGGAAAATACCCCTGCACAAAGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGCTGAAGGCCCAGCACCTCGTCGACCAGCAAACCGACGAAc                                                                                 <  2:1564185/150‑1 (MQ=255)
                                                                 tGAAGATCGCCCAGTTACGTTCACGCGGAAAATACCCCTGCACAAAGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAgg                                                                                                                                 >  2:2378331/1‑101 (MQ=255)
                                                                 tGAAGATCGCCCAGTTACGTTCACGCGGAAAATACCCCTGCACAAAGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAgg                                                                                                                                 <  1:2378331/101‑1 (MQ=255)
                                                                  gAAGATCGCCCAGTTACGTTCACGCGGAAAATACCCCTGCACAAAGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGCTGAAGGCCCAGCACCTCGTCGACCAGCAAACCGACGAACAg                                                                               <  1:341265/150‑1 (MQ=255)
                                                                  gAAGATCGCCCAGTTACGTTCACGCGGAAAATACCCCTGCACAAAGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGCTGAAGGCCCAGCACCTCGTCGACCAGCAAACAGACGAACAg                                                                               >  1:1752669/1‑150 (MQ=255)
                                                                     gATCGCCCAGTTACGTTCACGCGGAAAATACCCCTGCACAAAGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGCTGAAGGCCCAGCACCTCGTCGAcc                                                                                                >  1:493108/1‑130 (MQ=255)
                                                                     gATCGCCCAGTTACGTTCACGCGGAAAATACCCCTGCACAAAGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGCTGAAGGCCCAGCACCTCGTCGAcc                                                                                                <  2:493108/130‑1 (MQ=255)
                                                                     gATCGCCCAGTTACGTTCACGCGGAAAATACCCCTGCACAAAGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGCTGAAGGCCCAGCACCTCGTCGACCAGCAAACCGACGAACAGgtc                                                                            >  1:312909/1‑150 (MQ=255)
                                                                      aTCGCCCAGTTACGTTCACGCGGAAAATACCCCTGCACAAAGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGCTGAAGGCCCAGTACCTCGTCGACCAGCAAACCGACGAACAGgtcg                                                                           >  1:482858/1‑150 (MQ=255)
                                                                          cccAGTTACGTTCACGCGGAAAATACCCCTGCACAAGGCGTGCATCGGCATGGAGCAGCGGCTGCGGCGTCGCCAA‑CTGAACGGAAGCGCACGCATATTGCTGAAGGCCCAGCACCTCGTCGACCAGCAAACCGACGAACAGGTcgtcgt                                                                       >  1:2350486/1‑150 (MQ=255)
                                                                          cccAGTTACGTTCACGCGGAAAATACCCCTGCACAAAGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGCTGAAGGCCCAGCACCTCGTCGACCAGCAAACCGACGAACAGGTcgtcgt                                                                       >  1:2305185/1‑150 (MQ=255)
                                                                          cccAGTTACGTTCACGCGGAAAATACCCCTGCACAAAGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGCTGAAGGCCCAGCACCTCGTCGACCAGCAAACCGACGAACAGGTcgtcgt                                                                       <  2:312909/150‑1 (MQ=255)
                                                                           ccAGTTACGTTCACGCGGAAAATACCCCTGCACAAAGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGCTGAAGGCCCAGCACCTCGTCGACCAGCAAACCGACGAACAGGTCGTCGTg                                                                      >  2:1630085/1‑150 (MQ=255)
                                                                            cAGTTACGTTCACGCGGAAAATACCCCTGCACAAAGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGCTGAAGGCCCAGCACCTCGTCGACCAGCAAACCGACGAACAGGTCGTCGTGa                                                                     <  2:545061/150‑1 (MQ=255)
                                                                               ttACGTTCACGCGGAAAATACCCCTGCACAAAGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGCTGAAGGCCCAGCACCTCGTCGACCAGCAAACCGACGAACAGGTCGTCGTGATcc                                                                  <  1:1630085/150‑1 (MQ=255)
                                                                               ttACGTTCACGCGGAAAATACCCCTGCACAAAGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGCTGAAGGCCCAGCACCTCGTCGACCAGCAAACCGACCAACAGGTCGTCGTGATcc                                                                  <  1:94848/150‑1 (MQ=255)
                                                                                   gTTCACGCGGAAAATACCCCTGCACAAAGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGCTGAAGGCCCAGCACCTCGTCGACCAGCAAACCGACGAACAGGTCGTCGTGATccagca                                                              <  2:1163010/150‑1 (MQ=255)
                                                                                      cACGCGGAAAATACCCCTGCACAAAGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGCTGAAGGCCCAGCACCTCGTCGACCAGCAAACCGACGAACAGGTCGTCGTGATccagcacca                                                           <  2:2305185/150‑1 (MQ=255)
                                                                                         gcgGAAAATACCCCTGCACAAAGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGCTGAAGGCCCAGCACCTCGTCGACCAGCAAACCGACGAACAGGTCGTCGTGATccagcaccagca                                                        >  1:82344/1‑150 (MQ=255)
                                                                                               aaTACCCCTGCACAAAGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGCTGAAGGCCCAGCACCTCGTCGACCAGCAAACcga                                                                                      <  1:1525323/114‑1 (MQ=255)
                                                                                               aaTACCCCTGCACAAAGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGCTGAAGGCCCAGCACCTCGTCGACCAGCAAACcga                                                                                      >  2:1525323/1‑114 (MQ=255)
                                                                                               aaTACCCCTGCACAAAGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGCTGAAGGCCCAGCACCTCGTCGACCAGCAAACCGACGAACAGGTCGTCGTGATCCAGCACCAGCACTcgcc                                                  >  1:2055219/1‑150 (MQ=255)
                                                                                                    cccTGCACAAAGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGCTGGAGGCCCAGCACCTCGTCGACCAGCAAACCGACGAACAGGTCGTCGTGATCCAGCACCAGCACTCGCCgctgc                                             >  2:1058305/1‑150 (MQ=255)
                                                                                                    cccTGCACAAAGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGCTGAAGGCCCAGCACCTCGTCGACCAGCAAACCGACGAACAGGTCGTCGTGATCCAGCACCAGCACTCGCCgctgc                                             >  2:2249178/1‑150 (MQ=255)
                                                                                                      cTGCACAAAGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGCTGAAGGCCCAGCACCTCGTCGACCAGCAAACCGACGAACAGGTcgtcgt                                                                       >  1:990188/1‑122 (MQ=255)
                                                                                                      cTGCACAAAGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGCTGAAGGCCCAGCACCTCGTCGACCAGCAAACCGACGAACAGGTcgtcgt                                                                       <  2:990188/122‑1 (MQ=255)
                                                                                                      cTGCACAAAGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGCTGAAGGCCCAGCACCTCGTCGACCAGCAAACCGACGAACAGGTCGTCGTGATCCAGCACCAGCACTCGCCGCTGCtt                                           >  1:743987/1‑150 (MQ=255)
                                                                                                      cTGCACAAAGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGGCGAGGGCAAAGTGCTGAAGGCCCAGCACCTCGTCGACCAGCAAACCGACGAACAGGTCGTCGTGATCCAGCACCAGCACTCGCCGCTGCtt                                           <  1:1602069/150‑1 (MQ=255)
                                                                                                          acaAAGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGCTGAAGGCCCAGCACCTCGTCGACCAGCAAACCGACGAACAGGTCGTCGTGATCCAGCACCAGCACTCGCCGCTGCTTGccc                                       >  1:1311654/1‑150 (MQ=255)
                                                                                                          acaAAGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGCTGAAGGCCCAGCACCTCGTCGACCAGCAAACCGACGAACAGGTCGTCGTGATCCAGCACCAGCACTCGCCGCTGCTTGccc                                       >  1:1310555/1‑150 (MQ=255)
                                                                                                          acaAAGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGCTGAAGGCCCAGCACCTCGTCGACCAGCAAACCGACGAACAGGTCGTCGTGATCCAGCACCAGCACTCGCCGCTGCTTGcca                                       >  1:2327542/1‑149 (MQ=255)
                                                                                                            aaaGCGTGCAGCGGCAGGGATCAGCGTCTGCGGCGGCGACAG‑GCGCAGGCTGTCGAGGGCAAAGTGCTGAAGGCCCAGCACCTCGTCGACCAGCAAACCGACGAACAGGTCGTCGTGATCCAGCACCAGCACTCGCCGCTGCTTGCCCgg                                     >  1:328795/1‑150 (MQ=255)
                                                                                                            aaaGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGCTGAAGGCCCAGCACCTCGTCGACCAGCAAACCGACGAACAgg                                                                              >  1:2199933/1‑109 (MQ=255)
                                                                                                            aaaGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGCTGAAGGCCCAGCACCTCGTCGACCAGCAAACCGACGAACAgg                                                                              <  2:2199933/109‑1 (MQ=255)
                                                                                                            aaaGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGCTGAAGGCCCAGCACCTCGTCGACCAGCAAACCGACGAACAGGTCGTCGTGATCCAGCACCAGCACTCGCCGCTGCTTGCCCgg                                     >  1:1679514/1‑150 (MQ=255)
                                                                                                            aaaGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGCTGAAGGCCCAGCACCTCGTCGACCAGCAAACCGACGAACAGGTCGTCGTGATCCAGCACCAGCACTCGCCGCTGCTTGCCCgg                                     >  1:1561675/1‑150 (MQ=255)
                                                                                                             aaGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGCTGAAGGCCCAGCACCTCGTCGACCAGCAAACCGACGAACAGGTCGTCGTGATCCAGCACCAGCACTCGCCGCTGCTTGCCCGGt                                    <  2:82344/150‑1 (MQ=255)
                                                                                                                  tGCAGCGGTAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGCTGAACGCCCAGCACCTCGTCGACCAGCAAACCGACGAACAGGTCGTCGTGATCCAGCACCAGCACTCGCCGCTGCTTGCCCGGTgcggc                               >  2:1734812/1‑150 (MQ=255)
                                                                                                                  tGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGCTGAAGGCCCAGCACCTCGTCGACCAGCAAACCGACGAACAGGTCGTCGTGATCCAGCACCAGCACTCGCCGCTGCTTGCCCGGTgcggc                               <  2:2215770/150‑1 (MQ=255)
                                                                                                                  tGCAGCGGAAGGGATCAGTGGCTGCGGCGCCGCTAG‑CTTCAGGCTGTCGAGGGCAAAGTGTTGAAGGCCTAGCATCTTGTCGACCAGCAAACCGACGAACAGGTCGTCGTGATCCAGCACCAGCACTCGCCGCTGCTTGCCCGGTgcggc                               <  2:801993/150‑1 (MQ=255)
                                                                                                                     aGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGCTGAAGGCCCAGCACCTCGTCGACCAGCAAACCGACGAACAGGTCGTCGTGATCCAGCACCAGCACTCGCCGCTGCTTGCCCGGTGCGGCGTg                            <  2:1310555/150‑1 (MQ=255)
                                                                                                                     aGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGCTGAAGGCCCAGCACCTCGTCGACCAGCAAACCGACGAACAGGTCGTCGTGATCCAGCACCAGCACTCGCCGCTGCTTGCCCGGTGCGGCGTg                            <  2:1311654/150‑1 (MQ=255)
                                                                                                                        ggCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGCTGAAGGCCCAGCACCTCGTCGACCAGCAAACCGACGAACAGGTCGTCGTGATCCAGCACCAGCACTCGCCGCTGCTTGCCCGGTgcgggggtatt                         <  2:2055219/150‑8 (MQ=255)
                                                                                                                          cAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGGTGAAGGCCCAGCACCTCGTCGACCAGCAAACCGACGAACAGGTCGTCGTGATCCAGCACCAGCACTCGCCGCTGCTTGCCCGGTGCGGCGTGGCCaa                       <  1:2201525/150‑1 (MQ=255)
                                                                                                                              gATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGCTGAAGGCCCAGCACCTCGTCGACCAGCAAACCGACGAACAGGTCGTCGTGa                                                                     >  1:1018303/1‑100 (MQ=255)
                                                                                                                              gATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGCTGAAGGCCCAGCACCTCGTCGACCAGCAAACCGACGAACAGGTCGTCGTGa                                                                     <  2:1018303/100‑1 (MQ=255)
                                                                                                                                tCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGCTGAAGGCCCAGCACCTCGTCGACCAGCaaa                                                                                          >  1:819242/1‑77 (MQ=255)
                                                                                                                                tCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGCTGAAGGCCCAGCACCTCGTCGACCAGCaaa                                                                                          <  2:819242/77‑1 (MQ=255)
                                                                                                                                tCAGCGGCTGCGGCGGCGACAG‑CTGACGGCTGTCGAGGGCAAAGTGCTGAAGGCCCAGCACCTCGTCGACCAGCAAACCGACGAACAGGTCGTCGTGATCCAGCACCAGCACTCGCCGCTGCTTGCCCGGTGCGGCGTGGCCAAGGCCga                 >  2:25443/1‑150 (MQ=255)
                                                                                                                                 cAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGCTGAAGGCCCAGCACCTCGTCGACCAGCAAACCGACGAACAGGTCGTCGTGATCCAGCACCAGCACTCGCCGCTGCTTGCCCGGTGCGGCGTGGCCAAGGCCgaa                <  2:2029039/150‑1 (MQ=255)
                                                                                                                                 cAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGCTGAAGGCCCAGCACCTCGTCGACCAGCAAACCGACGAACAGGTCGTCGTGATCCAGCACCAGCACTCGCCGCTGCTTGCCCGGTGCGGCGTGGCCAAGGCCgaa                <  2:2029776/150‑1 (MQ=255)
                                                                                                                                   gCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGCTGAAgg                                                                                                                  <  1:2208588/50‑1 (MQ=255)
                                                                                                                                   gCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGCTGAAgg                                                                                                                  >  2:2208588/1‑50 (MQ=255)
                                                                                                                                   gCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAACTGCTGAAGGCCCAGCACCTCGTCGACCAGCAAACCGACGAACAGGTCGTCGTGATACAGCACCAGCACTCGCCGCTGCTTGCCCGGTGCGGCGTGGCCAAGGCCgaaga              >  1:2293145/1‑150 (MQ=255)
                                                                                                                                       cTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGCTGAAGGCCCAGCACCTCGTCGACCAGCAAAc                                                                                         <  2:299393/71‑1 (MQ=255)
                                                                                                                                       cTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGCTGAAGGCCCAGCACCTCGTCGACCAGCAAAc                                                                                         >  1:299393/1‑71 (MQ=255)
                                                                                                                                       cTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGCTGAAGGCCCAGCACCTCGTCGACCAGCAAACCGACGAACAGGTCGTCGTGATCCAGCACCAGCACTCGCCGCTGCTTGCCCGGTGCGGCGTGGCCAAGGCCGAAGAagca          >  1:2398929/1‑150 (MQ=255)
                                                                                                                                              ggcgACAG‑CTGCAGGCTGTCGAGGGCAAAGTGCTGAAGGCCCAGCACCTCGTCGACCAGCAAACCGACGAACAGGTCGTCGTGATCCAGCACCAGCACTCGCCGCTGCTTGCCCGGTGCGGCGTGGCCAAGGCCGAAGAAGCAGCTCAgg   <  2:2350486/150‑1 (MQ=255)
                                                                                                                                              ctcgACAGCCTGCA‑GCTGTCGAGGGCAAAGTGCTGAAGGCCCAGCACCTCGTCGACCAGCAAACCGACGAACAGGTCGTCGTGATCCAGCACCAGCACTCGCCGCTGCTTGCCCGGTGCGGCGTGGCCAAGGCCGAAGAAGCAGCTCAgg   >  2:1449826/3‑150 (MQ=255)
                                                                                                                                               gcgACAG‑CTGCAGGCTGTCGAGGGCAAAGTGCTGAAGGCCCAGCACCTCGTCGACCAGCAAACAGACGAACAGGTCGTCGTGATCCAGCACCAGCACTCGCCGCTGCTTGCCCGGTGCGGCGTGGCCAAGGCCGAAGAAGCAGCTCAGGt  <  2:1752669/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                               
TCACGTTCGCAGGTCCTATAACGCCACATCGAGAAAGCCCGGCGCCTGGGCCAGGGCAAAGGGGCTGAAGATCGCCCAGTTACGTTCACGCGGAAAATACCCCTGCACAAAGCGTGCAGCGGCAGGGATCAGCGGCTGCGGCGGCGACAG‑CTGCAGGCTGTCGAGGGCAAAGTGCTGAAGGCCCAGCACCTCGTCGACCAGCAAACCGACGAACAGGTCGTCGTGATCCAGCACCAGCACTCGCCGCTGCTTGCCCGGTGCGGCGTGGCCAAGGCCGAAGAAGCAGCTCAGGT  >  NC_002947/5689018‑5689310

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