Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,462,438 C→A 8.9% 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 8.9% 3190.8 / 43.4 1129A259S (GCA→TCA) ‡chiAperiplasmic endochitinase
Reads supporting (aligned to +/- strand):  ref base C (201/784);  new base A (38/60);  total (241/888)
Fisher's exact test for biased strand distribution p-value = 9.78e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 1133 total aligned reads displayed.

TTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCT  >  CP009273/3462372‑3462505
                                                                  |                                                                   
    cgcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCcgcg                                                                 <  1:48994724/67‑1 (MQ=255)
    cgcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCcgcg                                                                 <  1:42846243/67‑1 (MQ=255)
     gcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGcc                                                               <  1:55403986/68‑1 (MQ=255)
         gAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTAcc                                                           <  1:28034725/68‑1 (MQ=255)
         gAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTAcc                                                           <  1:38663250/68‑1 (MQ=255)
           aCATCACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTAccc                                                          <  1:32327017/67‑1 (MQ=255)
           aCATATCGACATGGGTTACTGTCATCGTTTTCTTTCGCTTCGACTGGTAAATTTGCCTCGCCTACcca                                                         <  1:10366484/68‑1 (MQ=255)
           aCATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTAccc                                                          <  1:26648330/67‑1 (MQ=255)
           aCATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACcca                                                         <  1:40922128/68‑1 (MQ=255)
                 cGCCATGGGTTACTGGCATCTTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGc                                                   <  1:45388399/68‑1 (MQ=255)
                  gCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCg                                                  <  1:55837320/68‑1 (MQ=255)
                  gCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGAAGCGCCTACCCACCAGGca                                                  >  1:11415174/1‑67 (MQ=25)
                  gCCATGGGTTACTGGCATCGTTTGCTTTCGCATCTCCAGGGCAATTTGAAGAACATACCCACCAGGca                                                  >  1:30815953/1‑67 (MQ=9)
                   ccATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGt                                                 >  1:31047114/1‑68 (MQ=255)
                    cATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGtt                                                <  1:42296577/68‑1 (MQ=255)
                    cATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGtt                                                <  1:12386752/68‑1 (MQ=255)
                      tGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt                                              <  1:7562836/68‑1 (MQ=255)
                      tGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt                                              <  1:4846607/68‑1 (MQ=255)
                        ggTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGACGCGGCTACCCACCAGGCGTTTTta                                            >  1:19545402/1‑68 (MQ=37)
                        ggTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGAAGCGCCTACCCACCAGGCGttatta                                            >  1:53028344/1‑68 (MQ=37)
                              tGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCTTTCTTATaaag                                       >  1:37664397/1‑66 (MQ=255)
                              tGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAAt                                       >  1:14396319/1‑67 (MQ=255)
                              tGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAAt                                       >  1:57337002/1‑67 (MQ=255)
                              tGGCATCGTTTTCTTTCGCATCGCCAGGGCAATTTGACGCGCGGTCCAAACAGGCGATCATATaaaa                                       >  1:23716930/1‑66 (MQ=255)
                               ggCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATg                                      >  1:21363135/1‑67 (MQ=255)
                               ggCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATg                                      >  1:11929861/1‑67 (MQ=255)
                               ggCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATg                                      >  1:35378035/1‑67 (MQ=255)
                               ggCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATg                                      >  1:24616995/1‑67 (MQ=255)
                                   tctttttCTTTGGCATCGCGTGGGCAATATTCAGAGCCTACCCACCAGGCGTTCTTATAAATgccgcc                                 <  1:25750769/65‑1 (MQ=18)
                                   tctttttCATTCGACTCGCCTTGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATgccgcc                                 <  1:33163823/65‑1 (MQ=37)
                                   tCGTTTTATATCGTACCCCCTGGGCAATTTTAAGAGCCTACCCACCAGGCGTTCTTATAAATgccgcc                                 <  1:4490291/68‑1 (MQ=12)
                                     tttttCTTTCGAATCGCCTGTGCAATTTGCCGCGCCTACCCCCCAGGCGTTCTTATAAATGCCGCCGt                               <  1:46866787/67‑1 (MQ=255)
                                     tttttCTTTCGAATCGCCTGGGCAATTTTGCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:38062777/67‑1 (MQ=38)
                                     ttcttctttgGCATCGCCTGGGCAATTTGCCGCCCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:47678596/65‑1 (MQ=38)
                                     gtcttctttcGCATCCCCTGGGCAATCTTACGAGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:49623021/65‑1 (MQ=21)
                                     gTTTTCTTTCGCATCGCCTGGGCAATTTTCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:31706925/68‑1 (MQ=255)
                                     gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:10037143/68‑1 (MQ=255)
                                     gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:35841678/68‑1 (MQ=255)
                                     gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:43208656/68‑1 (MQ=255)
                                     gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:44628340/68‑1 (MQ=255)
                                     gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:46426242/68‑1 (MQ=255)
                                     gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:2924106/68‑1 (MQ=255)
                                     gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:17175299/68‑1 (MQ=255)
                                     gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:56951187/68‑1 (MQ=255)
                                     gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:27244540/68‑1 (MQ=255)
                                     gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:12964606/68‑1 (MQ=255)
                                     gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:15515176/68‑1 (MQ=255)
                                     gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:16274553/68‑1 (MQ=255)
                                     gTTTTCTTTCGCAGCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:6223551/68‑1 (MQ=255)
                                     gTTTTCTTGCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:54211021/68‑1 (MQ=255)
                                     gTTTTCTATCGCATCGCCTGGGCAATTTGCAGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:24168574/68‑1 (MQ=255)
                                      ttttCTTTGGTATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:33441398/67‑1 (MQ=255)
                                      ttttCTTTCGCATCTCCTTGTCAATTTTAAGAGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:51094191/67‑1 (MQ=25)
                                      ttttCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCCTTATTATAAATGCCGCCGtt                              >  1:50192153/1‑68 (MQ=255)
                                      ttttCTTTCGCATCGCCTGGGCAATTTGCCGCCCCTACCCACCAGGCGTTATTATAAATGGAGCCGtt                              >  1:57953378/1‑68 (MQ=255)
                                      tcttcgttcgTATCCCCTGGGCAATATTTCGAGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:52654616/65‑1 (MQ=18)
                                        ttctttcGCTTCTACTGTGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTtat                            <  1:43779827/68‑1 (MQ=38)
                                        ttctttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTtat                            <  1:30444181/68‑1 (MQ=255)
                                        ttctttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTtat                            <  1:14837631/68‑1 (MQ=255)
                                        ttctttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTtat                            <  1:19098720/68‑1 (MQ=255)
                                        atctttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTtat                            <  1:13618519/67‑1 (MQ=255)
                                         tctttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTtata                           >  1:53822293/1‑68 (MQ=255)
                                          ctttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAg                          <  1:52352538/68‑1 (MQ=255)
                                           tttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGa                         >  1:20971921/1‑68 (MQ=255)
                                            ttggCACCGCGTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                        <  1:29624437/65‑1 (MQ=38)
                                            ttcGCATCGCCTGGGCAATTTTCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                        <  1:16843751/68‑1 (MQ=255)
                                            ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                        <  1:47986478/68‑1 (MQ=255)
                                            ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                        <  1:35027397/68‑1 (MQ=255)
                                            ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                        <  1:40422604/68‑1 (MQ=255)
                                            ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                        <  1:18276467/68‑1 (MQ=255)
                                            ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                        <  1:19885405/68‑1 (MQ=255)
                                              cGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATaa                      <  1:10621110/68‑1 (MQ=255)
                                              cGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATaa                      <  1:36502710/68‑1 (MQ=255)
                                               gCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAAc                     <  1:53523310/68‑1 (MQ=255)
                                                caccGCCTGGGCAATTTTCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAAc                     <  1:32782507/64‑1 (MQ=255)
                                                 aTCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                   <  1:4140378/68‑1 (MQ=255)
                                                 aTCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                   <  1:28792979/68‑1 (MQ=255)
                                                 aTCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                   >  1:2033150/1‑68 (MQ=255)
                                                   cGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTAc                 <  1:8259344/68‑1 (MQ=255)
                                                   cGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTAc                 <  1:34003226/68‑1 (MQ=255)
                                                   cGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTAc                 <  1:25202774/68‑1 (MQ=255)
                                                    gCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTAcc                <  1:46096401/68‑1 (MQ=255)
                                                      cggggCAATCTTCCGCGCCAACCCAACAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt              <  1:930458/66‑1 (MQ=255)
                                                      cTGGGCAATTTTCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt              <  1:8836555/68‑1 (MQ=255)
                                                      cTGGGCAATTTTACGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt              <  1:40054928/68‑1 (MQ=255)
                                                      cTGGGCAATTTGCCGCGGCTACCCACCAGGCGTTCTTTTAAATGCCGCCGTTATAGATAACGTACCAt              <  1:56369353/68‑1 (MQ=39)
                                                      cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt              <  1:51622033/68‑1 (MQ=255)
                                                      cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt              <  1:18638319/68‑1 (MQ=255)
                                                      cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAACGCCGCCGTTATAGATAACGTACCAt              <  1:50572445/68‑1 (MQ=255)
                                                      cTGGGCAATGTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt              <  1:27675000/68‑1 (MQ=255)
                                                        gggCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGt            <  1:23118455/68‑1 (MQ=255)
                                                          gCAATTTTCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTa           <  1:21753215/67‑1 (MQ=255)
                                                          gCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTag          >  1:26255152/1‑68 (MQ=255)
                                                          gCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTag          >  1:34529416/1‑68 (MQ=255)
                                                          gCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTa           <  1:38944293/67‑1 (MQ=255)
                                                          gCAATTTGCAGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTa           <  1:7202751/67‑1 (MQ=255)
                                                           cAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTag          <  1:22694223/67‑1 (MQ=255)
                                                            aaTTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTagaa        >  1:39504630/1‑68 (MQ=255)
                                                              tttGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTagaaga      <  1:22281132/68‑1 (MQ=255)
                                                              tttGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTagaaga      <  1:17710293/68‑1 (MQ=255)
                                                                  |                                                                   
TTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCT  >  CP009273/3462372‑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: