Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,168,459 C→T 25.9% S434S (AGC→AGT ynjE → predicted thiosulfate sulfur transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,168,4590CT25.9% 115.8 / 34.9 96S434S (AGC→AGTynjEpredicted thiosulfate sulfur transferase
Reads supporting (aligned to +/- strand):  ref base C (37/34);  new base T (8/17);  total (45/51)
Fisher's exact test for biased strand distribution p-value = 1.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCC  >  minE/1168398‑1168520
                                                             |                                                             
aCGAATGGAGCAGCGATCGAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAAc                                                        >  1:1395062/1‑69 (MQ=255)
aCGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAAc                                                        >  1:2269456/1‑69 (MQ=255)
aCGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAAc                                                        >  1:2984604/1‑69 (MQ=255)
aCGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAAc                                                        >  1:2988156/1‑69 (MQ=255)
aCGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAAc                                                        >  1:807647/1‑69 (MQ=255)
aCGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAAc                                                        >  1:1716861/1‑69 (MQ=255)
aCGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAAc                                                        >  1:2599239/1‑69 (MQ=255)
aCGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAAc                                                        >  1:1044434/1‑69 (MQ=255)
aCGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAAAAAc                                                        >  1:3041673/1‑69 (MQ=255)
       gAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCa                                                  >  1:3129455/1‑68 (MQ=255)
       gAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCa                                                  >  1:2069387/1‑68 (MQ=255)
       gAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCa                                                  >  1:2941720/1‑68 (MQ=255)
       gAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCa                                                  >  1:468379/1‑68 (MQ=255)
       gAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCa                                                  >  1:597315/1‑68 (MQ=255)
       gAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCa                                                  >  1:2125747/1‑68 (MQ=255)
       gAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCa                                                  >  1:1418274/1‑68 (MQ=255)
       gAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCa                                                  >  1:1153787/1‑68 (MQ=255)
         gcagcGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAAt                                                <  1:1995535/68‑1 (MQ=255)
            gcGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTAc                                             >  1:2331810/1‑68 (MQ=255)
            gcGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACt                                            >  1:3008120/1‑69 (MQ=255)
            gcGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACt                                            >  1:2163516/1‑69 (MQ=255)
            gcGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACt                                            >  1:2412216/1‑69 (MQ=255)
            gcGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACt                                            >  1:2435591/1‑69 (MQ=255)
            gcGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACt                                            >  1:2394353/1‑69 (MQ=255)
            gcGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACt                                            >  1:3155666/1‑69 (MQ=255)
            gcGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACt                                            >  1:3145821/1‑69 (MQ=255)
            gcGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACt                                            >  1:1881418/1‑69 (MQ=255)
            gcGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACt                                            >  1:1877416/1‑69 (MQ=255)
            gcGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACt                                            >  1:1839303/1‑69 (MQ=255)
            gcGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACt                                            >  1:354998/1‑69 (MQ=255)
            gcGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACt                                            >  1:1771773/1‑69 (MQ=255)
            gcGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACt                                            >  1:2682051/1‑69 (MQ=255)
            gcGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACt                                            >  1:837873/1‑69 (MQ=255)
            gcGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACt                                            >  1:917445/1‑69 (MQ=255)
            gcGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACt                                            >  1:2943838/1‑69 (MQ=255)
            gcGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAAATACt                                            >  1:1384250/1‑69 (MQ=255)
            gcGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATAGCAATTACt                                            >  1:185341/1‑69 (MQ=255)
                 ccAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATaa                                                         <  1:2458633/51‑1 (MQ=255)
                 ccAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATaa                                                         <  1:2294612/51‑1 (MQ=255)
                 ccAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGCAGTAAATAACATTGCAATTTCTGACGc                                       >  1:2638737/1‑69 (MQ=255)
                 ccAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGCAGTAAATAACATTGCAATTACTGACGc                                       >  1:1113496/1‑69 (MQ=255)
                 ccAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGCAGTAAATAACATTGCAATTACTGACGc                                       >  1:223873/1‑69 (MQ=255)
                 ccAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGCAGTAAATAACATTGCAATTACTGACGc                                       >  1:2040188/1‑69 (MQ=255)
                 ccAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGCAGTAAATAACATTGCAATTACTGACGc                                       >  1:2021985/1‑69 (MQ=255)
                 ccAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGCAGTAAATAACATTGCAATTACTGACGc                                       >  1:498655/1‑69 (MQ=255)
                 ccAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGCAGTAAATAACATTGCAATTACTGACGc                                       >  1:1580239/1‑69 (MQ=255)
                 ccAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGCAGTAAATAACATTGCAATTACTGACGc                                       >  1:1548961/1‑69 (MQ=255)
                    aaaaaTCCGGTAGCAACCGGTGAACGCGGCCCGTACAGTAGCAAATAACATTGc                                                   <  1:280391/54‑1 (MQ=255)
                    aaaaaTCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGc                                                   <  1:2013539/54‑1 (MQ=255)
                          ccGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGAc                              <  1:2942014/69‑1 (MQ=255)
                          ccGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGAc                              <  1:3172539/69‑1 (MQ=255)
                          ccGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGAc                              <  1:2751791/69‑1 (MQ=255)
                          ccGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGAc                              <  1:2586762/69‑1 (MQ=255)
                             gTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTc                           <  1:2643767/69‑1 (MQ=255)
                             gTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTc                           <  1:2694872/69‑1 (MQ=255)
                             gTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTc                           <  1:1372997/69‑1 (MQ=255)
                             gTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTc                           <  1:774708/69‑1 (MQ=255)
                              tAGCAACCGGTGAACGCGGCCCGGACAGCAGTAAATAACATTGCAATTACTGACGCT‑‑AGCGACTTCAGc                        <  1:3316742/69‑1 (MQ=255)
                                gCAACCGGTGAACGCGGCCCGGACAGCAGTAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGc                        <  1:2626279/69‑1 (MQ=255)
                                gCAACCGGTGAACGCGGCCCGGACAGCAGTAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGc                        <  1:845176/69‑1 (MQ=255)
                                gCAACCGGTGAACGCGGCCCGGACAGCAGTAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGc                        <  1:1405739/69‑1 (MQ=255)
                                gCAACCGGTGAACGCGGCCCGGACAGCAGTAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGc                        <  1:562118/69‑1 (MQ=255)
                                gCAACCGGTGAACGCGGCCCGGACAGCAGTAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGc                        <  1:1620072/69‑1 (MQ=255)
                                gCAACCGGTGAACGCGGCCCGGACAGCAGTAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGc                        <  1:1670838/69‑1 (MQ=255)
                                gCAACCGGTGAACGCGGCCCGGACAGCAGTAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGc                        <  1:413090/69‑1 (MQ=255)
                                gCAACCGGTGAACGCGGCCCGGACAGCAGTAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGc                        <  1:352447/69‑1 (MQ=255)
                                gCAACCGGTGAACGCGGCCCGGACAGCAGTAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGc                        <  1:1828302/69‑1 (MQ=255)
                                gCAACCGGTGAACGCGGCCCGGACAGCAGTAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGc                        <  1:3019911/69‑1 (MQ=255)
                                gCAACCGGTGAACGCGGCCCGGACAGCAGTAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGc                        <  1:2470075/69‑1 (MQ=255)
                                gCAACCGGTGAACGCGGCCCGGACAGCAGTAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGc                        <  1:2469935/69‑1 (MQ=255)
                                gCAACCGGTGAACGCGGCCCGGACAGCAGTAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGc                        <  1:2390263/69‑1 (MQ=255)
                                gCAACCGGTGAACGCGGCCCGGACAGCAGTAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGc                        <  1:2369145/69‑1 (MQ=255)
                                 cAACCGGTGAACGCGGCCCGGACAGCAGTAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGc                        <  1:845136/68‑1 (MQ=255)
                                 cAACCGGTGAACGCGGCCCGGACAGCAGTAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGc                        <  1:585009/68‑1 (MQ=255)
                                      ggTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGa                  <  1:1948019/69‑1 (MQ=255)
                                      ggTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGa                  <  1:1605668/69‑1 (MQ=255)
                                      ggTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGa                  <  1:591554/69‑1 (MQ=255)
                                      ggTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGa                  <  1:2155405/69‑1 (MQ=255)
                                       gTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGa                  <  1:976756/68‑1 (MQ=255)
                                       gTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGa                  <  1:2624896/68‑1 (MQ=255)
                                       gTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGa                  <  1:454838/68‑1 (MQ=255)
                                       gTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGa                  <  1:2142980/68‑1 (MQ=255)
                                         gAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTc                           <  1:1942522/57‑1 (MQ=255)
                                          aaCGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATAcc              <  1:295968/69‑1 (MQ=255)
                                          aaCGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATAcc              <  1:2121706/69‑1 (MQ=255)
                                          aaCGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATAcc              <  1:2843736/69‑1 (MQ=255)
                                          aaCGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATAcc              <  1:777815/69‑1 (MQ=255)
                                          aaCGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATAcc              <  1:920348/69‑1 (MQ=255)
                                          aaCGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGCAGCGACTTCAGCGTCAGATAcc              <  1:3142246/69‑1 (MQ=255)
                                            cgcgGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCa            <  1:1781067/69‑1 (MQ=255)
                                            cgcgGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCa            <  1:349310/69‑1 (MQ=255)
                                            cgcgGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCa            <  1:3169764/69‑1 (MQ=255)
                                            cgcgGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCa            <  1:1420412/69‑1 (MQ=255)
                                            cgcgGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCa            <  1:2727957/69‑1 (MQ=255)
                                            cgcgGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCa            <  1:2093501/69‑1 (MQ=255)
                                                      aCAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACAccc  >  1:1828891/1‑69 (MQ=255)
                                                             |                                                             
ACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCC  >  minE/1168398‑1168520

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