Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 675,311 T→A 36.5% V43E (GTG→GAG)  ycbG → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE675,3110TA36.5% 30.0 / 30.9 41V43E (GTG→GAG) ycbGconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (21/5);  new base A (13/2);  total (34/7)
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.98e-01

AATCACCCGTTACATAGAAGCCAGTGCCGCCCAGGAAGCCGTCGATGTGTTGCTCTCGCTGGAAAATGAACCCGTGCTGGTAAATGGCTGGATTGACAAGCATATGAATCCGGAAC  >  minE/675264‑675379
                                               |                                                                    
aaTCACCCGTTACATAGAAGCCAGTGCCGCCCAGGAAGCCGTCGATGTGTTGCTCTCGCTGGAAAATGa                                                 >  1:2954299/1‑69 (MQ=255)
aaTCACCCGTTACATAGAAGCCAGTGCCGCCCAGGAAGCCGTCGATGTGTTGCTCTCGCTGGAAAATGa                                                 >  1:1117631/1‑69 (MQ=255)
   cACCCGTTACATAGAAGCCAGTGCCGCCCAGGAAGCCGTCGATGTGTTGCTCTCGCTGGAAAATGAAcc                                              <  1:1959423/69‑1 (MQ=255)
   cACCCGTTACATAGAAGCCAGTGCCGCCCAGGAAGCCGTCGATGAGTTGCTCTCGCTGGAAAATGAAc                                               <  1:2752105/68‑1 (MQ=255)
   cACCCGTTACATAGAAGCCAGTGCCGCCCAGGAAGCCGTCGATGAGTTGCTCTCGCTGGAAAATGAAc                                               <  1:725039/68‑1 (MQ=255)
           aCATAGAAGCCAGTGCCGCCCAGGAAGCCGTCGAtgt                                                                      >  1:1537057/1‑37 (MQ=255)
           aCATAGAAGCCAGTGCCGCCCAGGAAGCCGTCGAtgt                                                                      >  1:431804/1‑37 (MQ=255)
           aCATAGAAGCCAGTGCCGCCCAGGAAGCCGTCGAtgt                                                                      >  1:1937965/1‑37 (MQ=255)
           aCATAGAAGCCAGTGCCGCCCAGGAAGCCGTCGAtgt                                                                      >  1:563682/1‑37 (MQ=255)
           aCATAGAAGCCAGTGCCGCCCAGGAAGCCGTCGAtgt                                                                      >  1:605219/1‑37 (MQ=255)
           aCATAGAAGCCAGTGCCGCCCAGGAAGCCGTCGAtgt                                                                      >  1:125132/1‑37 (MQ=255)
           aCATAGAAGCCAGTGCCGCCCAGGAAGCCGTCGAtgt                                                                      >  1:1243892/1‑37 (MQ=255)
           aCATAGAAGCCAGTGCCGCCCAGGAAGCCGTCGAtgt                                                                      >  1:1122241/1‑37 (MQ=255)
           aCATAGAAGCCAGTGCCGCCCAGGAAGCCGTCGATGTGTTGCTCTCGCTGGAAAATGAACCCGTGCTg                                       >  1:643072/1‑68 (MQ=255)
           aCATAGAAGCCAGTGCCGCCCAGGAAGCCGTCGATGTGTTGCTCTCGCTGGAAAATGAACCCGTGCTg                                       >  1:2303953/1‑68 (MQ=255)
           aCATAGAAGCCAGTGCCGCCCAGGAAGCCGTCGATGTGTTGCTCTCGCTGGAAAATGAACCCGTGCTg                                       >  1:1821889/1‑68 (MQ=255)
           aCATAGAAGCCAGTGCCGCCCAGGAAGCCGTCGATGTGTTGCTCTCGCTGGAAAATGAACCCGTGCTg                                       >  1:2307743/1‑68 (MQ=255)
           aCATAGAAGCCAGTGCCGCCCAGGAAGCCGTCGATGTGTTGCTCTCGCTGGAAAATGAACCCGTGCTg                                       >  1:449146/1‑68 (MQ=255)
                gAAGCCAGTGCCGCCCAGGAAGCCGTCGATGTGTTGCTCTCGCTGGAAAATGAACCCGTGCTGGTaaa                                  >  1:2777632/1‑68 (MQ=255)
                gAAGCCAGTGCCGCCCAGGAAGCCGTCGATGTGTTGCTCTCGCTGGAAAATGAACCCGTGCTGGTa                                    <  1:266241/66‑1 (MQ=255)
                gAAGCCAGTGCCGCCCAGGAAGCCGTCGATGTGTTGCTCTCGCTGGAAAATGAACCCGTGCTGGTa                                    >  1:2258797/1‑66 (MQ=255)
                gAAGCCAGTGCCGCCCAGGAAGCCGTCGATGTGTTGCTCTCGCTGGAAAATGAACCCGTGCTGGTa                                    <  1:258759/66‑1 (MQ=255)
                gAAGCCAGTGCCGCCCAGGAAGCCGTCGATGTGTTGCTCTCGCTGGAAAATGAACCCGTGCTGGTa                                    >  1:2373894/1‑66 (MQ=255)
                gAAGCCAGTGCCGCCCAGGAAGCCGTCGATGTGTTGCTCTCGCTGGAAAATGAACCCGTGCTGGTAAAt                                 >  1:2037386/1‑69 (MQ=255)
                gAAGCCAGTGCCGCCCAGGAAGCCGTCGATGTGTTGCTCTCGCTGGAAAATGAACCCGTGCTGGTAAAt                                 >  1:1723618/1‑69 (MQ=255)
                gAAGCCAGTGCCGCCCAGGAAGCCGTCGATGTGTTGCTCTCGCTGGAAAATGAACCCGTGCTGGTAAAt                                 >  1:3162124/1‑69 (MQ=255)
                      aGTGCCGCCCAGGAAGCCGTCGATGTGTTGCTCTCGCTGGAAAATGAACCCGTGCTGGTa                                    <  1:3098250/60‑1 (MQ=255)
                      aGTGCCGCCCAGGAAGCCGTCGATGTGTTGCTCTCGCTGGAAAATGAACCCGTGCTGGTAAATGGCTgg                           >  1:1517021/1‑69 (MQ=255)
                      aGTGCCGCCCAGGAAGCCGTCGATGTGTTGCTCTCGCTGGAAAATGAACCCGTGCTGGTAAATGGCTgg                           >  1:2651106/1‑69 (MQ=255)
                        tGCCGCCCAGGAAGCCGTCGATGAGTTGCTCTCGCTGGaaa                                                     >  1:2473856/1‑41 (MQ=255)
                        tGCCGCCCAGGAAGCCGTCGATGAGTTGCTCTCGCTGGAAAATGAACCCGTGCTGGTAAATGGTTGGa                          >  1:474073/1‑68 (MQ=255)
                        tGCCGCCCAGGAAGCCGTCGATGAGTTGCTCTCGCTGGAAAATGAACCCGTGCTGGTAAATGGTTGGa                          >  1:730203/1‑68 (MQ=255)
                        tGCCGCCCAGGAAGCCGTCGATGAGTTGCTCTCGCTGGAAAATGAACCCGTGCTGGTAAATGGTTGGa                          >  1:3211009/1‑68 (MQ=255)
                        tGCCGCCCAGGAAGCCGTCGATGAGTTGCTCTCGCTGGAAAATGAACCCGTGCTGGTAAATGGTTGGa                          >  1:2822904/1‑68 (MQ=255)
                        tGCCGCCCAGGAAGCCGTCGATGAGTTGCTCTCGCTGGAAAATGAACCCGTGCTGGTAAATGGTTGGa                          >  1:844925/1‑68 (MQ=255)
                        tGCCGCCCAGGAAGCCGTCGATGAGTTGCTCTCGCTGGAAAATGAACCCGTGCTGGTAAATGGTTGGa                          >  1:1078704/1‑68 (MQ=255)
                        tGCCGCCCAGGAAGCCGTCGATGAGTTGCTCTCGCTGGAAAATGAACCCGTGCTGGTAAATGGTTGGa                          >  1:103637/1‑68 (MQ=255)
                        tGCCGCCCAGGAAGCCGTCGATGAGTTGCTCTCGCTGGAAAATGAACCCGTGCTGGTAAATGGTTGGa                          >  1:1597083/1‑68 (MQ=255)
                        tGCCGCCCAGGAAGCCGTCGATGAGTTGCTCTCGCTGGAAAATGAACCCGTGCTGGTAAATGGTTGGa                          >  1:1681238/1‑68 (MQ=255)
                        tGCCGCCCAGGAAGCCGTCGATGAGTTGCTCTCGCTGGAAAATGAACCCGTGCTGGTAAATGGTTGGa                          >  1:1918862/1‑68 (MQ=255)
                        tGCCGCCCAGGAAGCCGTCGATGAGTTGCTCTCGCTGGAAAATGAACCCGTGCTGGTAAATGGTTGGa                          >  1:1999902/1‑68 (MQ=255)
                        tGCCGCCCAGGAAGCCGTCGATGAGATGCTCTCGCTCGAAAATGAACCCGTGCTgga                                     >  1:3222274/1‑56 (MQ=255)
                         gccgccCAGGAAGCCGTCGATGTGTTGCTCTCGCTGGAAAATGAACCCGTGCTGGTAAATGGCTGGAtt                        >  1:3055641/1‑69 (MQ=255)
                         gccgccCAGGAAGCCGTCGATGTGTTGCTCTCGCTGGAAAATGAACCCGTGCTGGTAAATGGCTGGAtt                        >  1:2456648/1‑69 (MQ=255)
                                 ggAAGCCGTCGATGTGTTGCTCTCGCTGGAAAATGAACCCGTGCTGGTAAATGGCTGGATTGACAAGCa                <  1:692404/69‑1 (MQ=255)
                                  gAAGCCGTCGATGTGTTGCTCTCGCTGGAAAATGAACCCGTGCTGGTAAATGGCTGGATTGCCAAGCat               >  1:2124708/1‑69 (MQ=255)
                                  gAAGCCGTCGATGTGTTGCTCTCGCTGGAAAATGAACCCGTGCTGGTAAATGGCTGGATTGACAAGCat               >  1:3098782/1‑69 (MQ=255)
                                   aaGCCGTCGATGTGTTGCTCTCGCTGGAAAATGAAcc                                              >  1:2687343/1‑37 (MQ=255)
                                   aaGCCGTCGATGTGTTGCTCTCGCTGGAAAATGAACCCGTGCTGGTAAATGGCTGGATTGACAAGCat               >  1:620002/1‑68 (MQ=255)
                                             tgtgTTGCTCTCGCTGGAAAATGAACCCGTGCTGGTAAATGGCTGGATTGACAAGCATATGAATCCgg     <  1:381073/68‑1 (MQ=255)
                                             tgtgTTGCTCTCGCTGGAAAATGAACCCGTGCTGGTAAATGGCTGGATTGACAAGCATATGAATCCgg     <  1:202299/68‑1 (MQ=255)
                                             tgtgTTGCTCTCGCTGGAAAATGAACCCGTGCTGGTAAATGGCTGGATTGACAAGCATATGAATCCgg     <  1:2014757/68‑1 (MQ=255)
                                             tgtgTTGCTCTCGCTGGAAAATGAACCCGTGCTGGTAAATGGCTGGATTGACAAGCATATGAATCCgg     <  1:317548/68‑1 (MQ=255)
                                             tgtgTTGCTCTCGCTGGAAAATGAACCCGTGCTGGTAAATGGCTGGATTGACAAGCATATGAATCCgg     <  1:518622/68‑1 (MQ=255)
                                             tgtgTTGCTCTCGCTGGAAAATGAACCCGTGCTGGTAAATGGCTGGATTGACAAGCATATGAATCCgg     <  1:2448032/68‑1 (MQ=255)
                                             tgtgTTGCTCTCGCTGGAAAATGAACCCGTGCTGGTAAATGGCTGGATTGACAAGCATATGAATCCgg     <  1:313497/68‑1 (MQ=255)
                                             tgtgTTGCTCTCGCTGGAAAATGAACCCGTGCTGGTAAATGGCTGGATTGACAAGCATATGAATCCgg     <  1:2658776/68‑1 (MQ=255)
                                             tgtgTTGCTCTCGCTGGAAAATGAACCCGTGCTGGTAAATGGCTGGATTGACAAGCATATGAATCCgg     <  1:661034/68‑1 (MQ=255)
                                             tgtgTTGCTCTCGCTGGAAAATGAACCCGTGCTGGTAAATGGCTGGATTGACAAGCATATGAATCCgg     <  1:2777484/68‑1 (MQ=255)
                                             tgtgTTGCTCTCGCCGGAAAATGAACCCGTGCTGGTAAATGGCTGGATTGACAAGCATATGAATCCgg     <  1:283657/68‑1 (MQ=255)
                                              gtgtTGCTCTCGCTGGAAAATGAACCCGTGCTGGTAAATGGCTGGATTGACAAGCATATGAATCCgg     <  1:1396616/67‑1 (MQ=255)
                                               tgtTGCTCTCGCTGGAAAATGAACCCGTGCTGGTAAATGGCTGGATTGACAAGCATATGAATCCGGAAc  >  1:1103392/1‑69 (MQ=255)
                                               |                                                                    
AATCACCCGTTACATAGAAGCCAGTGCCGCCCAGGAAGCCGTCGATGTGTTGCTCTCGCTGGAAAATGAACCCGTGCTGGTAAATGGCTGGATTGACAAGCATATGAATCCGGAAC  >  minE/675264‑675379

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