Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,462,438 C→A 10.3% A259S (GCA→TCA)  chiA ← periplasmic endochitinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,462,4380CA10.3% 2080.8 / 38.1 741A259S (GCA→TCA) chiAperiplasmic endochitinase
Reads supporting (aligned to +/- strand):  ref base C (165/476);  new base A (28/43);  total (194/547)
Fisher's exact test for biased strand distribution p-value = 1.68e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 742 total aligned reads displayed.

GCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCT  >  CP009273/3462374‑3462505
                                                                |                                                                   
   gcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGcc                                                               <  1:6289195/68‑1 (MQ=255)
   gcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGcc                                                               <  1:12620893/68‑1 (MQ=255)
   gcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGcc                                                               <  1:29434788/68‑1 (MQ=255)
     gcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTa                                                             <  1:16157979/68‑1 (MQ=255)
       gAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTAcc                                                           <  1:19579120/68‑1 (MQ=255)
         aCATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACcca                                                         <  1:25840486/68‑1 (MQ=255)
           aTAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACccacc                                                       >  1:26088294/1‑68 (MQ=255)
           aTAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACccacc                                                       >  1:14047147/1‑68 (MQ=255)
                gCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCg                                                  >  1:10845072/1‑68 (MQ=255)
                gCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCg                                                  >  1:11589575/1‑68 (MQ=255)
                gCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTAGCCGCGCCTACCCCCCAGGct                                                  >  1:21971488/1‑67 (MQ=255)
                 ccATGGGTTACTGGCATCGTTTTGTTTCGCATCGCCCGGTCATTTAGAAGAGCCTACCCACCAGGcat                                                 >  1:38743089/1‑66 (MQ=2)
                 ccATGGGTTACTGGCATCGTTTTCTTTCGCATCTCCAGGGCAATTAGAAGCGGCTACCCACCAGGcat                                                 >  1:28991020/1‑66 (MQ=17)
                 ccATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGt                                                 >  1:11075275/1‑68 (MQ=255)
                  cATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGtt                                                <  1:13539071/68‑1 (MQ=255)
                  cATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGtt                                                <  1:2273586/68‑1 (MQ=255)
                  cATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGt                                                 <  1:31397661/67‑1 (MQ=255)
                  cATGGGTTACTGGCATCGTTTTCTATCGCATCGCGTGGGCAATTTGACGCGCCTACCCACCAGGCGtt                                                <  1:30597060/68‑1 (MQ=38)
                    tGGGTTACTGGCATCTTTATCTTTGGCAACCCCTGGGCAATCTGCCGCGCCTACCCACCAGGCGTTCt                                              <  1:28253677/68‑1 (MQ=255)
                    tGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt                                              <  1:5307691/68‑1 (MQ=255)
                    tGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt                                              <  1:25143376/68‑1 (MQ=255)
                    tGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt                                              <  1:19111393/68‑1 (MQ=255)
                    tGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt                                              <  1:32531435/68‑1 (MQ=255)
                            tGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTTTTATATATg                                      >  1:12165164/1‑68 (MQ=255)
                            tGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAAt                                       >  1:28547016/1‑67 (MQ=255)
                            tGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAAt                                       >  1:31065054/1‑67 (MQ=255)
                            tGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAAt                                       >  1:8654197/1‑67 (MQ=255)
                            tGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTATTATAAAt                                       >  1:2664149/1‑67 (MQ=255)
                            tGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCAGCGCCTACCCACCAGGCATTTTtatataa                                       >  1:8351611/1‑64 (MQ=25)
                            tGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGAAGAGCATACCCACCAGGCGTTTTTATAAAtt                                      >  1:35083734/1‑67 (MQ=255)
                            tGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTAGCCGCGCATACACACCAGCCTTTCTTATAAAt                                       >  1:31755053/1‑67 (MQ=255)
                             ggCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATg                                      >  1:7836242/1‑67 (MQ=255)
                             ggCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATg                                      >  1:11286761/1‑67 (MQ=255)
                             ggCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCATTATTATaaaag                                      >  1:33864452/1‑65 (MQ=38)
                                aTCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATgccg                                   <  1:1871843/67‑1 (MQ=255)
                                 tCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATgccgcc                                 <  1:2094337/68‑1 (MQ=255)
                                 tCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATgccgcc                                 <  1:22316938/68‑1 (MQ=255)
                                 tCGTGTTCTTTCGCATCGCGTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATgccgcc                                 <  1:10615977/68‑1 (MQ=255)
                                  cGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGccgccg                                >  1:15293469/1‑68 (MQ=255)
                                   tttttCTTTCGCTTCGCCTGGTCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:36148019/67‑1 (MQ=39)
                                   tttttCTTTCGCATCGCGTGGGCAATTTGACGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:3833574/67‑1 (MQ=39)
                                   tttttCTATCGCAGCGCCTGGGCAATTTGTCGAGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:34737681/67‑1 (MQ=37)
                                   gtgttctttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:33361830/65‑1 (MQ=255)
                                   gtattctctcGAATCGCCTGGGCAATGTTTAGAGCCAACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:13800479/65‑1 (MQ=9)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:3430514/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:9646193/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:394781/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:4662741/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:28724494/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:5661261/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:8094326/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:21501629/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:17918144/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:15844747/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:15406308/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:20120409/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:2070716/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:24078595/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:2476500/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:24349788/68‑1 (MQ=255)
                                   gTTTTCTATCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:35610874/68‑1 (MQ=255)
                                   cggttctttcGCATCGCCTGGGCAATTTGGCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:20439212/65‑1 (MQ=255)
                                    ttttCTTTCGCATCTAATGGGCAATTTTACGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:29820523/67‑1 (MQ=37)
                                    ttttCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAAAGCCACCGtt                              >  1:26425221/1‑68 (MQ=255)
                                    ttttCTTTCGCACCGCGTGGGCAATTTGACGCGCCTACCCACCAGGAGTTCTTATAAATGCCGCCGt                               <  1:27426094/67‑1 (MQ=255)
                                    ttttCTATGGCACCGCGTGGGCAATATTCCGAGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGtt                              <  1:4975885/68‑1 (MQ=21)
                                     tatctttgGCACCGCCTTGGGAATTTGCCGACCCTACCCACCAGGCGTTCTTATATATGCCGCCGt                               <  1:36389048/64‑1 (MQ=21)
                                      ttctttgGCATCGCGTGGGCAATATTAAGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTtat                            <  1:19884689/68‑1 (MQ=25)
                                      ttctttcGCATCCCCTGGGCAATTTGACGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTtat                            <  1:7331419/68‑1 (MQ=255)
                                        ctttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAg                          <  1:6568916/68‑1 (MQ=255)
                                         tttcGCATCGCCTGGGCAATTTGCCGCGCCTCCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                        <  1:23177060/68‑1 (MQ=255)
                                          ttggCATCGCCTGGGCAATTTTCAGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                        <  1:37040022/65‑1 (MQ=38)
                                          ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                        <  1:28106184/68‑1 (MQ=255)
                                          ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                        <  1:33030437/68‑1 (MQ=255)
                                          ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                        >  1:17454180/1‑68 (MQ=255)
                                          ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                        <  1:14645564/68‑1 (MQ=255)
                                          ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                        <  1:13178979/68‑1 (MQ=255)
                                          ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                        <  1:8845090/68‑1 (MQ=255)
                                          atcGCATCGCCTGGGCAATTTTACGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                        <  1:25629383/67‑1 (MQ=39)
                                          atcGCACCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                        <  1:2786022/67‑1 (MQ=255)
                                           tcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATa                       <  1:7137616/68‑1 (MQ=255)
                                               aTCGCCTGGTCAATTTGGCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                   <  1:10046675/68‑1 (MQ=255)
                                               aTCGCCTGGGCAATTTGTCGAGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                   <  1:18100009/68‑1 (MQ=255)
                                               aTCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                   <  1:30143279/68‑1 (MQ=255)
                                               aTCCCCTGGGCAATATTCCGCGCCAACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                   <  1:37550503/68‑1 (MQ=37)
                                                tCGCGTGGGCAATTTGCCGAGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                   <  1:26983723/67‑1 (MQ=255)
                                                    cTGGGCAATTTTTCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt              <  1:21155737/68‑1 (MQ=255)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt              <  1:32880106/68‑1 (MQ=255)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt              <  1:18513355/68‑1 (MQ=255)
                                                       ggCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTa           >  1:21765452/1‑68 (MQ=255)
                                                        gCAATTTTCCGAGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTa           <  1:14438397/67‑1 (MQ=255)
                                                        gCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTag          >  1:16343604/1‑68 (MQ=255)
                                                        gCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTag          >  1:4380312/1‑68 (MQ=255)
                                                        gCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTa           <  1:37890588/67‑1 (MQ=255)
                                                        gCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTa           <  1:23491859/67‑1 (MQ=255)
                                                         cAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTag          <  1:12850436/67‑1 (MQ=255)
                                                         cAATTTGCCGCGCCTACCAACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTag          <  1:30360934/67‑1 (MQ=255)
                                                          aaTTTTTCGAGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTaga         <  1:16804826/67‑1 (MQ=38)
                                                             ttGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAGGTagaaga      >  1:9325194/1‑67 (MQ=255)
                                                               gCCGCGCCTACCCCCCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAgaagaa     >  1:6852836/1‑66 (MQ=255)
                                                                |                                                                   
GCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCT  >  CP009273/3462374‑3462505

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: