Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,462,676 A→T 7.9% A179A (GCT→GCA chiA ← periplasmic endochitinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,462,6760AT 7.9% 3264.0 / 74.6 990A179A (GCT→GCAchiAperiplasmic endochitinase
Reads supporting (aligned to +/- strand):  ref base A (536/358);  new base T (11/85);  total (547/443)
Fisher's exact test for biased strand distribution p-value = 5.94e-21
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 1000 total aligned reads displayed.

CCGGTGAAGGCTGTGGTGCTGAAGTACAGGACTGTGGATTACTGGTTTCTGAGATTTCCGTTGCCGTAGCGGCACGAACATAACGCCATGGGTTGCTGGCATCATTGCTTTTCGCATCACCCGGACAGTTAGAAG  >  CP009273/3462609‑3462743
                                                                   |                                                                   
     gAAGGCTGTGGTGCTGAAGTACAGGACTGTGGATTACTGGTTTCTGAGATTTCCGTTGCCGTAGCgg                                                                 <  1:52989882/67‑1 (MQ=255)
      aaGGCTGTGGTGCTGAAGTACAGGACTGTGGATTACTGGTTTCTGAGATTTCCGTTGCCGTAGCGGc                                                                >  1:31487584/1‑67 (MQ=255)
       aGGCTGTGGTGCTGAAGTACAGGACTGTGGATTACTGGTTTCTGAGATTTCCGTTGCCGTAGCGGCa                                                               <  1:19019947/67‑1 (MQ=255)
       aGGCTGTGGTGCTGAAGTACAGGACTGTGGATTACTGGTTTCTGAGATTTCCGTTGCCGTAGCGGCAc                                                              >  1:11881677/1‑68 (MQ=255)
       aGGATGTGGTGTTGAAGTACAGTAGTGTGGATTACTGGTTTCTGAGATTTCCGTTGCCGTAGCGGCa                                                               <  1:12168659/67‑1 (MQ=255)
         gCTGTGGTGCTGAAGTACAGGACTGTGGATTACTGGTTTCTGAGATTTCCGTTGCCGTAGCGGCACGa                                                            >  1:2455678/1‑68 (MQ=255)
         gCTGTGGTGCTGAAGTACAGGACTGTGGATTACTGGTTTCTGAGATTTCCGTTGCCGTAGCGGCACGa                                                            >  1:27080096/1‑68 (MQ=255)
         gCTGTGGTGCTGAAGTACAGGACTGTGGATTACTGGTTTCTGAGATTTCCGTTGCCGTAGCGGCACGa                                                            >  1:43888403/1‑68 (MQ=255)
          cTGTGGTGCTGAAGTACAGGACTGTGGATTACTGGTTTCTGAGATTTCCGTTGCCGTAGCGGCACGa                                                            >  1:19470376/1‑67 (MQ=255)
          cTGTGGTGCTGAAGTACAGGACTGTGGATTACTGGTTTCTGAGATTTCCGTTGCCGTAGCGGCACGa                                                            >  1:52568007/1‑67 (MQ=255)
          cTGTGGTGCTGAAGTACAGGACTGTGGATTACTGGTTTCTGAGATTTCCGTTGCCGTAGCGGCACGa                                                            <  1:12677412/67‑1 (MQ=255)
          cTGTGGTGCTGAAGTACAGGACTGTGGATTACTGGTTTCTGAGATTTCCGTTGCCGTAGCGGCACGa                                                            >  1:58031496/1‑67 (MQ=255)
           tgtgGTGCTGAAGTACAGGACTGTGGATTACTGGTTTCTGAGATTTCCGTTGCCGTAGCGGCACGaa                                                           <  1:27899114/67‑1 (MQ=255)
           tgtgGTGCTGAAGTACAGGACTGTGGATTACTGGTTTCTGAGATTTCCGTTGCCGTAGCGGCACGAAc                                                          <  1:42299418/68‑1 (MQ=255)
           tgtgGTGCTGAAGTACAGGACTGTGGATTACTGGTTTCTGAGATTTCCGTTGCCGTAGCGGCACGAAc                                                          <  1:31974972/68‑1 (MQ=255)
            gtggtgCTGAAGTACAGGACTGTGGATTACTGGTTTCTGAGATTTCCGTTGCCGTAGCGGCACGAACa                                                         <  1:40803177/68‑1 (MQ=255)
              ggtgCTGAAGTACAGGACTGTGGATTACTGGTTTCTGAGATTTCCGTTGCCGTAGCGGCACGAACATa                                                       >  1:17536716/1‑68 (MQ=255)
                 gCTGAAGTACAGGACTGTGGATTACTGGTTTCTGAGATTTCCGTTGCCGTAGCGGCACGAACATAAc                                                     >  1:19689659/1‑67 (MQ=255)
                 gCTGAAGTACAGGACTGTGGATTACTGGTTTCTGAGATTTCCGTTGCCGTAGCGGCACGAACATAACg                                                    <  1:772244/68‑1 (MQ=255)
                   tGAAGTACAGGACTGTGGATTACTGGTTTCTGAGATTTCCGTTGCCGTAGCGGCACGAACATAACGcc                                                  >  1:14475329/1‑68 (MQ=255)
                   tGAAGTACAGGACTGTGGATTACTGGTTTCTGAGATTTCCGTTGCCGTAGCGGCACGAACATAACGcc                                                  >  1:22983267/1‑68 (MQ=255)
                      aGTACAGGACTGTGGATTACTGGTTTCTGAGATTTCCGTTGCCGTAGCGGCACGAACATAACGCCAt                                                <  1:17016002/67‑1 (MQ=255)
                      aGTACAGGACTGTGGATTACTGGTTTCTGAGATTTCCGTTGCCGTAGCGGCACGAACATAACGCCATg                                               >  1:16630286/1‑68 (MQ=255)
                      aGTACAGGACTGTGGATTACTGGTTTCTGAGATTTCCGTTGCCGTAGCGGCACGAACATAACGCCATg                                               >  1:43393526/1‑68 (MQ=255)
                      aGTACAGGACTGTGGATTACTGGTTTCTGAGATTTCCGTTGCCGTAGCGGCACGAACATAACGCCATg                                               >  1:41998091/1‑68 (MQ=255)
                        tACAGGACTGTGGATTACTGGTTTCTGAGATTTCCGTTGCCGTAGCGGCACGAAc                                                          <  1:6439742/55‑1 (MQ=255)
                        tACAGGACTGTGGATTACTGGTTTCTGAGATTTCCGTTGCCGTAGCGGCACGAACATAACGCCATggg                                             >  1:55439442/1‑68 (MQ=255)
                        tACAGGACTGTGGATTACTGGTTTCTGAGATTTCCGTTGCCGTAGCGGCACGAACATAACGCCATggg                                             >  1:38661837/1‑68 (MQ=255)
                         aCAGGACTGTGGATTACTGGTTTCTGAGATTTCCGTTGCCGTAGCGGCACGAACATAACGCCATggg                                             <  1:22376188/67‑1 (MQ=255)
                         aCAGGACTGTGGATTACTGGTTTCTGAGATTTCCGTTGCCGTAGCGGCACGAACATAACGCCATggg                                             <  1:48344650/67‑1 (MQ=255)
                         aCAGGACTGTGGATTACTGGTTTCTGAGATTTCCGTTGCCGTAGCGGCACGAACATAACGCCATgg                                              <  1:36606998/66‑1 (MQ=255)
                           aGGACTGTGGATTACTGGTTTCTGAGATTTCCGTTGCCGTAGCGGCACGAACATAACGCCATGGGtt                                           >  1:41389119/1‑67 (MQ=255)
                           aGGACTGTGGATTACTGGTTTCTGAGATTTCCGTTGCCGTAGCGGCACGAACATAACGCCATGGGtt                                           >  1:26407686/1‑67 (MQ=255)
                           aGGACTGTGGATTACTGGTTTCTGAGATTTCCGTTGCCGTAGCGGCACGAACATAACGCCATGGGtt                                           >  1:21408825/1‑67 (MQ=255)
                           aGGACTGTGGATTACTGGTTTCTGAGATTTCCGTTGCCGTAGCGGCACGAACATAACGCCATGGGtt                                           >  1:39006148/1‑67 (MQ=255)
                            ggACTGTGGATTACTGGTTTCTGAGATTTCCGTTGCCGTAGCGGCACGAACATAACGCCATGGGtt                                           >  1:47438100/1‑66 (MQ=255)
                              aCTGTGGATTACTGGTTTCTGAGATTTCCGTTGCCGTAGCGGCACGAACATAACGCCATGGGTTGCTg                                       >  1:51234363/1‑68 (MQ=255)
                              aCTGTGGATTACTGGTTTCTGAGATTTCCGTTGCCGTAGCGGCACGAACATAACGCCATGGGTTGCTg                                       >  1:35924129/1‑68 (MQ=255)
                              aCTGTGGATTACTGGTTTCTGAGATTTCCGTTGCCGTAGCGGCACGAACATAACGCCATGGGTTGCTg                                       >  1:37095998/1‑68 (MQ=255)
                              aCTGTGGATTACTGGTTTCTGAGATTTCCGTTGCCGTAGCGGCACGAACATAACGCCATGGGTTGCTg                                       >  1:37702487/1‑68 (MQ=255)
                              aCTGTGGATTACTGGTTTCTGAGATTTCCGTTGCCGTAGCGGCACGAACATAACGCCATGGGTTGCTg                                       >  1:49093581/1‑68 (MQ=255)
                              aCTGTGGATTACTGGTTTCTGAGATTTCCGTTGCCGTAGCGGCACGAACATAACGCCATGGGTTGCTg                                       >  1:13589219/1‑68 (MQ=255)
                              aCTGTGGATTACTGGTTTCTGAGATTTCCGTTGCCGTAGCGGCACGAACATAACGCCATGGGTTGCTg                                       >  1:42721413/1‑68 (MQ=255)
                              aCTGTGGATTACTGGTTTCTGAGATTTCCGTTGCCGTAGCGGCACGAACATAACGCCATGGGTTGCTg                                       >  1:13955226/1‑68 (MQ=255)
                              aCTGTGGATTACTGGTTTCTGAGATTTCCGTTGCCGTAGCGGCACGAACATAACGCCATGGGTTGCTg                                       >  1:8455937/1‑68 (MQ=255)
                              aCTGTGGATTACTGGTTTCTGAGATTTCCGTTGCCGTAGCGGCACGAACATAACGCCATGGGTTGCTg                                       >  1:38348081/1‑68 (MQ=255)
                              aCTGTGGATTACTGGTTTCTGAGATTTCCGTTGCCGTAGCGGCACGAACATAACGCCATGGGTTGCTg                                       >  1:1544389/1‑68 (MQ=255)
                              aCTGTGGATTACTGGTTTCTGAGATTTCCGTTGCCGTAGCGGCACGAACATAACGCCATGGGTTGCTg                                       >  1:28575517/1‑68 (MQ=255)
                              aCTGTGGATTACTGGTTTCTGAGATTTCCGTTGCCGTAGCGGCACGAACATAACGCCATGGGTTGCTg                                       >  1:23686707/1‑68 (MQ=255)
                              aCTGTGGATTACTGGTTTCTGAGATTTCCGTTGCCGTAGCGGCACGAACATAACGCCAAGGGTTGCTg                                       >  1:45679924/1‑68 (MQ=255)
                                tgtgGATTACTGGTTTCTGAGATTTCCGTTGCCGTAGCGGCACGAACATAACGCCATGGGTTGCTGGc                                     <  1:5733787/68‑1 (MQ=255)
                                 gtgGATTACTGGTTTCTGAGATTTCCGTTGCCGTAGCGGCACGAACATAACGCCATGGGTTGCTGGca                                    <  1:4474143/68‑1 (MQ=255)
                                   ggATTACTGGTTTCTGAGATTTCCGTTGCCGTAGCGGCACGAACATAACGCCATGGGTTGCTGGcatc                                  >  1:24137186/1‑68 (MQ=255)
                                   ggATTACTGGTTTCTGAGATTTCCGTTGCCGTAGCGGCACGAACATAACGCCATGGGTTGCTGGcatc                                  >  1:34817691/1‑68 (MQ=255)
                                   ggATTACTGGTTTCTGAGATTTCCGTTGCCGTAGCGGCACGAACATAACGCCATGGGTTGCTGGcatc                                  >  1:34160456/1‑68 (MQ=255)
                                   ggATTACTGGTTTCTGAGATTTCCGTTGCCGTAGCGGCACGAACATAACGCCATGGGTTGCTGGcatc                                  >  1:30716076/1‑68 (MQ=255)
                                   ggATTACTGGTTTCTGAGATTTCCGTTGCCGTAGCGGCACGAACATAACGCCATGGGTTGCTGGcatc                                  >  1:47734944/1‑68 (MQ=255)
                                   ggATTACTGGTTTCTGAGATTTCCGTTGCCGTAGCGGCACGAACATAACGCCATGGGTTGCTGGcatc                                  >  1:11227854/1‑68 (MQ=255)
                                   ggATTACTGGTTTCTGAGATTTCCGTTGCCGTAGCGGCACGAACATAACGCCATGGGTTGCTGGcatc                                  >  1:50226862/1‑68 (MQ=255)
                                     aTTACTGGTTTCTGAGATTTCCGTTGCCGTAGCGGCACGAACATAACGCCATGGGTTGCTGGcatcat                                >  1:20465347/1‑68 (MQ=255)
                                        aCTGGTTTCTGAGATTTCCGTTGCCGTAGCGGCACGAACATAACGCCATGGGTTGCTGGCATCAtt                               >  1:54710712/1‑66 (MQ=255)
                                        aCTGGTTTCTGAGATTTCCGTTGCCGTAGCGGCACGAACATAACGCCATGGGTTGCTGGCATCAtt                               >  1:1158950/1‑66 (MQ=255)
                                        aCTGGTTTCTGAGATTTCCGTTGCCGTAGCGGCACGAACATAACGCCATGGGTTGCTGGCATCAtt                               >  1:30219328/1‑66 (MQ=255)
                                        aCTGGTTTCTGAGATTTCCGTTGCCGTAGCGGCACGAACATAACGCCATGGGTTGCTGGCATCAtt                               >  1:13242176/1‑66 (MQ=255)
                                        aCTGGTTTCTGAGATTTCCGTTGCCGTAGCGGCACGAACATAACGCCATGGGTTGCTGGCATCATTGc                             >  1:54028233/1‑68 (MQ=255)
                                           ggTTTCTGAGATTTCCGTTGCCGTAGCGGCACGAACATAACGCCATGGGTTGCTGGCATCATTTCttt                          >  1:48736292/1‑68 (MQ=255)
                                                  gagaTTTCCGTTGCCGTAGCGGCACGAACATAACGCCATGGGTTGCTGGCATCATTGCTTTTCGCAt                    >  1:9534268/1‑67 (MQ=255)
                                                   tgaTTTCTGTTGCCGTTGCCGCTCGAACATAACGCCATGGGTTGCTGGCATCATTGCTTTTCGCATCa                  <  1:25599155/67‑1 (MQ=14)
                                                   tgaTTTCTGTTGCCGTAGCCGCTCGAACATAACGCCATGGGTTGCTGGCATCATTGCTTTTCGCATc                   <  1:33351710/66‑1 (MQ=17)
                                                   tgaTTTCTGTTGCCGTAGCCGCACGAACATAACGCCATGGGTTGCTGGCATCATTGCTTTTCGCATCa                  <  1:33763280/67‑1 (MQ=18)
                                                   agaTTTCCGTTGCCGTAGCGGCACGAACATAACGCCATGGGTTGCTGGCATCATTGCTTTTCGCATc                   >  1:51590152/1‑67 (MQ=255)
                                                    gaTTTCGGTTGCCGTAGCGGCACGAACATAACGCCATGGGTTGCTGGCATCATTGCTTTTCGCATCa                  <  1:46543466/67‑1 (MQ=255)
                                                    gaTTTCCGTTGCCGTAGCGGCACGAACATAACGCCATGGGTTGCTGGCATCATTGCTTTTCGCATCa                  >  1:53477712/1‑67 (MQ=255)
                                                    gaTTTCCGTTGCCGTAGCGGCACGAACATAACGCCATGGGTTGCTGGCATCATTGCTTTTCGCATCa                  <  1:25136432/67‑1 (MQ=255)
                                                    gaTTTCCGTTGCCGTAGCGGCACGAACATAACGCCATGGGTTGCTGGCATCATTGCTTTTCGCATCa                  <  1:7245486/67‑1 (MQ=255)
                                                    gaTTTCCGTTGCCGTAGCGGCACGAACATAACGCCATGGGTTGCTGGCATCATTGCTTTTCGCACCa                  <  1:56422689/67‑1 (MQ=255)
                                                      tttCGGTTGCCGTAGCTGTACGAACATAACGGCATGGGTTGCTGGCATCATTGCTTTTCGCATCAcc                <  1:29541770/67‑1 (MQ=25)
                                                      tttCCGTTGCCGTAGCGGCACGAACATAACGCCATGGGTTGCTGGCATCATTGCTTTTCGCATCAcc                <  1:10333693/67‑1 (MQ=255)
                                                      tttCCGTTGCCGTAGCGGCACGAACATAACGCCATGGGTTGCTGGCATCATTGCTTTTCGCATCAcc                <  1:20105982/67‑1 (MQ=255)
                                                          cGTTGCCGTAGCGGCACGAACATAACGCCATGGGTTGCTGGCATCATTGCTTTTCGCATCACCCGGAc           >  1:403378/1‑68 (MQ=255)
                                                             tGCCGTTGCGGCTCGAACATAACGCCATGGGTTGCTGGCATCCTTGCTTTTCGCATCACCCGGACAg         <  1:56169098/67‑1 (MQ=21)
                                                             tGCCGTTGCGGCTCGAACATAACGCCATGGGTTGCTGGCATCATTGCTTTTCGCATCACCCGGACAg         <  1:20963060/67‑1 (MQ=255)
                                                             tGCCGTTGCGGCTCGAACATAACGCCATGGGTTGCTGGCATCATTGCTTTTCGCATCACCCGGACAg         <  1:56953779/67‑1 (MQ=255)
                                                             tGCCGTTGCGGCTCGAACATAACGCCATGGGTTGCTGGCATCATTGCTTTTCGCATCACCCGGACAGt        <  1:28202637/68‑1 (MQ=255)
                                                             tGCCGTTGCGGCACGAACATAACGCCATGGGTTGCTGGCATCATTGCTTTTCGCATCACCCGGACAg         <  1:32496894/67‑1 (MQ=255)
                                                             tGCCGTTGCCGCTCGAACATAACGCCATGGGTTGCTGGCATCATTGCTTTTCGCATCACCCGGACAg         <  1:9040391/67‑1 (MQ=18)
                                                             tGCCGTTGCCGCTCGAACATAACGCCATGGGTTGCTGGCATCATTGCTTTTCGCATCACCCGGACAg         <  1:39757376/67‑1 (MQ=18)
                                                             tGCCGTTGCCGCTCGAACATAACGCCATGGGTTGCTGGCATCATTGCTTTTCGCATCACCCGGACAg         <  1:14833035/67‑1 (MQ=18)
                                                             tGCCGTTGCCGCTCGAACATAACGCCATGGGTTGCTGGCATCATTGCTTTTCGCATCACCCGGACAg         <  1:21911719/67‑1 (MQ=18)
                                                             tGCCGTAGCGGCTCGAACATAACGCCATGGGTTGCTGGCATCATTGCTTTTCGCATCACCCGGACAg         <  1:35529457/67‑1 (MQ=255)
                                                             tGCCGTAGCGGCACGAACATAACGCCATGGGTTGCTGGCATCATTGCTTTTCGCATCACCCGGACAg         <  1:29121194/67‑1 (MQ=255)
                                                             tGCCGTAGCGGCACGAACATAACGCCATGGGTTGCTGGCATCATTGCTTTTCGCATCACCCGGACAg         <  1:17960207/67‑1 (MQ=255)
                                                             tGCCGTAGCGGCACGAACATAACGCCATGGGTTGCTGGCATCATTGCTTTTCGCATCACCCGGACAg         <  1:45262322/67‑1 (MQ=255)
                                                             tGCCGTAGCGGCACGAACATAACGCCATGGGTTGCTGGCATCATTGCTTTTCGCATCACCCGGACAg         <  1:49735617/67‑1 (MQ=255)
                                                             tGCCGTAGCGGCACGAACATAACGCCATGGGTTGCTGGCATCATTGCTTTTCGCATCACCCGGACAGt        <  1:1015013/68‑1 (MQ=255)
                                                             tGCCGTAGCGGCACGAACATAACGCCATGGGTTGCTGGCATCATTGCTTTTCGCATCACCCGGACAGt        <  1:18364759/68‑1 (MQ=255)
                                                              gCCGTAGCGGCACGAACATAACGCCATGGGTTGCTGGCATCATTGCTTTTCGCATCACCCGGACAGtt       <  1:55192407/68‑1 (MQ=255)
                                                              gCCGTAGCGGCACGAACATAACGCCATGGGTTGCTGGCATCATTGCTTTTCGCATCACCCGGACAGtt       >  1:50552757/1‑68 (MQ=255)
                                                                 gTAGCGGCACGAACATAACGCCATGGGTTGCTGGCATCATTGCTTTTCGCATCACCCGGACAGTTaga    <  1:15899460/68‑1 (MQ=255)
                                                                 gTAGCGGCACGAACATAACGCCATGGGTTGCTGGCATCATTGCTTTTCGCATCACCCGGACAGTTaga    <  1:46950989/68‑1 (MQ=255)
                                                                   |                                                                   
CCGGTGAAGGCTGTGGTGCTGAAGTACAGGACTGTGGATTACTGGTTTCTGAGATTTCCGTTGCCGTAGCGGCACGAACATAACGCCATGGGTTGCTGGCATCATTGCTTTTCGCATCACCCGGACAGTTAGAAG  >  CP009273/3462609‑3462743

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