Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,462,437 G→T 6.5% A259E (GCA→GAA) ‡ chiA ← periplasmic endochitinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,462,4370GT 6.5% 3114.5 / 31.5 980A259E (GCA→GAA) ‡chiAperiplasmic endochitinase
Reads supporting (aligned to +/- strand):  ref base G (256/629);  new base T (2/93);  total (258/722)
Fisher's exact test for biased strand distribution p-value = 6.13e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 985 total aligned reads displayed.

CGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACC  >  CP009273/3462370‑3462504
                                                                   |                                                                   
       gcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGcc                                                              <  1:3603181/68‑1 (MQ=255)
       gcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGcc                                                              <  1:38359140/68‑1 (MQ=255)
       gcgcgAACATAACCCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGcc                                                              <  1:46857095/68‑1 (MQ=255)
         gcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCt                                                             <  1:29222669/67‑1 (MQ=255)
           gAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTAcc                                                          <  1:22927858/68‑1 (MQ=255)
           gAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTAc                                                           >  1:48534868/1‑67 (MQ=255)
           aaaCATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTAcc                                                          <  1:7099063/67‑1 (MQ=255)
             tcATAACGCCATGGGTTACTGGCAGCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTAccc                                                         <  1:10428228/66‑1 (MQ=255)
             aCATCACTCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTAccc                                                         <  1:14142364/67‑1 (MQ=255)
             aCATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACcca                                                        <  1:40763216/68‑1 (MQ=255)
             aCATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACcca                                                        <  1:5367076/68‑1 (MQ=255)
             aCATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACcca                                                        <  1:5855776/68‑1 (MQ=255)
               aTAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGAAGCGCCTACccacc                                                      >  1:49259258/1‑68 (MQ=255)
               aTAACGCCATGGGTTACTGGCATCGTTTACTTTCGCATCGCCAGGCCAATTTGCCGAGCCTAcccccc                                                      >  1:30548293/1‑68 (MQ=25)
                tAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACccacca                                                     <  1:10677090/68‑1 (MQ=255)
                  aCGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGAAGCGGCTACCCACCAg                                                    >  1:39521951/1‑67 (MQ=25)
                   cGCCATGGGTTACTGGCATCGTTTTATTTCGCATAGCCTGGGCAATTTGGAGAACCAACCCAACAgga                                                  >  1:25424010/1‑67 (MQ=19)
                    gCCATGGGTTTCTGGCATCGTTTTCTTTCGCATCGCCAGGACAGTTTGAAGCGGCTACCCACCAGGca                                                 >  1:8701720/1‑67 (MQ=14)
                    gCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCg                                                 >  1:49883696/1‑68 (MQ=255)
                     ccATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGt                                                >  1:25011391/1‑68 (MQ=255)
                      cATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGtt                                               <  1:7890339/68‑1 (MQ=255)
                      cATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGt                                                <  1:47821590/67‑1 (MQ=255)
                        tGGGTTACTGGCATCTTTATCTCTCGCACCGCCTGGGCAATTTGCCGAGCCTACCCACCAGGCGTTCt                                             <  1:42042452/68‑1 (MQ=255)
                        tGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt                                             <  1:28438333/68‑1 (MQ=255)
                        tGGGTTACTGGCATCGTTATCTATCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt                                             <  1:18347540/68‑1 (MQ=255)
                        tGGGTTACTGGCATCGTCATCTATCGCAGCCCCTGGGCAATCTGACGCGCCAACCCACCAGGCGTTCt                                             <  1:33592933/68‑1 (MQ=21)
                         gggTTACTGGCTTAGTGTTCTATCGCAGCGCGTGGGCAATTTTCCGCGCCTACCCACCAGGCGTTCtt                                            <  1:38843202/68‑1 (MQ=25)
                          ggTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGAAGAGCCTACCCACCAGGCATTTTta                                           >  1:48902933/1‑68 (MQ=18)
                          ggTTACTGGCATCGTTTTCTTTCGCATCGCCGGGGCAATTAGAAGAGGCTACCCACCAGGCTTTTTta                                           >  1:48165372/1‑68 (MQ=14)
                          ggTTACCGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGttatta                                           >  1:23937740/1‑68 (MQ=255)
                             tACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATaa                                        <  1:43139079/68‑1 (MQ=255)
                              aCTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCTTTATtatata                                       >  1:10006482/1‑66 (MQ=38)
                                tGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTATTATaaaa                                      >  1:8241272/1‑66 (MQ=255)
                                 ggCATCGTCTTCTATAGCTTCGCCTGGGCAATTTTCAGAGCCAACCCACCAGGCGTTCTTATAAAt                                      <  1:16769808/66‑1 (MQ=12)
                                     tCGTTTTCTTTCGCATCGCCTGGGCAATTTTCCGCGCCTACCCACCAGGCGTTCTTATAAATgccgcc                                <  1:35715004/68‑1 (MQ=255)
                                     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)
                                       tttttCTTTCGCATCGCCTGGGGAATTTTACGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                              <  1:4038074/67‑1 (MQ=38)
                                       tttttCTATCGTACCGCCTGGGCAATTTGACGAGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                              <  1:40015013/67‑1 (MQ=37)
                                       gtcttctttcGCAACGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                              <  1:45044283/65‑1 (MQ=255)
                                       gTTTTCTTTCGGATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                              <  1:26286795/68‑1 (MQ=255)
                                       gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                              <  1:9531451/68‑1 (MQ=255)
                                       gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                              <  1:3515583/68‑1 (MQ=255)
                                       gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                              <  1:3454604/68‑1 (MQ=255)
                                       gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                              <  1:7549203/68‑1 (MQ=255)
                                       gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                              <  1:51209532/68‑1 (MQ=255)
                                       gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                              <  1:42402915/68‑1 (MQ=255)
                                       gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                              <  1:42839142/68‑1 (MQ=255)
                                       gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                              <  1:47244896/68‑1 (MQ=255)
                                       gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                              <  1:4354798/68‑1 (MQ=255)
                                       gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                              <  1:1533473/68‑1 (MQ=255)
                                       gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                              <  1:22585951/68‑1 (MQ=255)
                                       gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                              <  1:24297630/68‑1 (MQ=255)
                                       gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                              <  1:16999613/68‑1 (MQ=255)
                                       gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                              <  1:28076333/68‑1 (MQ=255)
                                       gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                              <  1:16402679/68‑1 (MQ=255)
                                       gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                              <  1:2991699/68‑1 (MQ=255)
                                       gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                              <  1:1348244/68‑1 (MQ=255)
                                       gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                              <  1:12836852/68‑1 (MQ=255)
                                       gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                              <  1:30899750/68‑1 (MQ=255)
                                       gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                              <  1:31873960/68‑1 (MQ=255)
                                       gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                              <  1:11580668/68‑1 (MQ=255)
                                       gTTATCTTTCGCATCGCCTGGGCAATTTTACGCGCCAACCCACCAGGCGTTCTTATAAATGCCGCCGt                              <  1:17580224/68‑1 (MQ=37)
                                          ttctttcGCATCGCGTGGGCAATTTGGCGAGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTtat                           <  1:19263393/68‑1 (MQ=38)
                                          ttctttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTtat                           <  1:37388600/68‑1 (MQ=255)
                                          ttctttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTtat                           <  1:31381728/68‑1 (MQ=255)
                                          ttctttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTtat                           <  1:43960686/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)
                                              ttcGCATCGGCTGGGCAATTTGCCGCGCCTTCCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                       <  1:36759635/68‑1 (MQ=255)
                                              ttcGCATCGCCTGGGCAATTTGGCGCGCCAACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                       <  1:27747961/68‑1 (MQ=255)
                                              ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTGTTATAAATCCCGCCTTTATAGAt                       >  1:44445396/1‑68 (MQ=255)
                                              ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAGATGCCGCCGTTATAGAt                       <  1:6508585/68‑1 (MQ=255)
                                              ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                       >  1:26907884/1‑68 (MQ=255)
                                              ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                       <  1:12415152/68‑1 (MQ=255)
                                              ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                       <  1:21574594/68‑1 (MQ=255)
                                              ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                       <  1:29596762/68‑1 (MQ=255)
                                              atcGCATCGCCTGGGCAATCTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                       <  1:22186034/67‑1 (MQ=255)
                                                 gCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAAc                    <  1:11119744/68‑1 (MQ=255)
                                                 gCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAAc                    <  1:27104150/68‑1 (MQ=255)
                                                  cATCGCGTGGGCAATTTTCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAAc                    <  1:50341294/67‑1 (MQ=255)
                                                   aTCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                  <  1:37843626/68‑1 (MQ=255)
                                                   aTCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                  <  1:15879/68‑1 (MQ=255)
                                                   aTCGCCTGGGCAATTTGACGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                  <  1:32423380/68‑1 (MQ=255)
                                                       ccTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCa              >  1:21055324/1‑68 (MQ=255)
                                                        cTGGGCAATTTTCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt             <  1:46341279/68‑1 (MQ=255)
                                                        cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt             <  1:20648492/68‑1 (MQ=255)
                                                        cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt             <  1:18649446/68‑1 (MQ=255)
                                                        cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt             <  1:36373408/68‑1 (MQ=255)
                                                            gCAATTTTCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCTTTATAGATAACGTACCATGTa          <  1:41575765/67‑1 (MQ=255)
                                                            gCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTag         >  1:12160667/1‑68 (MQ=255)
                                                            gCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTag         >  1:2595895/1‑68 (MQ=255)
                                                            gCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTa          <  1:28650512/67‑1 (MQ=255)
                                                            gCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTa          <  1:27394911/67‑1 (MQ=255)
                                                              aaTTTGCCGCGCCTACCCACCAGGCGTTGTTATATAGCCCCCCGTTAAAGAGGACCTACCAGGTacta       >  1:41185370/1‑65 (MQ=9)
                                                                 ttGCCGCGCCTACCCACCAGGCGTTGTTATAAATGACGCCGTAATAGATAACGTCCCATGAagaaga     >  1:45477384/1‑67 (MQ=255)
                                                                  tGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAAc   <  1:50929972/68‑1 (MQ=255)
                                                                   |                                                                   
CGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACC  >  CP009273/3462370‑3462504

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: