Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,259,087 C→T 51.5% V613I (GTT→ATT)  torZ ← trimethylamine N‑oxide reductase system III, catalytic subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,259,0870CT51.5% ‑3.3 / 59.3 66V613I (GTT→ATT) torZtrimethylamine N‑oxide reductase system III, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base C (22/10);  new base T (15/19);  total (37/29)
Fisher's exact test for biased strand distribution p-value = 5.16e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

TCAATTTTGCCGCTTGGCGTACCCAGCGCATTTTTCACCGGATCGGCGCGGAAATCACCATAACGAACGTACTGTTCGTTCTTCTCGCTGTGGCGCATTTCGATCAGTTTATTTTGCTGCCAGAAGGCA  >  minE/1259020‑1259148
                                                                   |                                                             
tCAATTTTGCCGCTTGGCGTACCCAGCGCATTTTTCACCGGATCGGCGCGGAAATCACCATaacgaacg                                                              <  1:2247200/69‑1 (MQ=255)
tCAATTTTGCCGCTTGGCGTACCCAGCGCATTTTTCACCGGATCGGCGCGGAAATCACCATaacgaacg                                                              <  1:1073316/69‑1 (MQ=255)
tCAATTTTGCCGCTTGGCGTACCCAGCGCATTTTTCACCGGATCGGCGCGGAAATCACCATaacgaacg                                                              <  1:110280/69‑1 (MQ=255)
tCAATTTTGCCGCTTGGCGTACCCAGCGCATTTTTCACCGGATCGGCGCGGAAATCACCATaacgaacg                                                              <  1:898785/69‑1 (MQ=255)
           gCTTGGTGTACCCAGCGCATTTTTCACCGGATCGGCGCGGAAATCACCATAACGAATGTACTgttcgtt                                                   >  1:1743041/1‑69 (MQ=255)
           gCTTGGTGTACCCAGCGCATTTTTCACCGGATCGGCGCGGAAATCACCATAACGAATGTACTgttcgtt                                                   >  1:2540112/1‑69 (MQ=255)
           gCTTGGTGTACCCAGCGCATTTTTCACCGGATCGGCGCGGAAATCACCATAACGAATGTACTgttcgtt                                                   >  1:1525375/1‑69 (MQ=255)
           gCTTGGTGTACCCAGCGCATTTTTCACCGGATCGGCGCGGAAATCACCATAACGAATGTACTgttcgtt                                                   >  1:39912/1‑69 (MQ=255)
           gCTTGGTGTACCCAGCGCATTTTTCACCGGATCGGCGCGGAAATCACCATAACGAATGTACTgttcgtt                                                   >  1:2540795/1‑69 (MQ=255)
           gCTTGGTGTACCCAGCGCATTTTTCACCGGATCGGCGCGGAAATCACCATAACGAATGTACTgttcgtt                                                   >  1:2785550/1‑69 (MQ=255)
           gCTTGGTGTACCCAGCGCATTTTTCACCGGATCGGCGCGGAAATCACCATAACGAATGTACTgttcgtt                                                   >  1:2561814/1‑69 (MQ=255)
           gCTTGGTGTACCCAGCGCATTTTTCACCGGATCGGCGCGGAAATCACCATAACGAATGTACTgttcgtt                                                   >  1:239270/1‑69 (MQ=255)
           gCTTGGTGTACCCAGCGCATTTTTCACCGGATCGGCGCGGAAATCACCATAACGAATGTACTgttcgtt                                                   >  1:2021248/1‑69 (MQ=255)
           gCTTGGTGTACCCAGCGCATTTTTCACCGGATCGGCGCGGAAATCACCATAACGAATGTACTgttcgtt                                                   >  1:1197792/1‑69 (MQ=255)
           gCTTGGTGTACCCAGCGCATTTTTCACCGGATCGGCGCGGAAATCACCATAACGAATGTACTgttcgtt                                                   >  1:2837589/1‑69 (MQ=255)
           gCTTGGTGTACCCAGCGCATTTTTCACCGGATCGGCGCGGAAATCACCATAACGAATGTACTgttcgtt                                                   >  1:816802/1‑69 (MQ=255)
           gCTTGGTGTACCCAGCGCATTTTTCACCGGATCGGCGCGGAAATCACCATAACGAATGTACTgttcgtt                                                   >  1:2566909/1‑69 (MQ=255)
           gCTTGGTGTACCCAGCGCATTTTTCACCGGATCGGCGCGGAAATCACCATAACGAATGTACTgttcgtt                                                   >  1:2694245/1‑69 (MQ=255)
           gCTTGGTGTACCCAGCGCATTTTTCACCGGATCGGCGCGGAAATCACCATAACGAATGGACTgttcgtt                                                   >  1:721456/1‑69 (MQ=255)
            cTTGGCGTACCCAGCGCATTTTTCACCGGATCGGCGCGGAAATCACCATAACGAACGTACTgttcgttc                                                  <  1:2630839/69‑1 (MQ=255)
            cTTGGCGTACCCAGCGCATTTTTCACCGGATCGGCGCGGAAATCACCATAACGAACGTACTgttcgttc                                                  <  1:1592379/69‑1 (MQ=255)
            cTTGGCGTACCCAGCGCATTTTTCACCGGATCGGCGCGGAAATCACCATAACGAACGTACTgttcgttc                                                  <  1:1261402/69‑1 (MQ=255)
                                tttCACCGGATCGGCGCGGAAATCACCATAACGAACGTACTGTTCGTTCTTCTCGCTGTGGCGCAttt                               >  1:1163259/1‑68 (MQ=255)
                                tttCACCGGATCGGCGCGGAAATCACCATAACGAACGTACTGTTCGTTCTTCTCGCTGTGGCGCATTTc                              >  1:2707620/1‑69 (MQ=255)
                                tttCACCGGATCGGCGCGGAAATCACCATAACGAACGTACTGTTCGTTCTTCTCGCTGTGGCGCATTTc                              >  1:779789/1‑69 (MQ=255)
                                tttCACCGGATCGGCGCGGAAATCACCATAACGAACGTACTGTTCGTTCTTCTCGCTGTGGCGCATTTc                              >  1:1413419/1‑69 (MQ=255)
                                tttCACCGGATCGGCGCGGAAATCACCATAACGAACGTACTGTTCGTTCTTCTCGCTGTGGCGCATTTc                              >  1:2395457/1‑69 (MQ=255)
                                tttCACCGGATCGGCGCGGAAATCACCATAACGAACGTACTGTTCGTTCTTCTCGCTGTGGCGCATTTc                              >  1:12444/1‑69 (MQ=255)
                                tttCACCGGATCGGCGCGGAAATCACCATAACGAACGTACTGTTCGTTCTTCTCGCTGTGGCGCATTTc                              >  1:2040808/1‑69 (MQ=255)
                                  tCACCGGATCGGCGCGGAAATCACCATAACGAACGTACTGTTCGTTCTTCTCGCTGTGGCGCATTTCg                             <  1:2612020/68‑1 (MQ=255)
                                        gATCGGCGCGGAAATCACCATAACGAACGTACTGTTCGTTCTTCTCGCTGTGGCGCATTTCGATCAGtt                      <  1:2301606/69‑1 (MQ=255)
                                          tCGGCGCGGAAATCACCATAACGAACGTACTGTTCGTTCTTCTCGCTGTGGCGCATTTCGATCAGTTTa                    >  1:3185817/1‑69 (MQ=255)
                                          tCGGCGCGGAAATCACCATAACGAACGTACTGTTCGTTCTTCTCGCTGTGGCGCATTTCGATCAGTTTa                    >  1:96087/1‑69 (MQ=255)
                                          tCGGCGCGGAAATCACCATAACGAACGTACTGTTCGTTCTTCTCGCTGTGGCGCATTTCGATCAGTTTa                    >  1:955412/1‑69 (MQ=255)
                                          tCGGCGCGGAAATCACCATAACGAACGTACTGTTCGTTCTTCTCGCTGTGGCGCATTTCGATCAGTTTa                    >  1:1069335/1‑69 (MQ=255)
                                          tCGGCGCGGAAATCACCATAACGAACGTACTGTTCGTTCTTCTCGCTGTGGCGCATTTCGATCAGTTTa                    >  1:1344535/1‑69 (MQ=255)
                                          tCGGCGCGGAAATCACCATAACGAACGTACTGTTCGTTCTTCTCGCTGTGGCGCATTTCGATCAGTTTa                    <  1:1424469/69‑1 (MQ=255)
                                          tCGGCGCGGAAATCACCATAACGAACGTACTGTTCGTTCTTCTCGCTGTGGCGCATTTCGATCAGTTTa                    <  1:1762605/69‑1 (MQ=255)
                                          tCGGCGCGGAAATCACCATAACGAACGTACTGTTCGTTCTTCTCGCTGTGGCGCATTTCGATCAGTTTa                    >  1:2539304/1‑69 (MQ=255)
                                          tCGGCGCGGAAATCACCATAACGAACGTACTGTTCGTTCTTCTCGCTGTGGCGCATTTCGATCAGTTTa                    <  1:1975907/69‑1 (MQ=255)
                                          tCGGCGCGGAAATCACCATAACGAACGTACTGTTCGTTCTTCTCGCTGTGGCGCATTTCGATCAGTTTa                    <  1:2460424/69‑1 (MQ=255)
                                          tCGGCGCGGAAATCACCATAACGAACGTACTGTTCGTTCTTCTCGCTGTGGCGCATTTCGATCAGTTTa                    >  1:2044770/1‑69 (MQ=255)
                                          tCGGCGCGGAAATCACCATAACGAACGTACTGTTCGTTCTTCTCGCTGTGGCGCATTTCGATCAGTTTa                    >  1:2897839/1‑69 (MQ=255)
                                          tCGGCGCGGAAATCACCATAACGAACGTACTGTTCGTTCTTCTCGCTGTGGCGCATTTCGATCAGTTTa                    >  1:2962352/1‑69 (MQ=255)
                                            ggCGCGGAAATCACCATAACGAACGTACTGTTCGTTCTTCTCGCTGTGGCGCATTTCGATCAGTTTAtt                  >  1:3258766/1‑69 (MQ=255)
                                            ggCGCGGAAATCACCATAACGAACGTACTGTTCGTTCTTCTCGCTGTGGCGCATTTCGATCAGTTTAtt                  >  1:2582671/1‑69 (MQ=255)
                                                   aaaTCACCATAACGAATGTACTGTTCGTTCTTCTCGCTgcggcg                                    <  1:369634/44‑1 (MQ=255)
                                                   aaaTCACCATAACGAATGTACTGTTCGTTCTTCTCGCTgcggcg                                    <  1:2784943/44‑1 (MQ=255)
                                                   aaaTCACCATAACGAATGTACTGTTCGTTCTTCTCGCTgcggcg                                    <  1:70588/44‑1 (MQ=255)
                                                   aaaTCACCATAACGAATGTACTGTTCGTTCTTCTCGCTgcggcg                                    <  1:341818/44‑1 (MQ=255)
                                                   aaaTCACCATAACGAATGTACTGTTCGTTCTTCTCGCTgcggcg                                    <  1:3141247/44‑1 (MQ=255)
                                                   aaaTCACCATAACGAATGTACTGTTCGTTCTTCTCGCTgcggcg                                    <  1:3210926/44‑1 (MQ=255)
                                                   aaaTCACCATAACGAATGTACTGTTCGTTCTTCTCGCTgcggcg                                    <  1:1047101/44‑1 (MQ=255)
                                                   aaaTCACCATAACGAATGTACTGTTCGTTCTTCTCGCTgcggcg                                    <  1:2523774/44‑1 (MQ=255)
                                                   aaaTCACCATAACGAATGTACTGTTCGTTCTTCTCGCTgcggcg                                    <  1:146820/44‑1 (MQ=255)
                                                   aaaTCACCATAACGAATGTACTGTTCGTTCTTCTCGCTgcggcg                                    <  1:2466853/44‑1 (MQ=255)
                                                   aaaTCACCATAACGAATGTACTGTTCGTTCTTCTCGCTgcggcg                                    <  1:2000979/44‑1 (MQ=255)
                                                   aaaTCACCATAACGAATGTACTGTTCGTTCTTCTCGCTgcggcg                                    <  1:2363116/44‑1 (MQ=255)
                                                   aaaTCACCATAACGAATGTACTGTTCGTTCTTCTCGCTgcggcg                                    <  1:2351227/44‑1 (MQ=255)
                                                   aaaTCACCATAACGAATGTACTGTTCGTTCTTCTCGCTgcggcg                                    <  1:229388/44‑1 (MQ=255)
                                                   aaaTCACCATAACGAATGTACTGTTCGTTCTTCTCGCTgcggcg                                    <  1:1105687/44‑1 (MQ=255)
                                                    aaTCACCATAACGAATGTACTGTTCGTTCTTCTCGCTgcggcg                                    <  1:2907213/43‑1 (MQ=255)
                                                    aaTCACCATAACGAATGTACTGTTCGTTCTTCTCGCTgcggcg                                    <  1:2135067/43‑1 (MQ=255)
                                                     aTCACCATAACGAATGTACTGTTCGTTCTTCTCGCTgcggcg                                    <  1:1571136/42‑1 (MQ=255)
                                                     aTCACCATAACGAATGTACTGTTCGTTCTTCTCGCTgcggcg                                    <  1:2317191/42‑1 (MQ=255)
                                                           aTAACGAACGTACTGTTCGTTCTTCTCGCTGTGGCGCATTTCGATCAGTTTATTTTGCTGCCAGAAGGc   >  1:3051617/1‑69 (MQ=255)
                                                           aTAACGAACGTACTGTTCGTTCTTCTCGCTGTGGCGCATTTCGATCAGTTTATTTTGCTGCCAGAAGGc   >  1:3225519/1‑69 (MQ=255)
                                                           aTAACGAACGTACTGTTCGTTCTTCTCGCTGTGGCGCATTTCGATCAGTTTATTTTGCTGCCAGAAGGc   >  1:2883303/1‑69 (MQ=255)
                                                           aTAACGAACGTACTGTTCGTTCTTCTCGCTGTGGCGCATTTCGATCAGTTTATTTTGCTGCCAGAAGGc   >  1:2517065/1‑69 (MQ=255)
                                                            tAACGAACGTACTGTTCGTTCTTCTCGCTGTGGCGCATTTCGATCAGTTTATTTTGCTGCCAGAAGGCa  <  1:1106456/69‑1 (MQ=255)
                                                                   cgTACTGTTCGTTCTTCTCGCTGTGGCGCATTTCGATCAGTTTATTTtgctgc           <  1:1018446/53‑1 (MQ=255)
                                                                   cgTACTGTTCGTTCTTCTCGCTGTGGCGCATTTCGATCAGTTTATTTtgctgc           <  1:872272/53‑1 (MQ=255)
                                                                   cgTACTGTTCGTTCTTCTCGCTGTGGCGCATTTCGATCAGTTTATTTtgctgc           <  1:1067751/53‑1 (MQ=255)
                                                                   cgTACTGTTCGTTCTTCTCGCTGTGGCGCATTTCGATCAGTTTATTTtgctgc           <  1:76336/53‑1 (MQ=255)
                                                                   cgTACTGTTCGTTCTTCTCGCTGTGGCGCATTTCGATCAGTTTATTTtgctgc           <  1:2603409/53‑1 (MQ=255)
                                                                   cgTACTGTTCGTTCTTCTCGCTGTGGCGCATTTCGATCAGTTTATTTtgctgc           <  1:1886305/53‑1 (MQ=255)
                                                                   cgTACTGTTCGTTCTTCTCGCTGTGGCGCATTTCGATCAGTTTATTTtgctgc           <  1:1896246/53‑1 (MQ=255)
                                                                   cgTACTGTTCGTTCTTCTCGCTGTGGCGCATTTCGATCAGTTTATTTtgctgc           <  1:2000824/53‑1 (MQ=255)
                                                                   cgTACTGTTCGTTCTTCTCGCTGTGGCGCATTTCGATCAGTTTATTTtgctgc           <  1:355038/53‑1 (MQ=255)
                                                                   cgTACTGTTCGTTCTTCTCGCTGTGGCGCATTTCGATCAGTTTATTTtgctgc           <  1:336176/53‑1 (MQ=255)
                                                                   cgTACTGTTCGTTCTTCTCGCTGTGGCGCATTTCGATCAGTTTATTTtgctgc           <  1:2111448/53‑1 (MQ=255)
                                                                   cgTACTGTTCGTTCTTCTCGCTGTGGCGCATTTCGATCAGTTTATTTtgctgc           <  1:2152040/53‑1 (MQ=255)
                                                                   cgTACTGTTCGTTCTTCTCGCTGTGGCGCATTTCGATCAGTTTATTTtgctgc           <  1:2818982/53‑1 (MQ=255)
                                                                   cgTACTGTTCGTTCTTCTCGCTGTGGCGCATTTCGATCAGTTTATTTtgctgc           <  1:2568762/53‑1 (MQ=255)
                                                                   |                                                             
TCAATTTTGCCGCTTGGCGTACCCAGCGCATTTTTCACCGGATCGGCGCGGAAATCACCATAACGAACGTACTGTTCGTTCTTCTCGCTGTGGCGCATTTCGATCAGTTTATTTTGCTGCCAGAAGGCA  >  minE/1259020‑1259148

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: