Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,903,226 T→C 53.2% S273S (TCA→TCG nrdD ← anaerobic ribonucleoside‑triphosphate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,903,2260TC53.2% 9.0 / 65.6 79S273S (TCA→TCGnrdDanaerobic ribonucleoside‑triphosphate reductase
Reads supporting (aligned to +/- strand):  ref base T (26/11);  new base C (28/14);  total (54/25)
Fisher's exact test for biased strand distribution p-value = 8.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.25e-01

CTGACCGTTGGCGGTGTGCAGGGTGTTTACTTCGTACTCCAGTGACTGGAAGGCATCGTAGCACTCTTTGATGGTTCGAGAGTTAGCGTAGCCTTCGGCGTCCGGGA  >  minE/2903184‑2903290
                                          |                                                                
cTGACCGTTGGCGGTGTGCAGGGTGTTTACTTCGTACTCCAGTGACTGGAAGGCATCGTAGCACTCttt                                        <  1:1025439/69‑1 (MQ=255)
cTGACCGTTGGCGGTGTGCAGGGTGTTTACTTCGTACTCCAGTGACTGGAAGGCATCGTAGCACTCttt                                        <  1:3278614/69‑1 (MQ=255)
cTGACCGTTGGCGGTGTGCAGGGTGTTTACTTCGTACTCCAGTGACTGGAAGGCATCGTAGCACTCttt                                        <  1:621343/69‑1 (MQ=255)
cTGACCGTTGGCGGTGTGCAGGGTGTTTACTTCGTACTCCAGTGACTGGAAGGCATCGTAGCACTCttt                                        <  1:684164/69‑1 (MQ=255)
cTGACCGTTGGCGGTGTGCAGGGTGTTTACTTCGTACTCCAGTGACTGGAAGGCATCGTAGCACTCttt                                        <  1:1573063/69‑1 (MQ=255)
cTGACCGTTGGCGGTGTGCAGGGTGTTTACTTCGTACTCCAGTGACTGGAAGGCATCGTAGCACTCttt                                        <  1:1467730/69‑1 (MQ=255)
 tGACCGTTGGCGGTGTGCAGGGTGTTTACTTCGTACTCCAGTGACTGGAAGGCATCGTAGCACTCTTTg                                       >  1:2551173/1‑69 (MQ=255)
 tGACCGTTGGCGGTGTGCAGGGTGTTTACTTCGTACTCCAGTGACTGGAAGGCATCGTAGCACTCTTTg                                       >  1:3294978/1‑69 (MQ=255)
 tGACCGTTGGCGGTGTGCAGGGTGTTTACTTCGTACTCCAGTGACTGGAAGGCATCGTAGCACTCTTTg                                       >  1:2768669/1‑69 (MQ=255)
 tGACCGTTGGCGGTGTGCAGGGTGTTTACTTCGTACTCCAGTGACTGGAAGGCATCGTAGCACTCTTTg                                       >  1:846557/1‑69 (MQ=255)
 tGACCGTTGGCGGTGTGCAGGGTGTTTACTTCGTACTCCAGTGACTGGAAGGCATCGTAGCACTCTTTg                                       >  1:81764/1‑69 (MQ=255)
 tGACCGTTGGCGGTGTGCAGGGTGTTTACTTCGTACTCCAGTGACTGGAAGGCATCGTAGCACTCTTTg                                       >  1:891439/1‑69 (MQ=255)
   aCCGTTGGCGGTGTGCAGGGTGTTTACTTCGTACTCCAGTGACTGGAAGGCATCGTAGCACTCTTTGAt                                     >  1:11794/1‑69 (MQ=255)
   aCCGTTGGCGGTGTGCAGGGTGTTTACTTCGTACTCCAGTGACTGGAAGGCATCGTAGCACTCTTTGAt                                     >  1:1623591/1‑69 (MQ=255)
   aCCGTTGGCGGTGTGCAGGGTGTTTACTTCGTACTCCAGTGACTGGAAGGCATCGTAGCACTCTTTGAt                                     >  1:2952362/1‑69 (MQ=255)
   aCCGTTGGCGGTGTGCAGGGTGTTTACTTCGTACTCCAGTGACTGGAAGGCATCGTAGCACTCTTTGAt                                     >  1:3225997/1‑69 (MQ=255)
     cGTTGGCGGTGTGCAGGGTGTTTACTTCGTACTCCAGTGACTGGAAGGCATCGTAGCACTCTTTGAt                                     <  1:3081422/67‑1 (MQ=255)
                      gtgtTTACTTCGTACTCCAGTGACTGGAAGGCATCGTAGCACTCTTTGATGGTTCGAGAGTTAGCGTAg                  >  1:2988975/1‑69 (MQ=255)
                       tgtTTACTTCGTACTCCAGTGACTGGAAGGCATCGTAGCACTCTTTGATGGTTCGAGAGTTAGCGTAGc                 >  1:1711875/1‑69 (MQ=255)
                        gtTTACTTCGTACTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGTTACGAGAGTTAGCGTAGcc                >  1:826740/1‑69 (MQ=255)
                        gtTTACTTCGTACTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGcc                >  1:2944615/1‑69 (MQ=255)
                        gtTTACTTCGTACTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGcc                >  1:3013902/1‑69 (MQ=255)
                        gtTTACTTCGTACTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGcc                >  1:2887426/1‑69 (MQ=255)
                        gtTTACTTCGTACTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGcc                >  1:2855595/1‑69 (MQ=255)
                        gtTTACTTCGTACTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGcc                >  1:3287095/1‑69 (MQ=255)
                        gtTTACTTCGTACTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGcc                >  1:266952/1‑69 (MQ=255)
                        gtTTACTTCGTACTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGcc                >  1:3305356/1‑69 (MQ=255)
                        gtTTACTTCGTACTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGcc                >  1:1219124/1‑69 (MQ=255)
                        gtTTACTTCGTACTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGcc                >  1:3383766/1‑69 (MQ=255)
                        gtTTACTTCGTACTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGcc                >  1:2428882/1‑69 (MQ=255)
                        gtTTACTTCGTACTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGcc                >  1:2230939/1‑69 (MQ=255)
                        gtTTACTTCGTACTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGcc                >  1:852117/1‑69 (MQ=255)
                        gtTTACTTCGTACTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGcc                >  1:1977072/1‑69 (MQ=255)
                        gtTTACTTCGTACTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGcc                >  1:422011/1‑69 (MQ=255)
                        gtTTACTTCGTACTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGcc                >  1:473954/1‑69 (MQ=255)
                        gtTTACTTCGTACTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGcc                >  1:851720/1‑69 (MQ=255)
                        gtTTACTTCGTACTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGcc                >  1:575154/1‑69 (MQ=255)
                        gtTTACTTCGTACTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGcc                >  1:924106/1‑69 (MQ=255)
                        gtTTACTTCGTACTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGcc                >  1:701712/1‑69 (MQ=255)
                        gtTTACTTCGTACTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGcc                >  1:2402322/1‑69 (MQ=255)
                        gtTTACTTCGTACTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGcc                >  1:791284/1‑69 (MQ=255)
                        gtTTACTTCGTACTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGcc                >  1:1470325/1‑69 (MQ=255)
                        gtTTACTTCGTACTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGcc                >  1:16131/1‑69 (MQ=255)
                        gtTTACTTCGTACTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGcc                >  1:11374/1‑69 (MQ=255)
                        gtTTACTTCGTACTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGc                 >  1:1282757/1‑68 (MQ=255)
                        gtTTACTTCGTACTACAGCGACTGGAAGGCATCGTAGCACTCTTTGATCGTACGAGAGTTAGCGTAGcc                >  1:723793/1‑69 (MQ=255)
                        gtTTACTTCGAACTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGcc                >  1:3147553/1‑69 (MQ=255)
                           tACTTCGTACTCCAGTGACTGGAAGGCATCGTAGCACTCttt                                        <  1:1028251/42‑1 (MQ=255)
                           tACTTCGTACTCCAGTGACTGGAAGGCATCGTAGCACTCttt                                        <  1:2576456/42‑1 (MQ=255)
                           tACTTCGTACTCCAGTGACTGGAAGGCATCGTAGCACTCttt                                        <  1:2499195/42‑1 (MQ=255)
                           tACTTCGTACTCCAGTGACTGGAAGGCATCGTAGCACTCTTTGATGGTTc                                >  1:2921489/1‑50 (MQ=255)
                           tACTTCGTACTCCAGTGACTGGAAGGCATCGTAGCACTCTTTGATGGTTCGAGAGTTAGCGTAGCCTTc             >  1:999483/1‑69 (MQ=255)
                                   aCTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGCCTTCGGCGTcc      <  1:597096/68‑1 (MQ=255)
                                   aCTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGCCTTCGGCGTcc      <  1:1349412/68‑1 (MQ=255)
                                   aCTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGCCTTCGGCGTcc      <  1:1659195/68‑1 (MQ=255)
                                   aCTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGCCTTCGGCGTcc      <  1:126753/68‑1 (MQ=255)
                                   aCTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGCCTTCGGCGTcc      <  1:1708827/68‑1 (MQ=255)
                                   aCTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGCCTTCGGCGTcc      <  1:651332/68‑1 (MQ=255)
                                   aCTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGCCTTCGGCGTcc      <  1:12670/68‑1 (MQ=255)
                                   aCTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGCCTTCGGCGTcc      <  1:1862862/68‑1 (MQ=255)
                                   aCTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGCCTTCGGCGTcc      <  1:1161713/68‑1 (MQ=255)
                                   aCTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGCCTTCGGCGTcc      <  1:483424/68‑1 (MQ=255)
                                   aCTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGCCTTCGGCGTcc      <  1:1038472/68‑1 (MQ=255)
                                   aCTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGCCTTCGGCGTcc      <  1:1937609/68‑1 (MQ=255)
                                   aCTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGCCTTCGGCGTcc      <  1:953349/68‑1 (MQ=255)
                                   aCTCCAGCGACTGGAAGGCATCGTAGCACTCTTTGATGGTACGAGAGTTAGCGTAGCCTTCGGCGTcc      <  1:2576157/68‑1 (MQ=255)
                                    cTCCAGTGACTTGAAGGCATCGTAGCACTCTTTGATGGTTCGAGAGTTAGCGTAGCCTTCGGCGTCCgg    >  1:2815910/1‑69 (MQ=255)
                                    cTCCAGTGACTGGAAGGCATCGTAGCACTCTTTGATGGTTCGAGAGTTAGCGTAGCCTTCGGCGTCCgg    >  1:863879/1‑69 (MQ=255)
                                    cTCCAGTGACTGGAAGGCATCGTAGCACTCTTTGATGGTTCGAGAGTTAGCGTAGCCTTCGGCGTCCgg    >  1:1821228/1‑69 (MQ=255)
                                    cTCCAGTGACTGGAAGGCATCGTAGCACTCTTTGATGGTTCGAGAGTTAGCGTAGCCTTCGGCGTCCgg    >  1:1028205/1‑69 (MQ=255)
                                    cTCCAGTGACTGGAAGGCATCGTAGCACTCTTTGATGGTTCGAGAGTTAGCGTAGCCTTCGGCGTCCgg    >  1:1376271/1‑69 (MQ=255)
                                    cTCCAGTGACTGGAAGGCATCGTAGCACTCTTTGATGGTTCGAGAGTTAGCGTAGCCTTCGGCGTCCgg    >  1:1434078/1‑69 (MQ=255)
                                    cTCCAGTGACTGGAAGGCATCGTAGCACTCTTTGATGGTTCGAGAGTTAGCGTAGCCTTCGGCGTCCgg    >  1:1683010/1‑69 (MQ=255)
                                    cTCCAGTGACTGGAAGGCATCGTAGCACTCTTTGATGGTTCGAGAGTTAGCGTAGCCTTCGGCGTCCgg    >  1:641670/1‑69 (MQ=255)
                                    cTCCAGTGACTGGAAGGCATCGTAGCACTCTTTGATGGTTCGAGAGTTAGCGTAGCCTTCGGCGTCCgg    >  1:1873364/1‑69 (MQ=255)
                                    cTCCAGTGACTGGAAGGCATCGTAGCACTCTTTGATGGTTCGAGAGTTAGCGTAGCCTTCGGCGTCCgg    >  1:1876762/1‑69 (MQ=255)
                                    cTCCAGTGACTGGAAGGCATCGTAGCACTCTTTGATGGTTCGAGAGTTAGCGTAGCCTTCGGCGTCCgg    >  1:2545112/1‑69 (MQ=255)
                                    cTCCAGTGACTGGAAGGCATCGTAGCACTCTTTGATGGTTCGAGAGTTAGCGTAGCCTTCGGCGTCCgg    >  1:2734640/1‑69 (MQ=255)
                                      ccAGTGACTGGAAGGCATCGTAGCACTCTTTGATGGTTCGAGAGTTAGCGTAGCCTTCGGCGTCCGGGa  <  1:494129/69‑1 (MQ=255)
                                          |                                                                
CTGACCGTTGGCGGTGTGCAGGGTGTTTACTTCGTACTCCAGTGACTGGAAGGCATCGTAGCACTCTTTGATGGTTCGAGAGTTAGCGTAGCCTTCGGCGTCCGGGA  >  minE/2903184‑2903290

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: