Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,462,441 C→A 5.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 5.4% 3408.5 / 27.8 1099A258S (GCG→TCG) ‡chiAperiplasmic endochitinase
Reads supporting (aligned to +/- strand):  ref base C (224/798);  new base A (16/61);  total (240/859)
Fisher's exact test for biased strand distribution p-value = 8.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 1170 total aligned reads displayed.

GCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCTT  >  CP009273/3462374‑3462506
                                                                   |                                                                 
  cgcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCcgcg                                                                  <  1:54239964/67‑1 (MQ=255)
   gcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGcc                                                                <  1:55403986/68‑1 (MQ=255)
       gAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTAcc                                                            <  1:47721134/68‑1 (MQ=255)
       gAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTAc                                                             <  1:507024/67‑1 (MQ=255)
         aCATATCGACATGGGTTACTGTCATCGTTTTCTTTCGCTTCGACTGGTAAATTTGCCTCGCCTACcca                                                          <  1:10366484/68‑1 (MQ=255)
         aCAAAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACcca                                                          <  1:8826981/68‑1 (MQ=255)
                gCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCATTTAGAAGCGCCTACCCACCAGGCg                                                   >  1:41710856/1‑68 (MQ=21)
                gCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCg                                                   >  1:51668884/1‑68 (MQ=255)
                gCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCg                                                   <  1:55837320/68‑1 (MQ=255)
                 ccATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGt                                                  >  1:38387612/1‑68 (MQ=255)
                 ccATGGGTTACTGGCATCGTTTTATTTCGCATCGCCAGGGCAATTTGAAGCGGCTACCCACCAGGCGt                                                  >  1:34622982/1‑68 (MQ=21)
                  cATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGtt                                                 <  1:42616276/68‑1 (MQ=255)
                  cATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGtt                                                 <  1:40775172/68‑1 (MQ=255)
                  cATGGGTTACTGGAATCTTAATCTCTGGCATCGCCTGGGCAATTTTCCGAGCCTACCCACCAGGCGtt                                                 <  1:33347301/68‑1 (MQ=21)
                   aTGGGTTACTGGCATCGTTTTCTTTCGCATCCCCCGGACAGTTTGAAGCGGCTACccacca                                                       >  1:32395148/1‑61 (MQ=11)
                    tGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt                                               <  1:7562836/68‑1 (MQ=255)
                    tGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt                                               <  1:25256972/68‑1 (MQ=255)
                    tGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt                                               <  1:43594780/68‑1 (MQ=255)
                    tGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGAGCCTACCCACCAGGCGTTCt                                               <  1:8150681/68‑1 (MQ=255)
                       gTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGtt                                                 <  1:4034186/63‑1 (MQ=255)
                        ttACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTtata                                           >  1:13026675/1‑68 (MQ=255)
                            tGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAAt                                        >  1:57337002/1‑67 (MQ=255)
                            tGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAAt                                        >  1:45226192/1‑67 (MQ=255)
                            tGGCATCGTTTTCTTTCGCATCGCCAGGGCAATTTGACGCGCGGTCCAAACAGGCGATCATATaaaa                                        >  1:23716930/1‑66 (MQ=255)
                             ggCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATg                                       >  1:11929861/1‑67 (MQ=255)
                             ggCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATg                                       >  1:24616995/1‑67 (MQ=255)
                                 tggTTTTCTTTATCATCTCATGCGCAATTTGTAGAGCCTACCCACCAGGCGTTCTTATAAATgccgcc                                  <  1:53539561/66‑1 (MQ=16)
                                 tctttttCTTTGGCATCGCGTGGGCAATATTCAGAGCCTACCCACCAGGCGTTCTTATAAATgccgcc                                  <  1:25750769/65‑1 (MQ=18)
                                 tagTTTTCTTTCGCATCGCCTGGGCAATTTGCCGAGCCTACCCACCAGGCGTTCTTATAAATgccgcc                                  <  1:49065890/66‑1 (MQ=255)
                                 tCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATgccgcc                                  <  1:17678471/68‑1 (MQ=255)
                                   tttttCTTTGGGATCGCGTGGGCAATTTTCAGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:719451/67‑1 (MQ=25)
                                   tttttCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:2100588/67‑1 (MQ=255)
                                   ttcttctatcGCATCCCCTGGGCAATTTTACGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:45768325/65‑1 (MQ=37)
                                   gtcttctttcGCATCGCCTGGGCAATTTTCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:31812624/65‑1 (MQ=255)
                                   gtattctttcGAAGCGCGTGGGCAATTTTAAGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:44781593/65‑1 (MQ=21)
                                   gTTTTCTTTGGCATCGCCTGGGCAATTTGTCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:46227724/68‑1 (MQ=255)
                                   gTTTTCTTTCGGATCGCCTGGGCAATTTTACGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:29256154/68‑1 (MQ=38)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:10037143/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:10670205/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:56951187/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:27244540/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:32993232/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:36370265/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:5031637/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:53107854/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAAATTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:6163987/68‑1 (MQ=255)
                                   gTTTTCTATGGCACCGCCTGGGCAATTTGACGAGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:17122375/68‑1 (MQ=37)
                                   gTTTTCTATCGCATCGCCTGGGCAATTTGCAGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:24168574/68‑1 (MQ=255)
                                    ttttCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:20459414/67‑1 (MQ=255)
                                    ttttCTTTCGCATCGCCTGGGCAATTTGCCGCCCCTACCCACCAGGCGTTATTATAAATGGAGCCGtt                               >  1:57953378/1‑68 (MQ=255)
                                    tttgCTATCGAATCCCCTGGGCAATTTTCAGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:37790296/67‑1 (MQ=25)
                                    ttatctttcGAATCGCGTGGGCAATTTTACGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:33620485/64‑1 (MQ=37)
                                      ttctttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTtat                             >  1:29854372/1‑68 (MQ=255)
                                      ttctttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTtat                             <  1:19098720/68‑1 (MQ=255)
                                      atctttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTtat                             <  1:13618519/67‑1 (MQ=255)
                                          ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         <  1:18276467/68‑1 (MQ=255)
                                          ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         <  1:21451365/68‑1 (MQ=255)
                                          ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         >  1:26653026/1‑68 (MQ=255)
                                          ttcGCATCGCCTGGGCAATATTCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         <  1:5269085/68‑1 (MQ=255)
                                             gCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAAc                      <  1:48601732/68‑1 (MQ=255)
                                               aTCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                    <  1:41207753/68‑1 (MQ=255)
                                               aTCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                    <  1:4809970/68‑1 (MQ=255)
                                               aTCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                    <  1:30828393/68‑1 (MQ=255)
                                               aTCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                    >  1:4297363/1‑68 (MQ=255)
                                               aTCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                    <  1:28792979/68‑1 (MQ=255)
                                               aTCGCCTGGGCAATTTGACGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                    <  1:37271522/68‑1 (MQ=255)
                                               aTCGCCTGGGCAATCTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                    <  1:11451397/68‑1 (MQ=255)
                                                  gCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTAcc                 <  1:49699714/68‑1 (MQ=255)
                                                  gCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTAcc                 <  1:38776377/68‑1 (MQ=255)
                                                   ccTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCa                >  1:30452546/1‑68 (MQ=255)
                                                    gtGGGCAATTTGACGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:39256349/67‑1 (MQ=255)
                                                    cggggCAATCTTCCGCGCCAACCCAACAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:930458/66‑1 (MQ=255)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:53838574/68‑1 (MQ=255)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:12420270/68‑1 (MQ=255)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:1440742/68‑1 (MQ=255)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:18638319/68‑1 (MQ=255)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:19590523/68‑1 (MQ=255)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:19912226/68‑1 (MQ=255)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:5640722/68‑1 (MQ=255)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:21872815/68‑1 (MQ=255)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:22292435/68‑1 (MQ=255)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:26319158/68‑1 (MQ=255)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:28091395/68‑1 (MQ=255)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               >  1:44305544/1‑68 (MQ=255)
                                                    cTGGGCAATGTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:27675000/68‑1 (MQ=255)
                                                       ggCAATTTTTCGAGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTAcc                 <  1:35254088/63‑1 (MQ=37)
                                                        gCAATTTTCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTa            <  1:35640119/67‑1 (MQ=255)
                                                        gCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTag           >  1:34190481/1‑68 (MQ=255)
                                                        gCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTag           >  1:16776205/1‑68 (MQ=255)
                                                         cAATTTGCCGCGCCTACCCACCAGGCGTTGTTATAAAGGCCGCCGTTATAGAAAACGTACCATGTaga          >  1:47054646/1‑68 (MQ=255)
                                                         cAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTag           <  1:22694223/67‑1 (MQ=255)
                                                          aaTTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTagaa         >  1:31091256/1‑68 (MQ=255)
                                                            tttGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTagaaga       <  1:39660696/68‑1 (MQ=255)
                                                             ttGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTCCCATGTagaaga       >  1:1622455/1‑67 (MQ=255)
                                                              tGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAgaagaa      >  1:52372089/1‑67 (MQ=255)
                                                                ccGCGCCAACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCt   <  1:15680135/68‑1 (MQ=255)
                                                                 cgcgCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCtt  <  1:51130770/68‑1 (MQ=255)
                                                                 cgcgCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCtt  <  1:14974608/68‑1 (MQ=255)
                                                                 cgcgCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCtt  <  1:23169936/68‑1 (MQ=255)
                                                                  gcgcCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACGATGTAGAAGCACCtt  >  1:46481191/1‑67 (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: