Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,462,441 C→A 6.4% 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 6.4% 1569.0 / 10.9 509A258S (GCG→TCG) ‡chiAperiplasmic endochitinase
Reads supporting (aligned to +/- strand):  ref base C (103/365);  new base A (7/34);  total (110/399)
Fisher's exact test for biased strand distribution p-value = 5.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 552 total aligned reads displayed.

GCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCTT  >  CP009273/3462374‑3462506
                                                                   |                                                                 
  cgcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCcgcg                                                                  <  1:10412117/67‑1 (MQ=255)
  cgcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCcgcg                                                                  <  1:7455889/67‑1 (MQ=255)
  cgcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCcgcg                                                                  <  1:4625347/67‑1 (MQ=255)
  cgcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCcgcg                                                                  <  1:22345232/67‑1 (MQ=255)
  cgcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCcgcg                                                                  <  1:18239123/67‑1 (MQ=255)
  cgcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCcgcg                                                                  <  1:12586938/67‑1 (MQ=255)
  cgcgcgAACATAACGCCATGGGTTACTGGCATCGTTGAGTTTCGCATCGCCAGGGCAGTTAGCagcg                                                                  >  1:961521/1‑67 (MQ=14)
   gcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGcc                                                                <  1:8077486/68‑1 (MQ=255)
         aCATTACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTAccc                                                           <  1:2254613/67‑1 (MQ=255)
         aCATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACcca                                                          <  1:5967845/68‑1 (MQ=255)
           aTAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGAAGCGGCTACccacc                                                        >  1:13286486/1‑68 (MQ=21)
               cGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGc                                                    >  1:9894150/1‑68 (MQ=255)
                 ccATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGt                                                  >  1:12928377/1‑68 (MQ=255)
                 ccATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGt                                                  >  1:23623716/1‑68 (MQ=255)
                  cATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGtt                                                 <  1:20175780/68‑1 (MQ=255)
                  cATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGt                                                  <  1:25065170/67‑1 (MQ=255)
                  cATGGGTTACTGGCAGCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGtt                                                 <  1:15261314/68‑1 (MQ=255)
                  cATGGGTTACTGGATTATTAATCTTTCGCATCGCCTGGGCAATTTTACGAGCCTACCCACCAGGCGtt                                                 <  1:14414928/68‑1 (MQ=255)
                    tGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt                                               <  1:7741850/68‑1 (MQ=255)
                    tGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt                                               <  1:29266390/68‑1 (MQ=255)
                    tGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt                                               <  1:26164055/68‑1 (MQ=255)
                    tGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt                                               <  1:14838843/68‑1 (MQ=255)
                    tGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt                                               <  1:2562956/68‑1 (MQ=255)
                    tGGGTTAATGGAATCGTGTTCTATCGCATCGCGTGGGCAATTTTTCGCGCCTACCCACCAGGCGTTCt                                               <  1:13529821/68‑1 (MQ=25)
                     gggTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCtt                                              >  1:9001010/1‑68 (MQ=255)
                            tGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAAt                                        >  1:15897027/1‑67 (MQ=255)
                            tGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAAt                                        >  1:15604915/1‑67 (MQ=255)
                            tGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAAt                                        >  1:19754591/1‑67 (MQ=255)
                            tGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCATTATTATaaaa                                        >  1:17535774/1‑66 (MQ=38)
                                 tCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATgccgcc                                  <  1:23169873/68‑1 (MQ=255)
                                 tCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATgccgcc                                  <  1:24223156/68‑1 (MQ=255)
                                  cGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCTTTCTTATAAATgccgc                                   >  1:25151475/1‑66 (MQ=255)
                                   tttttCTTTGGCATCGCCTGGGCAATTTGAAGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:24441119/67‑1 (MQ=38)
                                   tttttCTTTCGGATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:10055629/67‑1 (MQ=255)
                                   ttcttctttgGCATCGCCTGGGCAATTTGGCGCGCCTACCCACCATGCGTTCTTATAAATGCCGCCGt                                <  1:12146970/65‑1 (MQ=255)
                                   gtattctttcGCATCGCCTGGGCAATGTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:9291024/65‑1 (MQ=255)
                                   gTTTTCTTTGGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:428271/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTTACGAGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:273834/68‑1 (MQ=38)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:8827535/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:5358049/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:26596261/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:24714334/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:28963340/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:21331620/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:22239048/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:11809685/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:16728576/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:23446862/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:23927037/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCCCCTGGGCAATTTGTCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:26792555/68‑1 (MQ=255)
                                     tttctttcGTATCGCCTGTGCAATTTTCCGCGCCTACCCACCAGTCGTTCTTATATATGCCGCCGt                                <  1:29743871/66‑1 (MQ=255)
                                      ttctttcGCATCGCCTGGGCAATTTTCAGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTtat                             <  1:15985591/68‑1 (MQ=255)
                                      ttctttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTtat                             <  1:20908438/68‑1 (MQ=255)
                                      ttctttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTtat                             <  1:29409720/68‑1 (MQ=255)
                                      ttctttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTtat                             <  1:11655681/68‑1 (MQ=255)
                                      tgctttaGCATCTCATGTGAAATATGAAGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTtat                             <  1:27796376/66‑1 (MQ=19)
                                        ctttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAg                           <  1:833029/68‑1 (MQ=255)
                                          ttggCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         <  1:18576024/65‑1 (MQ=255)
                                          ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         <  1:14173203/68‑1 (MQ=255)
                                          ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         <  1:3905965/68‑1 (MQ=255)
                                          ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         <  1:8635470/68‑1 (MQ=255)
                                          ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         <  1:6272783/68‑1 (MQ=255)
                                          ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         <  1:12298397/68‑1 (MQ=255)
                                          ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         >  1:9449972/1‑68 (MQ=255)
                                          ttcGCATCCCCTGGGCAATTTTAAGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         <  1:28360134/68‑1 (MQ=37)
                                          ttcGAATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         <  1:7356755/68‑1 (MQ=255)
                                           tcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATa                        <  1:22870178/68‑1 (MQ=255)
                                           tcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCCCCGTTAAAGATa                        >  1:19263351/1‑68 (MQ=255)
                                               aTCGCGTGGGCAATTTTCAGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                    <  1:11011551/68‑1 (MQ=37)
                                               aTCGCCTTGGGAATCTTTCGCGCCTACCCAACATTCGTGATTATAAATGCCGCCGTTATAGATAACGt                    <  1:22670596/68‑1 (MQ=255)
                                               aTCGCCTGGGCAATTTTCCGAGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                    <  1:11176886/68‑1 (MQ=255)
                                               aTCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTGTTATATATGCCGCCGTGATAGAGAACgg                    >  1:10648484/1‑67 (MQ=255)
                                               aTCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACg                     <  1:2463135/67‑1 (MQ=255)
                                               aTCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACg                     <  1:15812016/67‑1 (MQ=255)
                                               aTCGCCTGGGAAATTTGACGCGCCTACCCACCAGGCGTTCTTATAAATGCAGCCGTTATAGATAACGt                    <  1:16074242/68‑1 (MQ=39)
                                               aTCCCCTGGGCAATTTTAAGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                    <  1:17914202/68‑1 (MQ=37)
                                                 ccccTGGGCAATTTTTCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTAc                  <  1:16328275/66‑1 (MQ=38)
                                                  gCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTAcc                 <  1:3606016/68‑1 (MQ=255)
                                                  gCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTAcc                 <  1:28110552/68‑1 (MQ=255)
                                                    cTGTGCAATTTGCAGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:27656481/68‑1 (MQ=255)
                                                    cTGGGCAATTTTCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:5479794/68‑1 (MQ=255)
                                                    cTGGGCAATTTGGCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:28900695/68‑1 (MQ=255)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:3425509/68‑1 (MQ=255)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:15039622/68‑1 (MQ=255)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:18947100/68‑1 (MQ=255)
                                                    cTGGGCAATGTGTCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:6959216/68‑1 (MQ=255)
                                                      gggCAATTTTAAGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGt             <  1:17198201/68‑1 (MQ=38)
                                                        gCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTag           >  1:13814786/1‑68 (MQ=255)
                                                        gCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTa            <  1:7998434/67‑1 (MQ=255)
                                                        gCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTa            <  1:2596538/67‑1 (MQ=255)
                                                        gCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTa            <  1:4881026/67‑1 (MQ=255)
                                                        gCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTa            <  1:14009191/67‑1 (MQ=255)
                                                        gCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTa            <  1:26975928/67‑1 (MQ=255)
                                                        gCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTAACATGTag           >  1:19473187/1‑68 (MQ=255)
                                                         cAATTTTTCGCGCCTACCCACCAGTCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTat           <  1:672680/67‑2 (MQ=255)
                                                         cAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTaga          >  1:12421366/1‑68 (MQ=255)
                                                         cAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTaga          >  1:20753350/1‑68 (MQ=255)
                                                              tGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTagaagac      >  1:24872273/1‑66 (MQ=255)
                                                                 cgcgCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCtt  <  1:11421355/68‑1 (MQ=255)
                                                                 cgcgCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCt   <  1:20495842/67‑1 (MQ=255)
                                                                   |                                                                 
GCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCTT  >  CP009273/3462374‑3462506

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: