Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,462,439 C→A 8.6% 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 8.6% 1909.1 / 33.3 635A258A (GCG→GCTchiAperiplasmic endochitinase
Reads supporting (aligned to +/- strand):  ref base C (133/433);  new base A (39/29);  total (172/463)
Fisher's exact test for biased strand distribution p-value = 2.63e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 680 total aligned reads displayed.

GCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCTT  >  CP009273/3462374‑3462506
                                                                 |                                                                   
  cgcgcgAACATTAGGCCAGGGGTTACTGGCATCGGTTTCTTTCGCATCGCCTGGGCAATTTGCcgcg                                                                  <  1:13207631/67‑1 (MQ=255)
  cgcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCcgcg                                                                  <  1:27113237/67‑1 (MQ=255)
  cgcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCcgcg                                                                  <  1:6566147/67‑1 (MQ=255)
  cgcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCcgcg                                                                  <  1:8031879/67‑1 (MQ=255)
     gcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCt                                                               <  1:11048134/67‑1 (MQ=255)
       gAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCCTCGCCTGGGCAATTTGCCGCGCCTAcc                                                            <  1:14565325/68‑1 (MQ=255)
         aCATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACcca                                                          <  1:29971915/68‑1 (MQ=255)
            tAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACccacca                                                       <  1:28684018/68‑1 (MQ=255)
               cGCCATGGGTTACTGGCATCTTTTGGTTTCGCATCCCCAGGGCAGTTTGAAGCGCCTACCCACCAGGc                                                    >  1:25761197/1‑68 (MQ=14)
               cGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGc                                                    >  1:19979566/1‑68 (MQ=255)
                gCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCg                                                   >  1:29644569/1‑68 (MQ=255)
                gCCATGGGTTACTGGCATCGTTTAGTTTCGCATCGCCAGGACAATTAGAAGAGCCTACCCACCAGGca                                                   >  1:8457003/1‑67 (MQ=2)
                 ccATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGt                                                  >  1:20756407/1‑68 (MQ=255)
                 ccATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGt                                                  >  1:8226326/1‑68 (MQ=255)
                   aTGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGt                                                  >  1:32618263/1‑66 (MQ=255)
                   aTGGGTTACTGGCATCGTTTTCTTTCGCATCCCCAGGCCAATTTCACGCGCCTACCCAACAGGCttta                                                >  1:21061537/1‑67 (MQ=255)
                    ttggTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt                                               <  1:4088899/66‑1 (MQ=255)
                    tGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt                                               <  1:1774426/68‑1 (MQ=255)
                      ggTTACTGGCATCGTTTTCTTTCGCATCGCCAGGGCAATTTGCAGCGCCTACCCACCAGGCTTTTTta                                             >  1:12709119/1‑68 (MQ=37)
                      ggTTACTGGCATCGTTTTCTTTCGCATAGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTta                                             >  1:10025706/1‑68 (MQ=255)
                       gTTACTGGCATCGTTTTCTTTCGCATCGCCAGGGCAATTTGCAGCGGCTACCCACCAGGCATTTTtat                                            >  1:12056178/1‑68 (MQ=21)
                            tGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAAt                                        >  1:19641002/1‑67 (MQ=255)
                            tGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAAt                                        >  1:30589398/1‑67 (MQ=255)
                            tGGCATCGTTTTCTTTCGCATCGCCTGGGCAACTTGCCGCGCCTACCCACCAGGCGTTATTAaaaaa                                        >  1:18869099/1‑66 (MQ=255)
                             ggCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTGTGATAAATga                                      >  1:15449138/1‑67 (MQ=255)
                             ggCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATg                                       >  1:28999922/1‑67 (MQ=255)
                             ggCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATg                                       >  1:7565019/1‑67 (MQ=255)
                                 tCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATgccgcc                                  <  1:9494935/68‑1 (MQ=255)
                                 tCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATgccgcc                                  <  1:33462504/68‑1 (MQ=255)
                                 tCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATgccgcc                                  <  1:21452925/68‑1 (MQ=255)
                                 tCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATgccgcc                                  <  1:22811422/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:24149501/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:30795799/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:6202219/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:11623807/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:29254633/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:1142627/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:1031439/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:22335864/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:12395397/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:1289032/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:20252975/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:33806666/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:7158312/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:5057406/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:5800319/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:16340146/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:15733293/68‑1 (MQ=255)
                                   gTTTTCTTTCGAATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:8850320/68‑1 (MQ=255)
                                   gTTTTCTCTGGCACCGCCTGGGCAATTTTTCGCGCCAACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:34302818/68‑1 (MQ=21)
                                   gTTTTCTATCGCAGCCCGTGGGCAATTTTACGCGCCTACCCACCCGGCGTTCTTATAAATGCCGCCGt                                <  1:26601627/68‑1 (MQ=25)
                                    ttttCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt                                               <  1:339878/52‑1 (MQ=255)
                                    ttatctttcGCATCGCCTGGGCAATTTGCAGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:21999696/64‑1 (MQ=255)
                                    tcttctttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:28488975/65‑1 (MQ=255)
                                      ttctatcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTtat                             <  1:28199852/68‑1 (MQ=255)
                                      ttctatcGCATCCCGTGGGCAATTTGTAGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTtat                             <  1:5956808/68‑1 (MQ=25)
                                        ctttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAg                           <  1:7767307/68‑1 (MQ=255)
                                        ctttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAg                           <  1:4384979/68‑1 (MQ=255)
                                          ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         <  1:27390707/68‑1 (MQ=255)
                                          ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         <  1:9815106/68‑1 (MQ=255)
                                          ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         <  1:117995/68‑1 (MQ=255)
                                          ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         <  1:13608316/68‑1 (MQ=255)
                                          ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         >  1:6903075/1‑68 (MQ=255)
                                          ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         <  1:4721481/68‑1 (MQ=255)
                                          ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         >  1:32212007/1‑68 (MQ=255)
                                          ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         <  1:31185408/68‑1 (MQ=255)
                                          ttcGCAACGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         <  1:27913084/68‑1 (MQ=255)
                                          atggCATCGCCTGGGCAATATTCCGCGCCTACCCACCACGCGTTCTTATAAAAGCCGCCGTTATAGAt                         <  1:27698789/65‑1 (MQ=255)
                                            cGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATaa                       <  1:26271266/68‑1 (MQ=255)
                                             gCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAAc                      <  1:33180728/68‑1 (MQ=255)
                                               aTCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                    <  1:10651248/68‑1 (MQ=255)
                                               aTCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTATTATAAATTGAGCCGTTTTAAATAACg                     >  1:14951583/1‑67 (MQ=255)
                                                  gCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTAcc                 <  1:9270800/68‑1 (MQ=255)
                                                  gCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTAcc                 <  1:17031932/68‑1 (MQ=255)
                                                    cTGGGCAATTTTACGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:15957060/68‑1 (MQ=255)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:31050750/68‑1 (MQ=255)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:31641056/68‑1 (MQ=255)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:19237778/68‑1 (MQ=255)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:1403590/68‑1 (MQ=255)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:17442046/68‑1 (MQ=255)
                                                    cTGGGCAATATGACGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:14291501/68‑1 (MQ=255)
                                                      gggCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGt             <  1:31355495/68‑1 (MQ=255)
                                                        gCAATTTTAAGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTa            <  1:16794811/67‑1 (MQ=38)
                                                        gCAATTTGCCGCGCCTACCCACCAGGCGTTGTTATAAAGGGCGCCGTTATAGAGGACGTACCATGTat           >  1:18270012/1‑67 (MQ=37)
                                                        gCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGTCGTTATAGATAACGTACCATGTag           >  1:23192675/1‑68 (MQ=255)
                                                        gCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTa            <  1:24807348/67‑1 (MQ=255)
                                                        gCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTa            <  1:25232398/67‑1 (MQ=255)
                                                        gCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTa            <  1:32067395/67‑1 (MQ=255)
                                                        gCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTa            <  1:24563896/67‑1 (MQ=255)
                                                        gCAATTTGCCGCGCCTACCCACAAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTag           >  1:23886813/1‑68 (MQ=255)
                                                         cAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTaga          >  1:31865732/1‑68 (MQ=255)
                                                             ttGCCGCGCCTACCCACCAGGCGTTGTTATAAATGCCGCCGTTATAGATAACGTACCATGTagaaga       >  1:22593224/1‑67 (MQ=255)
                                                              tgacgcgCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAAc     <  1:9038029/65‑1 (MQ=255)
                                                              tGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGGTATAGATAACGTACCATGTAgaagaa      >  1:13870198/1‑67 (MQ=255)
                                                               gCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAgaagaa      >  1:3512142/1‑66 (MQ=255)
                                                                 cgcgCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCtt  <  1:25505971/68‑1 (MQ=255)
                                                                 cgcgCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCtt  <  1:23500391/68‑1 (MQ=255)
                                                                 cgcgCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCtt  <  1:30232600/68‑1 (MQ=255)
                                                                 cgcgCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCtt  <  1:14457035/68‑1 (MQ=255)
                                                                 cgcgCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCtt  <  1:5373830/68‑1 (MQ=255)
                                                                 |                                                                   
GCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCTT  >  CP009273/3462374‑3462506

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: