Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 597,552 C→T 19.7% intergenic (+133/‑106) serS → / → dmsA seryl‑tRNA synthetase, also charges selenocysteinyl‑tRNA with serine/dimethyl sulfoxide reductase, anaerobic, subunit A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE597,5520CT19.7% 151.0 / 29.0 66intergenic (+133/‑106)serS/dmsAseryl‑tRNA synthetase, also charges selenocysteinyl‑tRNA with serine/dimethyl sulfoxide reductase, anaerobic, subunit A
Reads supporting (aligned to +/- strand):  ref base C (33/20);  new base T (8/5);  total (41/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

TAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAACTGCTATAGAAATAATTACACAATACGGTTTGTTAC  >  minE/597488‑597622
                                                                |                                                                      
taatTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAat                                                                    >  1:382251/1‑69 (MQ=255)
 aatTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAattat                                                                 >  1:1249065/1‑71 (MQ=255)
   ttACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTTAAATTATTc                                                               <  1:702437/71‑1 (MQ=255)
       tcctcACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATg                                                           >  1:1005625/1‑71 (MQ=255)
         ctcACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGta                                                         >  1:1565488/1‑71 (MQ=255)
           cACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGtata                                                       >  1:1702382/1‑71 (MQ=255)
           cACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGtata                                                       >  1:1400031/1‑71 (MQ=255)
           cACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGtata                                                       >  1:1402483/1‑71 (MQ=255)
           cACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGtata                                                       >  1:1553718/1‑71 (MQ=255)
           cACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGtat                                                        >  1:847557/1‑70 (MQ=255)
            aCTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCa                                                                       <  1:658257/54‑1 (MQ=255)
            aCTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGta                                                         <  1:474940/68‑1 (MQ=255)
              ttACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCaaa                                                                     >  1:247508/1‑54 (MQ=255)
                acacGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGc                                                  >  1:2060016/1‑71 (MQ=255)
                acacGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGc                                                  >  1:414889/1‑71 (MQ=255)
                  acGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCt                                                <  1:557249/71‑1 (MQ=255)
                   cGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTTAAATTATTCGATGTATACAAGCCta                                               <  1:1750073/71‑1 (MQ=255)
                   cGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGcc                                                 >  1:2169223/1‑69 (MQ=255)
                   cGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCta                                               >  1:504270/1‑71 (MQ=255)
                   cGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCta                                               <  1:926970/71‑1 (MQ=255)
                   cGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCta                                               >  1:2171745/1‑71 (MQ=255)
                    gTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCta                                               <  1:1632700/70‑1 (MQ=255)
                     tAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCtata                                             >  1:855856/1‑71 (MQ=255)
                      aaTACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCtata                                             >  1:1037203/1‑70 (MQ=255)
                        tactacTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCg                                                              <  1:403372/51‑1 (MQ=255)
                        tactacTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCg                                                              <  1:959570/51‑1 (MQ=255)
                          ctactTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTc                                                               >  1:1655664/1‑48 (MQ=255)
                           tactTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCtata                                             >  1:1225869/1‑65 (MQ=255)
                           tactTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCtata                                             >  1:2085039/1‑65 (MQ=255)
                              tttCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCg                                        <  1:1094580/67‑1 (MQ=255)
                                tCGAGTGAAAATCTACCTATCTCTTTGATTTTTAAATTATTCGATGTATACAAGCCTATATAGCGAACTGc                                  >  1:150065/1‑71 (MQ=255)
                                tCGAGTGAAAATCTACCTATCTCTTTGATTTTTAAATTATTCGATGTATACAAGCCTATATAGCGAACTGc                                  >  1:287919/1‑71 (MQ=255)
                                tCGAGTGAAAATCTACCTATCTCTTTGATTTTTAAATTATTCGATGTATACAAGCCTATATAGCGAACTGc                                  >  1:365557/1‑71 (MQ=255)
                                tCGAGTGAAAATCTACCTATCTCTTTGATTTTCaaa                                                                     >  1:1152158/1‑36 (MQ=255)
                                tCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGtata                                                       >  1:1966021/1‑50 (MQ=255)
                                tCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAACTGc                                  >  1:634282/1‑71 (MQ=255)
                                tCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAACTGc                                  >  1:364404/1‑71 (MQ=255)
                                     tGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATg                                                           <  1:422594/41‑1 (MQ=255)
                                      gAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAACTGCt                                 <  1:1120637/66‑1 (MQ=255)
                                       aaaaTCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCtata                                             >  1:1031363/1‑53 (MQ=255)
                                       aaaaTCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCg                                        <  1:326385/58‑1 (MQ=255)
                                       aaaaTCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCg                                        <  1:1424708/58‑1 (MQ=255)
                                       aaaaTCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCg                                        <  1:727307/58‑1 (MQ=255)
                                       aaaaTCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAACTGCTATAGa                            <  1:1498631/70‑1 (MQ=255)
                                         aaTCTACCTATCTCTTAGATTTTTAAATTATTCGATGTATACAAGCCTATATAGCGAACTGCTATAg                             >  1:1729619/1‑67 (MQ=255)
                                            cTACCTATCTCTTTGATTTTTAAATTATTCGATGTATACAAGCCta                                               >  1:1384285/1‑46 (MQ=255)
                                            cTACCTATCTCTTTGATTTTTAAATTATTCGATGTATAAAAGCCta                                               >  1:1509938/1‑46 (MQ=255)
                                             taacTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAACTGCTATAGAaat                         <  1:1839182/64‑1 (MQ=255)
                                             tACCTATCTCTTTGATTTTTAAATTATTCGATGTATACAAGcc                                                 <  1:2134689/43‑1 (MQ=255)
                                             tACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGaa                                      >  1:80426/1‑54 (MQ=255)
                                             tACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGaa                                      >  1:1596873/1‑54 (MQ=255)
                                              aCCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCtata                                             <  1:1448138/46‑1 (MQ=255)
                                               ccTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAACTGCTATAGAAATAATTa                    >  1:331628/1‑70 (MQ=255)
                                               ccTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAACTGCTATAGAAATAATTa                    <  1:1789532/70‑1 (MQ=255)
                                                  aTCTCTTTTATTTTCAAATTATTCGATGTATACAAGCCta                                               <  1:1564497/40‑1 (MQ=255)
                                                  aTCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAAc                                     >  1:147264/1‑50 (MQ=255)
                                                  aTCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAAc                                     >  1:1050988/1‑50 (MQ=255)
                                                     tctTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAACTGCTATAGAAATAATTACACAATa             >  1:2096977/1‑71 (MQ=255)
                                                     tctTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAACTGCTATAGAAATAATTACACAATa             >  1:863494/1‑71 (MQ=255)
                                                       tttGATTTTTAAATTATTCGATGTATACAAGCCTATATAGCGAACTGCTATAGAaataa                       >  1:999994/1‑59 (MQ=255)
                                                       tttGATTTTTAAATTATTCGATGTATACAAGCATATATAGCGAACTGCTATAGAaataa                       >  1:2117874/1‑59 (MQ=255)
                                                       tttGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCg                                        <  1:2135562/42‑1 (MQ=255)
                                                       tttGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCAAACTGCTATAGAAATAATTACACAATa             <  1:610250/69‑1 (MQ=255)
                                                         tGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAACTGCTATAGAAATAATTACACAATACGGt         >  1:1159569/1‑71 (MQ=255)
                                                            ttttTAAATTATTCGATGTATACAAGCCTATATAGCGAACTGCtata                              <  1:446039/47‑1 (MQ=255)
                                                            ttttTAAATTATTCGATGTATACAAGCCTATATAGCGAACTGCtata                              <  1:1712063/47‑1 (MQ=255)
                                                                cAAATTATTCGATGTATACAAGCCTATATAGCGAACTGCTATAGAAATAATTACACAATACGGTTTGTTAc  >  1:1779971/1‑71 (MQ=255)
                                                                |                                                                      
TAATTACTCCTCACTTACACGTAATACTACTTTCGAGTGAAAATCTACCTATCTCTTTGATTTTCAAATTATTCGATGTATACAAGCCTATATAGCGAACTGCTATAGAAATAATTACACAATACGGTTTGTTAC  >  minE/597488‑597622

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: