Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,462,438 C→A 9.8% 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 9.8% 2762.9 / 47.9 983A259S (GCA→TCA) ‡chiAperiplasmic endochitinase
Reads supporting (aligned to +/- strand):  ref base C (221/631);  new base A (39/60);  total (261/722)
Fisher's exact test for biased strand distribution p-value = 5.99e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 984 total aligned reads displayed.

GCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCT  >  CP009273/3462374‑3462505
                                                                |                                                                   
     gcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCt                                                              <  1:48938882/67‑1 (MQ=255)
       gAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTAcc                                                           <  1:22927858/68‑1 (MQ=255)
         tcATAACGCCATGGGTTACTGGCAGCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTAccc                                                          <  1:10428228/66‑1 (MQ=255)
         aCATCACTCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTAccc                                                          <  1:14142364/67‑1 (MQ=255)
           aTAACGCCATGGGTTACTGGCATCGTTTACTTTCGCATCGCCAGGCCAATTTGCCGAGCCTAcccccc                                                       >  1:30548293/1‑68 (MQ=25)
            tAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACccacca                                                      <  1:10677090/68‑1 (MQ=255)
              aCGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAg                                                     >  1:46639509/1‑67 (MQ=255)
                gCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCAGCGCCTACCCACCAGGCg                                                  >  1:25411174/1‑68 (MQ=255)
                 ccATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGt                                                 >  1:5034379/1‑68 (MQ=255)
                 ccATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGt                                                 >  1:16392032/1‑68 (MQ=255)
                 ccATGGGTTACTGGCATCGTTTTATTTCGCATCGCCAGGCCAATTTGAAGCGGCTACCCACCAGGcat                                                 >  1:8249430/1‑66 (MQ=17)
                  cATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGtt                                                <  1:40797089/68‑1 (MQ=255)
                  cATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGtt                                                <  1:31906331/68‑1 (MQ=255)
                   aTGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTc                                               <  1:35163275/68‑1 (MQ=255)
                    tGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt                                              <  1:19241276/68‑1 (MQ=255)
                    tGGGTTACTGGCATCGTTATCTATCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt                                              <  1:18347540/68‑1 (MQ=255)
                    tGGGTTACTGGCATCGTCATCTATCGCAGCCCCTGGGCAATCTGACGCGCCAACCCACCAGGCGTTCt                                              <  1:33592933/68‑1 (MQ=21)
                    tGGGTTACTGGCATCGTATTCTATCGCATCGCCTGGGCAATTTTCCGCGCCTACCCACCAGGCGTTCt                                              <  1:7962534/68‑1 (MQ=255)
                     gggTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCATTTTGCCGCGCCTACCCACCAGGCtttttt                                             >  1:18607345/1‑68 (MQ=38)
                      ggTTACTGGCATCGTTTTCTTTCGCATCGCCAGGGCACTTTGCCGCGCCTACCCACCAGGCGttatta                                            >  1:29268966/1‑68 (MQ=38)
                      ggTTACCGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGttatta                                            >  1:23937740/1‑68 (MQ=255)
                       gTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGAAGCGCCTACCCACCAGGCATTTTtat                                           >  1:2865745/1‑68 (MQ=21)
                       gTTACTGGCATCGTTTTCTTTCGCATCGCCTGGCCAATTTGCAGCGCCTACCCACCAGGCATTTTtat                                           >  1:42454306/1‑68 (MQ=25)
                        ttACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTtata                                          >  1:42876126/1‑68 (MQ=255)
                        ttACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGAAGCGGCTACCCACCAGGCGTTTTtata                                          >  1:48198432/1‑68 (MQ=21)
                          aCTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCTTTATtatata                                        >  1:10006482/1‑66 (MQ=38)
                            tGGCATCGTTTTCTTTCGCATCGCCTGGGCAGTTTGCCGCGCCTACCCACCAGGCGTTATTATaaaa                                       >  1:47860772/1‑66 (MQ=38)
                            tGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTATTATaaaac                                      >  1:3638786/1‑66 (MQ=255)
                             ggCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATg                                      >  1:7168998/1‑67 (MQ=255)
                             ggCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTATTATaaagc                                      >  1:3604352/1‑65 (MQ=255)
                             ggCATCGTTTTATTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTATTATAAATg                                      >  1:31429157/1‑67 (MQ=255)
                             ggCATCGTCTTCTATAGCTTCGCCTGGGCAATTTTCAGAGCCAACCCACCAGGCGTTCTTATAAAt                                       <  1:16769808/66‑1 (MQ=12)
                                 tCGTTTTCTTTCGCATCGCCTGGGCAATTTTCCGCGCCTACCCACCAGGCGTTCTTATAAATgccgcc                                 <  1:23326747/68‑1 (MQ=255)
                                 tCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATgccgcc                                 <  1:1806493/68‑1 (MQ=255)
                                 tCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATgccgcc                                 <  1:20216482/68‑1 (MQ=255)
                                  cGTTTTCTTTCGCATCGCCTGTGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGccgccg                                >  1:14643378/1‑68 (MQ=255)
                                   ttgttctttgGCATCGCGTGGGCAATATTCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:35726304/65‑1 (MQ=37)
                                   gTTTTCTTTCGGATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:26286795/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAATTGCCGCCGt                               <  1:27165307/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:9531451/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:34592415/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:39669502/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:45046612/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:47492945/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:48593580/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:6659790/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:7608963/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:2991699/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:16999613/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:22585951/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:24297630/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:24913810/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:1533473/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:1348244/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:12836852/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:11580668/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:28076333/68‑1 (MQ=255)
                                   gTTTTCTTTCGCAGCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:45633387/68‑1 (MQ=255)
                                   gTTTGCTATCGCATCCCCTGGGCAATCTTACGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:50982306/68‑1 (MQ=25)
                                   gTTATCTTTCGCATCGCCTGGGCAATTTTACGCGCCAACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:17580224/68‑1 (MQ=37)
                                    ttttATTTCCCAGCGCGTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:871365/67‑1 (MQ=38)
                                    tgttctttcGCATCGCCTGGGCAATTTGACGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:44613628/65‑1 (MQ=255)
                                     gtgctttcGCATCCCCTGGGCAATTTTCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:30917557/63‑1 (MQ=38)
                                      ttctttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTtat                            <  1:50053351/68‑1 (MQ=255)
                                      ttctttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTtat                            <  1:22167681/68‑1 (MQ=255)
                                       tctttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTtata                           >  1:33013456/1‑68 (MQ=255)
                                        ctttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAg                          >  1:19550438/1‑68 (MQ=255)
                                         tttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGa                         >  1:32758633/1‑68 (MQ=255)
                                          ttcGCATCGCCTGGGCAATTTGGCGCGCCAACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                        <  1:27747961/68‑1 (MQ=255)
                                          ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                        <  1:43612318/68‑1 (MQ=255)
                                          ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                        <  1:21574594/68‑1 (MQ=255)
                                          ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                        <  1:40431147/68‑1 (MQ=255)
                                          ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                        <  1:38455685/68‑1 (MQ=255)
                                          ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                        <  1:12415152/68‑1 (MQ=255)
                                          ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                        >  1:26907884/1‑68 (MQ=255)
                                          ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                        <  1:29596762/68‑1 (MQ=255)
                                          atcGGATCGCCTGGGCAATTTGTCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                        <  1:49350643/67‑1 (MQ=39)
                                             gCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAAc                     <  1:11119744/68‑1 (MQ=255)
                                             gCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAAc                     <  1:28468673/68‑1 (MQ=255)
                                               aTCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                   <  1:5954832/68‑1 (MQ=255)
                                               aTCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                   <  1:15879/68‑1 (MQ=255)
                                               aTCGCCTGGGCAATTTGACGAGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                   <  1:40110387/68‑1 (MQ=255)
                                                  gCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTAcc                <  1:5533599/68‑1 (MQ=255)
                                                  gCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTAcc                <  1:41231653/68‑1 (MQ=255)
                                                  ccgtGGGCAATTTTCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTAcc                <  1:44141331/65‑1 (MQ=255)
                                                   ccTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCa               >  1:21055324/1‑68 (MQ=255)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt              <  1:32470636/68‑1 (MQ=255)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt              <  1:3741311/68‑1 (MQ=255)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt              <  1:37995687/68‑1 (MQ=255)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt              <  1:43202201/68‑1 (MQ=255)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt              <  1:20648492/68‑1 (MQ=255)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt              <  1:46861510/68‑1 (MQ=255)
                                                    cTGGGCAATTTGCCGCGCCTACCAACAAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt              <  1:51243598/68‑1 (MQ=255)
                                                      gggCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGt            <  1:41608230/68‑1 (MQ=255)
                                                        gCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTag          >  1:2595895/1‑68 (MQ=255)
                                                        gCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTag          >  1:12160667/1‑68 (MQ=255)
                                                        gCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTa           <  1:27394911/67‑1 (MQ=255)
                                                        gCAATCTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTa           <  1:36777839/67‑1 (MQ=255)
                                                          aaTTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTaga         <  1:42062824/67‑1 (MQ=255)
                                                            tttGACGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTagaaga      <  1:38944586/68‑1 (MQ=255)
                                                                |                                                                   
GCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCT  >  CP009273/3462374‑3462505

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: