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.2% 180.4 / NA 57A378V (GCG→GTG) proPproline/glycine betaine transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (31/25);  minor base G (0/1);  total (31/26)
Fisher's exact test for biased strand distribution p-value = 4.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGCCCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGCGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTA  >  NC_000913/4331499‑4331769
                                                                                                                                       |                                                                                                                                       
ttGCCCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGAt                                                                                                                                      >  2:965575/1‑139 (MQ=255)
                  cGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCg                                                                                                                    >  2:720513/1‑139 (MQ=255)
                               tATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATTCTGGCGGTGATCCTTAACTGCATTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCAGTTACAGCGCGCTg                                                                                                       >  2:901022/1‑139 (MQ=255)
                               tATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTg                                                                                                       >  1:1107869/1‑139 (MQ=255)
                               tATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTg                                                                                                       >  1:540595/1‑139 (MQ=255)
                                   cTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcgg                                                                                                   <  1:720513/139‑1 (MQ=255)
                                    tGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcggc                                                                                                  <  1:78528/139‑1 (MQ=255)
                                     gATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcggcg                                                                                                 <  2:1174119/139‑1 (MQ=255)
                                              tAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCg                                                                                                                               >  2:729072/1‑100 (MQ=255)
                                              tAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGGGATGTTCCCg                                                                                                                               <  1:729072/100‑1 (MQ=255)
                                                aCGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGc                                                                                                                           >  1:249351/1‑102 (MQ=255)
                                                aCGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGc                                                                                                                           <  2:249351/102‑1 (MQ=255)
                                                aCGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcggcg                                                                                                 <  2:1409497/128‑1 (MQ=255)
                                                aCGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcggcg                                                                                                 >  1:1409497/1‑128 (MQ=255)
                                                aCGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAatat                                                                                      >  2:428357/1‑139 (MQ=255)
                                                aCGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAatat                                                                                      >  2:1014690/1‑139 (MQ=255)
                                                      tCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCATCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGt                                                                                >  1:858139/1‑139 (MQ=255)
                                                        ggCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGc                                                                              >  2:77283/1‑139 (MQ=255)
                                                                   tGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGt                                                                   >  2:1277029/1‑139 (MQ=255)
                                                                    gCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTg                                                                               >  2:944786/1‑126 (MQ=255)
                                                                    gCCGGGTTACTGATGCTGGCGGTGATCCTTAAATGCTTTACGGGCGTTATTGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTg                                                                               <  1:944786/126‑1 (MQ=255)
                                                                        ggTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCt                                                                 >  2:1224521/1‑136 (MQ=255)
                                                                        ggTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCt                                                                 <  1:1224521/136‑1 (MQ=255)
                                                                            aCTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCa                                                          <  1:965575/139‑1 (MQ=255)
                                                                                  gCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTg                                                    >  2:314419/1‑139 (MQ=255)
                                                                                   cTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTgg                                                   >  2:308650/1‑139 (MQ=255)
                                                                                   cTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGATTGCCGGTCTGACGCCAACGCTgg                                                   >  2:818920/1‑139 (MQ=255)
                                                                                     ggCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGc                                                                                                                           <  2:878800/65‑1 (MQ=255)
                                                                                     ggCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGc                                                                                                                           >  1:878800/1‑65 (MQ=255)
                                                                                     ggCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTggcg                                                 >  2:1063368/1‑139 (MQ=255)
                                                                                      gCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTggcgg                                                >  2:103890/1‑139 (MQ=255)
                                                                                       cGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGt                                                                   >  1:208746/1‑119 (MQ=255)
                                                                                       cGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGt                                                                   <  2:208746/119‑1 (MQ=255)
                                                                                                       cTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGc                               >  1:1426397/1‑139 (MQ=255)
                                                                                                       cTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGc                               >  2:346289/1‑139 (MQ=255)
                                                                                                        tttACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcgg                                                                                                   >  2:1356514/1‑70 (MQ=255)
                                                                                                        tttACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcgg                                                                                                   <  1:1356514/70‑1 (MQ=255)
                                                                                                        tttACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCt                              >  2:1300321/1‑139 (MQ=255)
                                                                                                         ttACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGc                                                                                                                           >  1:522072/1‑45 (MQ=255)
                                                                                                         ttACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGc                                                                                                                           <  2:522072/45‑1 (MQ=255)
                                                                                                            cGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCa                          <  1:1277029/139‑1 (MQ=255)
                                                                                                            cGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCa                          <  1:77283/139‑1 (MQ=255)
                                                                                                                 gTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGt                                      <  1:41330/122‑1 (MQ=255)
                                                                                                                 gTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGt                                      >  2:41330/1‑122 (MQ=255)
                                                                                                                    aTGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCtga                  <  1:1014690/139‑1 (MQ=255)
                                                                                                                    aTGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCtga                  <  1:818920/139‑1 (MQ=255)
                                                                                                                          tCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTgatgatac            <  2:858139/139‑3 (MQ=255)
                                                                                                                          tCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGAtgc            <  1:1300321/139‑1 (MQ=255)
                                                                                                                          tCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGAtgc            <  1:346289/139‑1 (MQ=255)
                                                                                                                           cTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCt                                                                 <  1:1293903/85‑1 (MQ=255)
                                                                                                                           cTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCt                                                                 >  2:1293903/1‑85 (MQ=255)
                                                                                                                           cTACCTTGCCAGTGATGTTCCCGACGCATATCAGTTACAGCGCGCTGGCGGCGGCATTTACTATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCGGGCTGGTCGAAAGCTCGCAGAATCTGATGAtgcc           <  1:901022/139‑1 (MQ=255)
                                                                                                                              ccTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGAtgcctgc        <  2:1407513/139‑1 (MQ=255)
                                                                                                                                 tGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTa     >  1:987067/1‑139 (MQ=255)
                                                                                                                                 tGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTa     >  1:1293673/1‑139 (MQ=255)
                                                                                                                                  gCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTAt    <  1:103890/139‑1 (MQ=255)
                                                                                                                                    cAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTa  <  1:314419/139‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                       
TTGCCCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGCGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTA  >  NC_000913/4331499‑4331769

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