Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,462,439 C→A 6.1% 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 6.1% 1965.7 / 27.4 620A258A (GCG→GCTchiAperiplasmic endochitinase
Reads supporting (aligned to +/- strand):  ref base C (112/459);  new base A (22/25);  total (134/486)
Fisher's exact test for biased strand distribution p-value = 6.30e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 663 total aligned reads displayed.

GCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCTT  >  CP009273/3462374‑3462506
                                                                 |                                                                   
gccgCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCcgc                                                                   <  1:25448606/68‑1 (MQ=255)
gccgCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCcgc                                                                   <  1:25747062/68‑1 (MQ=255)
  cgcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCcgcg                                                                  <  1:24277165/67‑1 (MQ=255)
  cgcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCcgcg                                                                  <  1:18757396/67‑1 (MQ=255)
  cgcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCcgcg                                                                  <  1:10052349/67‑1 (MQ=255)
  cgcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCcgcg                                                                  <  1:14853810/67‑1 (MQ=255)
  cgcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCgcgc                                                                 <  1:12781124/68‑1 (MQ=255)
   gcgcgAACAAAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGcc                                                                <  1:6866584/68‑1 (MQ=255)
     gcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCt                                                               <  1:29898235/67‑1 (MQ=255)
     gcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTa                                                              <  1:8800942/68‑1 (MQ=255)
       gAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATGGCCTGGGCAATTTGCCGCGCCTAcc                                                            <  1:13808514/68‑1 (MQ=255)
       gAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTAcc                                                            <  1:16010226/68‑1 (MQ=255)
         aCATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACcca                                                          <  1:1993961/68‑1 (MQ=255)
         aCATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACcca                                                          <  1:9252825/68‑1 (MQ=255)
          cATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACccac                                                         <  1:24652255/68‑1 (MQ=255)
                gCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCg                                                   >  1:14570628/1‑68 (MQ=255)
                 ccATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGAAGCGCCTACCCACCAGGcat                                                  >  1:31975469/1‑66 (MQ=25)
                 ccATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGAAGAACCTACCCACCAGGctt                                                  >  1:30402593/1‑66 (MQ=21)
                 ccATGGGTTACTGGCATCGTTTTATTTCGCATCGCCGGGGCATTTTGAAGAGCCTACCCACCAGGcat                                                  >  1:31627675/1‑66 (MQ=17)
                  cATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGtt                                                 <  1:15245407/68‑1 (MQ=255)
                  cATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGtt                                                 <  1:28018225/68‑1 (MQ=255)
                   aTGGGTTACTGGCATCTTTTTCTATCGCATCCCCTGGGCAATTTTCAGCGCCTACCCACCAGGCGTTc                                                <  1:17932136/68‑1 (MQ=37)
                    tGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt                                               <  1:2013213/68‑1 (MQ=255)
                      ggTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTta                                             >  1:21664335/1‑68 (MQ=255)
                      ggTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTta                                             >  1:4647688/1‑68 (MQ=255)
                       gTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTtat                                            >  1:32999042/1‑68 (MQ=255)
                        ttACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTtata                                           >  1:27552967/1‑68 (MQ=255)
                            tGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCTTTGTtatataa                                        >  1:20757932/1‑64 (MQ=38)
                            tGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAAt                                        >  1:24962960/1‑67 (MQ=255)
                            tGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTATTATATATg                                       >  1:33102075/1‑68 (MQ=255)
                            tGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCAGCGCCTACCCACCAGGCTTTGTCATAAATg                                       >  1:8471198/1‑68 (MQ=38)
                             ggCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATg                                       >  1:29693334/1‑67 (MQ=255)
                             ggCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATg                                       >  1:730651/1‑67 (MQ=255)
                             ggCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATg                                       >  1:8011995/1‑67 (MQ=255)
                             ggCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATg                                       >  1:13150189/1‑67 (MQ=255)
                              gCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATgcc                                     >  1:6565258/1‑68 (MQ=255)
                                 tCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATgccgcc                                  <  1:29119728/68‑1 (MQ=255)
                                 tCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATgccgcc                                  <  1:9558497/68‑1 (MQ=255)
                                 tCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATgccgcc                                  <  1:16099846/68‑1 (MQ=255)
                                   tttttCTTTGGTTTCGCGTGGGGAATTTTACGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:30648631/67‑1 (MQ=25)
                                   tttttCTTTCGCTTCGCCTGGGCAATTTTCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:22300094/67‑1 (MQ=39)
                                   ttgttctttgGCAGCGCCTGGGCAATTTGCAGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:513440/65‑1 (MQ=38)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTTCCGCGCCTACCCACCCGGCGTTCTTATAAATGCCGCCGt                                <  1:23961943/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:26517905/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:15635224/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:30183277/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:12451273/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:6110517/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:18316807/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:19176157/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:19399829/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:19679732/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:23681496/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:12111114/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:3391756/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:10498821/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCAACAGTCGTTCTTATAAATGCCGCCGt                                <  1:15760389/68‑1 (MQ=255)
                                   gTTTTCTTTCGAATCGCCTGGGCAATTTGCAGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:11418141/68‑1 (MQ=255)
                                   gTTTTCTATCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:10786013/68‑1 (MQ=255)
                                   cttttCTATGGCTTCCCCTGGGCAATTTTACGTGCCAACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:20932550/67‑1 (MQ=16)
                                      ttctttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTtat                             <  1:6348071/68‑1 (MQ=255)
                                        ctttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAg                           >  1:15345130/1‑68 (MQ=255)
                                          ttcGCATCGCCTTGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         <  1:32361558/68‑1 (MQ=255)
                                          ttcGCATCGCCTGGGCAATTTTCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         <  1:34069876/68‑1 (MQ=255)
                                          ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         <  1:28768919/68‑1 (MQ=255)
                                          ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         <  1:31062910/68‑1 (MQ=255)
                                          ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         <  1:20447585/68‑1 (MQ=255)
                                          ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         >  1:3791597/1‑68 (MQ=255)
                                          ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         <  1:23391796/68‑1 (MQ=255)
                                          atggCATCGCCTGGGCAATTTTCAGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         <  1:22777524/65‑1 (MQ=38)
                                             gCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAAc                      <  1:28635661/68‑1 (MQ=255)
                                               tccGCCTGGGCAATTTGTCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                    <  1:34581976/66‑1 (MQ=255)
                                               aTCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                    <  1:5453668/68‑1 (MQ=255)
                                               aTCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                    <  1:34290826/68‑1 (MQ=255)
                                               aTCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                    <  1:22009166/68‑1 (MQ=255)
                                               aTCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                    <  1:14166374/68‑1 (MQ=255)
                                               aTCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                    <  1:17169136/68‑1 (MQ=255)
                                               aTCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                    >  1:8612788/1‑68 (MQ=255)
                                                    gtGGGCAATTTTACGAGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:3562407/67‑1 (MQ=37)
                                                    cTGGGCAATTTTCCGCGCCAACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:28427826/68‑1 (MQ=255)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:7574932/68‑1 (MQ=255)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:8972489/68‑1 (MQ=255)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:32751416/68‑1 (MQ=255)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:4112283/68‑1 (MQ=255)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:16302957/68‑1 (MQ=255)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:17534116/68‑1 (MQ=255)
                                                    cTGGGCAATTTGACGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:16627779/68‑1 (MQ=255)
                                                       ggCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:34729179/48‑1 (MQ=255)
                                                        gCAATTTGCCGCGCCTACCCACCAGGCGTTGTTATAAATGCCGCCGTTATAGATAACGTACCATGTag           >  1:26239642/1‑68 (MQ=255)
                                                        gCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTag           >  1:1177615/1‑68 (MQ=255)
                                                        gCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTa            <  1:30931827/67‑1 (MQ=255)
                                                        gCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTa            <  1:13474483/67‑1 (MQ=255)
                                                        gCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTa            <  1:9371939/67‑1 (MQ=255)
                                                        gCAATTTGCCGCGCCTACCCACCAAGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTa            <  1:10303852/67‑1 (MQ=255)
                                                        gCAATTTGCCGAGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTa            <  1:33593864/67‑1 (MQ=255)
                                                         cAATTAGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTaga          <  1:11040415/68‑1 (MQ=255)
                                                          aaTTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTagaa         >  1:21253309/1‑68 (MQ=255)
                                                              tGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAgaagaa      >  1:783681/1‑67 (MQ=255)
                                                                 cgcgCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCtt  <  1:5829642/68‑1 (MQ=255)
                                                                 cgcgCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCtt  <  1:27050614/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 < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: