Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 603,393 T→A 8.6% intergenic (+58/+23) pheP → / ← ybdG phenylalanine:H(+) symporter PheP/miniconductance mechanosensitive channel YbdG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913603,3930TA 8.6% 156.1 / 3.9 57intergenic (+58/+23)pheP/ybdGphenylalanine:H(+) symporter PheP/miniconductance mechanosensitive channel YbdG
Reads supporting (aligned to +/- strand):  ref base T (28/23);  new base A (2/3);  total (30/27)
Fisher's exact test for biased strand distribution p-value = 6.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

GGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAATATATCAGCCTGAATGCT  >  NC_000913/603251‑603535
                                                                                                                                              |                                                                                                                                               
ggatgATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAACCCGCTCTCATCCAt                                                                                                                                                 >  1:172929/1‑143 (MQ=255)
  atgatATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACCAAACGCTGTCAACCATCCGAAACAGACCTGATCCATTc                                                                                                                                               >  2:914974/1‑143 (MQ=255)
    gatATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCCTCCGAAACCGCTCTCATCCAATCGa                                                                                                                                             >  1:376890/1‑143 (MQ=255)
    gatATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGa                                                                                                                                             >  1:327014/1‑143 (MQ=255)
       atGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATga                                                                                                                                          <  1:140024/143‑1 (MQ=255)
          cgcTTGTCGGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTTTCATCGAATGGATGagag                                                                                                                                       <  1:57550/143‑1 (MQ=255)
           gcTTGGCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGTCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGAGAGc                                                                                                                                      <  2:532292/143‑1 (MQ=255)
                 tatctaTCTTGCTGCCGGTGTGGATTGTATTCCGGTTGACGGCAGTCAAACCGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAGCGATGTCAACCATCCGAAACCGCTCTCATCGAATGGCTGAGAGCGGtttc                                                                                                                                <  1:837854/138‑2 (MQ=255)
                 cAGCCATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGtttt                                                                                                                                <  1:648711/143‑1 (MQ=255)
                     gacccTGCTGCGGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTa                                                                                                                            <  2:350398/140‑1 (MQ=255)
                       cccggctgcCGGTGCGGATTGTATTCCTGTTTACGGCATTTACAACGCCGCGTCGGAAATAAGGCATGCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAAt                                                                                                                          <  2:311637/139‑1 (MQ=255)
                         ctgctgCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATCCGATGAGAGCGGTTTTTTTACTTa                                                                                                                        >  2:372103/1‑143 (MQ=255)
                                    ggggATTGTATTCCTGTCTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATg                                                                                                             <  1:54694/141‑1 (MQ=255)
                                       gATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAATGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAGTTACTGCTTAAATGCAc                                                                                                          >  1:367190/1‑143 (MQ=255)
                                       gATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTACATGCAc                                                                                                          >  1:804736/1‑143 (MQ=255)
                                               tCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTCCTGCTTAAATGCACCCGCCaga                                                                                                  >  1:546082/1‑143 (MQ=255)
                                                          gCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCACCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATc                                                                                       >  1:535344/1‑143 (MQ=255)
                                                                  aaCGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCCTCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGACgg                                                                              >  2:814389/1‑143 (MQ=255)
                                                                          gagagAAAGAAGGCATTGAGGGTATGTACGACAAGGCGTTGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAAGGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGAc                                                                      <  1:1051799/139‑1 (MQ=255)
                                                                           tCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGCGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGCCGACt                                                                     >  2:1018017/1‑143 (MQ=255)
                                                                             gcgaaTGAGGCATTCACGCTACATCCGACAATACGATGTCAACGCTCCGAAACCGCTCTCGTCCGCGCGGGGCGAGGAAAATAATTGATCACTGCTTAAATGCACCCGCTGGAGAGCGGATAGCAGGGCCGGTTGGCGACtga                                                                   <  2:224592/140‑1 (MQ=255)
                                                                              gacttacggCATTGACGCTCCACTGCACAAATCGCTGTCACCCCTCCGAATCCGCTCTCATCCATTCGATTCTAGCGGTTTTTTTAATTACTGCTTACATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACtgat                                                                  <  2:104615/136‑1 (MQ=255)
                                                                                aaTAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGAACCCGCCAGAGAGCGAATATCATTCCCGGTTTGCGACtgacga                                                                >  2:1032226/1‑143 (MQ=255)
                                                                                 aTAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGaa                                                               <  1:897119/143‑1 (MQ=255)
                                                                                    aGGCATTCACGCTACATCCGACAAAACGATGTGAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGCTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATGATTGCCGGTTGGCGACTGATGAAggc                                                            >  1:679704/1‑141 (MQ=255)
                                                                                     tgccttCTCGCTACTTACGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTGAAATGCACCCGCCAGAGCGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCg                                                           <  2:363331/139‑1 (MQ=255)
                                                                                     ggCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGc                                                                                  >  2:251351/1‑121 (MQ=255)
                                                                                     ggCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGc                                                                                  <  1:251351/121‑1 (MQ=255)
                                                                                       cATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCa                                                         <  1:995762/143‑1 (MQ=255)
                                                                                        aTTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGccg                                                        >  1:787660/1‑141 (MQ=255)
                                                                                            aCGCTACATCCGACAAAACGCTGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCa                                                    >  1:609617/1‑143 (MQ=255)
                                                                                               cTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAAc                                                 <  1:654656/143‑1 (MQ=255)
                                                                                                    tCCTACAAAACGATGTCAACCATCCGAACCTGCTCTCATCCATTCGATGAGAGCGGTTTTTTTATTTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTc                                            <  2:102007/143‑1 (MQ=255)
                                                                                                    tCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTc                                            <  1:432041/143‑1 (MQ=255)
                                                                                                    tCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATCAAGCCGCAGACCCAActcctc                                            >  2:962269/1‑143 (MQ=255)
                                                                                                    tCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCGTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTc                                            <  2:172929/143‑1 (MQ=255)
                                                                                                     ccGACACAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTCTTTTAATTACTCCTCAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCg                                           <  1:682924/143‑1 (MQ=255)
                                                                                                          aaaaCGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAAt                                      >  1:178398/1‑143 (MQ=255)
                                                                                                           aaaCGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATTATTGCCGGTTGGCGCCTAAACATTACTCAGACCAAACTCTTCGACAATc                                     >  1:393288/1‑143 (MQ=255)
                                                                                                             aCGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGGTTTTATAATTACTGCTTAAATGCACACGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGc                                   >  1:126100/1‑143 (MQ=255)
                                                                                                                 tGTCAACCATCCAAAACCGCTATGATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGACAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAat                               >  2:131080/1‑143 (MQ=255)
                                                                                                                     aaCCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGCCCGCGCATGTCTGGCGCGGTTGGCGACTGATGAAGTCGCAGACCCACCTCTTCGACAATCGCAACTAtgt                           >  1:242923/1‑143 (MQ=255)
                                                                                                                     aaCCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATAtgt                           <  2:656770/143‑1 (MQ=255)
                                                                                                                       ccATCCGAAACCGCTCTCATCCCTTCGCTGCGACCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGa                         <  1:823094/143‑1 (MQ=255)
                                                                                                                         aTCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCCGATGGCATCGGATTAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATc                       >  1:1054051/1‑143 (MQ=255)
                                                                                                                            cgacgcggCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGaa                    <  2:1020262/136‑1 (MQ=255)
                                                                                                                             gAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGACGTCGCTGACCAAACTCTTCGACAATCGCAAAGTTGTTATCGatt                   >  1:196695/1‑141 (MQ=255)
                                                                                                                             gAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAat                   >  1:428327/1‑143 (MQ=255)
                                                                                                                                aCCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACGGCTTAAATGCACCCGCCAGAGGGAGAATAACAAAACCGGTTGGCGACTGATTAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAatata                <  1:245947/143‑1 (MQ=255)
                                                                                                                                aCCGCTCTCATCCATTCGATGAGAGCGGTATGTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAatata                <  1:670997/143‑1 (MQ=255)
                                                                                                                                  cGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAATATATc              >  1:790579/1‑143 (MQ=255)
                                                                                                                                    ctctCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTACCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCATCGACAATAGCAAATATGTGATCGAATATATCac            >  2:993071/1‑142 (MQ=255)
                                                                                                                                    ctctCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGCCAATCGCAACTATGTGATCGAATATAGCAg            >  1:995610/1‑143 (MQ=255)
                                                                                                                                         aTCCATTCGATGAGAGCGGTTTTTTTAGTTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAATATATCAGCCTGa       <  2:367190/143‑1 (MQ=255)
                                                                                                                                         aTCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAATATATCAGCCTGa       >  1:159483/1‑143 (MQ=255)
                                                                                                                                           caaTTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAATATATCAGCCTGAAt     <  1:56639/141‑1 (MQ=255)
                                                                                                                                            cATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGGCAATCGCAAATATGTCATCGAATATATCAGCCTAAAtt    >  2:545387/1‑142 (MQ=255)
                                                                                                                                             aTTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCACAGAGCGAATATCATTTCCGGTTGGCGACTGAAGAAATAGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAATATATCAGCCTGAATGc   >  1:364146/1‑143 (MQ=255)
                                                                                                                                              ttCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAATATATCAGCCTGAATGCt  >  1:341053/1‑143 (MQ=255)
                                                                                                                                              |                                                                                                                                               
GGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAATATATCAGCCTGAATGCT  >  NC_000913/603251‑603535

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