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% 1920.8 / 12.6 636A259S (GCA→TCA) chiAperiplasmic endochitinase
Reads supporting (aligned to +/- strand):  ref base C (135/432);  new base A (19/25);  total (156/480)
Fisher's exact test for biased strand distribution p-value = 6.57e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 638 total aligned reads displayed.

GCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCT  >  CP009273/3462374‑3462505
                                                                |                                                                   
  cgcgcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCcgcg                                                                 <  1:27996404/67‑1 (MQ=255)
     gcgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTa                                                             <  1:12851194/68‑1 (MQ=255)
       gAACTTTACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTAcc                                                           <  1:31995527/68‑1 (MQ=255)
       gAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTAcc                                                           <  1:25230450/68‑1 (MQ=255)
              aCGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAgg                                                    <  1:26752052/68‑1 (MQ=255)
               cGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGc                                                   <  1:23489715/68‑1 (MQ=255)
               cGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCAGGGCAATTTGCCGCGCCTACCCACCAGGc                                                   >  1:76517/1‑68 (MQ=255)
                gCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCAGCGCCTACCCACCAGGCg                                                  >  1:17330302/1‑68 (MQ=255)
                 ccATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGCCATTTAGCAGCGCCTACCCACCAGGCGt                                                 >  1:11665102/1‑68 (MQ=25)
                  cATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGtt                                                <  1:10001066/68‑1 (MQ=255)
                  cATGGGTTACTGGATTCGTTTGCTATCGCATCGCCTGGGCAATTTGCAGCGCCTACCCACCAGGCGt                                                 <  1:12531594/67‑1 (MQ=255)
                    tGGGTTATTGGATTCTTGTTCTTTGGCATCGCCTGGGCAATTTTCAGCGCCTACCCACCAGGCGTTCt                                              <  1:3322080/68‑1 (MQ=21)
                    tGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt                                              <  1:2887900/68‑1 (MQ=255)
                    tGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt                                              >  1:21877675/1‑68 (MQ=255)
                    tGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt                                              <  1:30462873/68‑1 (MQ=255)
                      ggTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGACGCGCCTACCCACCAGGCTTTTTta                                            >  1:8518406/1‑68 (MQ=37)
                      ggTTACTGGCATCGTTTTCTTTCGCATCGCCAGGGCAGTTTGAAGAGGCTACCCACCAGGCATTTTta                                            >  1:19921427/1‑68 (MQ=11)
                       gTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTta                                            <  1:29122450/67‑1 (MQ=255)
                       gTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGAAGAGCCTACCCACCAGGCGTTCTta                                            <  1:29382904/67‑1 (MQ=37)
                       gTTACTGGACTCGTCTTCTATCGCACCCCCTGGGCAATCTTACGCGCCTACCCACCAGGCGTTCt                                              <  1:20773077/65‑1 (MQ=16)
                       gTTAATGGCTTCTTCTTCTATGGCACCGCGTGGGCAATCTTCCGCCCCTACCCACCAGGCGTTCTta                                            <  1:11015848/67‑1 (MQ=9)
                            tGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCATTTTtatataa                                       >  1:15453224/1‑64 (MQ=37)
                            tGGCATCATTTTCTTTCGAATCCCCTGGCCAATTAGACGCGCCTACCCACCAGCCGTTATTATaaaa                                       >  1:8840627/1‑66 (MQ=19)
                             ggCATCGTTTTCTTTCGCATCGCCTGGGCAATTTTCCGCGCCTACCCACCAGGCGTTCTTATAAATg                                      >  1:18454442/1‑67 (MQ=255)
                               cATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCCTTGTCATAAATgaa                                    >  1:20470530/1‑65 (MQ=255)
                                 tCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACTAGGCGTTCTTATAAATgccgcc                                 <  1:22207824/68‑1 (MQ=255)
                                 tCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACTAGGCGTTCTTATAAATgccgcc                                 <  1:10546552/68‑1 (MQ=255)
                                 tCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATgccgcc                                 <  1:586247/68‑1 (MQ=255)
                                 tCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATgccgcc                                 <  1:15676710/68‑1 (MQ=255)
                                 tCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATgccgcc                                 <  1:8283570/68‑1 (MQ=255)
                                   ttattctttgGCATCGCCTGGGCAATTTGACGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:13882485/65‑1 (MQ=39)
                                   gtcttctttcGCACCGCCTGGGCAATCTTCCGCGCCAACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:6130477/65‑1 (MQ=25)
                                   gTTTTCTTTGGCATCGCCTGGGCAATTTTCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:10852835/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:25489524/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:5014183/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:25956075/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:2635292/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:5561995/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:3088967/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:26967947/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:27244357/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:6755889/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:28366708/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:747706/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:3533883/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:22961996/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:29649442/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:23918547/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:14445239/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:21161635/68‑1 (MQ=255)
                                   gTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:20153600/68‑1 (MQ=255)
                                    ttttCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGtt                              >  1:29890189/1‑68 (MQ=255)
                                    ttttCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:1178640/67‑1 (MQ=255)
                                     tttctttgGGATCGCCTGGGCAATTTTACGCGCCAACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:933380/66‑1 (MQ=25)
                                     cttctatcGCATCCCCTGGGCAATTTTTCGAGCCAACCCACCAGGCGTTCTTATAAATGCCGCCGt                               <  1:19098733/65‑1 (MQ=18)
                                      ttctttcGCATCGCCTGGGCAATTTTCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTtat                            <  1:1148660/68‑1 (MQ=255)
                                      ttctttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTtat                            <  1:6021613/68‑1 (MQ=255)
                                      ttctttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTtat                            <  1:24083725/68‑1 (MQ=255)
                                      ttctttcGCATCGCCCGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTtat                            <  1:14223554/68‑1 (MQ=255)
                                        ctttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAg                          <  1:19602261/68‑1 (MQ=255)
                                          ttggCATCGCCTGGGCAATTTGACGCGCCAACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                        <  1:27626967/65‑1 (MQ=38)
                                          ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                        <  1:24817734/68‑1 (MQ=255)
                                          ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                        <  1:5271665/68‑1 (MQ=255)
                                          ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                        <  1:12192063/68‑1 (MQ=255)
                                          ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                        <  1:13098933/68‑1 (MQ=255)
                                          ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                        <  1:21032175/68‑1 (MQ=255)
                                          ttcGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                        <  1:18141347/68‑1 (MQ=255)
                                          ctggCACCGCGTGGGCAATCTTTCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                        <  1:9088349/65‑1 (MQ=18)
                                               ctCGCCTTGTCATTTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                   <  1:9613129/67‑1 (MQ=38)
                                               aTCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                   <  1:19288523/68‑1 (MQ=255)
                                               aTCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                   <  1:3125241/68‑1 (MQ=255)
                                               aTCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                   <  1:23105550/68‑1 (MQ=255)
                                               aTCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                   <  1:15232692/68‑1 (MQ=255)
                                               aTCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                   <  1:4558896/68‑1 (MQ=255)
                                               aTCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                   <  1:17718557/68‑1 (MQ=255)
                                               aTCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                   <  1:16650429/68‑1 (MQ=255)
                                               aTCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGt                   <  1:18783924/68‑1 (MQ=255)
                                                 cGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTAc                 <  1:30802494/68‑1 (MQ=255)
                                                  gcgtGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTAcc                <  1:11992727/65‑1 (MQ=255)
                                                  gCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTAcc                <  1:7869720/68‑1 (MQ=255)
                                                    cTGGGCTATTTGACGCGCCTACCCACCAGGCGTTCTTATAAATTCCGCCGTTATAGATAACGTACCAt              <  1:9426965/68‑1 (MQ=39)
                                                    cTGGGCAATTTTTAGAGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt              <  1:7082299/68‑1 (MQ=25)
                                                    cTGGGCAATTTTAAGAGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt              <  1:31767942/68‑1 (MQ=37)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAATTGCCGCCGTTATAGATAACGTACCAt              <  1:17559209/68‑1 (MQ=255)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt              <  1:719298/68‑1 (MQ=255)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt              <  1:10220884/68‑1 (MQ=255)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt              <  1:5669989/68‑1 (MQ=255)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt              <  1:4691996/68‑1 (MQ=255)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt              <  1:2453631/68‑1 (MQ=255)
                                                    cTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCAt              <  1:132254/68‑1 (MQ=255)
                                                        gCAATTTTCAGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTa           <  1:28514339/67‑1 (MQ=255)
                                                        gCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTa           <  1:13481593/67‑1 (MQ=255)
                                                         cAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTaga         >  1:4078131/1‑68 (MQ=255)
                                                         cAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTag          <  1:16031364/67‑1 (MQ=255)
                                                            tttGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTagaaga      <  1:26093313/68‑1 (MQ=255)
                                                             ttGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTagaaga      >  1:21582487/1‑67 (MQ=255)
                                                             ttGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTagaaga      >  1:15050715/1‑67 (MQ=255)
                                                             ttGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTagaaca      >  1:2134275/1‑65 (MQ=255)
                                                              tGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAgaagaa     >  1:31539120/1‑67 (MQ=255)
                                                               gCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAAcc   <  1:25011364/68‑1 (MQ=255)
                                                                |                                                                   
GCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCT  >  CP009273/3462374‑3462505

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: