Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,462,441 C→A 7.9% A258S (GCG→TCG) ‡ chiA ← periplasmic endochitinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,462,4410CA 7.9% 2199.8 / 27.4 724A258S (GCG→TCG) ‡chiAperiplasmic endochitinase
Reads supporting (aligned to +/- strand):  ref base C (170/483);  new base A (21/48);  total (191/533)
Fisher's exact test for biased strand distribution p-value = 4.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 769 total aligned reads displayed.

GCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCTTC  >  CP009273/3462374‑3462507
                                                                   |                                                                  
gccgCACGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCcgc                                                                    <  1:25710390/68‑1 (MQ=255)
  cgcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCcgcg                                                                   <  1:24145404/67‑1 (MQ=255)
  cgcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCcgcg                                                                   <  1:37217879/67‑1 (MQ=255)
   gcgcgAACATCACGCCCTGGGGTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGcc                                                                 <  1:36530582/68‑1 (MQ=255)
   gcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGcc                                                                 <  1:12620893/68‑1 (MQ=255)
       gAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTAcc                                                             <  1:4578451/68‑1 (MQ=255)
         aCATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTAccc                                                            <  1:29904048/67‑1 (MQ=255)
         aCATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACcca                                                           <  1:28349196/68‑1 (MQ=255)
           aTAACGCCATGGGTTACTGGCATCGTTTTTTTTCGCATCGCCAGGCCAGTTTGAAGCGGCTACccacc                                                         >  1:37589932/1‑68 (MQ=14)
           aTAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACccacc                                                         >  1:14047147/1‑68 (MQ=255)
             aaCGCCATTGGTTACTTGCATCGTTATCATTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAg                                                       <  1:25893898/68‑1 (MQ=255)
              aCGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGAAGAGGATACCCACCAg                                                       >  1:5874705/1‑67 (MQ=18)
               cGCCATGGGTTTCTGGCATCGTTTTTTTTCGCATCGCCAGGGCATTTAGAAGAGCCTACCCACCAGGc                                                     >  1:18718156/1‑68 (MQ=11)
                gCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCATTTTGAAGCGCCTACCCACCAGGca                                                    >  1:8392091/1‑67 (MQ=25)
                gCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCg                                                    >  1:11589575/1‑68 (MQ=255)
                gCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCg                                                    >  1:10845072/1‑68 (MQ=255)
                gCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGCCAATTTGAAGCGGCTACCCACCAGGca                                                    >  1:31812568/1‑67 (MQ=21)
                gCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCAGGGCAATTAGCAGAGGCTACCCACCAGGca                                                    >  1:16234289/1‑67 (MQ=17)
                gCCATGGGTTACTGGCATCGTTTTATTTCGCATCGCCAGGGCAGTTTGAAGAGGCTACCCACCAGGca                                                    >  1:22771553/1‑67 (MQ=14)
                 ccATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCCCCAGGCGt                                                   >  1:30479928/1‑68 (MQ=255)
                 ccATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGt                                                   >  1:14673639/1‑68 (MQ=255)
                 ccATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGt                                                   >  1:32974343/1‑68 (MQ=255)
                 ccATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGt                                                   >  1:11075275/1‑68 (MQ=255)
                 ccATGGGTTACTGGCATCGTTTTCTTTCGCATATCCTGGGCAATTTGAAGAGGCTACCCACCAGGcat                                                   >  1:18444040/1‑66 (MQ=17)
                 ccATGGGTTACTGGCATCGTTTTATTTCGCATCGCCAGGGCAATTAGCAGCGCCTACCCACCAGGCGt                                                   >  1:2095768/1‑68 (MQ=25)
                  cATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGtt                                                  <  1:35861585/68‑1 (MQ=255)
                  cATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGtt                                                  <  1:13539071/68‑1 (MQ=255)
                   aTGGGTTATTGGCATCGTCTTCTTTGGTATCGCGTGGGCAATTTGTCGCGCCTACCCACCAGGCGTTc                                                 <  1:26725728/68‑1 (MQ=37)
                   aTGGGTTACTGGCCTCGTATTCTATCGCATCGCGTGGGCAATTTTCCGCGCCTACCCACCAGGCGTTc                                                 <  1:9314707/68‑1 (MQ=255)
                    tGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTTCCGCGCCTACCCACCAGGCGTTCt                                                <  1:24353067/68‑1 (MQ=255)
                    tGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt                                                <  1:25269928/68‑1 (MQ=255)
                    tGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt                                                <  1:28588950/68‑1 (MQ=255)
                    tGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt                                                <  1:5307691/68‑1 (MQ=255)
                       gTTACTGGCATCGTTTGCTATCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTta                                              <  1:21965000/67‑1 (MQ=255)
                            tGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAAt                                         >  1:12168663/1‑67 (MQ=255)
                            tGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAAt                                         >  1:16530539/1‑67 (MQ=255)
                            tGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGCCGTTATTATAAAt                                         >  1:16831679/1‑67 (MQ=255)
                             ggCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATg                                        >  1:7836242/1‑67 (MQ=255)
                             ggCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATg                                        >  1:11286761/1‑67 (MQ=255)
                                 tCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATgccgcc                                   <  1:22316938/68‑1 (MQ=255)
                                 tCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATgccgcc                                   <  1:21440059/68‑1 (MQ=255)
                                 tCGTTTTCTATCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATgccgcc                                   <  1:34843438/68‑1 (MQ=255)
                                   tttttCTTTCGCACCGCGTGGGCAATTTGTAGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                 <  1:38659571/67‑1 (MQ=25)
                                   tttttCTATGGCATCGCGTGGGCAATTTGTAGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                 <  1:34002757/67‑1 (MQ=25)
                                   ttcttctttcGCATCGCCTGGGCAATTTGTCGAGCCAACCCACCAGGCGTTCTTATAAATGCCGCCGt                                 <  1:33652983/65‑1 (MQ=37)
                                   ttcttctttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                 <  1:2704099/65‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                 <  1:8674/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                 <  1:8094326/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                 <  1:7124508/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                 <  1:9646193/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                 <  1:394781/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                 <  1:20120409/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                 <  1:26472793/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                 <  1:2070716/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                 <  1:13207227/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                 <  1:19036365/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                 <  1:18088175/68‑1 (MQ=255)
                                   gTTTTCTTTAGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                 <  1:33047261/68‑1 (MQ=255)
                                   cggttctttcGCATCGCCTGGGCAATTTGGCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                 <  1:20439212/65‑1 (MQ=255)
                                    ttttCTTTCGCATCGCCTGGGCAATTTTCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                 <  1:31597200/67‑1 (MQ=255)
                                    ttttCTATGGCACCGCGTGGGCAATATTCCGAGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGtt                                <  1:4975885/68‑1 (MQ=21)
                                      ttctttgGCATCGCGTGGGCAATATTAAGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTtat                              <  1:19884689/68‑1 (MQ=25)
                                      ttctttgGCATCCCCTGGGCAATTTTCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTtat                              <  1:28140609/68‑1 (MQ=38)
                                      ttctttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTtat                              <  1:36243301/68‑1 (MQ=255)
                                      ttctttcGCATCCCCTGGGCAATTTGACGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTtat                              <  1:7331419/68‑1 (MQ=255)
                                        ctttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAg                            >  1:23578729/1‑68 (MQ=255)
                                        ctttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAg                            <  1:6568916/68‑1 (MQ=255)
                                         tttcGCATCGCCTGGGCAATTTGCCGCGCCTCCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                          <  1:23177060/68‑1 (MQ=255)
                                          ttggCACCGCCTGGGCAATTTTCAGAGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                          <  1:2750292/65‑1 (MQ=25)
                                          ttcGGATCGCGTGGGCAATTTGCAGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                          <  1:19575247/68‑1 (MQ=38)
                                          ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                          <  1:15920271/68‑1 (MQ=255)
                                          ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                          <  1:13813835/68‑1 (MQ=255)
                                          ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                          <  1:8896830/68‑1 (MQ=255)
                                            cGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATaa                        <  1:27938596/68‑1 (MQ=255)
                                               aTCGCCTGGTCAATTTGGCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                     <  1:10046675/68‑1 (MQ=255)
                                               aTCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                     <  1:29236153/68‑1 (MQ=255)
                                                  gcgtGGGCAATTTGGAGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTAcc                  <  1:28810506/65‑1 (MQ=38)
                                                  gCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTAcc                  <  1:17916429/68‑1 (MQ=255)
                                                   ccTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCa                 >  1:32605851/1‑68 (MQ=255)
                                                    cTGGGCAATTTTTCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt                <  1:21155737/68‑1 (MQ=255)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt                <  1:38383557/68‑1 (MQ=255)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt                <  1:29693831/68‑1 (MQ=255)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt                <  1:24795072/68‑1 (MQ=255)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt                <  1:30681748/68‑1 (MQ=255)
                                                      gggCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGt              <  1:34470082/68‑1 (MQ=255)
                                                       ggCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTa             >  1:21765452/1‑68 (MQ=255)
                                                       ggCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTa             >  1:35487931/1‑68 (MQ=255)
                                                        gCAATTTTCCGAGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTa             <  1:14438397/67‑1 (MQ=255)
                                                        gCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTag            >  1:15476812/1‑68 (MQ=255)
                                                        gCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTag            >  1:4380312/1‑68 (MQ=255)
                                                        gCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTa             <  1:31138960/67‑1 (MQ=255)
                                                         cAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTag            <  1:12850436/67‑1 (MQ=255)
                                                            tttGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTagaaga        <  1:26256022/68‑1 (MQ=255)
                                                             ttGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTagaaga        >  1:15324129/1‑67 (MQ=255)
                                                               gCCGCGCCTACCCCCCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAgaagaa       >  1:6852836/1‑66 (MQ=255)
                                                                 cgcgCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCtt   <  1:17304107/68‑1 (MQ=255)
                                                                 cgcgCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCtt   <  1:30190783/68‑1 (MQ=255)
                                                                 cgcgCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCtt   <  1:631703/68‑1 (MQ=255)
                                                                 cgcgCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCtt   <  1:37926764/68‑1 (MQ=255)
                                                                 cgcgCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCtt   <  1:10590909/68‑1 (MQ=255)
                                                                   |                                                                  
GCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCTTC  >  CP009273/3462374‑3462507

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