Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 417,885 T→G 10.6% 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,8850TG10.6% 2745.6 / 67.5 893intergenic (+70/‑89)proY/malZproline‑specific permease/maltodextrin glucosidase
Reads supporting (aligned to +/- strand):  ref base T (31/755);  new base G (99/7);  total (131/762)
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 1008 total aligned reads displayed.

CTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGAT  >  CP009273/417818‑417952
                                                                   |                                                                   
catGCGTGCTGTAAGCCCGTTCAGGAGGAAATCCAGATCCGGCTATTCCGAAAGTTATCCGCCCCCGt                                                                     <  1:43489986/66‑1 (MQ=255)
 ttGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGtc                                                                    <  1:37381129/68‑1 (MQ=255)
   gCGTTCTGTAAGCCGGTTCAGGATGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGtcct                                                                  <  1:55661241/68‑1 (MQ=255)
    cGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGgccg                                                                  >  1:45307886/1‑66 (MQ=255)
     gTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGtccgga                                                                >  1:36288994/1‑65 (MQ=255)
     gTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGgccgcc                                                                >  1:50519893/1‑68 (MQ=255)
      ttCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGgcctcc                                                                >  1:25582052/1‑67 (MQ=255)
        cTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGtcctcctcc                                                             <  1:20699639/68‑1 (MQ=255)
        cTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGtcctcctcc                                                             <  1:19970879/68‑1 (MQ=255)
        cTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGtcctcctcc                                                             <  1:38148892/68‑1 (MQ=255)
         tGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTccc                                                            <  1:11726517/68‑1 (MQ=255)
          gTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTcccc                                                           <  1:40114034/68‑1 (MQ=255)
          gTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTcccc                                                           <  1:33713262/68‑1 (MQ=255)
          gTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTcccc                                                           <  1:10079048/68‑1 (MQ=255)
          gTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTcccc                                                           <  1:15450881/68‑1 (MQ=255)
            aaGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCa                                                         <  1:29822866/68‑1 (MQ=255)
            aaGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCa                                                         <  1:214638/68‑1 (MQ=255)
               ccGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGGCGCCGGCGGCCTCGCCCaac                                                       >  1:3028927/1‑66 (MQ=255)
                   ttCAGGAGGAAATCCTGATCCGGCTATTCCGAATGTTATCCGCCCCCGTCCTCCTCCCCCAAATATc                                                   <  1:38714305/67‑1 (MQ=255)
                   ttCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATc                                                   <  1:41227963/67‑1 (MQ=255)
                   ttCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATc                                                   <  1:52993478/67‑1 (MQ=255)
                   ttCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATc                                                   <  1:49955549/67‑1 (MQ=255)
                   ttCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATc                                                   <  1:8027834/67‑1 (MQ=255)
                    tCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATcgt                                                 >  1:9430773/1‑66 (MQ=255)
                    tCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCt                                                 <  1:14116463/68‑1 (MQ=255)
                    tCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCt                                                 <  1:22729890/68‑1 (MQ=255)
                     cAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCtt                                                <  1:8661348/68‑1 (MQ=255)
                         aggaAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCa                                              <  1:13783985/66‑1 (MQ=255)
                         aggaAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCGCCGCCCCCAAATATCCGTCAga                                            >  1:16959077/1‑68 (MQ=255)
                          ggaAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAga                                            <  1:45206401/67‑1 (MQ=255)
                          ggaAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAga                                            <  1:53740781/67‑1 (MQ=255)
                          ggaAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAga                                            <  1:15131460/67‑1 (MQ=255)
                          ggaAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAga                                            <  1:7526536/67‑1 (MQ=255)
                          ggaAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAga                                            <  1:23478955/67‑1 (MQ=255)
                             aaTCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAgatgat                                        <  1:48768185/68‑1 (MQ=255)
                              aTCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGatgatg                                       <  1:25109359/68‑1 (MQ=255)
                              aTCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGatgatg                                       <  1:3567704/68‑1 (MQ=255)
                              aTCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGatgatg                                       <  1:4774542/68‑1 (MQ=255)
                              aTCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCCGCCTCCTCCCCCAAAATTCGTTGAGatcatg                                       >  1:29173208/1‑68 (MQ=255)
                               tCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGAtgatga                                      <  1:19095189/68‑1 (MQ=255)
                               tCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGAtgatga                                      <  1:12153919/68‑1 (MQ=255)
                                 atGATCCGGCTATTCCTAAAGTTATCCGCCCCCGTGCTCCTCCCCCAAATATCCTTCATATGATGAGt                                    <  1:27261140/67‑1 (MQ=255)
                                  tGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTg                                   <  1:31799864/68‑1 (MQ=255)
                                  tGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCGCCGCCCCCAAATATAATTCACATGATGGGTg                                   >  1:4664266/1‑68 (MQ=255)
                                  tGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCGCCGCCCCCCCATATCAATAAGAGGATGCGTg                                   >  1:174468/1‑68 (MQ=255)
                                    aTCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCGCCTCCCCCAAATATCCTTCATATGATGAGTGAt                                 >  1:39126416/1‑68 (MQ=255)
                                     tCCGGCTATTCCGAAAGTTATCCGGGCCCGGCCTCCGCCCCCAAATATGATTCAGATGGTGAGTggt                                 >  1:13171783/1‑65 (MQ=255)
                                     tCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAACTATCCTTCAGATGATGAGTGATc                                <  1:16140773/68‑1 (MQ=255)
                                     tCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGAt                                 <  1:48460886/67‑1 (MQ=255)
                                     tCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGAt                                 <  1:1661127/67‑1 (MQ=255)
                                     tCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGAt                                 <  1:12464372/67‑1 (MQ=255)
                                     tCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGAt                                 <  1:26565236/67‑1 (MQ=255)
                                     tCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGAt                                 <  1:6998270/67‑1 (MQ=255)
                                     tCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGAt                                 <  1:21096992/67‑1 (MQ=255)
                                     tCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGAt                                 <  1:34668528/67‑1 (MQ=255)
                                     tCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGAt                                 <  1:33950974/67‑1 (MQ=255)
                                     tCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGAt                                 <  1:56986691/67‑1 (MQ=255)
                                     tCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGAt                                 <  1:33270041/67‑1 (MQ=255)
                                     tCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGAt                                 <  1:32940371/67‑1 (MQ=255)
                                     tCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGAt                                 <  1:32054470/67‑1 (MQ=255)
                                     tCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGAt                                 <  1:47268866/67‑1 (MQ=255)
                                     tCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGAt                                 <  1:56212171/67‑1 (MQ=255)
                                     tCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCGCCGCCCCCAAATATCCTTCAGATGATGAGTggt                                 >  1:20321226/1‑65 (MQ=255)
                                          cTATTCCGAATGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGc                            <  1:54314014/67‑1 (MQ=255)
                                          cTATTCCGAAAGTTTTCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGc                            <  1:54704066/67‑1 (MQ=255)
                                          cTATTCCGAAAGTTATCCGCCGCCGTCCTCCGCCCCCAAATATCCTTCAGATGATGAGTGATGCTGCa                           >  1:53414509/1‑68 (MQ=255)
                                          cTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGc                            <  1:44079842/67‑1 (MQ=255)
                                          cTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGc                            <  1:44672662/67‑1 (MQ=255)
                                          cTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGc                            <  1:5180342/67‑1 (MQ=255)
                                          cTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGc                            <  1:6332806/67‑1 (MQ=255)
                                          cTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGc                            <  1:49453685/67‑1 (MQ=255)
                                          cTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGc                            <  1:4613581/67‑1 (MQ=255)
                                          cTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGc                            <  1:36706029/67‑1 (MQ=255)
                                          cTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGc                            <  1:2695166/67‑1 (MQ=255)
                                          cTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGc                            <  1:26204422/67‑1 (MQ=255)
                                           tATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCa                           <  1:51195561/67‑1 (MQ=255)
                                              tCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAg                       <  1:35199897/68‑1 (MQ=255)
                                              tCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAg                       <  1:41922539/68‑1 (MQ=255)
                                              tCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAg                       <  1:28503425/68‑1 (MQ=255)
                                              tCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAg                       <  1:22206081/68‑1 (MQ=255)
                                              tCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATACTTCAGATGATGAGTGAGCCGGCATTAg                       >  1:39628390/1‑68 (MQ=255)
                                               ccGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAgg                      <  1:42840452/68‑1 (MQ=255)
                                               ccGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAgg                      <  1:11126981/68‑1 (MQ=255)
                                               ccGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAgg                      <  1:14657407/68‑1 (MQ=255)
                                               ccGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAgg                      <  1:17872335/68‑1 (MQ=255)
                                               ccGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAgg                      <  1:57929868/68‑1 (MQ=255)
                                               ccGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAgg                      <  1:24569277/68‑1 (MQ=255)
                                               ccGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAgg                      <  1:45729025/68‑1 (MQ=255)
                                                cGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGc                     <  1:55189499/68‑1 (MQ=255)
                                                   aaGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTAt                  <  1:24161496/68‑1 (MQ=255)
                                                   aaGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTAt                  <  1:30938652/68‑1 (MQ=255)
                                                       tATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGc               <  1:23187903/67‑1 (MQ=255)
                                                          ccGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAgg           <  1:52480364/68‑1 (MQ=255)
                                                           cGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGt          <  1:56613947/68‑1 (MQ=255)
                                                           cGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGt          <  1:2786395/68‑1 (MQ=255)
                                                           cGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGt          <  1:37726261/68‑1 (MQ=255)
                                                             cccccGTCCTCCGCCCCCAAATATCCTTCAGGTGATTATTGATCCTGGAGTTGGCTATGGCAAGGTGa        >  1:19341/1‑68 (MQ=255)
                                                              ccccGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGa        <  1:10465482/67‑1 (MQ=255)
                                                                  gTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGa   <  1:18531656/68‑1 (MQ=255)
                                                                  gTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGa   <  1:40563904/68‑1 (MQ=255)
                                                                   |                                                                   
CTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGAT  >  CP009273/417818‑417952

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: