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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,462,4380CA 7.0% 1952.6 / 12.3 668A259S (GCA→TCA) chiAperiplasmic endochitinase
Reads supporting (aligned to +/- strand):  ref base C (163/434);  new base A (19/26);  total (182/486)
Fisher's exact test for biased strand distribution p-value = 3.93e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 668 total aligned reads displayed.

GCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAAC  >  CP009273/3462374‑3462503
                                                                |                                                                 
  cgcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCcgcg                                                               >  1:29606939/1‑67 (MQ=255)
  cgcgcgAACATAACGCCATGGGTTACTGGCATCGTTTGGTTTCGCATCGCCCGGGCATTTTGCagcg                                                               >  1:7981703/1‑67 (MQ=21)
   gctcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGcc                                                             <  1:5824018/65‑1 (MQ=255)
   gcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGcc                                                             <  1:35205627/68‑1 (MQ=255)
   gcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGcc                                                             <  1:9432161/68‑1 (MQ=255)
     gcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCt                                                            <  1:25295996/67‑1 (MQ=255)
       gAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTAcc                                                         <  1:7507576/68‑1 (MQ=255)
           aTAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACccacc                                                     >  1:11447996/1‑68 (MQ=255)
               cGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGc                                                 >  1:31065911/1‑68 (MQ=255)
               cGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGc                                                 >  1:11136091/1‑68 (MQ=255)
                gCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCg                                                >  1:26811613/1‑68 (MQ=255)
                 ccATGGGTTTCTGGCATCATTTATTTTAGCATCGCCAGGACAATTAGCCGCGCCTACCCACCAGGcat                                               >  1:479495/1‑66 (MQ=9)
                 ccATGGGTTACTGGCATCTTTTTGTTTCGCATCGCCAGGACAATTTGAAGAGGCTACCCACCAGGcat                                               >  1:14117970/1‑66 (MQ=2)
                 ccATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGt                                               >  1:9680845/1‑68 (MQ=255)
                 ccATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGt                                               >  1:5159550/1‑68 (MQ=255)
                 ccATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGt                                               >  1:23511871/1‑68 (MQ=255)
                 ccATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGAAGCGGCTACCCACCAGGctt                                               >  1:25766150/1‑66 (MQ=21)
                 ccATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGAAGCCGCTACCCACCAGGcat                                               >  1:6257018/1‑66 (MQ=21)
                 ccATGGGTTACTGGCATCGTTTTATTTCGCATCCCCTGGGCACTTAGAAGAGCCTACCCACCAGGcat                                               >  1:33102408/1‑66 (MQ=14)
                  cATGGGTTAGTGGCTTAGTTTGCTTTCGCATCGCCTGGGCAATTTGTAGAGCCTACCCACCAGGCGtt                                              <  1:30759659/68‑1 (MQ=21)
                  cATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGtt                                              <  1:32191488/68‑1 (MQ=255)
                    tGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt                                            <  1:20542815/68‑1 (MQ=255)
                    tGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt                                            <  1:36733763/68‑1 (MQ=255)
                    tGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt                                            <  1:14827793/68‑1 (MQ=255)
                    tGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt                                            <  1:7192959/68‑1 (MQ=255)
                    tGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt                                            <  1:8562002/68‑1 (MQ=255)
                      ggTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTta                                          >  1:30530878/1‑68 (MQ=255)
                      ggTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTta                                          >  1:28695939/1‑68 (MQ=255)
                          aCTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATaaa                                      >  1:28394421/1‑68 (MQ=255)
                          aCTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATaaa                                      >  1:12038103/1‑68 (MQ=255)
                            tGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCTGTCTTATaaag                                     >  1:26097332/1‑66 (MQ=255)
                            tGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAAt                                     >  1:2642235/1‑67 (MQ=255)
                            tGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAAt                                     >  1:19762640/1‑67 (MQ=255)
                             ggCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTGTTATAAATg                                    >  1:1515614/1‑67 (MQ=255)
                             ggCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATg                                    >  1:9254246/1‑67 (MQ=255)
                                 tgttttgCTTTGGCATCGCGTGGGCAATTTGAAGAGCCTACCCACCAGGCGTTCTTATAAATgcctcc                               <  1:34492224/65‑1 (MQ=18)
                                   tttttCTTTCGTATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                             <  1:24106414/67‑1 (MQ=255)
                                   tttttCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                             <  1:8211466/67‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                             <  1:34036057/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                             <  1:18901218/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                             <  1:15553525/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                             <  1:18500200/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                             <  1:27609131/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                             <  1:10803698/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                             <  1:32836934/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                             <  1:16984727/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                             <  1:17971133/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                             <  1:29361226/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                             <  1:17588822/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                             <  1:24576269/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                             <  1:23945544/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                             <  1:4954022/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                             <  1:6605606/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                             <  1:6723627/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                             <  1:884309/68‑1 (MQ=255)
                                    ttttCTTTCGGATCGCCTGGGCAATTTTCCGAGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                             <  1:35984628/67‑1 (MQ=38)
                                    ttttCTTTCGCAGCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                             <  1:21253265/67‑1 (MQ=255)
                                    cattctttgGCAACGCGTGGGCAATCTTTCGCGCCAACCCACCAGGCGTTCTTATAAATGCCGCCGt                             <  1:13667630/65‑1 (MQ=17)
                                      ttctttcGCATCGCCTGGGCAATTTTCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTtat                          <  1:4000685/68‑1 (MQ=255)
                                      ttctttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTtat                          <  1:3712809/68‑1 (MQ=255)
                                      ttctttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTtat                          <  1:16554568/68‑1 (MQ=255)
                                      ttctttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTtat                          <  1:24977056/68‑1 (MQ=255)
                                      ttctttcGCAGCTCCTGTGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTtat                          <  1:1904292/68‑1 (MQ=39)
                                        ctttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCCCAAGGCGTTCTTATAAATGCCGCCGTTATAg                        >  1:17373057/1‑68 (MQ=255)
                                          ttcGCATCGCGTGGGCAATTTGTCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                      <  1:23072720/68‑1 (MQ=255)
                                          ttcGCATCGCCTGGGCAATTTGCCGCGCGTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                      <  1:34777298/68‑1 (MQ=255)
                                          ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                      <  1:22248797/68‑1 (MQ=255)
                                          ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                      <  1:32475803/68‑1 (MQ=255)
                                          ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                      <  1:16025596/68‑1 (MQ=255)
                                          ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                      >  1:19505306/1‑68 (MQ=255)
                                          ttcGCACCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                      <  1:1469489/68‑1 (MQ=255)
                                          atcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                      <  1:13933003/67‑1 (MQ=255)
                                          atcGCATCCCCTGGGCAATCTTACGCGCCTACCCACCATGCGTTCTTATAAATGCCGCCGTTATAGAt                      <  1:14538157/67‑1 (MQ=255)
                                            cGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATaa                    <  1:20679592/68‑1 (MQ=255)
                                             gCATCGTCTGTGAAATTTGCAGAGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAAc                   <  1:28070191/68‑1 (MQ=37)
                                               aTCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACg                  <  1:12712230/67‑1 (MQ=255)
                                               aTCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                 <  1:33770725/68‑1 (MQ=255)
                                               aTCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                 <  1:4514854/68‑1 (MQ=255)
                                                  gcgtGGGCAATTTTTCGCGCCAACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTAcc              <  1:35476548/65‑1 (MQ=25)
                                                  gCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTAcc              <  1:21683859/68‑1 (MQ=255)
                                                  gCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTAcc              <  1:18358789/68‑1 (MQ=255)
                                                    gtGGGCAATTTTTAGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt            <  1:2274995/67‑1 (MQ=37)
                                                    cTGGGCAATTTTGCGCGCCAACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt            <  1:21044972/68‑1 (MQ=38)
                                                    cTGGGCAATTTTAAGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt            <  1:13482009/68‑1 (MQ=38)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt            <  1:22118115/68‑1 (MQ=255)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt            <  1:17467201/68‑1 (MQ=255)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt            <  1:4378966/68‑1 (MQ=255)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt            <  1:12351336/68‑1 (MQ=255)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCACGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt            <  1:12968464/68‑1 (MQ=255)
                                                        gCAATTTTTCGAGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTa         <  1:31604221/67‑1 (MQ=38)
                                                        gCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTag        >  1:26598581/1‑68 (MQ=255)
                                                        gCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTa         <  1:33349767/67‑1 (MQ=255)
                                                        gCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTa         <  1:10309639/67‑1 (MQ=255)
                                                        gCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAAAGGCGCCGTTATAGATGACGTACCAGGTag        >  1:27827527/1‑68 (MQ=37)
                                                        gCAATTTGACGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTa         <  1:30226164/67‑1 (MQ=255)
                                                        gCAATTTGACGAGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTa         <  1:27197705/67‑1 (MQ=255)
                                                         cAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTaga       >  1:5598740/1‑68 (MQ=255)
                                                         cAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTaga       >  1:20212788/1‑68 (MQ=255)
                                                         cAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTag        <  1:13189381/67‑1 (MQ=255)
                                                          aaTTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTagaa      >  1:21448802/1‑68 (MQ=255)
                                                                |                                                                 
GCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAAC  >  CP009273/3462374‑3462503

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: