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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,168,4560TC24.2% 126.3 / 35.9 103S433S (AGT→AGCynjEpredicted thiosulfate sulfur transferase
Reads supporting (aligned to +/- strand):  ref base T (44/34);  new base C (8/17);  total (52/51)
Fisher's exact test for biased strand distribution p-value = 4.03e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 103 total aligned reads displayed.

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