Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,462,439 C→A 7.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 7.4% 3329.8 / 56.2 1097A258A (GCG→GCT) ‡chiAperiplasmic endochitinase
Reads supporting (aligned to +/- strand):  ref base C (192/803);  new base A (47/53);  total (239/858)
Fisher's exact test for biased strand distribution p-value = 3.54e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 1167 total aligned reads displayed.

TTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCTT  >  CP009273/3462372‑3462506
                                                                   |                                                                   
    cgcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCcgcg                                                                  <  1:14444209/67‑1 (MQ=255)
    cgcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCcgcg                                                                  <  1:54239964/67‑1 (MQ=255)
           aCATATCGACATGGGTTACTGTCATCGTTTTCTTTCGCTTCGACTGGTAAATTTGCCTCGCCTACcca                                                          <  1:10366484/68‑1 (MQ=255)
           aCATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACcca                                                          <  1:4647950/68‑1 (MQ=255)
           aCAAAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACcca                                                          <  1:8826981/68‑1 (MQ=255)
             aTAACGCCATGGGTTACTGGCATCGTTTTATTTCGCATCGCCGGGGCAATTTGAAGCGCCTACccacc                                                        >  1:30878442/1‑68 (MQ=21)
                  gCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCg                                                   >  1:51668884/1‑68 (MQ=255)
                  gCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGAAGAGGCTACCCACCAGGca                                                   >  1:49064124/1‑67 (MQ=18)
                   ccATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGt                                                  >  1:13710415/1‑68 (MQ=255)
                    cATGGGTTACTGGCATCTTTTTCTTTCGCATCGCCTGGGCAATATTCCGCGCCTACCCACCAGGCGtt                                                 <  1:44282707/68‑1 (MQ=255)
                    cATGGGTTACTGGAATCGTTTTCTTTGGCATCGCCTGGGCAATTTTCCGCGCCTACCCACCAGGCGtt                                                 <  1:34609406/68‑1 (MQ=255)
                     aTGGGTTACTGGCATCGTTTTCTTTCGCATCCCCCGGACAGTTTGAAGCGGCTACccacca                                                       >  1:32395148/1‑61 (MQ=11)
                      tGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt                                               <  1:7165119/68‑1 (MQ=255)
                      tGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt                                               <  1:25256972/68‑1 (MQ=255)
                      tGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGAGCCTACCCACCAGGCGTTCt                                               <  1:8150681/68‑1 (MQ=255)
                      tGGGTTACTGGACTTGTTTTCTTTCGGACCGCGTGGGCAATATTAAGAGCCTACCCACCAGGCGTTCt                                               <  1:52669045/68‑1 (MQ=9)
                        ggTTACTGGCATCGTTTTCTTTCGCATCGCCAGGGCATTTAGAAGCGCCTACCCACCAGGCATTTTta                                             >  1:18354367/1‑68 (MQ=17)
                         gTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGACGCGCCTACCCACCAGGCGTTCTta                                             <  1:22316366/67‑1 (MQ=255)
                              tGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAAt                                        >  1:39656175/1‑67 (MQ=255)
                              tGGCATCGTTTTCTTTCGCATCGCCTGGCCAATTTGACGCGCCTACCCACCAGGCCTTATTATAACTg                                       >  1:41152673/1‑68 (MQ=37)
                               ggCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATg                                       >  1:24616995/1‑67 (MQ=255)
                               ggCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTATTATAAAtc                                       >  1:50671611/1‑66 (MQ=255)
                               ggCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCATTTTTATATATg                                       >  1:38757880/1‑67 (MQ=37)
                                   tctttttCTTTGGCATCGCGTGGGCAATATTCAGAGCCTACCCACCAGGCGTTCTTATAAATgccgcc                                  <  1:25750769/65‑1 (MQ=18)
                                   tCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATgccgcc                                  <  1:33298227/68‑1 (MQ=255)
                                   tCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATgccgcc                                  <  1:47720772/68‑1 (MQ=255)
                                   tCGTTTTCTTTCGCATCGCCGGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATgccgcc                                  <  1:7490656/68‑1 (MQ=255)
                                    cGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTATTATAAAAcccgcca                                 >  1:38254048/1‑67 (MQ=255)
                                     tttttCTTTCGGAACGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:19129138/67‑1 (MQ=39)
                                     gtcttctttcGCATCGCCTGGGCAATTTTCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:31812624/65‑1 (MQ=255)
                                     gTTTTCTTTCGGATCGCCTGGGCAATTTTACGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:29256154/68‑1 (MQ=38)
                                     gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:23875140/68‑1 (MQ=255)
                                     gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:10037143/68‑1 (MQ=255)
                                     gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:36236601/68‑1 (MQ=255)
                                     gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:56947564/68‑1 (MQ=255)
                                     gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:32993232/68‑1 (MQ=255)
                                     gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:37101649/68‑1 (MQ=255)
                                     gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:31202208/68‑1 (MQ=255)
                                     gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:55820528/68‑1 (MQ=255)
                                     gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:27244540/68‑1 (MQ=255)
                                     gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:5031637/68‑1 (MQ=255)
                                     gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:53107854/68‑1 (MQ=255)
                                     gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:40772907/68‑1 (MQ=255)
                                     gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:45161435/68‑1 (MQ=255)
                                     gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:10670205/68‑1 (MQ=255)
                                     gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCCTATAAATGCCGCCGt                                <  1:51128014/68‑1 (MQ=255)
                                     gTTTTCTATGGCACCGCCTGGGCAATTTGACGAGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:17122375/68‑1 (MQ=37)
                                     gTTTTCTATCGCATCGCCTGGGCAATTTGCAGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:24168574/68‑1 (MQ=255)
                                      ttttCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:20459414/67‑1 (MQ=255)
                                      tttgCTATCGAATCCCCTGGGCAATTTTCAGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:37790296/67‑1 (MQ=25)
                                      ttatctttcGAATCGCGTGGGCAATTTTACGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:33620485/64‑1 (MQ=37)
                                      tattctttcGCATCCCCTGGGCAATCTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:6134203/65‑1 (MQ=38)
                                        ttctttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTtat                             <  1:19616294/68‑1 (MQ=255)
                                        ttctttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTtat                             >  1:29854372/1‑68 (MQ=255)
                                         tctttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTtata                            >  1:53822293/1‑68 (MQ=255)
                                          ctttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAg                           <  1:52352538/68‑1 (MQ=255)
                                            ttcGCATCGCCTGGGCAATTTTCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         <  1:16843751/68‑1 (MQ=255)
                                            ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         <  1:45743770/68‑1 (MQ=255)
                                            ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         <  1:21451365/68‑1 (MQ=255)
                                            ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         <  1:42940157/68‑1 (MQ=255)
                                            ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         <  1:13076494/68‑1 (MQ=255)
                                            ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         >  1:26653026/1‑68 (MQ=255)
                                            ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTATTATAAATGCCGCCGTTATAGAt                         >  1:21012745/1‑68 (MQ=255)
                                               gCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAAc                      <  1:53523310/68‑1 (MQ=255)
                                               gCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAAc                      <  1:48601732/68‑1 (MQ=255)
                                                 aTCGCCTGGGGAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                    <  1:42586145/68‑1 (MQ=255)
                                                 aTCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                    <  1:4809970/68‑1 (MQ=255)
                                                 aTCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                    <  1:28792979/68‑1 (MQ=255)
                                                 aTCGCCTGGGCAATCTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                    <  1:11451397/68‑1 (MQ=255)
                                                    gCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTAcc                 <  1:49699714/68‑1 (MQ=255)
                                                    gCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTAcc                 <  1:12059493/68‑1 (MQ=255)
                                                     ccTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCa                >  1:30452546/1‑68 (MQ=255)
                                                     caggggCAATTTGTAGAGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTAcc                 <  1:35186118/64‑1 (MQ=37)
                                                      gtGGGCAATTTGACGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:39256349/67‑1 (MQ=255)
                                                      cggggCAATCTTCCGCGCCAACCCAACAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:930458/66‑1 (MQ=255)
                                                      cTGGGCAATTTTCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:43535631/68‑1 (MQ=255)
                                                      cTGGGCAATTTGCCGCGGCTACCCACCAGGCGTTCTTTTAAATGCCGCCGTTATAGATAACGTACCAt               <  1:56369353/68‑1 (MQ=39)
                                                      cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:28091395/68‑1 (MQ=255)
                                                      cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:18638319/68‑1 (MQ=255)
                                                      cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:19912226/68‑1 (MQ=255)
                                                      cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:21872815/68‑1 (MQ=255)
                                                      cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:1625850/68‑1 (MQ=255)
                                                      cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:26319158/68‑1 (MQ=255)
                                                      cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:4169794/68‑1 (MQ=255)
                                                      cTGGGCAATGTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:27675000/68‑1 (MQ=255)
                                                          gCAATTTTCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTa            <  1:35640119/67‑1 (MQ=255)
                                                          gCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTag           >  1:34190481/1‑68 (MQ=255)
                                                          gCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTa            <  1:12429514/67‑1 (MQ=255)
                                                          gCAATTTGCAGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTa            <  1:15070524/67‑1 (MQ=255)
                                                          gCAATCTGACGAGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTa            <  1:57936611/67‑1 (MQ=38)
                                                           cAATTTGCCGCGCCTACCCACCAGGCGTTGTTATAAAGGCCGCCGTTATAGAAAACGTACCATGTaga          >  1:47054646/1‑68 (MQ=255)
                                                           cAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTaga          >  1:46170195/1‑68 (MQ=255)
                                                           cAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTag           <  1:22694223/67‑1 (MQ=255)
                                                              tttGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTagaaga       <  1:17710293/68‑1 (MQ=255)
                                                              tttGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTagaaga       <  1:40338263/68‑1 (MQ=255)
                                                               ttGCCGCGCCTACCCACCAGGCGTTCTTATAACTGCCGCAGTTATAGATAACGTACCATGTagaaga       >  1:57322091/1‑67 (MQ=255)
                                                                  ccGCGCCAACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCt   <  1:15680135/68‑1 (MQ=255)
                                                                   cgcgCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCtt  <  1:55327594/68‑1 (MQ=255)
                                                                   cgcgCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCtt  <  1:23169936/68‑1 (MQ=255)
                                                                   cgcgCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCtt  <  1:44777010/68‑1 (MQ=255)
                                                                   |                                                                   
TTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCTT  >  CP009273/3462372‑3462506

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: