Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,462,400 A→G 5.1% S271S (AGT→AGC chiA ← periplasmic endochitinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,462,4000AG 5.1% 1327.1 / 18.5 410S271S (AGT→AGCchiAperiplasmic endochitinase
Reads supporting (aligned to +/- strand):  ref base A (208/170);  new base G (16/5);  total (233/177)
Fisher's exact test for biased strand distribution p-value = 7.07e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 411 total aligned reads displayed.

AGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAA  >  CP009273/3462336‑3462467
                                                                |                                                                   
aGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTg                                                                  >  1:25954846/1‑68 (MQ=255)
aGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTg                                                                  >  1:1802915/1‑68 (MQ=255)
aGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTg                                                                  >  1:21911274/1‑68 (MQ=255)
aGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTg                                                                  >  1:28294659/1‑68 (MQ=255)
aGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTg                                                                  >  1:29737413/1‑68 (MQ=255)
aGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTg                                                                  >  1:5657845/1‑68 (MQ=255)
 gCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTg                                                                  <  1:23158618/67‑1 (MQ=255)
   cAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGc                                                                <  1:28816569/67‑1 (MQ=255)
     gggTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATc                                                             >  1:2566293/1‑68 (MQ=255)
      ggTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTGCTGGCATCg                                                            >  1:20248862/1‑68 (MQ=255)
      ggTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATct                                                            >  1:19878522/1‑67 (MQ=255)
      ggTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCg                                                            >  1:7998514/1‑68 (MQ=255)
       gTTACCATACTGGCTGATTGCTGTTGCCGTAGCCGCGCGAACATAACGCCATGGGTTACTGGCATc                                                             <  1:11563393/66‑1 (MQ=255)
        ttACCTTAATGGCTTATTTCTGTGGCCGTTGCCGCTCGATCATCACTACATGGGTTACTGGCATCGt                                                           <  1:23505630/67‑1 (MQ=255)
        ttACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGt                                                           <  1:9154961/67‑1 (MQ=255)
        ttACCATACTGCCTGATTTCTGATGCCGTTGCCGCGCGAACCAAACGCCATGGGTTACTGGCATCGt                                                           <  1:25067155/67‑1 (MQ=25)
            cATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTCCTGGCATCAtttaat                                                      >  1:1211704/1‑65 (MQ=21)
            cATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTtttct                                                      >  1:25338119/1‑68 (MQ=255)
            cATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTtttct                                                      >  1:10270099/1‑68 (MQ=255)
            cATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTtttct                                                      >  1:29553622/1‑68 (MQ=255)
            cATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTtttct                                                      >  1:18349005/1‑68 (MQ=255)
            cATACTGGCTGATTTCTGTTGCCGTTGCAGATCGAACATTACGCCATGGGTTACTGGCATCGTtttct                                                      <  1:8554351/68‑1 (MQ=255)
             aTACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTGCt                                                                   >  1:34146649/1‑54 (MQ=255)
              tACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTTCTGGCATCTTtttctt                                                     >  1:32932139/1‑67 (MQ=255)
              tACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTtttctt                                                     >  1:26808711/1‑67 (MQ=255)
              tACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTtttctt                                                     >  1:3985030/1‑67 (MQ=255)
              tACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTtttctt                                                     >  1:34618231/1‑67 (MQ=255)
              tACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTtttctt                                                     >  1:33312954/1‑67 (MQ=255)
              tACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTtttctt                                                     >  1:13681983/1‑67 (MQ=255)
              tACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTtttctt                                                     >  1:32697763/1‑67 (MQ=255)
              tACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTtttctt                                                     >  1:27550758/1‑67 (MQ=255)
              tACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTtttctt                                                     >  1:18861535/1‑67 (MQ=255)
              tACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTtttctt                                                     >  1:27814598/1‑67 (MQ=255)
              tACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTtttctt                                                     >  1:31804813/1‑67 (MQ=255)
              tACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTtttctt                                                     >  1:17427132/1‑67 (MQ=255)
              tACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTttgctt                                                     >  1:14249739/1‑67 (MQ=255)
               aCTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTGCTGGCATCATttattttc                                                   >  1:3148522/1‑68 (MQ=17)
               aCTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTtttgtttc                                                   >  1:17739665/1‑68 (MQ=255)
               aCTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTtttctttc                                                   >  1:30488237/1‑68 (MQ=255)
               aCTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTtttctttc                                                   >  1:21133982/1‑68 (MQ=255)
               aCTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTtttctttc                                                   >  1:5329613/1‑68 (MQ=255)
               aCTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTtgtgtttc                                                   >  1:33931829/1‑68 (MQ=255)
               aCTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTtgactttc                                                   >  1:15183443/1‑68 (MQ=255)
               aCTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCATtttatttc                                                   >  1:8831605/1‑68 (MQ=255)
                cTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCg                                                  <  1:3571552/68‑1 (MQ=255)
                 tGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGc                                                 >  1:168049/1‑68 (MQ=255)
                    cTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATc                                              >  1:12871894/1‑68 (MQ=255)
                     tGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCg                                             <  1:10471873/68‑1 (MQ=255)
                     tGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCg                                             <  1:7784369/68‑1 (MQ=255)
                        tttCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGcc                                           <  1:13182930/67‑1 (MQ=255)
                        tttCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGcc                                           <  1:21262764/67‑1 (MQ=255)
                         ttCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTGCTGGCATCGTTTATTTTCGCATCTccag                                         >  1:5057697/1‑66 (MQ=17)
                         ttCTGTTGCCGTTGCCGCGAGAACATAACGCCATTGGTTACTCGCATCCTTTGCTTTCGCATCGCCTg                                         >  1:31224716/1‑68 (MQ=25)
                          tCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTgg                                        <  1:19424260/68‑1 (MQ=255)
                            tGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGc                                      <  1:8228728/68‑1 (MQ=255)
                               tgccgttgccgCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTGTTTCGCATCGCCGGGACagt                                   >  1:8372165/1‑66 (MQ=18)
                               tgccgttgccgCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCGGGACagt                                   >  1:16438761/1‑66 (MQ=21)
                                  cgttgccgCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAAttt                                 >  1:33076798/1‑67 (MQ=255)
                                  cgttgccgCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAAttt                                 >  1:4542457/1‑67 (MQ=255)
                                      gccgCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCcgc                            <  1:7090951/68‑1 (MQ=255)
                                        cgcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCcgcg                           >  1:632721/1‑67 (MQ=255)
                                        cgcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCcgcg                           <  1:10052349/67‑1 (MQ=255)
                                        cgcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCcgcg                           <  1:11078874/67‑1 (MQ=255)
                                        cgcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCcgcg                           <  1:14853810/67‑1 (MQ=255)
                                        cgcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCcgcg                           <  1:24277165/67‑1 (MQ=255)
                                        cgcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCgcgc                          <  1:12781124/68‑1 (MQ=255)
                                         gcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGcc                         <  1:9389550/68‑1 (MQ=255)
                                         gcgcgAACAAAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGcc                         <  1:6866584/68‑1 (MQ=255)
                                           gcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCt                        <  1:29898235/67‑1 (MQ=255)
                                           gcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTa                       <  1:4195480/68‑1 (MQ=255)
                                             gAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTAcc                     <  1:16010226/68‑1 (MQ=255)
                                             gAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTAcc                     <  1:6628098/68‑1 (MQ=255)
                                               aCATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTAccc                    <  1:16086172/67‑1 (MQ=255)
                                               aCATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACcca                   <  1:9252825/68‑1 (MQ=255)
                                               aCATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACcca                   <  1:19223147/68‑1 (MQ=255)
                                               aCATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACcca                   <  1:6048869/68‑1 (MQ=255)
                                               aCAAAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACcca                   <  1:32217253/68‑1 (MQ=255)
                                                cATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCACGCCTGGGCAATTTGCCGCGCCTACccacc                 <  1:34310409/68‑1 (MQ=255)
                                                     cGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGc             >  1:20462195/1‑68 (MQ=255)
                                                     cGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGc             >  1:1076469/1‑68 (MQ=255)
                                                     cGCCATGGGTTACTGGCATCATTTTCTTTCGCATCGCCCGGACAATTTGAAGCGCCTACCCACCAGGc             >  1:1533387/1‑68 (MQ=17)
                                                      gCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCg            >  1:14570628/1‑68 (MQ=255)
                                                      gCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCg            <  1:222118/68‑1 (MQ=255)
                                                      gCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCg            >  1:7407625/1‑68 (MQ=255)
                                                      gCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCAGAGCCTACCCACCAGGCg            >  1:24805809/1‑68 (MQ=255)
                                                      gCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGAAGCGCCTACCCACCAGGca            >  1:30855654/1‑67 (MQ=25)
                                                      gCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCAGGGCAATTTGCAGCGCCTACCCACCAGGct            >  1:17883980/1‑67 (MQ=37)
                                                       ccATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGAAGAACCTACCCACCAGGctt           >  1:30402593/1‑66 (MQ=21)
                                                        cATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGtt          <  1:26333590/68‑1 (MQ=255)
                                                          tGGGTTACTGTCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt        <  1:22523779/68‑1 (MQ=255)
                                                          tGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt        <  1:14352648/68‑1 (MQ=255)
                                                          tGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt        <  1:2013213/68‑1 (MQ=255)
                                                          tGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt        <  1:15520910/68‑1 (MQ=255)
                                                          tGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt        <  1:13855050/68‑1 (MQ=255)
                                                          tGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt        <  1:28363790/68‑1 (MQ=255)
                                                          tGGGTTACTGGCATCGTTTTCTATCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt        <  1:26594891/68‑1 (MQ=255)
                                                            ggTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTta      >  1:21664335/1‑68 (MQ=255)
                                                            ggTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTta      >  1:4647688/1‑68 (MQ=255)
                                                            ggTTACTGGCATCGTTTTCTTTCGCATCGCCGGGGCAATTTGCCGCGCCTACCCACCAGGCCTTTTta      >  1:2081372/1‑68 (MQ=38)
                                                                aCTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATaaa  >  1:19595008/1‑68 (MQ=255)
                                                                |                                                                   
AGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAA  >  CP009273/3462336‑3462467

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: