Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 712,019 G→A 79.0% intergenic (‑132/+8) fldA ← / ← ybfE flavodoxin 1/LexA regulated protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb712,0190GA79.0% 74.5 / 17.7 38intergenic (‑132/+8)fldA/ybfEflavodoxin 1/LexA regulated protein
Reads supporting (aligned to +/- strand):  ref base G (3/5);  new base A (19/11);  total (22/16)
Fisher's exact test for biased strand distribution p-value = 2.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

TCATAGCAGAACAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGACTGCACACTGTGCTACATGAAAGTGGAAATTTAAACGATGCCCTGACTACGCAGCGCCGCCAGTTGCTGCATTAACATCTCTTCGAT  >  W3110S.gb/711952‑712083
                                                                   |                                                                
tCATAGCAGAACAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGACTGCACACTGTGCTACATGAAAATgg                                                               <  1:217279/71‑1 (MQ=255)
 cATAGCAGAACAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGACTGCACACTGTGCTACATGAAAATgg                                                               >  1:505616/1‑70 (MQ=255)
 cATAGCAGAACAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGACTGCACACTGTGCTACATGAAAATgg                                                               >  1:1167309/1‑70 (MQ=255)
 cATAGCAGAACAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGACTGCACACTGTGCTACATGAAAATgg                                                               >  1:1529908/1‑70 (MQ=255)
          aCAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGACTGCACACTGTGCTACATGAAAGTGGAAATTTAAAc                                                     <  1:2194163/71‑1 (MQ=255)
          aCAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGACTGCACACTGTGCTACATGAAAGTGGAAATTTAAAc                                                     <  1:184600/71‑1 (MQ=255)
          aCAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGACTGCACACTGTGCTACATGAAAATGGAAATTTAAAc                                                     <  1:2028753/71‑1 (MQ=255)
           cAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGACTGCACACTGTGCTACATGAAAATGGAAATTTAAAc                                                     >  1:1294556/1‑70 (MQ=255)
                    ttGCGGAATCAGCAAACGAGCAGGACTGCACACTGTGCTACATGAAAATGGAAATTTAAACGATGCCCTGa                                           >  1:2117734/1‑71 (MQ=255)
                    ttGCGGAATCAGCAAACGAGCAGGACTGCACACTGTGCTACATGAAAATGGAAATTTAAACGATGCCCTGa                                           >  1:2390207/1‑71 (MQ=255)
                            tCAGCAAACGAGCAGGACTGCACACTGTGCTACATGAAAATGGAAATTTAAAc                                                     <  1:1115693/53‑1 (MQ=255)
                            tCAGCAAACGAGCAGGACTGCACACTGTGCTACATGAAAATGGAAATTTAAACGATGCCCTGAc                                          <  1:2257942/64‑1 (MQ=255)
                                cAAACGAGCAGGACTGCACACTGTGCTACATGAAAATGGAAATTTAAACGa                                                   >  1:1017879/1‑51 (MQ=255)
                                  aaCGAGCAGGACTGCACACTGTGCTACATGAAAATGGAAATTTAAACGATGCCCTGGCTACGCAGcgccgc                             >  1:2554089/1‑71 (MQ=255)
                                  aaCGAGCAGGACTGCACACTGTGCTACATGAAAATGGAAATTTAAACGATGCCCTGACTACGCAGcgccgc                             >  1:1354146/1‑71 (MQ=255)
                                    cGAGCAGGACTGCACACTGTGCTACATGAAAATGGAAATTTAAACGATGCCCTGACTACGCAg                                   >  1:962228/1‑63 (MQ=255)
                                    cGAGCAGGACTGCACACTGTGCTACATGAAAATGGAAATTTAAACGATGCCCTGACTACGCAGcgccgc                             >  1:860957/1‑69 (MQ=255)
                                       gCAGGACTGCACACTGTGCTACATGAAAGTGGAAATTTAAACGATGCCCTGACTACGCAGCGCCGCCAGtt                        <  1:2486520/71‑1 (MQ=255)
                                           gACTGCACACTGTGCTACATGAAAATGGAAATTTAAACGATGCCCTGACTACGc                                     <  1:2189394/54‑1 (MQ=255)
                                            aCTGCACACTGTGCTACATGAAAATGGAAATTTaaa                                                      >  1:1782862/1‑36 (MQ=255)
                                            aCTGCACACTGTGCTACATGAAAATGGAAATTTAAACGATGCCCTGACTACGc                                     <  1:1096556/53‑1 (MQ=255)
                                             cTGCACACTGTGCTACATGAAAATGGAAATTTAAACGATGCCCTGACTACGCAGCGCCGCCAGTTGCTGCa                  >  1:1650150/1‑71 (MQ=255)
                                              tGCACACTGTGCTACATGAAAGTGGAAATTTAAACGATGCCCTGACTACGCAGCgccgcc                            <  1:178729/60‑1 (MQ=255)
                                              tGCACACTGTGCTACATGAAAGTGGAAATTTAAACGATGCCCTGACTACGCAGCgccgcc                            <  1:1327376/60‑1 (MQ=255)
                                                cacaCTGTGCTACATGAAAATGGAAATTTAAACGa                                                   >  1:1165886/1‑35 (MQ=255)
                                                cacaCTGTGCTACATGAAAATGGAAATTTAAACGATGcc                                               >  1:371769/1‑39 (MQ=255)
                                                cacaCTGTGCTACATGAAAATGGAAATTTAAACGATGCCCTGa                                           >  1:2441153/1‑43 (MQ=255)
                                                cacaCTGTGCTACATGAAAATGGAAATTTAAACGATGCCCTGACTACGCAGCGCCGCCAGTTGCTGCAtt                >  1:1720196/1‑70 (MQ=255)
                                                cacaCTGTGCTACATGAAAATGGAAATTTAAACGATGCCCTGACTACGCAGCGCCGCCAGTTGCTGCATTa               <  1:1960265/71‑1 (MQ=255)
                                                cacaCTGTGCTACATGAAAATGGAAATTTAAACGATGCCCTGACTACGCAGCGCCGCCAGTTGCTGCATTa               <  1:2031369/71‑1 (MQ=255)
                                                 acacTGTGCTACATGAAAATGGAAATTTAAACGATGCCCTGACTACGCAGCGCCGCCAGTTGCTGCATTa               <  1:2556103/70‑1 (MQ=255)
                                                 acacTGTGCTACATGAAAATGGAAATTTAAACGATGCCCTGACTACGCAGCGCCGCCAGTTGCTGCATTa               <  1:1040670/70‑1 (MQ=255)
                                                   acTGTGCTACATGAAAATGGAAATTTAAACGATGCCCTGACTACGCAGcgccgc                             <  1:1126199/54‑1 (MQ=255)
                                                      gtgCTACATGAAAATGGAAATTTAAACGATGCCCTGACTACGCAGCGCCGCCAGTTGCTGCATTAACAtct         >  1:542720/1‑71 (MQ=255)
                                                            cATGAAAGTGGAAATTTAAACGATGCCCTGACTACGCAGCGCCGCCAGTTGCTGCATTAACATCTCTTCGa   >  1:1759983/1‑71 (MQ=255)
                                                            cATGAAAATGGAAATTTAAACGATGCCCTGACTACGCAGCGCCGCCAGTTGCTGCATTAACATCTCTTCGa   >  1:1059947/1‑71 (MQ=255)
                                                               gAAAGTGGAAATTTAAACGATGCCCTGACTACGCAGCGCCGCCAGTTGCTGCATTAACATCTCTTCGAt  >  1:1436257/1‑69 (MQ=255)
                                                                aaaGTGGAAATTTAAACGATGCCCTGACTACGCAGCGCCGCCAGTTGCTGCa                  >  1:89229/1‑52 (MQ=255)
                                                                   |                                                                
TCATAGCAGAACAGGCAGTCTTGCGGAATCAGCAAACGAGCAGGACTGCACACTGTGCTACATGAAAGTGGAAATTTAAACGATGCCCTGACTACGCAGCGCCGCCAGTTGCTGCATTAACATCTCTTCGAT  >  W3110S.gb/711952‑712083

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: