Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,462,439 C→A 5.8% A258A (GCG→GCT chiA ← periplasmic endochitinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,462,4390CA 5.8% 2516.3 / 30.7 801A258A (GCG→GCTchiAperiplasmic endochitinase
Reads supporting (aligned to +/- strand):  ref base C (159/581);  new base A (31/29);  total (190/611)
Fisher's exact test for biased strand distribution p-value = 9.89e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 862 total aligned reads displayed.

TGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCTT  >  CP009273/3462373‑3462506
                                                                  |                                                                   
   cgcgcgAACATACCGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCcgcg                                                                  <  1:39812829/67‑1 (MQ=255)
   cgcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTGTTTCGCATCACCTGGGCAATTAGCagcg                                                                  >  1:42287845/1‑67 (MQ=21)
   cgcgcgAACAAAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCcgcg                                                                  <  1:37714444/67‑1 (MQ=255)
    gcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGcc                                                                <  1:38273580/68‑1 (MQ=255)
      gcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCt                                                               <  1:10998518/67‑1 (MQ=255)
      gcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCt                                                               <  1:24946520/67‑1 (MQ=255)
        gAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTAcc                                                            <  1:13616027/68‑1 (MQ=255)
        gAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTAc                                                             <  1:30588747/67‑1 (MQ=255)
          aCATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACcca                                                          <  1:22897059/68‑1 (MQ=255)
          aCATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACcca                                                          <  1:4687726/68‑1 (MQ=255)
          aCATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACcca                                                          <  1:32930026/68‑1 (MQ=255)
            aTAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACccacc                                                        >  1:20201797/1‑68 (MQ=255)
                cGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGc                                                    >  1:34837270/1‑68 (MQ=255)
                cGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGc                                                    <  1:41703316/68‑1 (MQ=255)
                 gCCATGGTTTACTGGCATCGTTTGATTTCGCATCGCCAGGGCAATTAGAAGCGCCTACCCACCAGGca                                                   >  1:10301958/1‑67 (MQ=17)
                 gCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCg                                                   <  1:3896951/68‑1 (MQ=255)
                 gCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCg                                                   >  1:9555903/1‑68 (MQ=255)
                  ccATGGGTTACTGGCATCGTTTTGTTTCGCATCGCCAGGGCAGTTTGAAGAGGCTACCCACCAGGcat                                                  >  1:8350868/1‑66 (MQ=11)
                  ccATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGt                                                  >  1:27677378/1‑68 (MQ=255)
                  ccATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGt                                                  >  1:6308565/1‑68 (MQ=255)
                   cATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGtt                                                 <  1:42687189/68‑1 (MQ=255)
                   cATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGtt                                                 <  1:22621148/68‑1 (MQ=255)
                    aTGGGTTATTGGCTTTTTGTTCTTTGGCACCGCCTGGGCAATTTTACGAGCCTACCCACCAGGCGTTc                                                <  1:7102032/68‑1 (MQ=16)
                    aTGGGTTACTGGCATCGTTTTATTTCGCATCGCCAGGGCAGTTTGAAGCGCCTACCCACCAGGCAttt                                                >  1:14417410/1‑67 (MQ=17)
                     tGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt                                               <  1:17145457/68‑1 (MQ=255)
                     tGGGTTACTGGCATCGTTTTCTATCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt                                               <  1:11383011/68‑1 (MQ=255)
                      gggTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCATtttt                                              >  1:40730244/1‑68 (MQ=255)
                      gggTTACTGGCATCGTTTTCTTTCGCAGCCCCAGGGCAATTTGCAGCGCCTACCCACCAGGCATtttt                                              >  1:18538208/1‑68 (MQ=21)
                       ggTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTta                                             >  1:30919231/1‑68 (MQ=255)
                       ggTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGAAGAGCATACCCACCAGGCTttatta                                             >  1:15408329/1‑68 (MQ=255)
                       ggTTACTGGCATCGTTTTCTTTCGCATCGCCAGGGCAGTTTGAAGCGCCTACCCACCAGGCATTTTta                                             >  1:7357036/1‑68 (MQ=18)
                         ttACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTtata                                           >  1:20951671/1‑68 (MQ=255)
                         ttACTGGCATCGTTTTCTTTCGCATCGCCAGGGCAATTTGCAGCGCCTACCCACCAGGCGTTTTtata                                           >  1:26839461/1‑68 (MQ=37)
                             tGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAAt                                        >  1:28130018/1‑67 (MQ=255)
                              ggCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCTTTATTATaaaac                                       >  1:40153391/1‑65 (MQ=255)
                              ggCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTATTATaaaag                                       >  1:14210224/1‑65 (MQ=255)
                               gCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATgcc                                     >  1:3000587/1‑68 (MQ=255)
                               gCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTATTatatatcac                                     >  1:40331108/1‑65 (MQ=38)
                                    tttttCTTTCGCACCCCCTGGGCAATTTGACGAGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:29248790/67‑1 (MQ=37)
                                    tttttCTATGGGAACGCGTGGGCAATTTTCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:23459814/67‑1 (MQ=25)
                                    ttcttcgttcgTATCGCCTGGGCAATTTTACGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:19350995/65‑1 (MQ=37)
                                    gtcttctttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:36877048/65‑1 (MQ=255)
                                    gtcttctatcGCATCGCCTGGGCAATTTTCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:4136646/65‑1 (MQ=38)
                                    gTTTTCTTTCGCCTCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:33617865/68‑1 (MQ=255)
                                    gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:24644605/68‑1 (MQ=255)
                                    gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:29782519/68‑1 (MQ=255)
                                    gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:13444782/68‑1 (MQ=255)
                                    gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:36471143/68‑1 (MQ=255)
                                    gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:18055683/68‑1 (MQ=255)
                                    gTTTTCTATCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:2301526/68‑1 (MQ=255)
                                    gTTTGCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:5084866/68‑1 (MQ=255)
                                    gTTTGCATTCGCATCTCCTGGGCAATTTGGAGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:1796863/68‑1 (MQ=37)
                                    gTTATCTATCGTAGCGCCTGGGCAATTTTTCGAGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:19913832/68‑1 (MQ=21)
                                      tttctatcGCATCGCGTGGGCAATATTACGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:31275907/66‑1 (MQ=37)
                                       ttctttcGCATCGCCTGTGGAATTTGCCGCGCCTACACACCTGGCGTTCTTATAAATGCCTCCGTtat                             <  1:16732319/68‑1 (MQ=255)
                                       ttctttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTtat                             <  1:35770741/68‑1 (MQ=255)
                                       ttctttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTtat                             <  1:8731936/68‑1 (MQ=255)
                                       ttctttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTtat                             <  1:32479602/68‑1 (MQ=255)
                                        tctttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTtata                            <  1:35430186/68‑1 (MQ=255)
                                           ttcGTATCGCCTGGGCAATTTGGCGAGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         <  1:21721994/68‑1 (MQ=38)
                                           ttcGCATCGCCTGTTCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         <  1:27234863/68‑1 (MQ=255)
                                           ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         <  1:12616018/68‑1 (MQ=255)
                                           ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         <  1:28933565/68‑1 (MQ=255)
                                           ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         <  1:17760606/68‑1 (MQ=255)
                                           ttcGCATCGAGTGGGCAATGTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         <  1:13801841/68‑1 (MQ=39)
                                           ttcGAATCGCCTGGGCAATATGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         <  1:26008221/68‑1 (MQ=255)
                                              gCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTGTTATAAATGGCGCCGTTATAGATaa                       >  1:20504227/1‑67 (MQ=255)
                                              gCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAAc                      <  1:21157368/68‑1 (MQ=255)
                                                aTCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                    <  1:12854404/68‑1 (MQ=255)
                                                aTCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                    <  1:14890645/68‑1 (MQ=255)
                                                aTCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                    <  1:770509/68‑1 (MQ=255)
                                                 tCGCCTGGGCAATTTTCAGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                    <  1:24389883/67‑1 (MQ=255)
                                                   gcgtGGGCAATTTTCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTAcc                 <  1:28661599/65‑1 (MQ=255)
                                                   cccTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTAcc                 <  1:38530075/67‑1 (MQ=255)
                                                    ccTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCa                >  1:44139/1‑68 (MQ=255)
                                                     gtGGGCAATTTTCAGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:36108478/67‑1 (MQ=38)
                                                     cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:2401489/68‑1 (MQ=255)
                                                     cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:21970698/68‑1 (MQ=255)
                                                     cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:33262776/68‑1 (MQ=255)
                                                     cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:2248343/68‑1 (MQ=255)
                                                       gggCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGt             <  1:7992303/68‑1 (MQ=255)
                                                       gggCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGt             <  1:17467133/68‑1 (MQ=255)
                                                         gCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATTACGTACCAGGTag           >  1:1000192/1‑68 (MQ=255)
                                                         gCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTag           >  1:34224284/1‑68 (MQ=255)
                                                         gCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTa            <  1:25456332/67‑1 (MQ=255)
                                                         gCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTa            <  1:6709990/67‑1 (MQ=255)
                                                         gCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTa            <  1:5519838/67‑1 (MQ=255)
                                                         gCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTa            <  1:37554396/67‑1 (MQ=255)
                                                         gCAATTTGCCGCGCCTACCCACCAGGCATTATGATAAATGAAGCCGTGAAAGATAACGTACGATTTat           >  1:22197000/1‑67 (MQ=255)
                                                          cAATTTTCAGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTag           <  1:39350352/67‑1 (MQ=255)
                                                          cAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTaga          >  1:37204710/1‑68 (MQ=255)
                                                          cAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTag           <  1:16224118/67‑1 (MQ=255)
                                                          cAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTag           <  1:32150624/67‑1 (MQ=255)
                                                          cAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTag           <  1:35011397/67‑1 (MQ=255)
                                                           aaTTTGCCGCGCCTACCCTCCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTaga          <  1:11854855/67‑1 (MQ=255)
                                                           aaTTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTagaa         >  1:42994651/1‑68 (MQ=255)
                                                               tGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAAc     <  1:20801344/68‑1 (MQ=255)
                                                                tccGCGCCTACCCACCATGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAAcc    <  1:31691013/67‑1 (MQ=255)
                                                                  cgcgCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCtt  <  1:34486270/68‑1 (MQ=255)
                                                                  cgcgCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCt   <  1:9268003/67‑1 (MQ=255)
                                                                  |                                                                   
TGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCTT  >  CP009273/3462373‑3462506

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: