Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 194,465 C→T 5.2% G124G (GGC→GGT nlpE → lipoprotein involved with copper homeostasis and adhesion

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE194,4650CT 5.2% 302.0 / 4.2 96G124G (GGC→GGTnlpElipoprotein involved with copper homeostasis and adhesion
Reads supporting (aligned to +/- strand):  ref base C (40/51);  new base T (2/3);  total (42/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GACAGCAAAGGTGAAAAGTCATATTATCGGGCGAAAGGCGATGCGCTGGAGATGCTCGATCGTGAAGGCAATCCGATTGAATCGCAGTTCAACTATACGCTGGAAGCGGCACAATCCAGTTTACCTATGACGCCGATGA  >  minE/194397‑194535
                                                                    |                                                                      
gACAGCAAAGGTGAAAAGTCATATTATCGGGCGAAAGGCGATGCGCTGGAGATGCTCGATCGTGAAGGCaa                                                                      >  1:4456785/1‑71 (MQ=255)
   aGCAAAGGTGAAAAGTCATATTATCGGGCGAAAGGCGATGCGCTGGAGATGCTCGATCGTGAAGGCAATcc                                                                   >  1:2782839/1‑71 (MQ=255)
        aGGTGAAAAGTCATATTATCGGGCGAAAGGCGATGCGCTGGAGATGCTCGATCGTGAAGGCAATCCGAtt                                                               >  1:462101/1‑70 (MQ=255)
        aGGTGAAAAGTCATATTATCGGGCGAAAGGCGATGCGCTGGAGATGCTCGATCGTGAAGGCAATCCGATTg                                                              <  1:1244676/71‑1 (MQ=255)
        aGGTGAAAAGTCATATTATCGGGCGAAAGGCGATGCGCTGGAGATGCTCGATCGTGAAGGCAATCCGATTg                                                              <  1:3297494/71‑1 (MQ=255)
        aGGTGAAAAGTCATATTATCGGGCGAAAGGCGATGCGCTGGAGATGCTCGATCGTGAAGGCAATCCGATTg                                                              <  1:4545705/71‑1 (MQ=255)
        aGGTGAAAAGTCATATTATCGGGCGAAAGGCGATGCGCTGGAGATGCTCGATCGTGAAGGCAATCCGATTg                                                              <  1:4288525/71‑1 (MQ=255)
            gAAAAGTCATATTATCGGGCGAAAGGCGATGCGCTGGAGATGCTCGATCGTGAAGGCAATCCGa                                                                 >  1:3897572/1‑64 (MQ=255)
            gAAAAGTCATATTATCGGGCGAAAGGCGATGCGCTGGAGATGCTCGATCGTGAAGGCAATCCGATTGAATc                                                          >  1:3793905/1‑71 (MQ=255)
             aaaaGTCATATTATCGGGCGAAAGGCGATGCGCTGGAGATGCTCGATCGTGAAGGCAATCCGATTGAATCg                                                         >  1:3499623/1‑71 (MQ=255)
             aaaaGTCATATTATCGGGCGAAAGGCGATGCGCTGGAGATGCTCGATCGTGAAGGCAATCCGATTGAATCg                                                         <  1:4175524/71‑1 (MQ=255)
              aaaGTCATATTATCGGGCGAAAGGCGATGCGCTGGAGATGCTCGATCGTGAAGGCAATCCGATTGAATCg                                                         >  1:932418/1‑70 (MQ=255)
              aaaGTCATATTATCGGGCGAAAGGCGATGCGCTGGAGATGCTCGATCGTGAAGGCAATCCGATTGAATCg                                                         >  1:3678459/1‑70 (MQ=255)
              aaaGTCATATTATCGGGCGAAAGGCGATGCGCTGGAGATGCTCGATCGTGAAGGCAATCCGATTGAATCg                                                         <  1:3576062/70‑1 (MQ=255)
                aGTCATATTATCGGGCGAAAGGCGATGCGCTGGAGATGCTCCATCGTGAAGGCAATCCGATTGAATCGCAg                                                      >  1:4284955/1‑71 (MQ=255)
                  tCATATTATCGGGCGAAAGGCGATGCGCTGGAGATGCTCGATCGTGAAGGTAATCCg                                                                  <  1:604811/57‑1 (MQ=255)
                  tCATATTATCGGGCGAAAGGCGATGCGCTGGAGATGCTCGATCGTGAAGGCAATCCGATTGa                                                             <  1:1777231/62‑1 (MQ=255)
                  tCATATTATCGGGCGAAAGGCGATGCGCTGGAGATGCTCGATCGTGAAGGCAATCCGATTGAATCGCAGtt                                                    >  1:3228437/1‑71 (MQ=255)
                  tCATATTATCGGGCGAAAGGCGATGCGCTGGAGATGCTCGATCGTGAAGGCAATCCGATTGAATCGCAGtt                                                    <  1:4648619/71‑1 (MQ=255)
                   aggattatCGGGCGAAAGGCGATGCGCTGGAGATGCTCGATCGTGAAGGCAATCCGATTGAATCGCAGTTc                                                   <  1:2244822/68‑1 (MQ=255)
                    atatTATCGGGCGAAATGCGATGCGCTGGAGATGCTCGATCGTGAAGGCAATCCGATTGAATCGCAGTTc                                                   <  1:4695594/70‑1 (MQ=255)
                    atatTATCGGGCGAAAGGCGATGCGCTGGAGATGCTCGATCGTGAAGTCAATCCGATTGAATCGCAGTTCa                                                  <  1:3642481/71‑1 (MQ=255)
                    atatTATCGGGCGAAAGGCGATGCGCTGGAGATGCTCGATCGTGAAGGCAATCCGATTGAAt                                                           >  1:2527190/1‑62 (MQ=255)
                    atatTATCGGGCGAAAGGCGATGCGCTGGAGATGCTCGATCGTGAAGGCAATCCGATTGAATCGCAGTTc                                                   <  1:1011317/70‑1 (MQ=255)
                    atatTATCGGGCGAAAGGCGATGCGCTGGAGATGCTCGATCGTGAAGGCAATCCGATTGAATCGCAGTTCa                                                  <  1:3581380/71‑1 (MQ=255)
                    atatTATCGGGCGAAAGGCGATGCGCTGGAGATGCTCGATCGTGAAGGCAATCCGATTGAATCGCAGTTCa                                                  <  1:2975481/71‑1 (MQ=255)
                     tattatCGGGCGAAAGGCGATGCGCTGGAGATGCTCGATCGTGAAGGCAATCCGATTGAATCGCAGTTCa                                                  <  1:2719977/70‑1 (MQ=255)
                      attatCGGGCGAAAGGCGATGCGCTGGAGATGCTCGATCGTGAAGGCAATCCGa                                                                 >  1:221517/1‑54 (MQ=255)
                       ttatCGGGCGAAAGGCGATGCGCTGGAGATGCTCGATCGTGAAGGCAATCCGATTGAATCGCAGTTCAAc                                                >  1:3368762/1‑70 (MQ=255)
                       ttatCGGGCGAAAGGCGATGCGCTGGAGATGCTCGATCGTGAAGGCAATCCGATTGAATCGCAGTTCAAc                                                >  1:1577817/1‑70 (MQ=255)
                         atCGGGCGAAAGGCGATGCGCTGGAGATGCTCGATCGTGAAGGCAATCCGAt                                                                >  1:2404855/1‑52 (MQ=255)
                         atCGGGCGAAAGGCGATGCGCTGGAGATGCTCGATCGTGAAGGCAATCCGATTGAATCGCAGTTCAACta                                              >  1:4584889/1‑70 (MQ=255)
                         atCGGGCGAAAGGCGATGCGCTGGAGATGCTCGATCGTGAAGGCAATCCGATTGAATCGCAGTTCAACta                                              >  1:1469736/1‑70 (MQ=255)
                         atCGGGCGAAAGGCGATGCGCTGGAGATGCTCGATCGTGAAGCCAATCCGATTGAATCGCAGTTCAACta                                              <  1:704614/70‑1 (MQ=255)
                          tCGGGCGAAAGGCGATGCGCTGGAGATGCTCGATCGTGAAGGCAATCCGATTGAATCGCAGTTCa                                                  <  1:4573381/65‑1 (MQ=255)
                          tCGGGCGAAAGGCGATGCGCTGGAGATGCTCGATCGTGAAGGCAATCCGATTGAATCGCAGTTCa                                                  <  1:1278495/65‑1 (MQ=255)
                          tCGGGCGAAAGGCGATGCGCTGGAGATGCTCGATCGTGAAGGCAATCCGATTGAATCGCAGTTCAACtat                                             <  1:4179999/70‑1 (MQ=255)
                            gggCGAAAGGCGATGCGCTGGAGATGCTCGATCGTGAAGGCAATCCGATTGAATCGCAGTTCAACTATAc                                           <  1:738025/70‑1 (MQ=255)
                                gAAAGGCGATGCGCTGGAGATGCTCGATCGTGAAGGCAATCCGATTGAATCGCAGTTCAACTATACGCTgg                                      <  1:2908823/71‑1 (MQ=255)
                                   aGGCGATGCGCTGGAGATGCTCGATCGTGAAGGCAATCCGATTGAATCg                                                         >  1:4021914/1‑49 (MQ=255)
                                   aGGCGATGCGCTGGAGATGCTCGATCGTGAAGGCAATCCGATTGAATCGCAGTTCAACTATACGCTGGAAg                                   <  1:4565119/71‑1 (MQ=255)
                                     gCGATGCGCTGGAGATGCTCGATCGTGAAGGCAATCCGATTGAATCGCa                                                       >  1:1654710/1‑49 (MQ=255)
                                     gCGATGCGCTGGAGATGCTCGATCGTGAAGGCAATCCGATTGAATCGCa                                                       >  1:649391/1‑49 (MQ=255)
                                     gCGATGCGCTGGAGATGCTCGATCGTGAAGGCAATCCGATTGAATCGCAGTTCAACTATACGCTGGAAg                                   >  1:66758/1‑69 (MQ=255)
                                         tGCGCTGGAGATGCTCGATCGTGAAGGCAATCCGATTGAATCGCAGTTc                                                   <  1:398715/49‑1 (MQ=255)
                                         tGCGCTGGAGATGCTCGATCGTGAAGGCAATCCGATTGAATCGCAGTTCAACTATACGCTGGAAGCGGca                              >  1:451964/1‑70 (MQ=255)
                                          gcgcTGGAGATGCTCGATCGTGAAGGCAATCCGATTGAATCGCa                                                       >  1:1111353/1‑44 (MQ=255)
                                          gcgcTGGAGATGCTCGATCGTGAAGGCAATCCGATTGAATCGCAGTTCAACTATACGCTGGAAGCGGcac                             >  1:3954650/1‑70 (MQ=255)
                                           cgcTGGAGATGCTCGATCGTGAAGGCAATCCGATTGAATCGCAGTTCAACTATACGCTGGAAGCGGCACaa                           >  1:669168/1‑71 (MQ=255)
                                           cgcTGGAGATGCTCGATCGTGAAGGCAATCCGATTGAATCGCAGTTCAACTATACGCTGGAAGCGGCACaa                           <  1:2364849/71‑1 (MQ=255)
                                            gcTGGAGATGCTCGATCGTGAAGGCAATCCGATTGAATCGCAGTTCAACTATACGCTGGAAGCGGCACAAt                          <  1:1008027/71‑1 (MQ=255)
                                             cTGGAGATGCTCGATCGTGAAGGCAATCCGATTGAATCg                                                         <  1:3426459/39‑1 (MQ=255)
                                             cTGGAGATGCTCGATCGTGAAGGCAATCCGATTGAATCGCAGTTCAACTATACGCTGGaa                                    <  1:478022/60‑1 (MQ=255)
                                              tGGAGATGCTCGATCGTGAAGGCAATCCGATTGAATCGCAGTTCAAc                                                <  1:3953316/47‑1 (MQ=255)
                                              tGGAGATGCTCGATCGTGAAGGCAATCCGATTGAATCGCAGTTCAACTATAc                                           <  1:3525800/52‑1 (MQ=255)
                                               ggAGATGCTCGATCGTGAAGGCAATCCGATTGAATCGCAGTTCAACTATACGCTGGAAGCGGCACAATc                         >  1:1207774/1‑69 (MQ=255)
                                                gagaTGCTCGATCGTGAAGGCAATCCGATTGAATCGCAGTTCAACTATACGCTGGAAGCGGCACAATCCAg                      <  1:3656789/71‑1 (MQ=255)
                                                  gaTGCTCGATCGTGAAGGCAATCCGATTGAATCGCAg                                                      >  1:2308885/1‑37 (MQ=255)
                                                  gaTGCTCGATCGTGAAGGCAATCCGATTGAATCGCAg                                                      >  1:39843/1‑37 (MQ=255)
                                                   aTGCTCGATCGTGAAGGCAATCCGATTGAATCGCAGTTCa                                                  <  1:99361/40‑1 (MQ=255)
                                                   aTGCTCGATCGTGAAGGCAATCCGATTGAATCGCAGTTCa                                                  <  1:2322318/40‑1 (MQ=255)
                                                   aTGCTCGATCGTGAAGGCAATCCGATTGAATCGCAGTTCa                                                  <  1:2843740/40‑1 (MQ=255)
                                                   aTGCTCGATCGTGAAGGCAATCCGATTGAATCGCAGTTCAACTATACGCTGGAAGCGGca                              <  1:106694/60‑1 (MQ=255)
                                                   aTGCTCGATCGTGAAGGCAATCCGATTGAATCGCAGTTCAACTATACGCTGGAAGCGGca                              <  1:176665/60‑1 (MQ=255)
                                                   aTGCTCGATCGTGAAGGCAATCCGATTGAATCGCAGTTCAACTATACGCTGGAAGCGGCACAATCCAGttt                   <  1:3965547/71‑1 (MQ=255)
                                                   aTGCTCGATCGTGAAGGCAATCCGATTGAATCGCAGTTCAACTATACGCTGGAAGCGGCACAATCCAGttt                   <  1:2139750/71‑1 (MQ=255)
                                                   aTGCTCGATCGTGAAGGCAATCCGATTGAATCGCAGTTCAACTATACGCTGGAAGCGGCACAATCCAGtt                    >  1:1655820/1‑70 (MQ=255)
                                                    tGCTCGATCGTGAAGGCAATCCGATTGAATCGCAGTTc                                                   <  1:1063174/38‑1 (MQ=255)
                                                     gCTCGATCGTGAAGGCAATCCGATTGAATCGCAGt                                                     >  1:2993355/1‑35 (MQ=255)
                                                     gCTCGATCGTGAAGGCAATCCGATTGAATCGCAGTTCa                                                  <  1:3553370/38‑1 (MQ=255)
                                                     gCTCGATCGTGAAGGCAATCCGATTGAATCGCAGTTCAACTATACGCTGGa                                     >  1:3694785/1‑51 (MQ=255)
                                                     gCTCGATCGTGAAGGCAATCCGATTGAATCGCAGTTCAACTATACGCTGGa                                     >  1:592710/1‑51 (MQ=255)
                                                     gCTCGATCGTGAAGGCAATCCGATTGAATCGCAGTTCAACTATACGCTGGAAGCGGCACAATCCAGTTTAc                 >  1:1505101/1‑71 (MQ=255)
                                                      cTCGATCGTGAAGGCAATCCGATTGAATCGCAGTTCAACt                                               <  1:821580/40‑1 (MQ=255)
                                                      cTCGATCGTGAAGGCAATCCGATTGAATCGCAGTTCAACTATACGCTGGAAGCGGCACAATCCAGTTTAcc                <  1:3808556/71‑1 (MQ=255)
                                                      cTCGATCGTGAAGGCAATCCGATTGAATCGCAGTTCAACTATACGCTGGAAGCGGCACAATCCAGTTTAcc                >  1:3386455/1‑71 (MQ=255)
                                                        cGATCGTGAAGGCAATCCGATTGAATCGCAGTTCAACTATACGCTGGAAGCGGCACAATCCAg                      <  1:16879/63‑1 (MQ=255)
                                                         gATCGTGAAGGTAATCCGATTGAATCGCAGTTCAACTACACGCTGGaacc                                  >  1:2304392/1‑48 (MQ=38)
                                                         gATCGTGAAGGCAATCCGATTGAATCGCAGTTCAACTATACGCTGGAAGCGGCACAATCCAGTTTACCt               <  1:2899788/69‑1 (MQ=255)
                                                         gATCGTCAAGGTAATCCGATTGAATCGCAGTTCAACTACACGCTGGaacc                                  >  1:7713/1‑48 (MQ=255)
                                                          ttCGTGAAGGCAATCCGATTGAATCGCAGTTCAACTATACGCTGGAAGCGGCACAATCCAGTTTACCTa              <  1:1510806/68‑1 (MQ=255)
                                                          aTCGTGAAGGCAATCCGATTGAATCGCAGTTCAACTATACGCTGGAAGCGGCACAATCCAGTTTACCTAt             <  1:3578378/70‑1 (MQ=255)
                                                           tCGTGAAGGCAATCCGATTGAATCGCAGTTCAACtata                                            >  1:4332511/1‑38 (MQ=255)
                                                           tCGTGAAGGCAATCCGATTGAATCGCAGTTCAACtata                                            >  1:1148640/1‑38 (MQ=255)
                                                           tCGTGAAGGCAATCCGATTGAATCGCAGTTCAACTATACGCTGGAAGCGGCACAATCCAGTTTACCTATGa           >  1:3613669/1‑71 (MQ=255)
                                                            cGTGAAGGCACTCCGATTGAATCGCAGTTCAACTATACGCTGGAAGCGGCACaa                           >  1:2655993/1‑54 (MQ=255)
                                                             gTGAAGGCAATCCGATTGAATCGCAGTTCAACTATACGCTgg                                      >  1:2137827/1‑42 (MQ=255)
                                                             gTGAAGGCAATCCGATTGAATCGCAGTTCAACTATACGCTgg                                      >  1:4400753/1‑42 (MQ=255)
                                                             gTGAAGGCAATCCGATTGAATCGCAGTTCAACTATACGCTGGAAGCGGcac                             <  1:2963225/51‑1 (MQ=255)
                                                              tGAAGGCAATCCGATTGAATCGCAGTTCAACTATACGCTGGAAGc                                  <  1:1632245/45‑1 (MQ=255)
                                                              tGAAGGCAATCCGATTGAATCGCAGTTCAACTATACGCTGGAAGCGGCACAATCCAGttt                   <  1:3956167/60‑1 (MQ=255)
                                                               gAAGGTAATCCGATTGAATCGCAGTTCAACTACACGCt                                        <  1:2622143/38‑1 (MQ=38)
                                                                aaGGTAATCCGATTGAATCGCAGTTCAACTACACGCt                                        <  1:4123140/37‑1 (MQ=255)
                                                                aaGGCAATCCGATTGAATCGCAGTTCAACTATACGCTGGAAGCGGCACAATCCAGTTTACCTATGACGcc       <  1:4034735/70‑1 (MQ=255)
                                                                aaGGCAATCCGATTGAATCGCAGTTCAACTATACGCTGGAAGCGGCACAATCCAGTTTACCTATGACGcc       <  1:615370/70‑1 (MQ=255)
                                                                   gCAATCCGATTGAATCGCAGTTCAACTATACGCTGGAAGCGGCACAATCCAGTTTACCTATGACGCCGATg   <  1:2746143/71‑1 (MQ=255)
                                                                    cAATCCGATTGAATCGCAGTTCAACTATACGCTGGAAGCGGCACAATCCAGTTTACCTATGACGCCGATGa  <  1:1973698/71‑1 (MQ=255)
                                                                    |                                                                      
GACAGCAAAGGTGAAAAGTCATATTATCGGGCGAAAGGCGATGCGCTGGAGATGCTCGATCGTGAAGGCAATCCGATTGAATCGCAGTTCAACTATACGCTGGAAGCGGCACAATCCAGTTTACCTATGACGCCGATGA  >  minE/194397‑194535

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: