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% 1416.1 / 23.6 455S271S (AGT→AGCchiAperiplasmic endochitinase
Reads supporting (aligned to +/- strand):  ref base A (240/179);  new base G (19/6);  total (270/185)
Fisher's exact test for biased strand distribution p-value = 9.31e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 459 total aligned reads displayed.

AGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAA  >  CP009273/3462336‑3462467
                                                                |                                                                   
aGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTa                                                                     <  1:36130559/65‑1 (MQ=255)
aGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTGCTg                                                                  >  1:10187359/1‑68 (MQ=255)
aGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTg                                                                  >  1:27556388/1‑68 (MQ=255)
aGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTg                                                                  >  1:556394/1‑68 (MQ=255)
aGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTg                                                                  >  1:25107684/1‑68 (MQ=255)
aGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTg                                                                  >  1:33778623/1‑68 (MQ=255)
aGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTg                                                                  >  1:13995286/1‑68 (MQ=255)
   cAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTgg                                                                 <  1:31726045/66‑1 (MQ=255)
     gggTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATc                                                             >  1:7181589/1‑68 (MQ=255)
      ggTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCg                                                            >  1:26472025/1‑68 (MQ=255)
       gTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTCCTGGCATcat                                                           >  1:29945842/1‑66 (MQ=255)
        ttACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGt                                                           <  1:14244381/67‑1 (MQ=255)
        ttACCAAACTGACTGAATTCTGTTGCCGTTGCCGCGCGCACATAACGCCATGGGTTACTGGCATCGt                                                           <  1:35925156/67‑1 (MQ=25)
          aCCATACTGGCTGATTTCTGTTGCCGTTGCCGCTCGAACATAACGCCATGGGTTGCTGGCATCAttgc                                                        <  1:18297425/68‑3 (MQ=17)
           ccATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTGCTGGCATCGtttga                                                       >  1:16366969/1‑66 (MQ=255)
           ccATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGtttgg                                                       >  1:29621296/1‑66 (MQ=255)
            cATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACTTAACGCCATGGGTTACTGGCATCGTtttct                                                      <  1:21511331/68‑1 (MQ=255)
            cATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTGCt                                                                   >  1:22232685/1‑55 (MQ=255)
            cATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTGCTGGCATCAtttatt                                                      >  1:23941952/1‑65 (MQ=18)
            cATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGttttgt                                                      >  1:4018660/1‑66 (MQ=255)
            cATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGttttat                                                      >  1:10724997/1‑66 (MQ=255)
            cATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTtttct                                                      <  1:4844261/68‑1 (MQ=255)
            cATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTtttct                                                      >  1:33102996/1‑68 (MQ=255)
            cATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTtttct                                                      >  1:21815700/1‑68 (MQ=255)
            cATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTtttct                                                      >  1:3765092/1‑68 (MQ=255)
            cATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTtttct                                                      >  1:17792323/1‑68 (MQ=255)
            cATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTtttct                                                      <  1:14688904/68‑1 (MQ=255)
            cATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCAtttatt                                                      >  1:24407210/1‑65 (MQ=255)
            cATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCAtttagt                                                      >  1:13441591/1‑65 (MQ=255)
            cATACTGGCTGATTTCTGTTGCCGTAGCCGCTCGAACATAACGCCATGGGTTGCTGGc                                                                <  1:23620157/58‑1 (MQ=17)
              tACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTtttgtt                                                     >  1:26945244/1‑67 (MQ=255)
              tACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTtttctt                                                     >  1:12013521/1‑67 (MQ=255)
              tACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTtttctt                                                     >  1:15132780/1‑67 (MQ=255)
              tACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTtttctt                                                     >  1:22988354/1‑67 (MQ=255)
              tACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTtttctt                                                     >  1:31360542/1‑67 (MQ=255)
               aCTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTGCTGGCATCATtttgtttc                                                   >  1:17929904/1‑68 (MQ=18)
               aCTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTtttctttc                                                   >  1:924950/1‑68 (MQ=255)
               aCTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCATttattttc                                                   >  1:28429600/1‑68 (MQ=18)
               aCTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCAAGGGTTACTGGCATCGTtttctttc                                                   >  1:11145895/1‑68 (MQ=255)
                cTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATCACGCCATGGGTTACTGGCATCGTTTTCTTTCg                                                  <  1:15612995/68‑1 (MQ=255)
                cTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCg                                                  <  1:15808325/68‑1 (MQ=255)
                  ggCTGATTTCTGTTGCAGTTGCCGCGCGAACATATCGCCATGGGTTACTGGCATCTTTTCGTTTCGCa                                                >  1:36693627/1‑68 (MQ=255)
                    cTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCTTTTGCTTTCGctta                                              >  1:27234559/1‑65 (MQ=21)
                     tGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATc                                              <  1:36495436/67‑1 (MQ=255)
                     tGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCg                                             <  1:32068928/68‑1 (MQ=255)
                     tGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCg                                             <  1:12852694/68‑1 (MQ=255)
                     tGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCg                                             <  1:20939742/68‑1 (MQ=255)
                     tGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCg                                             <  1:20488638/68‑1 (MQ=255)
                     tGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGCTTACTGGCATCGTTTTCTTTCGCATCg                                             <  1:30129355/68‑1 (MQ=255)
                         ttCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTGCTTTCGCATCTccgg                                         >  1:19450996/1‑66 (MQ=21)
                            tGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGc                                      <  1:34504725/68‑1 (MQ=255)
                            tGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGc                                      <  1:8713899/68‑1 (MQ=255)
                             gttgccgttgccGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCa                                     >  1:32722504/1‑68 (MQ=255)
                             gttgccgttgccGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCa                                     >  1:29257014/1‑68 (MQ=255)
                             gttgccgttgccGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCGGGGCa                                     >  1:3296620/1‑68 (MQ=255)
                               agccgtcgccgCTCGACCATAACGCCATGGGGTACTGGCATCGGTTTCTTACGCATCGCCTGGGCAAt                                   <  1:13817896/67‑1 (MQ=25)
                                  cgttgccgCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCGGGCCAGttt                                 >  1:33370271/1‑67 (MQ=25)
                                      gccgCGCGAACATAACGCCTTGGGGTTATGGCATGCTAGTCGATCGCAACGCGTGGGCGATTTGCcgc                            >  1:2272906/1‑68 (MQ=255)
                                      gccgCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCcgc                            <  1:2607973/68‑1 (MQ=255)
                                      gccgCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCcgc                            <  1:17064445/68‑1 (MQ=255)
                                      gccgCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCcgc                            >  1:25823855/1‑68 (MQ=255)
                                        cgcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCcgcg                           <  1:21344599/67‑1 (MQ=255)
                                        cgcgcgAACATAACGCCATGGGTTACTGGCAGCGTTTTCTTTCGCATCGCCTGGGCATTTTGCcgcga                          >  1:32321025/1‑67 (MQ=38)
                                         gctcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGcc                         <  1:5824018/65‑1 (MQ=255)
                                         gcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGcc                         <  1:9432161/68‑1 (MQ=255)
                                         gcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGcc                         <  1:12188939/68‑1 (MQ=255)
                                         gcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGcc                         <  1:33641963/68‑1 (MQ=255)
                                         gcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGcc                         <  1:4235677/68‑1 (MQ=255)
                                           gcgAACTTCACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCt                        <  1:35808537/67‑1 (MQ=255)
                                           gcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTa                       <  1:16590580/68‑1 (MQ=255)
                                             gAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTAcc                     <  1:7838223/68‑1 (MQ=255)
                                             gAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTAcc                     <  1:30735574/68‑1 (MQ=255)
                                             gAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTAcc                     <  1:7507576/68‑1 (MQ=255)
                                               aCATAACGCGAGGGGTTACCGGCAGCGTTTTCTTTCGCATCGCCTGGGCAATTTGACGCGCCTACcca                   <  1:3704717/68‑1 (MQ=255)
                                                     cGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAgg              >  1:6702083/1‑67 (MQ=255)
                                                     cGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGc             >  1:25518158/1‑68 (MQ=255)
                                                     cGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGc             <  1:21076343/68‑1 (MQ=255)
                                                     cGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGACGCGCCTACCCACCAGGc             >  1:23228213/1‑68 (MQ=255)
                                                       ccATGGGTTACTGGCATCTTTTTGTTTCGCATCGCCAGGGCAGTTAGAAGAGCCTACCCACCAGGcat           >  1:19738623/1‑66 (MQ=2)
                                                       ccATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGt           >  1:5159550/1‑68 (MQ=255)
                                                       ccATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGt           >  1:7308619/1‑68 (MQ=255)
                                                       ccATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGt           >  1:9680845/1‑68 (MQ=255)
                                                       ccATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGAAGCCGCTACCCACCAGGcat           >  1:6257018/1‑66 (MQ=21)
                                                       ccATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCGGGGCAGTTTGCCGCGCCTACCCACCAGGcat           >  1:30825514/1‑66 (MQ=37)
                                                       ccATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCAGGGCAATTTGCCGCGGCTACCCACCAGGctt           >  1:8615341/1‑66 (MQ=37)
                                                       ccATGGGTTACTGGCATCGTTTACTTTCGCATCGCCAGGCCAATTTGAAGCGCCTACCCACCAGGcat           >  1:17397708/1‑66 (MQ=18)
                                                        cATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGtt          <  1:28041466/68‑1 (MQ=255)
                                                        cATGGGTTACTGGCATCGTTTTCTTTCGCATAGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGtt          >  1:24862091/1‑68 (MQ=255)
                                                        cATGGGTTACTGGACTCGTCTGCTATCGCATCGCCTGGGCAATTTGACGCGCCTACCCACCAGGCGtt          <  1:8332080/68‑1 (MQ=255)
                                                        cATGGGTTACTGGAATCTTTTTCTTTCGCATCGCCTGGGCAATTTTCCGCGCCTACCCACCAGGCGtt          <  1:18918890/68‑1 (MQ=255)
                                                        cATGGGTTAATGGACTGTTTATCTCTCGCATCGCCTGGGCAATTTGCCGCGCCAACCCACCAGGCGt           <  1:30622625/67‑1 (MQ=255)
                                                         aTGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCAGCGCCTACCCACCAGGCCttt         >  1:11404475/1‑67 (MQ=38)
                                                          tGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt        <  1:14827793/68‑1 (MQ=255)
                                                          cgggTTACTGGCAGCGTTTTCTTTCGCATCTCGTGCGCAATTTTCAGCGCCAACCCACCAGGCGTTCt        <  1:19975427/67‑1 (MQ=21)
                                                            ggTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTta      >  1:28695939/1‑68 (MQ=255)
                                                            ggTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTta      >  1:34765705/1‑68 (MQ=255)
                                                            ggTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTta      >  1:13249729/1‑68 (MQ=255)
                                                            ggTTACTGGCATCGTTTTCTTTCGCATCGCCGGGGCAATTTGCCGCGCCTACCCACCAGGCGTTTTta      >  1:28901206/1‑68 (MQ=255)
                                                             gTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTta      <  1:30427847/67‑1 (MQ=255)
                                                             gTTACTGGCATCGTTATCTATCGCATCGCCTGGGCAATATGCCGCGCCTACCCACCAGGCGTTCTta      <  1:28567644/67‑1 (MQ=255)
                                                                |                                                                   
AGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAA  >  CP009273/3462336‑3462467

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