Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,462,441 C→A 7.5% A258S (GCG→TCG) ‡ chiA ← periplasmic endochitinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,462,4410CA 7.5% 1903.2 / 20.5 619A258S (GCG→TCG) ‡chiAperiplasmic endochitinase
Reads supporting (aligned to +/- strand):  ref base C (139/421);  new base A (15/42);  total (154/465)
Fisher's exact test for biased strand distribution p-value = 8.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 666 total aligned reads displayed.

GCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCTT  >  CP009273/3462374‑3462506
                                                                   |                                                                 
  cgcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCcgcg                                                                  <  1:27996404/67‑1 (MQ=255)
   gcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGcc                                                                <  1:1399468/68‑1 (MQ=255)
       gAACTTTACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTAcc                                                            <  1:31995527/68‑1 (MQ=255)
       gAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTAcc                                                            <  1:25230450/68‑1 (MQ=255)
         aCATCACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACcca                                                          <  1:26943726/68‑1 (MQ=255)
         aCATCACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACcca                                                          <  1:7192780/68‑1 (MQ=255)
               cGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGc                                                    <  1:23489715/68‑1 (MQ=255)
               cGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGAAGAGGCTACCCACCAGGc                                                    >  1:13271445/1‑68 (MQ=18)
                gCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCAGCGCCTACCCACCAGGCg                                                   >  1:17330302/1‑68 (MQ=255)
                gCCATGGGTTACTGGCATCGTTTTATTTCGCATCGCCGGGGCATTTAGCAGAGCCTACCCACCAGGca                                                   >  1:20176022/1‑67 (MQ=17)
                 ccATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGAAGAGCCTACCCACCAGGcat                                                  >  1:22503181/1‑66 (MQ=21)
                 ccATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGCCATTTAGCAGCGCCTACCCACCAGGCGt                                                  >  1:11665102/1‑68 (MQ=25)
                  cATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGtt                                                 <  1:10001066/68‑1 (MQ=255)
                    tGGGTTATTGGATTCTTGTTCTTTGGCATCGCCTGGGCAATTTTCAGCGCCTACCCACCAGGCGTTCt                                               <  1:3322080/68‑1 (MQ=21)
                    tGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt                                               <  1:18136577/68‑1 (MQ=255)
                    tGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt                                               <  1:15054183/68‑1 (MQ=255)
                      ggTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCATTTTta                                             >  1:15244004/1‑68 (MQ=255)
                      ggTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGACGCGCCTACCCACCAGGCTTTTTta                                             >  1:8518406/1‑68 (MQ=37)
                      ggTTACTGGCATCGTTTTCTTTCGCATCGCCAGGGCAGTTTGAAGAGGCTACCCACCAGGCATTTTta                                             >  1:19921427/1‑68 (MQ=11)
                       gTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTtat                                            >  1:23188473/1‑68 (MQ=255)
                       gTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTta                                             <  1:29122450/67‑1 (MQ=255)
                       gTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTta                                             <  1:13539938/67‑1 (MQ=255)
                       gTTACTGGACTCGTCTTCTATCGCACCCCCTGGGCAATCTTACGCGCCTACCCACCAGGCGTTCt                                               <  1:20773077/65‑1 (MQ=16)
                            tGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATaaagg                                       >  1:5989913/1‑66 (MQ=255)
                            tGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCATTATtatataa                                        >  1:6977819/1‑64 (MQ=38)
                            tGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGAAGAGCCTACCCACCAGGCGTTTTAATAAAtt                                       >  1:11819109/1‑67 (MQ=21)
                             ggCATCGTTTTCTTTCGCATCGCCTGGGCAATTTTCCGCGCCTACCCACCAGGCGTTCTTATAAATg                                       >  1:18454442/1‑67 (MQ=255)
                             ggCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATg                                       >  1:1926506/1‑67 (MQ=255)
                               cATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCCTTGTCATAAATgaa                                     >  1:20470530/1‑65 (MQ=255)
                                aTCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATgccg                                    >  1:22964317/1‑67 (MQ=255)
                                 tCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATgccgcc                                  <  1:5809012/68‑1 (MQ=255)
                                 tCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATgccgcc                                  <  1:664823/68‑1 (MQ=255)
                                 tCGATTCCTTTCGGATCGCCTGGGCAATGTTCAGAGCCTACCCACCAGGCGTTCTTATAAATgccgcc                                  <  1:8202742/68‑1 (MQ=21)
                                  cGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGccgccg                                 <  1:24180154/68‑1 (MQ=255)
                                   tttttCTTTCGCATCGCCTGGGCAATTTTCCGCGGCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:14237339/67‑1 (MQ=39)
                                   tttttCTTTCGCATCCCCTGGGCAATATTTCGAGCCAACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:14520491/67‑1 (MQ=18)
                                   gTTTTCTTTGGCATCGCCTGGGCAATTTTCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:10852835/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:12635388/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:23918547/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:7572181/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:25563975/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:25956075/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:9425219/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:2635292/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:11164802/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:28366708/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:155344/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:15740988/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:21161635/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:31794164/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:3088967/68‑1 (MQ=255)
                                    ttttCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGtt                               >  1:29890189/1‑68 (MQ=255)
                                    ttttCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGtt                               >  1:29618352/1‑68 (MQ=255)
                                     tttctttgGGATCGCCTGGGCAATTTTACGCGCCAACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:933380/66‑1 (MQ=25)
                                      ttctttgCCATCGCCTGGGCAATTTGGCGCGCCTACCCACCACGCGTTCTTTTAAATGCCGCCGTtat                             <  1:27164276/68‑1 (MQ=255)
                                      ttctttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTtat                             <  1:21958922/68‑1 (MQ=255)
                                        ctttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAg                           <  1:19602261/68‑1 (MQ=255)
                                        ctttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAg                           <  1:6112715/68‑1 (MQ=255)
                                          ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         <  1:19059520/68‑1 (MQ=255)
                                          ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         <  1:13098933/68‑1 (MQ=255)
                                          ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         <  1:21032175/68‑1 (MQ=255)
                                          ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         <  1:16735033/68‑1 (MQ=255)
                                          ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         <  1:24817734/68‑1 (MQ=255)
                                          ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         <  1:12192063/68‑1 (MQ=255)
                                          ctggCACCGCGTGGGCAATCTTTCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         <  1:9088349/65‑1 (MQ=18)
                                           tcGCATCGCCTGGGCAATATTACGAGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATa                        <  1:8824147/68‑1 (MQ=37)
                                               ctCGCCTTGTCATTTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                    <  1:9613129/67‑1 (MQ=38)
                                               aTCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                    <  1:27591952/68‑1 (MQ=255)
                                               aTCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                    <  1:18783924/68‑1 (MQ=255)
                                               aTCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                    <  1:17718557/68‑1 (MQ=255)
                                               aTCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                    <  1:3125241/68‑1 (MQ=255)
                                               aTCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                    <  1:4558896/68‑1 (MQ=255)
                                               aTCCCCTGGGCAATTTTCAGCGCCAACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                    <  1:17614698/68‑1 (MQ=25)
                                                  gcgtGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTAcc                 <  1:11992727/65‑1 (MQ=255)
                                                   ccTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTATTATAAATCCCGCCGTTATAGAAAACGTACCa                >  1:5625465/1‑68 (MQ=255)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:10220884/68‑1 (MQ=255)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:5519798/68‑1 (MQ=255)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:4691996/68‑1 (MQ=255)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:1150743/68‑1 (MQ=255)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:1285308/68‑1 (MQ=255)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:30794016/68‑1 (MQ=255)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:2453631/68‑1 (MQ=255)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:28884389/68‑1 (MQ=255)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:26705087/68‑1 (MQ=255)
                                                        gCAATTTTCAGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTa            <  1:28514339/67‑1 (MQ=255)
                                                        gCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTa            <  1:4972912/67‑1 (MQ=255)
                                                         cAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTaga          >  1:4078131/1‑68 (MQ=255)
                                                            tttGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTagaaga       <  1:26093313/68‑1 (MQ=255)
                                                             ttGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTagaaga       >  1:10560726/1‑67 (MQ=255)
                                                             ttGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTagaaga       >  1:21582487/1‑67 (MQ=255)
                                                             ttGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTagaaga       >  1:3597633/1‑67 (MQ=255)
                                                             ttGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTagaaca       >  1:2134275/1‑65 (MQ=255)
                                                             ttGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCAGCCGTTTTAGATAACGTGCGATGTataagc       >  1:7861269/1‑66 (MQ=37)
                                                              tGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAgaagaa      >  1:31539120/1‑67 (MQ=255)
                                                               gCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAAcc    <  1:7410187/68‑1 (MQ=255)
                                                               gCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAAcc    <  1:25011364/68‑1 (MQ=255)
                                                                 cgcgCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCtt  <  1:5271058/68‑1 (MQ=255)
                                                                 cgcgCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCt   <  1:30167656/67‑1 (MQ=255)
                                                                 cgcgCCAACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCtt  <  1:16275459/68‑1 (MQ=255)
                                                                  gcgcCTACCCACCAGGCGTTTTTATAAAGGCCGCCGTTATAGAAAACCTACCATGTAGAAGATCCtt  >  1:29342992/1‑67 (MQ=37)
                                                                   |                                                                 
GCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCTT  >  CP009273/3462374‑3462506

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: