Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,288,799 A→T 8.0% intergenic (+108/‑490) argE → / → ppc acetylornithine deacetylase/phosphoenolpyruvate carboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,288,7990AT 8.0% 172.6 / 5.3 61intergenic (+108/‑490)argE/ppcacetylornithine deacetylase/phosphoenolpyruvate carboxylase
Reads supporting (aligned to +/- strand):  ref base A (32/24);  new base T (3/2);  total (35/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGCCAAACTCCAGTGCCGCAATAATGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAGTGACTCAAACGATGCCCAACCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCACTGTTGCCAAACTCCAG  >  minE/2288733‑2288869
                                                                  |                                                                      
ttGCCAAACTCCAGTGCCGCAATAATGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAGTGACTc                                                                     >  1:292032/1‑70 (MQ=40)
    cAAACTCCAGTGCCGCAATAATGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAGTGACTCAAAc                                                                 >  1:4678212/1‑70 (MQ=40)
     aaaCTCCAGTGCCGCAATAATGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAGTGACTCAAAc                                                                 <  1:4639973/69‑1 (MQ=40)
     aaaCTCCAGTGCCGCAATAATGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAGTGACTCAAAc                                                                 <  1:2084466/69‑1 (MQ=40)
      aaCTCCAGTGCCGCAATAATGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAGTGACTCAAACGAt                                                              >  1:41196/1‑71 (MQ=39)
       aCTCCAGTGCCGCAATAATGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAGTGACTCAAACGAt                                                              <  1:1206124/70‑1 (MQ=40)
                ccGCAATAATGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAGTGACTCAAACGATGCCCAACCg                                                     <  1:3957682/70‑1 (MQ=255)
                ccGCAATAATGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAGTGACTCAAACGATGCCCAACCg                                                     <  1:4730046/70‑1 (MQ=255)
                  gCAATAATGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAGTGACTCAAACGATGCCCAACCGTAg                                                  >  1:2271043/1‑71 (MQ=255)
                   cAATAATGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAGTGACTCAAACGATGCCCAACCGTAgg                                                 >  1:4089315/1‑71 (MQ=255)
                       aatGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAGTGa                                                                        >  1:3434698/1‑44 (MQ=255)
                       aatGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAGTGa                                                                        >  1:4025586/1‑44 (MQ=255)
                        atGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAGTGACTCAAACGATGCCCAACCGTa                                                   >  1:4086561/1‑64 (MQ=255)
                          gTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAGTGAc                                                                       >  1:1219942/1‑42 (MQ=255)
                           tCGGATGCGATGCTTGCGCATCTTATCCGACCTACAGTGACTCAAACGa                                                               <  1:3164547/49‑1 (MQ=255)
                           tCGGATGCGATGCTTGCGCATCTTATCCGACCTACAGTGACTCAAACGa                                                               <  1:828356/49‑1 (MQ=255)
                           tCGGATGCGATGCTTGCGCATCTTATCCGACCTACAGTGACTCAAACGa                                                               <  1:1415168/49‑1 (MQ=255)
                           tCGGATGCGATGCTTGCGCATCTTATCCGACCTACAGTGACTCAAACGATGCCCAACCGTAGGCCGGATa                                          <  1:165461/70‑1 (MQ=255)
                                tgcgatgcTTGCGCATCTTATCCGACCTACAGTGACTCAAACGATGCCCAACCGTAGGCCGGATAAGgcgc                                    >  1:1323680/1‑71 (MQ=255)
                                 gcgatgcTTGCGCATCTTATCCGACCTACAGTGTCTCa                                                                    <  1:723654/38‑1 (MQ=255)
                                 gcgatgcTTGCGCATCTTATCCGACCTACAGTGTCTCa                                                                    <  1:714655/38‑1 (MQ=255)
                                   gatgcTTGCGCATCTTATCCGACCTACAGTGTCTCaaa                                                                  >  1:3602547/1‑38 (MQ=255)
                                   gatgcTTGCGCATCTTATCCGACCTACAGTGTCTCaaa                                                                  >  1:1811896/1‑38 (MQ=255)
                                   gatgcTTGCGCATCTTATCCGACCTACAGTGACTCAAACGATGCCCAACCGTAgg                                                 <  1:2761899/55‑1 (MQ=255)
                                   gatgcTTGCGCATCTTATCAGACCTACAGTGTCTCaaa                                                                  >  1:1487945/1‑38 (MQ=17)
                                    atgcTTGCGCATCTTATCCGACCTACAGTGACTCAAACGATGCCCAACCGTa                                                   <  1:1393461/52‑1 (MQ=255)
                                    atgcTTGCGCATCTTATCCGACCTACAGTGACTCAAACGATGCCCAACCGTa                                                   <  1:303343/52‑1 (MQ=255)
                                    atgcTTGCGCATCTTATCCGACCTACAGTGACTCAAACGATGCCCAACCGTa                                                   <  1:507319/52‑1 (MQ=255)
                                    atgcTTGCGCATCTTATCCGACCTACAGTGACTCAAACGATGCCCAACCGTa                                                   <  1:4365436/52‑1 (MQ=255)
                                    atgcTTGCGCATCTTATCCGACCTACAGTGACTCAAACGATGCCCAACCGTa                                                   <  1:3256708/52‑1 (MQ=255)
                                    atgcTTGCGCATCTTATCCGACCTACAGTGACTCAAACGATGCCCAACCGTa                                                   <  1:746785/52‑1 (MQ=255)
                                     tgcTTGCGCATCTTATCCGACCTACAGTGACTCAAAc                                                                 <  1:3154464/37‑1 (MQ=255)
                                     tgcTTGCGCATCTTATCCGACCTACAGTGACTCAAAc                                                                 <  1:4485481/37‑1 (MQ=255)
                                     tgcTTGCGCATCTTATCCGACCTACAGTGACTCAAACGa                                                               >  1:470643/1‑39 (MQ=255)
                                     tgcTTGCGCATCTTATCCGACCTACAGTGACTCAAACGa                                                               >  1:3312112/1‑39 (MQ=255)
                                     tgcTTGCGCATCTTATCCGACCTACAGTGACTCAAACGa                                                               >  1:1193485/1‑39 (MQ=255)
                                     tgcTTGCGCATCTTATCCGACCTACAGTGACTCAAACGa                                                               >  1:2283154/1‑39 (MQ=255)
                                     tgcTTGCGCATCTTATCCGACCTACAGTGACTCAAACGATGCCCAACCGTAGGCCGGATAAGGCGCTcgcg                               >  1:852187/1‑71 (MQ=255)
                                       cTTGCGCATCTTATCCGACCTACAGTGACTCAAAc                                                                 <  1:3710366/35‑1 (MQ=40)
                                       cTTGCGCATCTTATCCGACCTACAGTGACTCAAAc                                                                 <  1:1920263/35‑1 (MQ=40)
                                       cTTGCGCATCTTATCCGACCTACAGTGACTCAAAc                                                                 <  1:273004/35‑1 (MQ=40)
                                       cTTGCGCATCTTATCCGACCTACAGTGACTCAAAc                                                                 <  1:4409404/35‑1 (MQ=40)
                                       cTTGCGCATCTTATCCGACCTACAGTGACTCAAAc                                                                 <  1:4542857/35‑1 (MQ=40)
                                         tGCGCATCTTATCCGACCTACAGTGACTCAAACGa                                                               >  1:3942981/1‑35 (MQ=40)
                                         tGCGCATCTTATCCGACCTACAGTGACTCAAACGa                                                               >  1:1656935/1‑35 (MQ=40)
                                         tGCGCATCTTATCCGACCTACAGTGACTCAAACGa                                                               >  1:2035556/1‑35 (MQ=40)
                                                 ttATCCGACCTACAGTGACTCAAACGATGCCCAACCGTAGGCCGGATAAGGCGCTCGCGCCGCATc                        >  1:1560600/1‑66 (MQ=40)
                                                  tATCCGACCTACAGTGACTCAAACGATGCCCAACCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCAc                  >  1:2516694/1‑71 (MQ=39)
                                                  tATCCGACCTACAGTGACTCAAACGATGCCCAACCGTAGGCCGGATAAGGCGCTAGCGCCGCATCCGGCAc                  >  1:3836201/1‑71 (MQ=255)
                                                   aTCCGACCTACAGTGACTCAAACGATGCCCAACCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCgg                     <  1:2839825/67‑1 (MQ=40)
                                                     ccGACCTACAGTGACTCAAACGATGCCCAACCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCACTGt               >  1:259484/1‑71 (MQ=39)
                                                      cGACCTACAGTGACTCAAACGATGCCCAACCGTAGGCCGGATAAGgcg                                     >  1:2774505/1‑48 (MQ=255)
                                                      cGACCTACAGTGACTCAAACGATGCCCAACCGTAGGCCGGATAAGgcg                                     >  1:564729/1‑48 (MQ=255)
                                                      cGACCTACAGTGACTCAAACGATGCCCAACCGTAGGCCGGATAAGgcg                                     >  1:84495/1‑48 (MQ=255)
                                                      cGACCTACAGTGACTCAAACGATGCCCAACCGTAGGCCGGATAAGgcg                                     >  1:3693483/1‑48 (MQ=255)
                                                        aCCTACAGTGACTCAAACGATGCCCAACCGTAGGCCGGATaa                                         >  1:2659163/1‑42 (MQ=255)
                                                        aCCTACAGTGACTCAAACGATGCCCAACCGTAGGCCGGATaa                                         >  1:1742635/1‑42 (MQ=255)
                                                        aCCTACAGTGACTCAAACGATGCCCAACCGTAGGCCGGATaa                                         >  1:3065324/1‑42 (MQ=255)
                                                          cTACAGTGACTCAAACGATGCCCAACCGTAGGCCGGATAAGGCGCTAGCGCCGCATCCGGCACt                 >  1:2981420/1‑64 (MQ=255)
                                                             cAGTGACTCAAACGATGCCCAACCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCACTGTTGCCAAAc       >  1:450586/1‑71 (MQ=34)
                                                             cAGTGACTCAAACGATGCCCAACCGTAGGCCGGATAAGGAGCTCGCGCCGCAt                         >  1:4431192/1‑53 (MQ=255)
                                                             cAGTGACTCAAACGATGCCCAACCGTAGGCCGCATAAGGAGCTCGcgccg                            >  1:4431190/1‑50 (MQ=11)
                                                               gTGACTCAAACGATGCCCAACCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCACTGTTGCCAAACt      >  1:3322902/1‑70 (MQ=35)
                                                                  aCTCAAACGATGCCCAACCGTAGGCCGGATAAGGCGCTCGcgccgc                           >  1:1284468/1‑46 (MQ=38)
                                                                  aCTCAAACGATGCCCAACCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCACTGtt              >  1:1377330/1‑59 (MQ=35)
                                                                  aCTCAAACGATGCCCAACCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCACTGTTGCCAAACTCCAg  <  1:934853/71‑1 (MQ=34)
                                                                  |                                                                      
TTGCCAAACTCCAGTGCCGCAATAATGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAGTGACTCAAACGATGCCCAACCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCACTGTTGCCAAACTCCAG  >  minE/2288733‑2288869

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: