Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,462,400 A→G 5.5% S271S (AGT→AGC chiA ← periplasmic endochitinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,462,4000AG 5.5% 1136.5 / 21.8 364S271S (AGT→AGCchiAperiplasmic endochitinase
Reads supporting (aligned to +/- strand):  ref base A (179/154);  new base G (14/6);  total (204/160)
Fisher's exact test for biased strand distribution p-value = 1.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 366 total aligned reads displayed.

AGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATA  >  CP009273/3462336‑3462465
                                                                |                                                                 
aGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTTCTg                                                                >  1:5755955/1‑68 (MQ=255)
aGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTTCTg                                                                >  1:25408629/1‑68 (MQ=255)
aGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTGCTg                                                                >  1:17902432/1‑68 (MQ=255)
aGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTg                                                                >  1:23388287/1‑68 (MQ=255)
aGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTg                                                                >  1:19702623/1‑68 (MQ=255)
aGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTg                                                                >  1:21937751/1‑68 (MQ=255)
  ccAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTTCTGGc                                                              >  1:1491034/1‑68 (MQ=255)
   cAGGTTTACCATACTGGCTGATTTCTGTTGCCGTTGCGGCTCGAACATAACGCCATGGGTTGCTGGCa                                                             <  1:24117455/68‑1 (MQ=18)
   cAGGGTTACCATACTGGGTGATTTCTGTTGCCGTAGCGGCACGAACATAACGCCATGGGTTGCTgg                                                               <  1:1713503/66‑1 (MQ=14)
   cAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTgg                                                               <  1:19468244/66‑1 (MQ=255)
     gggTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTGCTGGCATc                                                           >  1:7425879/1‑68 (MQ=255)
     gggTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATc                                                           >  1:2978249/1‑68 (MQ=255)
     gggTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATc                                                           >  1:3655122/1‑68 (MQ=255)
      ggTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATATCGCCATGGGTTGCTGGCATct                                                          >  1:20189429/1‑67 (MQ=25)
      ggTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTGCTGGCATca                                                          >  1:29178979/1‑67 (MQ=255)
       gTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCTCGAACATAACGCCATGGGTTGCTGGCATc                                                           <  1:27946846/66‑1 (MQ=255)
       gTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATcat                                                         >  1:27352109/1‑66 (MQ=255)
       gTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATcat                                                         >  1:12347632/1‑66 (MQ=255)
       gTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGt                                                         >  1:11577412/1‑68 (MQ=255)
       gTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGt                                                         >  1:12885362/1‑68 (MQ=255)
        ttACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGt                                                         <  1:5488306/67‑1 (MQ=255)
        ttACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGt                                                         <  1:1604561/67‑1 (MQ=255)
        ttACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGt                                                         <  1:18313532/67‑1 (MQ=255)
            cATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTTCTGGCATCAtttagt                                                    >  1:24169678/1‑65 (MQ=21)
            cATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGtttagt                                                    >  1:28793245/1‑65 (MQ=255)
            cATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTtttct                                                    <  1:9103105/68‑1 (MQ=255)
            cATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTtttct                                                    >  1:23248512/1‑68 (MQ=255)
            cATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTtttct                                                    >  1:21392931/1‑68 (MQ=255)
            cATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTtttct                                                    >  1:16953598/1‑68 (MQ=255)
            cATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCAtttagt                                                    >  1:8636020/1‑65 (MQ=255)
             aTACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTtttctt                                                   >  1:17792067/1‑68 (MQ=255)
             aTACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTtttctt                                                   >  1:28170488/1‑68 (MQ=255)
             aTACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCATtttatt                                                   >  1:15558124/1‑68 (MQ=255)
              tACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTtttctt                                                   >  1:24675962/1‑67 (MQ=255)
              tACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTtttctt                                                   >  1:2295304/1‑67 (MQ=255)
              tACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTtttctt                                                   >  1:25558723/1‑67 (MQ=255)
              tACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTtttctt                                                   >  1:6984928/1‑67 (MQ=255)
              tACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTtttctt                                                   >  1:10323139/1‑67 (MQ=255)
              tACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTtttctt                                                   >  1:1848202/1‑67 (MQ=255)
               aCTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTCCTGGCATCATttgttttc                                                 >  1:29341359/1‑68 (MQ=18)
               aCTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTtttctttc                                                 >  1:11139534/1‑68 (MQ=255)
               aCTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTtttctttc                                                 >  1:16561886/1‑68 (MQ=255)
               aCTGGCTGATTTCTGTTGCAGTTGCAGCGCGAACATAACGCCATGGGTTACTGGCATCATttagtttc                                                 >  1:21315412/1‑68 (MQ=18)
                  ggCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCa                                              >  1:24505838/1‑68 (MQ=255)
                    cTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATc                                            >  1:4434732/1‑68 (MQ=255)
                    cTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATc                                            >  1:22568046/1‑68 (MQ=255)
                     tGATTTCTGTTGCCGTTGCCGCGCGAACATCACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCg                                           <  1:19044116/68‑1 (MQ=255)
                     tGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCg                                           <  1:23632749/68‑1 (MQ=255)
                     tGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCg                                           <  1:9377947/68‑1 (MQ=255)
                      gATTTCTGTTGCCGTTGGCGCGCGAACATAACGCCTTGTGTTAGTGGCATGGTTTTCTTTCGCTTCgg                                          <  1:19948265/68‑2 (MQ=25)
                        tttCTGTTGCCGGTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCt                                        <  1:3962933/68‑1 (MQ=255)
                         ttCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTg                                       >  1:6424260/1‑68 (MQ=255)
                         ttCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTGCTTTCGCATCGccgg                                       >  1:26738460/1‑66 (MQ=255)
                          tCTGTTGCCGTTGCCGCGCGAACTTAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTgg                                      <  1:12273782/68‑1 (MQ=255)
                          tCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTgg                                      <  1:27754458/68‑1 (MQ=255)
                            tGTTGCCGTTGCCGCGCTAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGc                                    <  1:20670404/68‑1 (MQ=255)
                             gttgccgttgccGCGCGAACATAACGCCATGGGTTTCTGGCATCGTTGTATTTCGCATCGCCAGGGCa                                   >  1:764482/1‑68 (MQ=18)
                             gttgccgttgccGCGCGAACATAACGCCATGGGTTACTGGCATCATTTTATTTCGCATCCCCTGGGCa                                   >  1:10007608/1‑68 (MQ=21)
                                  cgttgccgCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTGTTTCGCATCGCCAGGCCATTTTg                              >  1:17394359/1‑68 (MQ=25)
                                  cgttgccgCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAAttt                               >  1:13169024/1‑67 (MQ=255)
                                  cgttgccgCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAAttt                               >  1:10785531/1‑67 (MQ=255)
                                      gccgCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCcgc                          >  1:28516140/1‑68 (MQ=255)
                                        cgcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCcgcg                         <  1:28936202/67‑1 (MQ=255)
                                        cgcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCcgcg                         <  1:11963678/67‑1 (MQ=255)
                                        cgcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCcgcg                         <  1:26108010/67‑1 (MQ=255)
                                        cgcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCcgcg                         <  1:22345232/67‑1 (MQ=255)
                                        cgcgcgAACATAACGCCATGGGTTACTGGCATCGTTGAGTTTCGCATCGCCAGGGCAGTTAGCagcg                         >  1:961521/1‑67 (MQ=14)
                                         gcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGcc                       <  1:8077486/68‑1 (MQ=255)
                                             gAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTAcc                   <  1:3440278/68‑1 (MQ=255)
                                               aCATAACGCCATGGGTTACTTGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACcca                 <  1:29492038/68‑1 (MQ=255)
                                               aCATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCTTCGCCTGGGCAATTTGCCGCGCCTACcca                 <  1:9793632/68‑1 (MQ=255)
                                               aCATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACcca                 <  1:7359860/68‑1 (MQ=255)
                                               aCATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACcca                 <  1:6535041/68‑1 (MQ=255)
                                               aCATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACcca                 <  1:8405765/68‑1 (MQ=255)
                                               aCATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACcca                 <  1:5967845/68‑1 (MQ=255)
                                                cATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCTCCAGGGCAATTTGAAGAGGCTACcca                 >  1:147033/1‑67 (MQ=17)
                                                     cGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGc           >  1:11767672/1‑68 (MQ=255)
                                                     cGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGc           >  1:27494505/1‑68 (MQ=255)
                                                     cGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGc           <  1:19129482/68‑1 (MQ=255)
                                                      gCCATGGGTTACTTGCATACTTTTATTACGAATCCCCAGGCCAATTTTCAGCGCCTACCCACCAGGCg          >  1:26856304/1‑68 (MQ=255)
                                                       ccATGGGTTACTGGCATCTTTTTATTTCGCATCTCCAGGGCAATTTGAAGAGGCTACCCACCAGGcat         >  1:14038184/1‑66 (MQ=2)
                                                       ccATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCAGGGCAGTTTGAAGAGCCTACCCACCAGGcat         >  1:8775246/1‑66 (MQ=17)
                                                        cATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGt         <  1:15919163/67‑1 (MQ=255)
                                                        cATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGt         <  1:26443359/67‑1 (MQ=255)
                                                        cATGGGTTACTGGCATCGTTTTCTATCGCACCGCCTGGGCAATTTGACGCGCCTACCCACCAGGCGtt        <  1:20819069/68‑1 (MQ=38)
                                                        cATGGGTTACTGGATTATTAATCTTTCGCATCGCCTGGGCAATTTTACGAGCCTACCCACCAGGCGtt        <  1:14414928/68‑1 (MQ=255)
                                                          tGGGTTACTGGCTTAGTTTTCTATCGCACCGCCTGGGCAATTTTACGAGCCTACCCACCAGGCGTTCt      <  1:25824748/68‑1 (MQ=25)
                                                          tGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt      <  1:26164055/68‑1 (MQ=255)
                                                          tGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt      <  1:5354971/68‑1 (MQ=255)
                                                          tGGGTTACTGGCATCGGTTTCTTTCGCAGCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt      <  1:4808459/68‑1 (MQ=255)
                                                          tGGGTTACTGGACTAGTCTTCTATCGCAGCGCCTGGGCAATATTCCGCGCCTACCCACCAGGCGTTCt      <  1:24966193/68‑1 (MQ=255)
                                                          tGGGTTAATGGAATCGTGTTCTATCGCATCGCGTGGGCAATTTTTCGCGCCTACCCACCAGGCGTTCt      <  1:13529821/68‑1 (MQ=25)
                                                           gggTTACTGGCATCGTTTTCTTTCGCATCGCCGGGGCAATTTGAAGCCGCTACCCACCAGGCATtttt     >  1:28240160/1‑68 (MQ=18)
                                                            ggTTACTGGCATCGTTTTCTTTCGCATCGCCTggg                                     >  1:21014811/1‑35 (MQ=255)
                                                            ggTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTta    >  1:18522684/1‑68 (MQ=255)
                                                            ggTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTta    >  1:12574418/1‑68 (MQ=255)
                                                             gTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTta    <  1:15238721/67‑1 (MQ=255)
                                                             gTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTta    <  1:6218236/67‑1 (MQ=255)
                                                             gTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTta    <  1:15091272/67‑1 (MQ=255)
                                                             gTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTta    <  1:18852931/67‑1 (MQ=255)
                                                                |                                                                 
AGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATA  >  CP009273/3462336‑3462465

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