Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,168,168 C→A 26.8% G337G (GGC→GGA ynjE → predicted thiosulfate sulfur transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,168,1680CA26.8% 98.1 / 27.9 63G337G (GGC→GGAynjEpredicted thiosulfate sulfur transferase
Reads supporting (aligned to +/- strand):  ref base C (31/15);  new base A (8/9);  total (39/24)
Fisher's exact test for biased strand distribution p-value = 1.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAATTTATCGGTACGACCAGCGGTTACAGCTATATTAAACCAAAAGGTGAAATAGCCGGCGCACGTTGGGGACACGCTGGTAGCGACTCGACGCATATGGAAGATTTTCATAACCCGGAT  >  minE/1168109‑1168228
                                                           |                                                            
gAATTTATCGGTACGACCAGCGGTTACAGCTATATTAAACCAAAAGGTGAAATAGCCGGAGCACGTTg                                                      >  1:248057/1‑68 (MQ=255)
gAATTTATCGGTACGACCAGCGGTTACAGCTATATTAAACCAAAAGGTGAAATAGCCGGAGCACGTTg                                                      >  1:48195/1‑68 (MQ=255)
gAATTTATCGGTACGACCAGCGGTTACAGCTATATTAAACCAAAAGGTGAAATAGCCGGAGCACGTTg                                                      >  1:766919/1‑68 (MQ=255)
gAATTTATCGGTACGACCAGCGGTTACAGCTATATTAAACCAAAAGGTGAAATAGCCGGAGCACGTTg                                                      >  1:2914988/1‑68 (MQ=255)
gAATTTATCGGTACGACCAGCGGTTACAGCTATATTAAACCAAAAGGTGAAATAGCCGGAGCACGTTg                                                      >  1:3121370/1‑68 (MQ=255)
gAATTTATCGGTACGACCAGCGGTTACAGCTATATTAAACCAAAAGGTGAAATAGCCGGAGCACGTTg                                                      >  1:2964481/1‑68 (MQ=255)
gAATTTATCGGTACGACCAGCGGTTACAGCTATATTAAACCAAAAGGTGAAATAGCCGGAGCACGTTg                                                      >  1:992032/1‑68 (MQ=255)
gAATTTATCGGTACGACCAGCGGTTACAGCTATATTAAACCAAAAGGTGAAATAGCCGGAGCACGTTg                                                      >  1:2193413/1‑68 (MQ=255)
   tttATCGGTACGACCAGCGGTTACAGCTATATTAAACCAAAAGGTGAAATAGCCGGCGCACGTTGGGGa                                                  <  1:2974167/69‑1 (MQ=255)
   tttATCGGTACGACCAGCGGTTACAGCTATATTAAACCAAAAGGTGAAATAGCCGGCGCACGTTGGGGa                                                  <  1:3156853/69‑1 (MQ=255)
   tttATCGGTACGACCAGCGGTTACAGCTATATTAAACCAAAAGGTGAAATAGCCGGCGCACGTTGGGGa                                                  <  1:1554845/69‑1 (MQ=255)
   tttATCGGTACGACCAGCGGTTACAGCTATATTAAACCAAAAGGTGAAATAGCCGGCGCACGTTGGGGa                                                  <  1:1595493/69‑1 (MQ=255)
    ttATCGGTACGACCAGCGGTTACAGCTATATTAAACCAAAAGGTGAAATAGCCGGCGCACGTTGGGGa                                                  <  1:1914999/68‑1 (MQ=255)
     tATCGGTACGACCAGCGGTTACAGCTATATTAAACCAAAAGGTGAAATAGCCGGCGCACGTTGGGGa                                                  <  1:2326307/67‑1 (MQ=255)
           tACGACCAGCGGTTACAGCTATATTAAACCAAAAGGTGAAATAGCCGGCGCACGTTGGGGACACGCTgg                                          <  1:2139249/69‑1 (MQ=255)
           tACGACCAGCGGTTACAGCTATATTAAACCAAAAGGTGAAATAGCCGGCGCACGTTGGGGACACGCTgg                                          <  1:194949/69‑1 (MQ=255)
                      gTTACAGCTATATTAAACCAAAAGGTGAAATAGCCGGCGCACGTTGGGGACACGCTGGTAGCGACTCGa                               >  1:1678319/1‑69 (MQ=255)
                        tACAGCTATATTAAACCAAAAGGTGAAATAGCCGGCGCACGTTGGGGACACGCTGGTAGCGACTCGAc                              >  1:2975780/1‑68 (MQ=255)
                        tACAGCTATATTAAACCAAAAGGTGAAATAGCCGGCGCACGTTGGGGACACGCTGGTAGCGACTCGACg                             >  1:433365/1‑69 (MQ=255)
                        tACAGCTATATTAAACCAAAAGGTGAAATAGCCGGCGCACGTTGGGGACACGCTGGTAGCGACTCGACg                             >  1:3282161/1‑69 (MQ=255)
                        tACAGCTATATTAAACCAAAAGGTGAAATAGCCGGCGCACGTTGGGGACACGCTGGTAGCGACTCGACg                             >  1:1033203/1‑69 (MQ=255)
                        tACAGCTATATTAAACCAAAAGGTGAAATAGCCGGCGCACGTTGGGGACACGCTGGTAGCGACTCGACg                             >  1:520185/1‑69 (MQ=255)
                        tACAGCTATATTAAACCAAAAGGTGAAATAGCCGGCGCACGTTGGGGACACGCTGGTAGCGACTCGACg                             >  1:2041079/1‑69 (MQ=255)
                           aGCTATATTAAACCAAAAGGTGAAATAGCCGGCGCACGTTGGGGACACGCTGGTAGCGACTCGACGCa                           <  1:2318080/68‑1 (MQ=255)
                           aGCTATATTAAACCAAAAGGTGAAATAGCCGGCGCACGTTGGGGACACGCTGGTAGCGACTCGACGCa                           <  1:3062652/68‑1 (MQ=255)
                           aGCTATATTAAACCAAAAGGTGAAATAGCCGGCGCACGTTGGGGACACGCTGATAGCGACTCGACGCa                           <  1:2461581/68‑1 (MQ=255)
                              tataTTAAACCAAAAGGTGAAATAGCCGGCGCACGTTGTGTACACGCTGGTAGCGACTCGACGCa                           <  1:2185239/65‑1 (MQ=255)
                                tatTAAACCAAAAGGTGAAATAGCCGGCGCACGTTGGGGACACGCTGGTAGCGACTCGACGCatat                        >  1:1434172/1‑66 (MQ=255)
                                tatTAAACCAAAAGGTGAAATAGCCGGCGCACGTTGGGGACACGCTGGTAGCGACTCGACGCatat                        >  1:3127838/1‑66 (MQ=255)
                                tatTAAACCAAAAGGTGAAATAGCCGGCGCACGTTGGGGACACGCTGGTAGCGACTCGACGCatat                        >  1:2859808/1‑66 (MQ=255)
                                tatTAAACCAAAAGGTGAAATAGCCGGCGCACGTTGGGGACACGCTGGTAGCGACTCGACGCatat                        >  1:933559/1‑66 (MQ=255)
                                  ttAAACCAAAAGGTGAAATAGCCGGCGCACGTTGGGGACACGCTGGTAGCGACTCGACGCata                         >  1:557511/1‑63 (MQ=255)
                                  ttAAACCAAAAGGTGAAATAGCCGGCGCACGTTGGGGACACGCTGGTAGCGACTCGACGCata                         >  1:2448851/1‑63 (MQ=255)
                                  ttAAACCAAAAGGTGAAATAGCCGGCGCACGTTGGGGACACGCTGGTAGCGACTCGACGCata                         >  1:1355710/1‑63 (MQ=255)
                                  ttAAACCAAAAGGTGAAATAGCCGGCGCACGTTGGGGACACGCTGGTAGCGACTCGACGCata                         >  1:3123319/1‑63 (MQ=255)
                                  ttAAACCAAAAGGTGAAATAGCCGGCGCACGTTGGGGACACGCTGGTAGCGACTCGACGCata                         >  1:975856/1‑63 (MQ=255)
                                  ttAAACCAAAAGGTGAAATAGCCGGCGCACGTTGGGGACACGCTGGTAGCGACTCGACGCata                         >  1:2955112/1‑63 (MQ=255)
                                  ttAAACCAAAAGGTGAAATAGCCGGCGCACGTTGGGGACACGCTGGTAGCGACTCGACGCata                         >  1:125803/1‑63 (MQ=255)
                                  ttAAACCAAAAGGTGAAATAGCCGGCGCACGTTGGGGACACGCTGGTAGCGACTCGACGCata                         >  1:1165098/1‑63 (MQ=255)
                                  ttAAACCAAAAGGTGAAATAGCCGGCGCACGTTGGGGACACGCTGGTAGCGACTCGACGCata                         >  1:2270825/1‑63 (MQ=255)
                                  ttAAACCAAAAGGTGAAATAGCCGGCGCACGTTGGGGACACGCTGGTAGCGACTCGACGCat                          >  1:1117236/1‑62 (MQ=255)
                                  ttAAACCAAAAGGTGAAATAGCCGGCGCACGTTGGGGACACGCTGGTAGCGACTCGACGCATATGGaa                    >  1:3431953/1‑68 (MQ=255)
                                  ttAAACCAAAAGGTGAAATAGCCGGCGCACGTTGGGGACACGCTGGTAGCGACTCGACGCATATGGAAg                   >  1:2360073/1‑69 (MQ=255)
                                  ttAAACCAAAAGGTGAAATAGCCGGCGCACGTTGGGGACACGCTGGTAGCGACTCGACGCATATGGAAg                   >  1:1143446/1‑69 (MQ=255)
                                  ttAAACCAAAAGGTGAAATAGCCGGCGCACGTTGGGGACACGCTGGTAGCGACTCGACGCATATGGAAg                   >  1:1262377/1‑69 (MQ=255)
                                  ttAAACCAAAAGGTGAAATAGCCGGCGCACGTTGGGGACACGCTGGTAGCGACTCGACGCATATGGAAg                   >  1:1304222/1‑69 (MQ=255)
                                  ttAAACCAAAAGGTGAAATAGCCGGCGCACGTTGGGGACACGCTGGTAGCGACTCGACGCATATGGAAg                   >  1:1780803/1‑69 (MQ=255)
                                  ttAAACCAAAAGGTGAAATAGCCGGCGCACGTTGGGGACACGCTGGTAGCGACTCGACGCATATGGAAg                   >  1:2613304/1‑69 (MQ=255)
                                  ttAAACCAAAAGGTGAAATAGCCGGCGCACGTTGGGGACACGCTGGTAGCGACTCGACGCATATGGAAg                   >  1:2872349/1‑69 (MQ=255)
                                  ttAAACCAAAAGGTGAAATAGCAGGCGCACGTTGGGGACACGCTGGTAGCGACTCGACGCata                         >  1:3029810/1‑63 (MQ=255)
                                  atAAACCAAAAGGTGAAATAGCCGGCGCACGTTGGGGACACGCTGGTAGCGACTCGACGCATATGGAAg                   >  1:3293010/2‑69 (MQ=255)
                                    aaaCCAAAAGGTGAAATAGCCGGCGCACGTTGGGGACACGCTGGTAGCGACTCGACGCATATGGAAGa                  <  1:3103393/68‑1 (MQ=255)
                                       ccAAAAGGTGAAATAGCCGGAGCACGTTGGGGACACGCAGGTAGCGACTCGACGCATATGGAAGAttt               <  1:2447107/68‑1 (MQ=255)
                                       ccAAAAGGTGAAATAGCCGGAGCACGTTGGGGACACGCAGGTAGCGACTCGACGCATATGGAAGAttt               <  1:240250/68‑1 (MQ=255)
                                       ccAAAAGGTGAAATAGCCGGAGCACGTTGGGGACACGCAGGTAGCGACTCGACGCATATGGAAGAttt               <  1:1459738/68‑1 (MQ=255)
                                       ccAAAAGGTGAAATAGCCGGAGCACGTTGGGGACACGCAGGTAGCGACTCGACGCATATGGAAGAttt               <  1:725243/68‑1 (MQ=255)
                                       ccAAAAGGTGAAATAGCCGGAGCACGTTGGGGACACGCAGGTAGCGACTCGACGCATATGGAAGAttt               <  1:818206/68‑1 (MQ=255)
                                       ccAAAAGGTGAAATAGCCGGAGCACGTTGGGGACACGCAGGTAGCGACTCGACGCATATGGAAGAttt               <  1:833095/68‑1 (MQ=255)
                                       ccAAAAGGTGAAATAGCCGGAGCACGTTGGGGACACGCAGGTAGCGACTCGACGCATATGGAAGAttt               <  1:293612/68‑1 (MQ=255)
                                       ccAAAAGGTGAAATAGCCGGAGCACGTTGGGGACACGCAGGTAGCGACTCGACGCATATGGAAGAttt               <  1:3013491/68‑1 (MQ=255)
                                                    tAGCCGGCGCACGTTGGGGACACGCTGGTAGCGACTCGACGCATATGGAAGATTTTCATAACCCGGAt  <  1:721680/68‑1 (MQ=255)
                                                    tAGCCGGCGCACGTTGGGGACACGCTGGTAGCGACTCGACGCATATGGAAGATTTTCATAACCCGGAt  <  1:141261/68‑1 (MQ=255)
                                                       ccGGAGCACGTTGGGGACACGCAGGTAGCGACTCGACGCATATGGAAGAttt               <  1:2796902/52‑1 (MQ=255)
                                                           |                                                            
GAATTTATCGGTACGACCAGCGGTTACAGCTATATTAAACCAAAAGGTGAAATAGCCGGCGCACGTTGGGGACACGCTGGTAGCGACTCGACGCATATGGAAGATTTTCATAACCCGGAT  >  minE/1168109‑1168228

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: