Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,462,438 C→A 6.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 6.9% 1946.0 / 16.5 656A259S (GCA→TCA) chiAperiplasmic endochitinase
Reads supporting (aligned to +/- strand):  ref base C (154/434);  new base A (21/29);  total (175/481)
Fisher's exact test for biased strand distribution p-value = 2.06e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 658 total aligned reads displayed.

TTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCT  >  CP009273/3462372‑3462505
                                                                  |                                                                   
ttgccgCGCGAACATAACGCCATGGGTTTCTGGCATCCTTGTTTTTCGCATCGCCGGGGCAGTTAGca                                                                    >  1:9296104/1‑67 (MQ=14)
    cgcgcgAACCTAGCGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCcgcg                                                                 <  1:16307677/67‑1 (MQ=255)
    cgcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCcgcg                                                                 <  1:15538075/67‑1 (MQ=255)
    cgcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCcgcg                                                                 <  1:16983254/67‑1 (MQ=255)
    cgcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCcgcg                                                                 <  1:8031879/67‑1 (MQ=255)
    cgcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCcgcg                                                                 <  1:24028498/67‑1 (MQ=255)
     gcgcgAACATCACGCGATCGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGTGCAATTTGCCGCGcc                                                               <  1:30513354/68‑1 (MQ=255)
     gcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGcc                                                               <  1:12774297/68‑1 (MQ=255)
         gAACATAACTCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTAcc                                                           <  1:6540380/68‑1 (MQ=255)
         gAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTAcc                                                           <  1:27813206/68‑1 (MQ=255)
           aCATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACcca                                                         <  1:5886476/68‑1 (MQ=255)
                aCGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCAGGGCAATTTGAAGAGGCAACCCACCAg                                                     >  1:3394919/1‑67 (MQ=18)
                 cGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGc                                                   >  1:31037545/1‑68 (MQ=255)
                 cGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGc                                                   >  1:10926913/1‑68 (MQ=255)
                 cGCCATGGGTTACTGGACTCGTTAGCTATCGCACCGCCTGGGCAATTTGCCGAGCCTACCCACCAGGc                                                   <  1:14400859/68‑1 (MQ=255)
                   ccATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGt                                                 >  1:24396597/1‑68 (MQ=255)
                   ccATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGt                                                 >  1:16604688/1‑68 (MQ=255)
                   ccATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCAGCGCCTACCCACCAGGCGt                                                 >  1:2801389/1‑68 (MQ=255)
                    cATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGtt                                                <  1:13425093/68‑1 (MQ=255)
                    cATGGGTTACTGGCATCGTATTCTATGGCATCGCCTGGGCAATATGCCGCGCCTACCCACCAGGCGt                                                 <  1:18626023/67‑1 (MQ=255)
                     aTGGGTTACTGGCATCGTTTTATTTCGCATCGCCTGGGCAGTTTGCAGCGCCTACCCACCAGGCAttt                                               >  1:4477663/1‑67 (MQ=21)
                      tGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt                                              <  1:1774426/68‑1 (MQ=255)
                      tGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt                                              <  1:20178570/68‑1 (MQ=255)
                      tGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt                                              <  1:19603653/68‑1 (MQ=255)
                        ggTTACTGGCATCGTTTTGTTTCGCATCGCCAGGGCAGTTAGCAGAGCCTACCCACCAGGCTTTTTta                                            >  1:21905571/1‑68 (MQ=14)
                        ggTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTta                                            >  1:15402591/1‑68 (MQ=255)
                        ggTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTta                                            >  1:26928539/1‑68 (MQ=255)
                        ggTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTta                                            >  1:12624359/1‑68 (MQ=255)
                        ggTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCAttatta                                            >  1:21018624/1‑68 (MQ=255)
                        ggTTACTGGCATCGTTTTCTTTCGCATAGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTta                                            >  1:10025706/1‑68 (MQ=255)
                         gTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTta                                            <  1:11380813/67‑1 (MQ=255)
                         gTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTta                                            <  1:31166203/67‑1 (MQ=255)
                             atGGCATCGTGTGCTTTCGAATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAAt                                       <  1:25210919/67‑1 (MQ=38)
                              tGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTGTTATAAATg                                      >  1:561736/1‑68 (MQ=255)
                              tGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAAt                                       >  1:225434/1‑67 (MQ=255)
                               ggCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTGTTATAAATg                                      >  1:20626442/1‑67 (MQ=255)
                               ggCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTGTTATAAATg                                      >  1:11536138/1‑67 (MQ=255)
                               ggCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATg                                      >  1:12009586/1‑67 (MQ=255)
                               ggCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATg                                      >  1:7565019/1‑67 (MQ=255)
                               ggCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCATTATTATaaaac                                      >  1:28663236/1‑65 (MQ=38)
                               ggCATCGTTTTCTTTCGCATCGCCTGGGCAATTGCCGCGCCTACCCACCAGGCGTTCTTATAAATgcc                                    >  1:32180525/1‑68 (MQ=255)
                                 cATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATgcc                                    >  1:28843787/1‑67 (MQ=255)
                                  aTCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATgccg                                   <  1:31849012/67‑1 (MQ=255)
                                   tctttttCTTTCGCACCCCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATgccgcc                                 <  1:24788049/65‑1 (MQ=38)
                                   tCGTTTTCTTTGGCACCGCCTGGGCAATTTGAAGCGCCTACCCACCAGGCGTTCTTATAAATgccgcc                                 <  1:6153431/68‑1 (MQ=38)
                                   tCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATgccgcc                                 <  1:26595166/68‑1 (MQ=255)
                                   tCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATgccgcc                                 <  1:9494935/68‑1 (MQ=255)
                                     ttcttctttgGCATCGCGTGGGCAATTTGCCGAGCCAACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:2917007/65‑1 (MQ=25)
                                     gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:8427978/68‑1 (MQ=255)
                                     gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:28450017/68‑1 (MQ=255)
                                     gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:7065042/68‑1 (MQ=255)
                                     gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:5222807/68‑1 (MQ=255)
                                     gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:1031439/68‑1 (MQ=255)
                                     gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:33284797/68‑1 (MQ=255)
                                     gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:27355853/68‑1 (MQ=255)
                                     gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:18139842/68‑1 (MQ=255)
                                     gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:22786373/68‑1 (MQ=255)
                                     gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:23177160/68‑1 (MQ=255)
                                     gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:14134476/68‑1 (MQ=255)
                                     gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:13140580/68‑1 (MQ=255)
                                     gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:12371185/68‑1 (MQ=255)
                                     gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:11798745/68‑1 (MQ=255)
                                     gTTTTCTTTCGCAGCGCCTGGGCAATTTTCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:29825208/68‑1 (MQ=255)
                                      ttttCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:31551979/67‑1 (MQ=255)
                                        atctatcGCATCGCCTGGGCAATTTTAAGAGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTtat                            <  1:8831078/67‑1 (MQ=25)
                                          ctttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAg                          <  1:7767307/68‑1 (MQ=255)
                                          ctttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAg                          <  1:6899591/68‑1 (MQ=255)
                                            ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                        <  1:15878667/68‑1 (MQ=255)
                                            ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                        <  1:9815106/68‑1 (MQ=255)
                                            ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                        <  1:3689663/68‑1 (MQ=255)
                                            ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                        <  1:34048359/68‑1 (MQ=255)
                                            ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                        <  1:19180383/68‑1 (MQ=255)
                                            ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                        <  1:33530362/68‑1 (MQ=255)
                                            ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                        <  1:3192848/68‑1 (MQ=255)
                                               gCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAAc                     <  1:30786365/68‑1 (MQ=255)
                                                 aTCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                   <  1:10651248/68‑1 (MQ=255)
                                                 aTCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                   <  1:26247402/68‑1 (MQ=255)
                                                 aTCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                   <  1:4196041/68‑1 (MQ=255)
                                                 aTCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                   <  1:25706193/68‑1 (MQ=255)
                                                    gCCTGGGCAATTTTCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTAcc                <  1:8161352/68‑1 (MQ=255)
                                                    gCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTAcc                <  1:27587729/68‑1 (MQ=255)
                                                    gCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTAcc                <  1:3233727/68‑1 (MQ=255)
                                                    cccTGGGCAACTTTCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTAcc                <  1:2229314/67‑1 (MQ=39)
                                                      gtGGGCAATATGACGAGCCAACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt              <  1:4940337/67‑1 (MQ=37)
                                                      cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt              <  1:32817657/68‑1 (MQ=255)
                                                      cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt              <  1:29448136/68‑1 (MQ=255)
                                                      cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt              <  1:13958149/68‑1 (MQ=255)
                                                      cTGGGCAATTTGCAGAGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt              <  1:17183072/68‑1 (MQ=255)
                                                          gCAATTTGCCGCGCCTACCCCCCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTa           <  1:14540444/67‑1 (MQ=255)
                                                          gCAATTTGCCGCGCCTACCCACCAGGCGTTGTTATAAATGCCGCCGTTATAGAATACCTACCATGTag          >  1:9121617/1‑68 (MQ=255)
                                                          gCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGGCGCCGTTATAGATAACGTACCATGTag          >  1:25505531/1‑68 (MQ=255)
                                                          gCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTa           <  1:30104555/67‑1 (MQ=255)
                                                          gCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTa           <  1:1365924/67‑1 (MQ=255)
                                                          gCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTa           <  1:3130880/67‑1 (MQ=255)
                                                          gCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTa           <  1:19925237/67‑1 (MQ=255)
                                                          gCAATTTGCCGCGCCTACCCACAAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTag          >  1:23886813/1‑68 (MQ=255)
                                                            aaTTTGCCGCGCCTACCCACCAGGCGTTGTTATAAATGCCGCCGTGATAGATAACGTCCCATGTagta        >  1:23467994/1‑66 (MQ=255)
                                                              tttGGCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTagaaga      <  1:21389983/68‑1 (MQ=255)
                                                               ttGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTagaaga      >  1:34414358/1‑67 (MQ=255)
                                                                tGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTagaagac     >  1:14832962/1‑66 (MQ=255)
                                                                  |                                                                   
TTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCT  >  CP009273/3462372‑3462505

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: