Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,331,634 C→T A378V (GCG→GTG)  proP → proline/glycine betaine transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,331,6340CT98.5% 208.0 / NA 66A378V (GCG→GTG) proPproline/glycine betaine transporter
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (34/31);  total (34/32)
Fisher's exact test for biased strand distribution p-value = 4.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTGCCCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGCGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATT  >  NC_000913/4331498‑4331768
                                                                                                                                        |                                                                                                                                      
gtTGTCCTGTTTGTGTTGGCGATCCCGGCGGTTATTCTGATTAACAGTAACGGCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGa                                                                                                                                      <  2:487153/139‑1 (MQ=255)
 ttGCCCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGAt                                                                                                                                     >  2:913107/1‑139 (MQ=255)
 ttGCCCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGAt                                                                                                                                     <  1:47725/139‑1 (MQ=255)
  tGCCCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATg                                                                                                                                    <  1:248519/139‑1 (MQ=255)
    cccTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATAGCTTCTACCTTGCCAGTGATGtt                                                                                                                                  >  1:673740/1‑139 (MQ=255)
     ccTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTc                                                                                                                                 >  2:181288/1‑139 (MQ=255)
      cTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATAGCTTCTACCTTGCCAGTGATGTTcc                                                                                                                                <  2:673740/139‑1 (MQ=255)
          ttGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGAc                                                                                                                            <  2:515103/139‑1 (MQ=255)
            gtgtTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGc                                                                                                                          <  1:913107/139‑1 (MQ=255)
            gtgtTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGc                                                                                                                          <  1:763870/139‑1 (MQ=255)
                 ggCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTcc                                                                                                                                >  1:535545/1‑128 (MQ=255)
                 ggCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTcc                                                                                                                                <  2:535545/128‑1 (MQ=255)
                   cGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCg                                                                                                                   >  1:1146047/1‑139 (MQ=255)
                                    cTGATTAACAGTATCGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGCGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcgg                                                                                                  <  1:474128/139‑1 (MQ=255)
                                       aTTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcggcgg                                                                                               >  2:364088/1‑139 (MQ=255)
                                        ttAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcggcggc                                                                                              <  2:358580/139‑1 (MQ=255)
                                        ttAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcggcggc                                                                                              >  2:979651/1‑139 (MQ=255)
                                                 aCGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcggcg                                                                                                <  2:948449/128‑1 (MQ=255)
                                                 aCGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcggcg                                                                                                >  1:948449/1‑128 (MQ=255)
                                                    tCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTc                                                                                  >  2:1052856/1‑139 (MQ=255)
                                                     cATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCg                                                                                 >  1:382147/1‑139 (MQ=255)
                                                       tCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGt                                                                               >  1:1173352/1‑139 (MQ=255)
                                                         ggCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGc                                                                             >  1:1218284/1‑139 (MQ=255)
                                                             tGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTgtt                                                                         >  1:150572/1‑137 (MQ=255)
                                                                         ggTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGAc                                                             >  1:171339/1‑139 (MQ=255)
                                                                                    cTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGcc                                                                     >  1:4867/1‑120 (MQ=255)
                                                                                    cTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGcc                                                                     <  2:4867/120‑1 (MQ=255)
                                                                                      ggCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTggcg                                                <  1:181288/139‑1 (MQ=255)
                                                                                        cGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTggcggc                                              >  1:777305/1‑139 (MQ=255)
                                                                                         ggTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGcc                                             >  2:864515/1‑139 (MQ=255)
                                                                                             aTCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCctggc                                         >  1:583122/1‑139 (MQ=255)
                                                                                              tCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcggcg                                                                                                <  2:878297/83‑1 (MQ=255)
                                                                                              tCCTTAACTACTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcggcg                                                                                                >  1:878297/1‑83 (MQ=255)
                                                                                                       gCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCat                                                                                                                        <  1:172986/50‑1 (MQ=255)
                                                                                                       gCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCat                                                                                                                        >  2:172986/1‑50 (MQ=255)
                                                                                                       gCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAg                               >  1:1187015/1‑139 (MQ=255)
                                                                                                       gCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAg                               >  1:1286226/1‑139 (MQ=255)
                                                                                                        cTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGc                                                                                                                          <  1:741387/47‑1 (MQ=255)
                                                                                                        cTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGc                                                                                                                          >  2:741387/1‑47 (MQ=255)
                                                                                                        cTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcggcg                                                                                                >  1:1128181/1‑73 (MQ=255)
                                                                                                        cTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcggcg                                                                                                <  2:1128181/73‑1 (MQ=255)
                                                                                                        cTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGc                              >  2:392731/1‑139 (MQ=255)
                                                                                                        cTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGc                              >  2:281782/1‑139 (MQ=255)
                                                                                                          ttACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTc                            >  1:74890/1‑139 (MQ=255)
                                                                                                          ttACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTc                            >  1:549162/1‑139 (MQ=255)
                                                                                                           tACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCg                           <  2:1218284/139‑1 (MQ=255)
                                                                                                           tACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCg                           <  1:716540/139‑1 (MQ=255)
                                                                                                            aCGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGc                          <  2:1286226/139‑1 (MQ=255)
                                                                                                             cGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTccc                                                                                                                               <  2:1192042/37‑1 (MQ=255)
                                                                                                             cGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTccc                                                                                                                               >  1:1192042/1‑37 (MQ=255)
                                                                                                             cGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCa                         <  2:852330/139‑1 (MQ=255)
                                                                                                                  gTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGc                                                                                                                          <  1:692444/37‑1 (MQ=255)
                                                                                                                  gTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGc                                                                                                                          >  2:692444/1‑37 (MQ=255)
                                                                                                                  gTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCt                                                                <  1:736604/95‑1 (MQ=255)
                                                                                                                  gTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCt                                                                >  2:736604/1‑95 (MQ=255)
                                                                                                                      tGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCtgat                <  1:471458/139‑1 (MQ=255)
                                                                                                                           tCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTgatgattc           <  2:777305/139‑3 (MQ=255)
                                                                                                                           tCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGAtgc           <  2:1146047/139‑1 (MQ=255)
                                                                                                                            cTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGAtgcc          <  1:979651/139‑1 (MQ=255)
                                                                                                                            cTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGAtgcc          <  1:392731/139‑1 (MQ=255)
                                                                                                                            cTACCTTGCCAGTGATGTTCCCCACGGATTTCCGTTTCAGCGCGCCGGCGGCGGCATTTAATATTTTGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGAtgcc          <  2:150572/139‑1 (MQ=255)
                                                                                                                                  tGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGAtgcct         <  2:748887/134‑1 (MQ=255)
                                                                                                                                  tGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGAtgcct         >  1:748887/1‑134 (MQ=255)
                                                                                                                                   gCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTgg                                                  <  2:136256/92‑1 (MQ=255)
                                                                                                                                   gCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTgg                                                  >  1:136256/1‑92 (MQ=255)
                                                                                                                                   gCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTAt   <  1:281782/139‑1 (MQ=255)
                                                                                                                                    ccAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTAtt  >  1:1106811/1‑139 (MQ=255)
                                                                                                                                        |                                                                                                                                      
GTTGCCCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGCGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATT  >  NC_000913/4331498‑4331768

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: