Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,858,376 A→T 24.7% G402G (GGT→GGA cysI ← sulfite reductase, beta subunit, NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,858,3760AT24.7% 145.6 / 46.9 93G402G (GGT→GGAcysIsulfite reductase, beta subunit, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base A (35/35);  new base T (18/5);  total (53/40)
Fisher's exact test for biased strand distribution p-value = 2.77e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGATTCTGGTTCGCCGTAATGCGGAAATCGCCTTTGTGGATCTT  >  minE/1858313‑1858432
                                                               |                                                        
ggCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGc                                                     >  1:173308/1‑69 (MQ=255)
ggCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGc                                                     <  1:2482406/69‑1 (MQ=255)
ggCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGc                                                     <  1:1540381/69‑1 (MQ=255)
ggCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGc                                                     <  1:1531970/69‑1 (MQ=255)
ggCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGc                                                     <  1:364255/69‑1 (MQ=255)
ggCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGc                                                     <  1:2825282/69‑1 (MQ=255)
ggCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGc                                                     <  1:456912/69‑1 (MQ=255)
 gCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGc                                                     <  1:2183075/68‑1 (MQ=255)
   attcattAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCgat                                                  >  1:2713394/1‑69 (MQ=255)
   attcattAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCgat                                                  >  1:2030295/1‑69 (MQ=255)
      cattAACCCGCTCTCTTTGGCTATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCgatg                                                 >  1:329765/1‑67 (MQ=255)
      cattAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCgatg                                                 >  1:390553/1‑67 (MQ=255)
      cattAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCgatg                                                 >  1:2083996/1‑67 (MQ=255)
      cattAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCgatg                                                 >  1:1905422/1‑67 (MQ=255)
      cattAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCgatg                                                 >  1:2317484/1‑67 (MQ=255)
      cattAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCgatg                                                 >  1:519828/1‑67 (MQ=255)
          aaCCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCgatgat                                               >  1:1167329/1‑65 (MQ=255)
          aaCCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCgatgat                                               >  1:171466/1‑65 (MQ=255)
          aaCCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCgatgat                                               >  1:1670622/1‑65 (MQ=255)
          aaCCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCgatgat                                               >  1:158950/1‑65 (MQ=255)
          aaCCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCgatgat                                               >  1:338225/1‑65 (MQ=255)
          aaCCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCgatgat                                               >  1:905903/1‑65 (MQ=255)
          aaCCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCgatgat                                               >  1:1532232/1‑65 (MQ=255)
          aaCCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCgatgat                                               >  1:1485105/1‑65 (MQ=255)
          aaCCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCgatgat                                               >  1:81136/1‑65 (MQ=255)
          aaCCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCgat                                                  >  1:871216/1‑62 (MQ=255)
          aaCCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGa                                           <  1:220515/69‑1 (MQ=255)
          aaCCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGa                                           >  1:461433/1‑69 (MQ=255)
          aaCCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGa                                           >  1:2543748/1‑69 (MQ=255)
          aaCCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGa                                           <  1:565222/69‑1 (MQ=255)
          aaCCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGa                                           <  1:723433/69‑1 (MQ=255)
          aaCCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGa                                           <  1:1037969/69‑1 (MQ=255)
          aaCCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGa                                           <  1:966399/69‑1 (MQ=255)
          aaCCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGa                                           >  1:1167940/1‑69 (MQ=255)
          aaCCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGa                                           <  1:1209235/69‑1 (MQ=255)
          aaCCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGa                                           <  1:796760/69‑1 (MQ=255)
          aaCCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGa                                           >  1:1343071/1‑69 (MQ=255)
          aaCCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGa                                           >  1:2637952/1‑69 (MQ=255)
          aaCCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGa                                           <  1:2848806/69‑1 (MQ=255)
          aaCCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGa                                           <  1:3201970/69‑1 (MQ=255)
          aaCCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGa                                           <  1:1492565/69‑1 (MQ=255)
                       tGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGATTCTGGTTCGc                                <  1:3195613/67‑1 (MQ=255)
                         gCGATCTTGTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGATTCTGGTTCGCCGt                             <  1:1495303/68‑1 (MQ=255)
                         gCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGa                                           <  1:2004126/54‑1 (MQ=255)
                         gCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGATTCTGGTTCGCCg                              <  1:2000730/67‑1 (MQ=255)
                         gCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGATTCTGGTTCGCCGt                             <  1:2556841/68‑1 (MQ=255)
                         gCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGATTCTGGTTCGCCGt                             <  1:2926435/68‑1 (MQ=255)
                         gCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGATTCTGGTTCGCCGt                             <  1:1881407/68‑1 (MQ=255)
                         gCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGATTCTGGTTCGCCGt                             <  1:1430509/68‑1 (MQ=255)
                         gCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGCACACCGGCGATGATCAGATTCTGGTTCGCCGt                             <  1:457201/68‑1 (MQ=255)
                             tcttctCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGATTCTGGTTCGCCGt                             <  1:786386/64‑1 (MQ=255)
                             tcttctCGATCTTCGCTTTCTCGCTTTCCGGCACTCCGGCGATGATCAGATTc                                        >  1:1802188/1‑53 (MQ=255)
                             tcttctCGATCTTCGCTTTCTCGCTTTCCGGCACTCCGGCGATGATCAGATTCTGGTTCGCCGTAATg                         >  1:3366582/1‑68 (MQ=255)
                             tcttctCGATCTTCGCTTTCTCGCTTTCCGGCACTCCGGCGATGATCAGATTCTGGTTCGCCGTAATg                         >  1:988154/1‑68 (MQ=255)
                             tcttctCGATCTTCGCTTTCTCGCTTTCCGGCACTCCGGCGATGATCAGATTCTGGTTCGCCGTAATg                         >  1:567480/1‑68 (MQ=255)
                             tcttctCGATCTTCGCTTTCTCGCTTTCCGGCACTCCGGCGATGATCAGATTCTGGTTCGCCGTAATg                         >  1:3214832/1‑68 (MQ=255)
                             tcttctCGATCTTCGCTTTCTCGCTTTCCGGCACTCCGGCGATGATCAGATTCTGGTTCGCCGTAATg                         >  1:569793/1‑68 (MQ=255)
                             tcttctCGATCTTCGCTTTCTCGCTTTCCGGCACTCCGGCGATGATCAGATTCTGGTTCGCCGTAATg                         >  1:699227/1‑68 (MQ=255)
                             tcttctCGATCTTCGCTTTCTCGCTTTCCGGCACTCCGGCGATGATCAGATTCTGGTTCGCCGTAATg                         >  1:979240/1‑68 (MQ=255)
                             tcttctCGATCTTCGCTTTCTCGCTTTCCGGCACTCCGGCGATGATCAGATTCTGGTTCGCCGTAATg                         >  1:3137671/1‑68 (MQ=255)
                             tcttctCGATCTTCGCTTTCTCGCTTTCCGGCACTCCGGCGATGATCAGATTCTGGTTCGCCGTAATg                         >  1:2069225/1‑68 (MQ=255)
                             tcttctCGATCTTCGCTTTCTCGCTTTCCGGCACTCCGGCGATGATCAGATTCTGGTTCGCCGTAATg                         >  1:1102506/1‑68 (MQ=255)
                             tcttctCGATCTTCGCTTTCTCGCTTTCCGGCACTCCGGCGATGATCAGATTCTGGTTCGCCGTAATg                         >  1:1154398/1‑68 (MQ=255)
                             tcttctCGATCTTCGCTTTCTCGCTTTCCGGCACTCCGGCGATGATCAGATTCTGGTTCGCCGTAATg                         >  1:2756697/1‑68 (MQ=255)
                             tcttctCGATCTTCGCTTTCTCGCTTTCCGGCACTCCGGCGATGATCAGATTCTGGTTCGCCGTAATg                         >  1:2157302/1‑68 (MQ=255)
                             tcttctCGATCTTCGCTTTCTCGCTTTCCGGCACTCCGGCGATGATCAGATTCTGGTTCGCCGTAATg                         >  1:131053/1‑68 (MQ=255)
                             tcttctCGATCTTCGCTTTCTCGCTTTCCGGCACTCCGGCGATGATCAGATTCTGGTTCGCCGTAATg                         >  1:2569412/1‑68 (MQ=255)
                             tcttctCGATCTTCGCTTTCTCGCTTTCCGGCACTCCGGCGATGATCAGATTCTGGTTCGCCGTAATg                         >  1:1386006/1‑68 (MQ=255)
                             tcttctCGATCTTCGCTTTCTCGCTTTCCGGCACTCCGGCGATGATCAGATTCTGGTTCGCCGTAATg                         >  1:1553003/1‑68 (MQ=255)
                                          cgctttcttgctttcCGGTACACCGTCGATGATCAGATTCTGGTTCGCCGt                             <  1:2452664/51‑1 (MQ=39)
                                          cgctttctggctttcCGGTACACCGGCGATGATCAGATTCTGGTTCGCCGt                             >  1:2013790/1‑51 (MQ=255)
                                          cgctttctcgctttcCGGTACACCGGCGATGATCAGATTCTGGTTCGCCGt                             <  1:1414499/51‑1 (MQ=255)
                                          cgctttctcgctttcCGGTACACCGGCGATGATCAGATTCTGGTTCGCCGt                             >  1:1443016/1‑51 (MQ=255)
                                          cgctttctcgctttcCGGTACACCGGCGATGATCAGATTCTGGTTCGCCGt                             <  1:710493/51‑1 (MQ=255)
                                          cgctttctcgctttcCGGTACACCGGCGATGATCAGATTCTGGTTCGCCGt                             >  1:1254684/1‑51 (MQ=255)
                                          cgctttctcgctttcCGGTACACCGGCGATGATCAGATTCTGGTTCGCCGt                             >  1:1079050/1‑51 (MQ=255)
                                          cgctttctcgctttcCGGTACACCGGCGATGATCAGATTCTGGTTCGCCGt                             >  1:3071614/1‑51 (MQ=255)
                                          cgctttctcgctttcCGGTACACCGGCGATGATCAGATTCTGGTTCGCCGt                             >  1:1584710/1‑51 (MQ=255)
                                          cgctttctcgctttcCGGTACACCGGCGATGATCAGATTCTGGTTCGCCGt                             >  1:2187576/1‑51 (MQ=255)
                                          cgctttctcgctttcCGGTACACCGGCGATGATCAGATTCTGGTTCGCCGt                             >  1:2240227/1‑51 (MQ=255)
                                          cgctttctcgctttcCGGTACACCGGCGATGATCAGATTCTGGTTCGCCGt                             <  1:3283095/51‑1 (MQ=255)
                                          cgctttctcgctttcCGGTACACCGGCGATGATCAGATTCTGGTTCGCCGt                             <  1:2819092/51‑1 (MQ=255)
                                          cgctttctcgctttcCGGTACACCGGCGATGATCAGATTCTGGTTCGCCGt                             >  1:2865766/1‑51 (MQ=255)
                                          cgctttctcgctttcCGGTACACCGGCGATGATCAGATTCTGGTTCGCCGt                             >  1:3118677/1‑51 (MQ=255)
                                          cgctttctcgctttcCGGTACACCGGCGATGATCAGATTCTGGATCGCCGt                             <  1:1265128/51‑1 (MQ=255)
                                             tttctcgctttcCGGCACTCCGGCGATGATCAGATTCTGGTTCGCCGTAATGCGGAAATCACCTTtgt         <  1:394661/68‑1 (MQ=255)
                                             tttctcgctttcCGGCACTCCGGCGATGATCAGATTCTGGTTCGCCGTAATGCGGAAATCACCTTtgt         <  1:2307096/68‑1 (MQ=255)
                                             tttctcgctttcCGGCACTCCGGCGATGATCAGATTCTGGTTCGCCGTAATGCGGAAATCACCTTtgt         <  1:356476/68‑1 (MQ=255)
                                             tttctcgctttcCGGCACTCCGGCGATGATCAGATTCTGGTTCGCCGTAATGCGGAAATCACCTTtgt         <  1:788597/68‑1 (MQ=255)
                                             tttctcgctttcCGGCACTCCGGCGATGATCAGATTCTGGTTCGCCGTAATGCGGAAATCACCTTtgt         <  1:1056363/68‑1 (MQ=255)
                                              ttctcgctttcCGGTACACCGGCGATGATCAGATTCTGGTTCGCCGTAATGCGGAAATCGCCTTtgtg        >  1:2635344/1‑68 (MQ=255)
                                                     tttcCGGTACACCGGCGATGATCAGATTCTGGTTCGCCGTAATGCGGAAATCGCCTTTGTGGATCtt  <  1:70372/67‑1 (MQ=255)
                                                     tttcCGGTACACCGGCGATGATCAGATTCTGGTTCGCCGTAATGCGGAAATCGCCTTTGTGGATCtt  <  1:1250261/67‑1 (MQ=255)
                                                               |                                                        
GGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGATTCTGGTTCGCCGTAATGCGGAAATCGCCTTTGTGGATCTT  >  minE/1858313‑1858432

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: