Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 194,942 G→T 20.7% intergenic (+138/+62) nlpE → / ← yaeF lipoprotein involved with copper homeostasis and adhesion/predicted lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE194,9420GT20.7% 87.9 / 9.4 35intergenic (+138/+62)nlpE/yaeFlipoprotein involved with copper homeostasis and adhesion/predicted lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (12/15);  new base T (3/4);  total (16/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.83e-01

TTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGG  >  minE/194875‑195010
                                                                    |                                                                    
tttttGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGgcg                                                                   >  1:823979/1‑71 (MQ=255)
tttttGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGgc                                                                    >  1:4768027/1‑70 (MQ=255)
  tccgTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACAGTGCCAACTGCCTGATGCGACGCTCgcgc                                                                  >  1:3064589/4‑70 (MQ=11)
   ccgTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACAGTGCCAACTGCCTGATGCGACGCTTGCGCGt                                                                >  1:4934577/3‑71 (MQ=11)
    cgTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAACAGTGCCAACTGCCTGATGCGACGCTTGCGCGTc                                                               >  1:3941813/2‑71 (MQ=11)
      tAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCtt                                                             <  1:1143037/71‑1 (MQ=255)
      tAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCtt                                                             <  1:5137867/71‑1 (MQ=255)
        ggCCGGATAAGGCGTTTACGCCGCATCCGGCAACAGTGCCAACTGCCTGATGCGACGCTTGCGCGTCTTAt                                                           <  1:1916953/71‑1 (MQ=11)
         gCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTAt                                                           <  1:5507982/70‑1 (MQ=255)
         gCCGGATAAGGCGTTTACGCCGCATCCGGCAACAGTGCCAACTGCCTGATGCGACGCTTGCGCGTCTTATc                                                          <  1:4794087/71‑1 (MQ=11)
          ccGTATAAGGCGTTTACGCCGCATCCGGCAACAGTGCCAACTGCCTGATGCGACGCTTGCGCGTCTTAt                                                           <  1:1212443/69‑1 (MQ=11)
             gATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGc                                                      >  1:883738/1‑71 (MQ=255)
                   gCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTg                                                                      <  1:957765/49‑1 (MQ=255)
                   gCGTTTACGCCGCATCCGGCAACAGTGCCAACTGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACa                                                 >  1:1408254/1‑70 (MQ=1)
                     gTTTACGCCGCATCCGGCAATGGTGCCGACTGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACaaat                                              <  1:5002092/71‑3 (MQ=11)
                       ttACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCa                                                         >  1:5319103/1‑58 (MQ=255)
                                aTCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATAcc                                   >  1:4895975/1‑71 (MQ=255)
                                 tCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCt                                  <  1:5325089/71‑1 (MQ=255)
                                    ggCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTaa                               <  1:3638862/71‑1 (MQ=255)
                                    ggCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTaa                               <  1:3645414/71‑1 (MQ=255)
                                        aTGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCAt                                      <  1:4844560/60‑1 (MQ=255)
                                        aTGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCAt                                      <  1:103898/60‑1 (MQ=255)
                                        aTGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCAt                                      <  1:4839344/60‑1 (MQ=255)
                                        aTGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCg                           <  1:4997723/71‑1 (MQ=255)
                                           gTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGa                         >  1:1363025/1‑70 (MQ=255)
                                           gTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGa                         <  1:2437194/70‑1 (MQ=255)
                                                 aaCGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTa                                                   >  1:5123926/1‑38 (MQ=40)
                                                  aCGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCgg                          <  1:451821/62‑1 (MQ=255)
                                                       cTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCtt                                 <  1:631807/51‑1 (MQ=255)
                                                         gATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACtt           >  1:2156525/1‑71 (MQ=255)
                                                           tGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTa         <  1:538903/71‑1 (MQ=255)
                                                            gCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACttt          <  1:2208516/69‑1 (MQ=255)
                                                               aCGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCAc     >  1:3509702/1‑71 (MQ=255)
                                                                cGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCAcc    >  1:4606618/1‑71 (MQ=255)
                                                                 gCTGGCGCGTCTTATCCTGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATc       >  1:3536427/1‑67 (MQ=255)
                                                                  cTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCgg  >  1:4721254/1‑71 (MQ=255)
                                                                    |                                                                    
TTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGG  >  minE/194875‑195010

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: