Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,462,439 C→A 5.5% 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 5.5% 1578.7 / 14.7 509A258A (GCG→GCT) ‡chiAperiplasmic endochitinase
Reads supporting (aligned to +/- strand):  ref base C (98/375);  new base A (12/24);  total (110/399)
Fisher's exact test for biased strand distribution p-value = 9.18e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 553 total aligned reads displayed.

TTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCTT  >  CP009273/3462372‑3462506
                                                                   |                                                                   
   ccgCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCAGGGCAATTTGCcgc                                                                   >  1:14899296/1‑67 (MQ=255)
    cgcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCcgcg                                                                  <  1:10412117/67‑1 (MQ=255)
    cgcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCcgcg                                                                  <  1:7455889/67‑1 (MQ=255)
    cgcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCcgcg                                                                  <  1:3891892/67‑1 (MQ=255)
    cgcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCcgcg                                                                  <  1:12586938/67‑1 (MQ=255)
    cgcgcgAACATAACGCCATGGGTTACTGGCATCGTTGAGTTTCGCATCGCCAGGGCAGTTAGCagcg                                                                  >  1:961521/1‑67 (MQ=14)
     gcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGcc                                                                <  1:8077486/68‑1 (MQ=255)
         gAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTAcc                                                            <  1:25498932/68‑1 (MQ=255)
           aCATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTAccc                                                           <  1:15843572/67‑1 (MQ=255)
             aTAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGAAGCGGCTACccacc                                                        >  1:13286486/1‑68 (MQ=21)
                aCGCCATGGGTTACTGGCATCGTTTGCTTTCGCATCCCCAGGGCAGTTTGAAGCGGCTACCCACCAg                                                      >  1:23401749/1‑67 (MQ=14)
                 cGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGc                                                    >  1:9894150/1‑68 (MQ=255)
                 cGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGc                                                    <  1:19129482/68‑1 (MQ=255)
                  gCCATGGGTTACTTGCATACTTTTATTACGAATCCCCAGGCCAATTTTCAGCGCCTACCCACCAGGCg                                                   >  1:26856304/1‑68 (MQ=255)
                   ccATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGt                                                  >  1:12928377/1‑68 (MQ=255)
                    cATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGtt                                                 <  1:26507357/68‑1 (MQ=255)
                    cATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGtt                                                 <  1:23577752/68‑1 (MQ=255)
                    cATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGt                                                  <  1:26687197/67‑1 (MQ=255)
                    cATGGGTTACTGGATTATTAATCTTTCGCATCGCCTGGGCAATTTTACGAGCCTACCCACCAGGCGtt                                                 <  1:14414928/68‑1 (MQ=255)
                      tGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt                                               <  1:7741850/68‑1 (MQ=255)
                      tGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt                                               <  1:5354971/68‑1 (MQ=255)
                      tGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt                                               <  1:29266390/68‑1 (MQ=255)
                      tGGGTTAATGGAATCGTGTTCTATCGCATCGCGTGGGCAATTTTTCGCGCCTACCCACCAGGCGTTCt                                               <  1:13529821/68‑1 (MQ=25)
                       gggTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCtt                                              >  1:9001010/1‑68 (MQ=255)
                        ggTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTta                                             >  1:18522684/1‑68 (MQ=255)
                              tGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCCCCAGGCGTTATTATAAAt                                        >  1:24841540/1‑67 (MQ=255)
                              tGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAAt                                        >  1:15897027/1‑67 (MQ=255)
                              tGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAAt                                        >  1:22800055/1‑67 (MQ=255)
                              tGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAAt                                        >  1:15604915/1‑67 (MQ=255)
                               ggCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTGTTATAAATg                                       >  1:1943113/1‑67 (MQ=255)
                                   tCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATgccgcc                                  <  1:14173219/68‑1 (MQ=255)
                                     tttttCTTTCGGATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:10055629/67‑1 (MQ=255)
                                     ttcttctttgGCATCGCCTGGGCAATTTGGCGCGCCTACCCACCATGCGTTCTTATAAATGCCGCCGt                                <  1:12146970/65‑1 (MQ=255)
                                     gtattctttcGCATCGCCTGGGCAATGTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:9291024/65‑1 (MQ=255)
                                     gTTTTCTTTGGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:428271/68‑1 (MQ=255)
                                     gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:28963340/68‑1 (MQ=255)
                                     gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:23900046/68‑1 (MQ=255)
                                     gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:11809685/68‑1 (MQ=255)
                                     gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:23085514/68‑1 (MQ=255)
                                     gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:24702432/68‑1 (MQ=255)
                                     gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:25062915/68‑1 (MQ=255)
                                     gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:5817732/68‑1 (MQ=255)
                                     gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:2068270/68‑1 (MQ=255)
                                     gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:8827535/68‑1 (MQ=255)
                                     gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:4530949/68‑1 (MQ=255)
                                     gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:27793672/68‑1 (MQ=255)
                                     gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:16728576/68‑1 (MQ=255)
                                       tttctttcGTATCGCCTGTGCAATTTTCCGCGCCTACCCACCAGTCGTTCTTATATATGCCGCCGt                                <  1:29743871/66‑1 (MQ=255)
                                        ttctttcGCATCGCCTGGGCAATTTTCAGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTtat                             <  1:15985591/68‑1 (MQ=255)
                                        ttctttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTtat                             <  1:11655681/68‑1 (MQ=255)
                                        ttctttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTtat                             <  1:27571693/68‑1 (MQ=255)
                                        ttctttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTtat                             <  1:29409720/68‑1 (MQ=255)
                                        ttctttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTtat                             <  1:2455856/68‑1 (MQ=255)
                                          ctttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAg                           <  1:833029/68‑1 (MQ=255)
                                            ttggCATCGCGTGGGCAATTTTTCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         <  1:21275855/65‑1 (MQ=37)
                                            ttggCATCGCCTGGGCAATATTTAGAGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         <  1:19662668/65‑1 (MQ=21)
                                            ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTGTTATAAATGCCGCCGTTATAGAt                         >  1:3266911/1‑68 (MQ=255)
                                            ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         <  1:12298397/68‑1 (MQ=255)
                                            ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         <  1:6272783/68‑1 (MQ=255)
                                            ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         >  1:9449972/1‑68 (MQ=255)
                                            ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         <  1:8635470/68‑1 (MQ=255)
                                            ttcGCATCCCCTGGGCAATTTTAAGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         <  1:28360134/68‑1 (MQ=37)
                                            ttcGCAACCCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         <  1:24374808/68‑1 (MQ=255)
                                            ttcGAATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         <  1:7356755/68‑1 (MQ=255)
                                               gCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAAc                      <  1:25931266/68‑1 (MQ=255)
                                                 aTCGCGTGGGCAATTTTCAGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                    <  1:11011551/68‑1 (MQ=37)
                                                 aTCGCCTGGGCAATTTTCCGAGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                    <  1:11176886/68‑1 (MQ=255)
                                                 aTCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTGTTATATATGCCGCCGTGATAGAGAACgg                    >  1:10648484/1‑67 (MQ=255)
                                                 aTCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                    <  1:17908111/68‑1 (MQ=255)
                                                 aTCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                    <  1:22667619/68‑1 (MQ=255)
                                                 aTCGCCTGGGCAATTTGCCGAGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                    <  1:17049944/68‑1 (MQ=255)
                                                 aTCGCCTGGGAAATTTGACGCGCCTACCCACCAGGCGTTCTTATAAATGCAGCCGTTATAGATAACGt                    <  1:16074242/68‑1 (MQ=39)
                                                   ccccTGGGCAATTTTTCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTAc                  <  1:16328275/66‑1 (MQ=38)
                                                    gCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTAcc                 <  1:28110552/68‑1 (MQ=255)
                                                    gCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTAcc                 <  1:3606016/68‑1 (MQ=255)
                                                     ccTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCa                >  1:15337206/1‑68 (MQ=255)
                                                      cTGGGCACTTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:20374289/68‑1 (MQ=255)
                                                      cTGGGCAATTTTCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:5479794/68‑1 (MQ=255)
                                                      cTGGGCAATTTGGCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:28900695/68‑1 (MQ=255)
                                                      cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               >  1:17475684/1‑68 (MQ=255)
                                                      cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:18947100/68‑1 (MQ=255)
                                                      cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:2228568/68‑1 (MQ=255)
                                                      cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:22424282/68‑1 (MQ=255)
                                                      cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:25105613/68‑1 (MQ=255)
                                                      cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:15039622/68‑1 (MQ=255)
                                                      cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:2732993/68‑1 (MQ=255)
                                                      cTGGGCAATGTGTCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:6959216/68‑1 (MQ=255)
                                                          gCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTag           >  1:13814786/1‑68 (MQ=255)
                                                          gCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTa            <  1:4881026/67‑1 (MQ=255)
                                                          gCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTa            <  1:14009191/67‑1 (MQ=255)
                                                          gCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTa            <  1:7998434/67‑1 (MQ=255)
                                                           cAATTTTTCGCGCCTACCCACCAGTCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTat           <  1:672680/67‑2 (MQ=255)
                                                           cAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTaga          >  1:12421366/1‑68 (MQ=255)
                                                                   cgcgCCTACCCACCAGGCGTTCTTATATATCCCGCCGTTATAGATAACGTACATTGTAGAAGAACCtt  >  1:26146876/1‑68 (MQ=255)
                                                                   cgcgCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCtt  <  1:21908418/68‑1 (MQ=255)
                                                                   cgcgCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCtt  <  1:20828274/68‑1 (MQ=255)
                                                                   cgcgCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCtt  <  1:20144753/68‑1 (MQ=255)
                                                                   cgcgCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCtt  <  1:11421355/68‑1 (MQ=255)
                                                                   cgcgCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCtt  <  1:18141728/68‑1 (MQ=255)
                                                                   cgcgCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCt   <  1:24153417/67‑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 < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: