Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,462,439 C→A 6.4% 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.4% 1999.8 / 30.5 646A258A (GCG→GCT) ‡chiAperiplasmic endochitinase
Reads supporting (aligned to +/- strand):  ref base C (151/443);  new base A (23/28);  total (174/472)
Fisher's exact test for biased strand distribution p-value = 4.59e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 699 total aligned reads displayed.

GCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCTT  >  CP009273/3462374‑3462506
                                                                 |                                                                   
 ccgCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCcgc                                                                   >  1:14798122/1‑67 (MQ=255)
  cgcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCcgcg                                                                  <  1:461803/67‑1 (MQ=255)
  cgcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCcgcg                                                                  <  1:5893081/67‑1 (MQ=255)
  cgcgcgAACAAAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCcgcg                                                                  <  1:18511078/67‑1 (MQ=255)
   gcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGcc                                                                <  1:4235677/68‑1 (MQ=255)
       gAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTAcc                                                            <  1:29650903/68‑1 (MQ=255)
       gAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTAcc                                                            <  1:6380547/68‑1 (MQ=255)
         aCATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTAccc                                                           <  1:20229431/67‑1 (MQ=255)
         aCATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACcca                                                          <  1:7582390/68‑1 (MQ=255)
           aTAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACccacc                                                        >  1:11447996/1‑68 (MQ=255)
               cGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGGCTACCCACCAGGc                                                    >  1:22857951/1‑68 (MQ=255)
               cGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGc                                                    >  1:25518158/1‑68 (MQ=255)
               cGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGc                                                    >  1:11136091/1‑68 (MQ=255)
               cGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGACGCGCCTACCCACCAGGc                                                    >  1:23228213/1‑68 (MQ=255)
                gCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCAGAGGCTACCCACCAGGct                                                   >  1:30837841/1‑67 (MQ=21)
                gCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCGGGGCAATTTGAAGCGGCTACCCACCAGGca                                                   >  1:28166895/1‑67 (MQ=21)
                 ccATGGGTTACTGGCATCGTTTTATTTCGCATCGCCAGGGCATTTAGAAGCGCCTACCCACCAGGcat                                                  >  1:12984012/1‑66 (MQ=17)
                  cATGGGTTATTGGAATCGTCTGCTTTCGCATCGCCTGGGCAATCTGCCGCGCCTACCCACCAGGCGtt                                                 <  1:32918682/68‑1 (MQ=255)
                  cATGGGTTACTGGAATCTTTTTCTTTCGCATCGCCTGGGCAATTTTCCGCGCCTACCCACCAGGCGtt                                                 <  1:18918890/68‑1 (MQ=255)
                    tGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt                                               <  1:19601273/68‑1 (MQ=255)
                    tGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt                                               <  1:31728710/68‑1 (MQ=255)
                    tGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt                                               <  1:30564269/68‑1 (MQ=255)
                    tGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt                                               <  1:8562002/68‑1 (MQ=255)
                    tGGGTTACTGGCATCGTTTTCTTTCCCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt                                               <  1:34645576/68‑1 (MQ=255)
                     gggTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCtt                                              >  1:26134891/1‑68 (MQ=255)
                       gTTACTGGCATCGTTATCTATCGCATCGCCTGGGCAATATGCCGCGCCTACCCACCAGGCGTTCTta                                             <  1:28567644/67‑1 (MQ=255)
                            tGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTGTtatataa                                        >  1:19892890/1‑64 (MQ=255)
                            tGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAAt                                        >  1:32321996/1‑67 (MQ=255)
                            tGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCAACAGGAGTTCTTATAAAt                                        >  1:3638629/1‑67 (MQ=255)
                             ggCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCTTTTTTATATATg                                       >  1:15981663/1‑67 (MQ=38)
                             ggCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATg                                       >  1:30392586/1‑67 (MQ=255)
                              gCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATgcc                                     >  1:21059320/1‑68 (MQ=255)
                                 tCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATgccgcc                                  <  1:36943676/68‑1 (MQ=255)
                                 tCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATgccgcc                                  <  1:14678495/68‑1 (MQ=255)
                                 tCGTTTGCTTTGGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATgccgcc                                  <  1:13447169/68‑1 (MQ=255)
                                   gtgttctttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:24008649/65‑1 (MQ=255)
                                   gTTTTCTTTGGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:1210224/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:15553367/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:10803698/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:884309/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:13917485/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:6723627/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:35771728/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:31154280/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:4398441/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:17588822/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:20543883/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:17971133/68‑1 (MQ=255)
                                   gTTTTCTATCGCATCGCCTGGGCAATTTTCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:7208872/68‑1 (MQ=255)
                                    ttttATATCGCATCGCGTGGGCTATTTTCAGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:9745187/67‑1 (MQ=25)
                                    tgttctctgGCACCGCGTGGGCAATCTTTAGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:27053354/65‑1 (MQ=9)
                                    tcttctttcGTCTCCCCTGGGCAATTTTACGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:27685204/65‑1 (MQ=37)
                                    tcttctttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:21535127/65‑1 (MQ=255)
                                    tcttctttaGCATCTTCAGCGCAATTTGCCGAGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:13607419/65‑1 (MQ=25)
                                      ttctttgGCATCCCCTGGGCAATTTTTCGAGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTtat                             <  1:24237657/68‑1 (MQ=25)
                                      ttctttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTtat                             <  1:25894920/68‑1 (MQ=255)
                                      ttctttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTtat                             <  1:19134034/68‑1 (MQ=255)
                                      ttctttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTtat                             <  1:16554568/68‑1 (MQ=255)
                                      ttctttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTtat                             <  1:32690518/68‑1 (MQ=255)
                                       tctttgGCATCGCCTGGGCAATCTTAAGAGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTtata                            <  1:5118310/68‑1 (MQ=21)
                                       tctttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTtata                            <  1:14115714/68‑1 (MQ=255)
                                        ctttgGCATCCCCTGGGCAATTTTCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAg                           <  1:5636406/68‑1 (MQ=38)
                                        ctttcGCATCGCCTGTGCTATATGCCGCGCCTACCCACCATGCGTTCTTTTAAATGCCGCCGTTATAg                           <  1:12362725/68‑1 (MQ=255)
                                        ctttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAg                           <  1:4818725/68‑1 (MQ=255)
                                        ctatcGCATCGCCTGGGCAATTTGTCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAg                           <  1:2246153/65‑1 (MQ=255)
                                          ttcTCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         <  1:22133186/68‑1 (MQ=255)
                                          ttcGCCTCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         <  1:23643808/68‑1 (MQ=255)
                                          ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         <  1:21326948/68‑1 (MQ=255)
                                          ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         <  1:17328446/68‑1 (MQ=255)
                                          ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         <  1:1444792/68‑1 (MQ=255)
                                          ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         <  1:35285155/68‑1 (MQ=255)
                                          ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         <  1:34950574/68‑1 (MQ=255)
                                          ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         >  1:20759895/1‑68 (MQ=255)
                                          ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         <  1:2472893/68‑1 (MQ=255)
                                          ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         >  1:25095414/1‑68 (MQ=255)
                                          ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         <  1:9454301/68‑1 (MQ=255)
                                          ttcGCATCGCCTGGGCAACTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         <  1:21845322/68‑1 (MQ=255)
                                          atcGCATCGCCTGGGCAATTTTCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         <  1:3776193/67‑1 (MQ=255)
                                          atcGCATCGCCTGGGCAATTTGACGAGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         <  1:33545163/67‑1 (MQ=38)
                                             gCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAAc                      <  1:29432775/68‑1 (MQ=255)
                                               aTCGCGTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                    <  1:26424129/68‑1 (MQ=255)
                                               aTCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACg                     <  1:8211629/67‑1 (MQ=255)
                                               aTCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACg                     <  1:12712230/67‑1 (MQ=255)
                                                  gCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTAcc                 <  1:18358789/68‑1 (MQ=255)
                                                    gtGGGCAATTTGTAGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:29094463/67‑1 (MQ=38)
                                                    cTGGGCAATTTTCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:26652164/68‑1 (MQ=255)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:15065451/68‑1 (MQ=255)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:17467201/68‑1 (MQ=255)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:34218310/68‑1 (MQ=255)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:8012807/68‑1 (MQ=255)
                                                     tGGGCGATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATg              >  1:9284529/1‑68 (MQ=255)
                                                        gCAATTTGCCGCGCCTACCCACCAGGCGTTTTTATATACCGCGCCGTTATAGAGGACGTACCagg              >  1:33920953/1‑63 (MQ=14)
                                                        gCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTa            <  1:33166152/67‑1 (MQ=255)
                                                        gCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTa            <  1:10309639/67‑1 (MQ=255)
                                                          aaTTTGCCGCGCCTACCCACCAGGCGTTCTTATAAAAGCCGCCGGTATAGATTACGGACCTGGTagaa         >  1:13174781/1‑68 (MQ=37)
                                                             ttGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATGACGTACCATGTagaaga       >  1:27940233/1‑67 (MQ=255)
                                                                 cgcgCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCtt  <  1:16967903/68‑1 (MQ=255)
                                                                 cgcgCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCtt  <  1:6627887/68‑1 (MQ=255)
                                                                 cgcgCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCt   <  1:5311787/67‑1 (MQ=255)
                                                                 cgcgCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCt   <  1:27361727/67‑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 < 21 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: