Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,462,439 C→A 7.9% A258A (GCG→GCT) ‡ chiA ← periplasmic endochitinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,462,4390CA 7.9% 1892.1 / 47.4 618A258A (GCG→GCT) ‡chiAperiplasmic endochitinase
Reads supporting (aligned to +/- strand):  ref base C (129/432);  new base A (25/31);  total (154/464)
Fisher's exact test for biased strand distribution p-value = 9.61e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 663 total aligned reads displayed.

GCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCTT  >  CP009273/3462374‑3462506
                                                                 |                                                                   
gccgCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCcgc                                                                   <  1:26302370/68‑1 (MQ=255)
gccgCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCcgc                                                                   <  1:17714266/68‑1 (MQ=255)
 ccgCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCcgc                                                                   >  1:2799056/1‑67 (MQ=255)
   gcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGcc                                                                <  1:1399468/68‑1 (MQ=255)
     gcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTa                                                              <  1:12851194/68‑1 (MQ=255)
       gAACTTTACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTAcc                                                            <  1:31995527/68‑1 (MQ=255)
         aCATCACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACcca                                                          <  1:7192780/68‑1 (MQ=255)
               cGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGAAGAGGCTACCCACCAGGc                                                    >  1:13271445/1‑68 (MQ=18)
                gCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCAGCGCCTACCCACCAGGCg                                                   >  1:17330302/1‑68 (MQ=255)
                 ccATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGCCATTTAGCAGCGCCTACCCACCAGGCGt                                                  >  1:11665102/1‑68 (MQ=25)
                  cATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGtt                                                 <  1:10001066/68‑1 (MQ=255)
                  cATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGt                                                  <  1:24520544/67‑1 (MQ=255)
                    tGGGTTATTGGATTCTTGTTCTTTGGCATCGCCTGGGCAATTTTCAGCGCCTACCCACCAGGCGTTCt                                               <  1:3322080/68‑1 (MQ=21)
                    tGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt                                               <  1:18136577/68‑1 (MQ=255)
                    tGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt                                               <  1:2887900/68‑1 (MQ=255)
                    tGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt                                               <  1:21830079/68‑1 (MQ=255)
                    tGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt                                               <  1:15054183/68‑1 (MQ=255)
                      ggTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGACGCGCCTACCCACCAGGCTTTTTta                                             >  1:8518406/1‑68 (MQ=37)
                          aCTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATaaa                                         >  1:19899486/1‑68 (MQ=255)
                            tGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTTTGATAAAtt                                       >  1:5926361/1‑67 (MQ=39)
                            tGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCATTTTtatataa                                        >  1:15453224/1‑64 (MQ=37)
                            tGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCATTATTATaaaa                                        >  1:24082184/1‑66 (MQ=38)
                            tGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGAAGAGCCTACCCACCAGGCGTTTTAATAAAtt                                       >  1:11819109/1‑67 (MQ=21)
                             ggCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATaaa                                         >  1:27149646/1‑65 (MQ=255)
                             ggCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGATCTTATAAATg                                       >  1:21013209/1‑67 (MQ=255)
                              gCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGAGCCTACCCACCACGCGTTCATATAAATgcc                                     >  1:18783507/1‑68 (MQ=255)
                                 tCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACTAGGCGTTCTTATAAATgccgcc                                  <  1:22207824/68‑1 (MQ=255)
                                 tCGATTCCTTTCGGATCGCCTGGGCAATGTTCAGAGCCTACCCACCAGGCGTTCTTATAAATgccgcc                                  <  1:8202742/68‑1 (MQ=21)
                                  cGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGccgccg                                 <  1:18452805/68‑1 (MQ=255)
                                   tttttCTTTCGCATCGCCTGGGCAATTTTCCGCGGCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:14237339/67‑1 (MQ=39)
                                   ttcttctatcGTATCGCCTGGGCAATTTGACGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:23092128/65‑1 (MQ=38)
                                   gTTTTCTTTGGGATCGCGTGGGCAATTTGCCGAGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:24780082/68‑1 (MQ=37)
                                   gTTTTCTTTGGCATCGCCTGGGCAATTTTCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:10852835/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:21485225/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:25222593/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:11164802/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:22961996/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:3088967/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:9425219/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:3533883/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:15740988/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:7572181/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:6572651/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:12635388/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:20153600/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:19573816/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:4688621/68‑1 (MQ=255)
                                    ttttCTTTGGCAGCGCCTGGGCAATTTTTAGAGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:24936969/67‑1 (MQ=21)
                                    ttttCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGtt                               >  1:29618352/1‑68 (MQ=255)
                                    ttttCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGtt                               >  1:29890189/1‑68 (MQ=255)
                                     tttctttgGGATCGCCTGGGCAATTTTACGCGCCAACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:933380/66‑1 (MQ=25)
                                     tttctttcGCATCGCCTGGTCAATTTTCCGCGCGTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                <  1:26091954/66‑1 (MQ=39)
                                     tttctttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGtt                               <  1:25483978/67‑1 (MQ=255)
                                      ttctttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTtat                             <  1:23483640/68‑1 (MQ=255)
                                        ctttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAg                           <  1:6112715/68‑1 (MQ=255)
                                          ttcGCATCGCCTGTGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         <  1:23910951/68‑1 (MQ=255)
                                          ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         <  1:12192063/68‑1 (MQ=255)
                                          ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         <  1:13098933/68‑1 (MQ=255)
                                          ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCCTTCTTATAAATGCCGCCGTTAAAGAt                         >  1:575415/1‑68 (MQ=255)
                                          ctggCACCGCGTGGGCAATCTTTCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                         <  1:9088349/65‑1 (MQ=18)
                                           tcGCATCGCCTGGGCAATATTACGAGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATa                        <  1:8824147/68‑1 (MQ=37)
                                              taTCGCGTGGGCAATCTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAAc                      <  1:28350699/66‑1 (MQ=38)
                                               ctCGCCTTGTCATTTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                    <  1:9613129/67‑1 (MQ=38)
                                               aTCGCCTGGGCAATTTGCCGGGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                    <  1:7390712/68‑1 (MQ=255)
                                               aTCGCCTGGGCAATTTGCCGCGCCTACCCAGCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                    <  1:26928433/68‑1 (MQ=255)
                                               aTCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                    <  1:27591952/68‑1 (MQ=255)
                                               aTCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                    <  1:3125241/68‑1 (MQ=255)
                                               aTCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                    <  1:15232692/68‑1 (MQ=255)
                                               aTCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                    <  1:6965303/68‑1 (MQ=255)
                                               aTCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                    <  1:16650429/68‑1 (MQ=255)
                                               aTCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                    <  1:4558896/68‑1 (MQ=255)
                                               aTCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTGATAAATGCCGCCGTTATAGATAACGt                    >  1:2044912/1‑68 (MQ=255)
                                                  gcgtGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTAcc                 <  1:11992727/65‑1 (MQ=255)
                                                  gCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTAcc                 <  1:20756166/68‑1 (MQ=255)
                                                   ccTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCa                >  1:5396787/1‑68 (MQ=255)
                                                   ccTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTATTATAAATCCCGCCGTTATAGAAAACGTACCa                >  1:5625465/1‑68 (MQ=255)
                                                    cTGGGCAATTTTAAGAGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:31767942/68‑1 (MQ=37)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAATTGCCGCCGTTATAGATAACGTACCAt               <  1:17559209/68‑1 (MQ=255)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:1150743/68‑1 (MQ=255)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:21142542/68‑1 (MQ=255)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:25924053/68‑1 (MQ=255)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:26705087/68‑1 (MQ=255)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:30794016/68‑1 (MQ=255)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt               <  1:10220884/68‑1 (MQ=255)
                                                       ggCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTa            >  1:21305945/1‑68 (MQ=255)
                                                        gCAATTTTCCGAGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTa            <  1:14486213/67‑1 (MQ=255)
                                                        gCAATTTTCAGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTa            <  1:28514339/67‑1 (MQ=255)
                                                        gCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTa            <  1:4972912/67‑1 (MQ=255)
                                                        gCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTa            <  1:13481593/67‑1 (MQ=255)
                                                         cAATTTGCCGCGCCTACCGACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTag           <  1:19172051/67‑1 (MQ=255)
                                                         cAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATATCGTACCATGTaca          >  1:29339284/1‑66 (MQ=255)
                                                         cAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTaga          >  1:4078131/1‑68 (MQ=255)
                                                         cAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTag           <  1:16031364/67‑1 (MQ=255)
                                                             ttGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTagaaga       >  1:10560726/1‑67 (MQ=255)
                                                             ttGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCAGCCGTTTTAGATAACGTGCGATGTataagc       >  1:7861269/1‑66 (MQ=37)
                                                              tGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAgaagaa      >  1:31539120/1‑67 (MQ=255)
                                                                 cgcgCCTACCCACCGGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCtt  <  1:29119453/68‑1 (MQ=255)
                                                                 cgcgCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCtt  <  1:5271058/68‑1 (MQ=255)
                                                                 cgcgCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCt   >  1:19050585/1‑67 (MQ=255)
                                                                 cgcgCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCt   <  1:30167656/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 < 15 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: