Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 417,885 T→G 13.7% intergenic (+70/‑89) proY → / → malZ proline‑specific permease/maltodextrin glucosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273417,8850TG13.7% 1557.8 / 53.4 542intergenic (+70/‑89)proY/malZproline‑specific permease/maltodextrin glucosidase
Reads supporting (aligned to +/- strand):  ref base T (22/438);  new base G (78/3);  total (101/441)
Fisher's exact test for biased strand distribution p-value = 0.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 617 total aligned reads displayed.

CTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGAT  >  CP009273/417818‑417952
                                                                   |                                                                   
cTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGt                                                                     <  1:10416828/68‑1 (MQ=255)
cTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGt                                                                     <  1:19995833/68‑1 (MQ=255)
cTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGt                                                                     <  1:5314836/68‑1 (MQ=255)
   gCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGtcct                                                                  <  1:35584383/68‑1 (MQ=255)
     gTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGtcctcc                                                                <  1:19845033/68‑1 (MQ=255)
        cTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGtcctcctcc                                                             <  1:29457657/68‑1 (MQ=255)
        cTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGtcctcctcc                                                             <  1:20239309/68‑1 (MQ=255)
        cTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGtcctcctcc                                                             <  1:29066434/68‑1 (MQ=255)
        cTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGtcctcctcc                                                             <  1:37335772/68‑1 (MQ=255)
        cTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGtcctcctcc                                                             <  1:12937819/68‑1 (MQ=255)
        cTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGtcctcctcc                                                             <  1:7925002/68‑1 (MQ=255)
         tGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGtcctcctcc                                                             <  1:17428291/67‑1 (MQ=255)
         tGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTccc                                                            <  1:27501426/68‑1 (MQ=255)
          gTAAGGCGGTTCAGGAGGTATTGCTGATCCGGGTATTCCGAAAGTTATCCGCCGCCGTCCTCCTcgcc                                                           <  1:12060196/68‑1 (MQ=255)
          gTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTcccc                                                           <  1:23455909/68‑1 (MQ=255)
          gTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTcccc                                                           <  1:29894803/68‑1 (MQ=255)
            aaGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTccccc                                                          <  1:1425138/67‑1 (MQ=255)
            aaGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCa                                                         <  1:19284545/68‑1 (MQ=255)
            aaGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCa                                                         <  1:32292524/68‑1 (MQ=255)
               ccGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGGGGCCGGGCGCCGCCCCCaaag                                                      >  1:6457388/1‑67 (MQ=255)
               ccGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGGCCCCGGCCTCCTCCCCCAAat                                                      >  1:38091770/1‑68 (MQ=255)
               ccGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCGCGTCCTCGTCCCCCCAat                                                      >  1:13472193/1‑68 (MQ=255)
               ccGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGGGTGCTCCCCCACat                                                      >  1:18753973/1‑68 (MQ=255)
               ccGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGGCGCCGCCCCCGAat                                                      >  1:37583480/1‑68 (MQ=255)
               ccGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCGCCGCCCCCCAat                                                      >  1:2285930/1‑68 (MQ=255)
                 ggTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAatat                                                    <  1:17838566/68‑1 (MQ=255)
                   ttCTGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATc                                                   <  1:5830367/67‑1 (MQ=255)
                    tCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATcc                                                  <  1:27027507/67‑1 (MQ=255)
                    tCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCt                                                 <  1:22032457/68‑1 (MQ=255)
                    tCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCt                                                 <  1:28585700/68‑1 (MQ=255)
                          ggaAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAga                                            <  1:28281178/67‑1 (MQ=255)
                          ggaAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAga                                            <  1:28777730/67‑1 (MQ=255)
                          ggaAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAga                                            <  1:26439272/67‑1 (MQ=255)
                          ggaAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAga                                            <  1:34694081/67‑1 (MQ=255)
                          ggaAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAga                                            <  1:18442082/67‑1 (MQ=255)
                          ggaAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAga                                            <  1:17167924/67‑1 (MQ=255)
                             aaTCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAgatgat                                        <  1:12447092/68‑1 (MQ=255)
                             aaTCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAgatgat                                        <  1:7769987/68‑1 (MQ=255)
                             aaTCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAgatga                                         <  1:34074388/67‑1 (MQ=255)
                               tCCTGATCCGGCTATTCCGAAAGTTATCCGCGCCCGGCCGCCGCCCCCAAATATACTTCAGAtgatga                                      >  1:36046709/1‑68 (MQ=255)
                               tCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCGCCGCCCCCAAATATCCAGCAGAggagga                                      >  1:1859645/1‑68 (MQ=255)
                                ccTGATCCGGCTATTCCGAAAGTTATCGGCCGCCGGCCGCCGCCCCCCAAGAAACGTAAGatgatggg                                     >  1:25338422/1‑66 (MQ=255)
                                 cTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCTCCGCCCCCACAGATCCGGCAGCTGATGAGt                                    >  1:4540260/1‑68 (MQ=255)
                                     tCCGGCTATTCCGATAGTTATCCGCCCGCGTCCTCCTCGCCCAAATATCCTTCAGATGATGAGTGAt                                 <  1:875486/67‑1 (MQ=255)
                                     tCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGAt                                 <  1:16803410/67‑1 (MQ=255)
                                     tCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGAt                                 <  1:14067353/67‑1 (MQ=255)
                                     tCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGAt                                 <  1:31113956/67‑1 (MQ=255)
                                     tCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGAt                                 <  1:33209111/67‑1 (MQ=255)
                                     tCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGAt                                 <  1:32644017/67‑1 (MQ=255)
                                     tCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGAt                                 <  1:30564969/67‑1 (MQ=255)
                                     tCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGAt                                 <  1:15317982/67‑1 (MQ=255)
                                     tCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGAt                                 <  1:39073211/67‑1 (MQ=255)
                                     tCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGAt                                 <  1:11757857/67‑1 (MQ=255)
                                     tCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGAt                                 <  1:10988004/67‑1 (MQ=255)
                                     tCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGAt                                 <  1:19053849/67‑1 (MQ=255)
                                     tCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGAt                                 <  1:9229528/67‑1 (MQ=255)
                                     tCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGAt                                 <  1:9609343/67‑1 (MQ=255)
                                     tCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGAt                                 <  1:21743017/67‑1 (MQ=255)
                                     tCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGAt                                 <  1:21125092/67‑1 (MQ=255)
                                     tCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCGCCGCCCCCAAATATACTTCAGATGATGAGTGAt                                 >  1:6159437/1‑67 (MQ=255)
                                     tCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCGCCGCCCCCAAATATCCTTGAGATGATGTGTGAt                                 >  1:8188697/1‑67 (MQ=255)
                                     tCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCGCCGCCCCCAAATATCCTTCAGATGATTAGTGAt                                 >  1:25559589/1‑67 (MQ=255)
                                     tCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCACCGCCCCCAAAAATCCTTCAGATGATGAGGGAt                                 >  1:10001616/1‑67 (MQ=255)
                                     tCCGGCTATTCCGAAAGTTATCCGCCCCCCGGCTGCCCCGCCGATAATACATCAGATGATGAGTGAt                                 >  1:22328994/1‑67 (MQ=255)
                                     tCCGGCTATTCCGAAAGTTATCCCCCCCCGTCCGCCTCCCCCCAATAACCGGCAGTTCATGAGTCAt                                 >  1:2139047/1‑67 (MQ=255)
                                          cTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGc                            <  1:32373420/67‑1 (MQ=255)
                                          cTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGc                            <  1:22666305/67‑1 (MQ=255)
                                          cTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGc                            <  1:15704504/67‑1 (MQ=255)
                                          cTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGc                            <  1:24619942/67‑1 (MQ=255)
                                          cTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGc                            <  1:24850740/67‑1 (MQ=255)
                                          cTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGc                            <  1:34522656/67‑1 (MQ=255)
                                          cTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGc                            <  1:26218309/67‑1 (MQ=255)
                                          cTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGc                            <  1:27822142/67‑1 (MQ=255)
                                          cTATTCCGAAAGTTATCCGCCCCCGGCCTCCTCCCCCAAATATACTTCAGATGATGAGTGTTCCTGCa                           >  1:30855473/1‑68 (MQ=255)
                                          cTATTCCGAAAGTTAACCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGc                            <  1:38382577/67‑1 (MQ=255)
                                          cTATACTGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTCATCCTGc                            <  1:7159748/67‑1 (MQ=255)
                                           tATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATAATGAGTGATCCTGCa                           <  1:36641443/67‑1 (MQ=255)
                                              tCCGAAAGTTATCCGCCCCCGTCCTCCTCCTCGAAATATCCTTCAGATGATGCGTGATCCTGCATTAg                       <  1:7290187/68‑1 (MQ=255)
                                              tCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAg                       <  1:1700594/68‑1 (MQ=255)
                                              tCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAg                       <  1:11492688/68‑1 (MQ=255)
                                              tCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAg                       <  1:4889105/68‑1 (MQ=255)
                                              tCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAg                       <  1:24040377/68‑1 (MQ=255)
                                               ccGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAgg                      <  1:2668754/68‑1 (MQ=255)
                                               ccGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAgg                      <  1:15001147/68‑1 (MQ=255)
                                               ccGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAgg                      <  1:4194391/68‑1 (MQ=255)
                                               ccGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAgg                      <  1:3365157/68‑1 (MQ=255)
                                               ccGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAgg                      <  1:35279640/68‑1 (MQ=255)
                                               ccGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAgg                      <  1:23203115/68‑1 (MQ=255)
                                               ccGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAgg                      <  1:3026509/68‑1 (MQ=255)
                                                   aaGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTAt                  <  1:26061948/68‑1 (MQ=255)
                                                   aaGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTAt                  <  1:14529737/68‑1 (MQ=255)
                                                    aGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATg                 <  1:11122712/68‑1 (MQ=255)
                                                     gTTATCCGCCCCCGGCCGCCTCCCCCAAATAAAATTCAGATGAGGAGTGATACTGCATTAGGCTATgg                >  1:2079530/1‑68 (MQ=255)
                                                       tATCCGCCCCCTTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATgg                <  1:358553/66‑1 (MQ=255)
                                                          ccGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAgg           <  1:19589851/68‑1 (MQ=255)
                                                          ccGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAgg           <  1:16645231/68‑1 (MQ=255)
                                                           cGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGt          <  1:16212185/68‑1 (MQ=255)
                                                              ccccGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGa        <  1:31836913/67‑1 (MQ=255)
                                                                  gTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGa   <  1:13193666/68‑1 (MQ=255)
                                                                   tcctcctccCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGAt  >  1:3874848/1‑68 (MQ=255)
                                                                   |                                                                   
CTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGAT  >  CP009273/417818‑417952

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: