Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 194,492 T→C 5.6% Y133Y (TAT→TAC 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,4920TC 5.6% 192.6 / 2.9 71Y133Y (TAT→TACnlpElipoprotein involved with copper homeostasis and adhesion
Reads supporting (aligned to +/- strand):  ref base T (28/39);  new base C (2/2);  total (30/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.45e-01

TCGGGCGAAAGGCGATGCGCTGGAGATGCTCGATCGTGAAGGCAATCCGATTGAATCGCAGTTCAACTATACGCTGGAAGCGGCACAATCCAGTTTACCTATGACGCCGATGACCCTGCGGGGCATGTATTTTTATATG  >  minE/194423‑194561
                                                                     |                                                                     
tCGGGCGAAAGGCGATGCGCTGGAGATGCTCGATCGTGAAGGCAATCCGATTGAATCGCAGTTCAACtat                                                                       <  1:4179999/70‑1 (MQ=255)
  gggCGAAAGGCGATGCGCTGGAGATGCTCGATCGTGAAGGCAATCCGATTGAATCGCAGTTCAACTATAc                                                                     <  1:738025/70‑1 (MQ=255)
      gAAAGGCGATGCGCTGGAGATGCTCGATCGTGAAGGCAATCCGATTGAATCGCAGTTCAACTATACGCTgg                                                                <  1:2908823/71‑1 (MQ=255)
         aGGCGATGCGCTGGAGATGCTCGATCGTGAAGGCAATCCGATTGAATCGCAGTTCAACTATACGCTGGAAg                                                             <  1:4565119/71‑1 (MQ=255)
           gCGATGCGCTGGAGATGCTCGATCGTGAAGGCAATCCGATTGAATCGCAGTTCAACTATACGCTGGAAg                                                             >  1:66758/1‑69 (MQ=255)
               tGCGCTGGAGATGCTCGATCGTGAAGGCAATCCGATTGAATCGCAGTTCAACTATACGCTGGAAGCGGca                                                        >  1:451964/1‑70 (MQ=255)
                gcgcTGGAGATGCTCGATCGTGAAGGCAATCCGATTGAATCGCAGTTCAACTATACGCTGGAAGCGGcac                                                       >  1:3954650/1‑70 (MQ=255)
                 cgcTGGAGATGCTCGATCGTGAAGGCAATCCGATTGAATCGCAGTTCAACTATACGCTGGAAGCGGCACaa                                                     <  1:2364849/71‑1 (MQ=255)
                 cgcTGGAGATGCTCGATCGTGAAGGCAATCCGATTGAATCGCAGTTCAACTATACGCTGGAAGCGGCACaa                                                     >  1:669168/1‑71 (MQ=255)
                  gcTGGAGATGCTCGATCGTGAAGGCAATCCGATTGAATCGCAGTTCAACTATACGCTGGAAGCGGCACAAt                                                    <  1:1008027/71‑1 (MQ=255)
                   cTGGAGATGCTCGATCGTGAAGGCAATCCGATTGAATCGCAGTTCAACTATACGCTGGaa                                                              <  1:478022/60‑1 (MQ=255)
                    tGGAGATGCTCGATCGTGAAGGCAATCCGATTGAATCGCAGTTCAACTATAc                                                                     <  1:3525800/52‑1 (MQ=255)
                     ggAGATGCTCGATCGTGAAGGCAATCCGATTGAATCGCAGTTCAACTATACGCTGGAAGCGGCACAATc                                                   >  1:1207774/1‑69 (MQ=255)
                      gagaTGCTCGATCGTGAAGGCAATCCGATTGAATCGCAGTTCAACTATACGCTGGAAGCGGCACAATCCAg                                                <  1:3656789/71‑1 (MQ=255)
                         aTGCTCGATCGTGAAGGCAATCCGATTGAATCGCAGTTCAACTATACGCTGGAAGCGGca                                                        <  1:176665/60‑1 (MQ=255)
                         aTGCTCGATCGTGAAGGCAATCCGATTGAATCGCAGTTCAACTATACGCTGGAAGCGGca                                                        <  1:106694/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)
                           gCTCGATCGTGAAGGCAATCCGATTGAATCGCAGTTCAACTATACGCTGGa                                                               >  1:592710/1‑51 (MQ=255)
                           gCTCGATCGTGAAGGCAATCCGATTGAATCGCAGTTCAACTATACGCTGGa                                                               >  1:3694785/1‑51 (MQ=255)
                           gCTCGATCGTGAAGGCAATCCGATTGAATCGCAGTTCAACTATACGCTGGAAGCGGCACAATCCAGTTTAc                                           >  1:1505101/1‑71 (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:4400753/1‑42 (MQ=255)
                                   gTGAAGGCAATCCGATTGAATCGCAGTTCAACTATACGCTgg                                                                >  1:2137827/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:615370/70‑1 (MQ=255)
                                      aaGGCAATCCGATTGAATCGCAGTTCAACTATACGCTGGAAGCGGCACAATCCAGTTTACCTATGACGcc                                 <  1:4034735/70‑1 (MQ=255)
                                         gCAATCCGATTGAATCGCAGTTCAACTATACGCTGGAAGCGGCACAATCCAGTTTACCTATGACGCCGATg                             <  1:2746143/71‑1 (MQ=255)
                                          cAATCCGATTGAATCGCAGTTCAACTATACGCTGGAAGCGGCACAATCCAGTTTACCTATGACGCCGATGa                            <  1:1973698/71‑1 (MQ=255)
                                           aaTCCGATTGAATCGCAGTTCAACTATACGCTGGAAGCGGCACAATCCAGTTTACCTATGACGCCGATGAc                           >  1:3638310/1‑71 (MQ=255)
                                            aTCCGATTGAATCGCAGTTCAACTATACGCTGGAAGCGGCACAATCCAGttt                                             <  1:3479266/52‑1 (MQ=255)
                                             tCCGATTGAATCGCAGTTCAACTATACGCTGGAAGCGGCACAATCGAGTTTAcc                                          <  1:1805559/54‑1 (MQ=255)
                                             tCCGATTGAATCGCAGTTCAACTATACGCTGGAAGCGGCACAATCCAGTTTAcc                                          <  1:1392310/54‑1 (MQ=255)
                                              ccGATTGAATCGCAGTTCAACTATACGCTGGAAGCGGCACAATCCAGTTTAc                                           >  1:2799435/1‑52 (MQ=255)
                                               cGATTGAATCGCAGTTCAACTATACGCTGGAAGCGGCACAATCCAGtt                                              <  1:355442/48‑1 (MQ=255)
                                                  ttGAATCGCAGTTCAACTATACGCTGGAAGCGGCACAATCCAGTTTACCTATGa                                     >  1:343470/1‑54 (MQ=255)
                                                   tGAATCGCAGTTCAACTATACGCTGGAAGCGGCACAATCCAGTTTAcc                                          >  1:2705404/1‑48 (MQ=255)
                                                   tGAATCGCAGTTCAACTATACGCTGGAAGCGGCACAATCCAGTTTACCt                                         >  1:443049/1‑49 (MQ=255)
                                                   tGAATCGCAGTTCAACTATACGCTGGAAGCGGCACAATCCAGTTTACCTATGACGCCGATGACCCTGCggg                   >  1:2375255/1‑71 (MQ=255)
                                                   tGAATCGCAGTTCAACTATACGCTGGAAGCGGCACAATCCAGTTTACCTATGACGCCGATGACCCTGCggg                   >  1:665469/1‑71 (MQ=255)
                                                    gAATCGCAGTTCAACTATACGCTGGAAGCGGCACAATCCAGTTTAc                                           >  1:149636/1‑46 (MQ=255)
                                                       tCGCAGTTCAACTATACGCTGGAAGCGGCACAATCCa                                                 >  1:566808/1‑37 (MQ=255)
                                                       tCGCAGTTCAACTATACGCTGGAAGCGGCACAATCCAGTTTACCTATGACGCCGATGACCCTGCGGGGCa                >  1:2442344/1‑70 (MQ=255)
                                                       tCGCAGTTCAACTATACGCTGGAAGCGGCACAATCCAGTTTACCTATGACGCCGATGACCCTGCGGGGCAt               >  1:1303052/1‑71 (MQ=255)
                                                           aGTTCAACTATACGCTGGAAGCGGCACAATCCAGTTTACCTATGACGCCGATGACCCTGCGGGGCATGTa            <  1:538674/70‑1 (MQ=255)
                                                           aGTTCAACTATACGCTGGAAGCGGCACAATCCAGTTTACCTATGACGCCGATGACCCTGCGGGGCATGTAt           <  1:508135/71‑1 (MQ=255)
                                                             ttCAACTATACGCTGGAAGCGGCACAATCCAGttt                                             <  1:1898144/35‑1 (MQ=255)
                                                             ttCAACTATACGCTGGAAGCGGCACAATCCAGTTTACCTATGACGCCGATGAc                           <  1:866489/53‑1 (MQ=255)
                                                               cAACTATACGCTGGAAGCGGCACAATCCAGTTTACCTATGACGCCGATGACCCTGCGGGGCATGTAtttt        <  1:429904/70‑1 (MQ=255)
                                                               cAACTATACGCTGGAAGCGGCACAATCCAGTTTACCTATGACGCCGATGACCCTGCGGGGCATGTAtttt        >  1:680114/1‑70 (MQ=255)
                                                                aaCTATACGCTGGAAGCGGCACAATCCAGTTTACCTATGACGCCGATGAccc                         >  1:3513228/1‑52 (MQ=255)
                                                                aaCTATACGCTGGAAGCGGCACAATCCAGTTTACCTATGACGCCGATGAccc                         >  1:1442219/1‑52 (MQ=255)
                                                                aaCTATACGCTGGAAGCGGCACAATCCAGTTTACCTATGACGCCGATGACCCTGCGGGGCATGTAttttt       <  1:93189/70‑1 (MQ=255)
                                                                aaCTATACGCTGGAAGCGGCACAATCCAGTTTACCTATGACGCCGATGACCCTGCGGGGCATGTAttttt       <  1:183093/70‑1 (MQ=255)
                                                                aaCTATACGCTGGAAGCGGCACAATCCAGTTTACCTATGACGCCGATGACCCTGCGGGGCATGTATTTTta      <  1:1875210/71‑1 (MQ=255)
                                                                aaCTATACGCTGGAAGCGGCACAATCCAGTTTACCTATGACGCCGATGACCCTGCGGGGCATGTATTTTta      <  1:1999672/71‑1 (MQ=255)
                                                                aaCTATACGCTGGAAGCGGCACAATCCAGTTTACCTATGACGCCGATGACCCTGCGGGGCATGTATTTTta      <  1:876974/71‑1 (MQ=255)
                                                                  cTATACGCTGGAAGCGGCACAATCCAGTTTACCTATGACGCCGAt                              <  1:2498821/45‑1 (MQ=255)
                                                                    ataCGCTGGAAGCGGCACAATCCAGTTTACCTATGa                                     >  1:2128857/1‑36 (MQ=255)
                                                                    ataCGCTGGAAGCGGCACAATCCAGTTTACCTATGACGCCGATGACCCTGCg                     >  1:278598/1‑52 (MQ=255)
                                                                    ataCGCTGGAAGCGGCACAATCCAGTTTACCTATGACGCCGATGACCCTGCGGGGCATGTATTTTTATATg  >  1:3003810/1‑71 (MQ=255)
                                                                    ataCGCTGGAAGCGGCACAATCCAGTTTACCTATGACGCCGATGACCCTGCGGGGCATGTATTTTTATATg  >  1:77678/1‑71 (MQ=255)
                                                                     taCGCTGGAAGCGGCACAATCCAGTTTACCTATGCCGCCGATGACCCTGCGGGGCATGTATTTTTATATg  >  1:3021666/1‑70 (MQ=255)
                                                                     taCGCTGGAAGCGGCACAATCCAGTTTACCTATGACGCCGATGACCCTGc                      <  1:2954030/50‑1 (MQ=255)
                                                                     |                                                                     
TCGGGCGAAAGGCGATGCGCTGGAGATGCTCGATCGTGAAGGCAATCCGATTGAATCGCAGTTCAACTATACGCTGGAAGCGGCACAATCCAGTTTACCTATGACGCCGATGACCCTGCGGGGCATGTATTTTTATATG  >  minE/194423‑194561

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: