Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,820,972 A→G 14.0% A93A (GCA→GCG norV → flavorubredoxin oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,820,9720AG14.0% 160.7 / 12.5 92A93A (GCA→GCGnorVflavorubredoxin oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (50/29);  new base G (7/6);  total (57/35)
Fisher's exact test for biased strand distribution p-value = 5.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

TATCGATTACATCGTGATTAACCATGCAGAAGAGGACCACGCTGGGGCGCTGACCGAACTGATGGCACAAATTCCCGATACGCCGATCTACTGTACAGCCAACGCTATCGACTCGATAAATG  >  minE/1820906‑1821027
                                                                  |                                                       
tatCGATTACATCGTGATTAACCATGCAGAAGAGGACCACGCTGGGGCGCTGACCGAACTGATGGcaca                                                       <  1:1080810/69‑1 (MQ=255)
tatCGATTACATCGTGATTAACCATGCAGAAGAGGACCACGCTGGGGCGCTGACCGAACTGATGGcaca                                                       <  1:277184/69‑1 (MQ=255)
tatCGATTACATCGTGATTAACCATGCAGAAGAGGACCACGCTGGGGCGCTGACCGAACTGATGGcaca                                                       <  1:2674512/69‑1 (MQ=255)
tatCGATTACATCGTGATTAACCATGCAGAAGAGGACCACGCTGGGGCGCTGACCGAACTGATGGcaca                                                       <  1:2214527/69‑1 (MQ=255)
 atCGATTACATCGTGATTAACCATGCAGAAGAGGACCACGCTGGGGCGCTGACCGAACTGATGGcaca                                                       <  1:1925645/68‑1 (MQ=255)
 atCGATTACATCGTGATTAACCATGCAGAAGAGGACCACGCTGGGGCGCTGACCGAACTGATGGcaca                                                       <  1:1915022/68‑1 (MQ=255)
   cGATTACATCGTGATTAACCATGCAGAAGAGGACCACGCTGTGGCGCTGACCGAACTGATGGCACAAAt                                                    >  1:1186340/1‑69 (MQ=255)
   cGATTACATCGTGATTAACCATGCAGAAGAGGACCACGCTGGGGCGCTGACCGAACTGATGGCACAAAt                                                    >  1:2077958/1‑69 (MQ=255)
   cGATTACATCGTGATTAACCATGCAGAAGAGGACCACGCTGGGGCGCTGACCGAACTGATGGCACAAAt                                                    >  1:1768653/1‑69 (MQ=255)
    gATTACATCGTGATTACCCATGCAGAAGAGGACCACGCTGGGGCGCTGACCGAACTGATGGCACAAAtt                                                   >  1:532299/1‑69 (MQ=255)
    gATTACATCGTGATTAACCATGCAGAAGAGGACCACGCTGGGGCGCTGACCGAACTGATGGCACAAAtt                                                   >  1:1959833/1‑69 (MQ=255)
    gATTACATCGTGATTAACCATGCAGAAGAGGACCACGCTGGGGCGCTGACCGAACTGATGGCACAAAtt                                                   >  1:1892150/1‑69 (MQ=255)
    gATTACATCGTGATTAACCATGCAGAAGAGGACCACGCTGGGGCGCTGACCGAACTGATGGCACAAAtt                                                   >  1:3034729/1‑69 (MQ=255)
       tACATCGTGATTAACCATGCAGAAGAGGACCACGCTGGGGCGCTGACCGAACTGATGGCACAAATTccc                                                >  1:3343386/1‑69 (MQ=255)
       tACATCGTGATTAACCATGCAGAAGAGGACCACGCTGGGGCGCTGACCGAACTGATGGCACAAATTccc                                                >  1:666652/1‑69 (MQ=255)
          aTCGTGATTAACCATGCAGAAGAGGACCACGCTGGGGCGCTGACCGAACTGATGGCACAAATTCCCGa                                              <  1:940991/68‑1 (MQ=255)
                aTTAACCATGCTGAAGAGGACCACGCCGGGGCGCTGACCGAACTGATGGCGCAAATTCCCGATACGcc                                        >  1:636028/1‑68 (MQ=255)
                aTTAACCATGCTGAAGAGGACCACGCCGGGGCGCTGACCGAACTGATGGCGCAAATTCCCGATACGCCg                                       >  1:2388422/1‑69 (MQ=255)
                aTTAACCATGCTGAAGAGGACCACGCCGGGGCGCTGACCGAACTGATGGCGCAAATTCCCGATACGCCg                                       >  1:1607372/1‑69 (MQ=255)
                aTTAACCATGCTGAAGAGGACCACGCCGGGGCGCTGACCGAACTGATGGCGCAAATTCCCGATACGCCg                                       >  1:3100442/1‑69 (MQ=255)
                aTTAACCATGCTGAAGAGGACCACGCCGGGGCGCTGACCGAACTGATGGCGCAAATTCCCGATACGCCg                                       >  1:2819752/1‑69 (MQ=255)
                aTTAACCATGCTGAAGAGGACCACGCCGGGGCGCTGACCGAACTGATGGCGCAAATTCCCGATACGCCg                                       >  1:3234135/1‑69 (MQ=255)
                aTTAACCATGCAGAAGAGGACCACGCTGGGGCGCTGACCGAACTGATGGCACAAATTCCCGATACGCCg                                       <  1:2235644/69‑1 (MQ=255)
                aTTAACCATGCAGAAGAGGACCACGCTGGGGCGCTGACCGAACTGATGGCACAAATTCCCGATACGCCg                                       <  1:595418/69‑1 (MQ=255)
                aTTAACCATGCAGAAGAGGACCACGCTGGGGCGCTGACCGAACTGATGGCACAAATTCCCGATACGCCg                                       <  1:1332025/69‑1 (MQ=255)
                aTTAACCATGCAGAAGAGGACCACGCTGGGGCGCTGACCGAACTGATGGCACAAATTCCCGATACGCCg                                       <  1:1359635/69‑1 (MQ=255)
                 ttAACCATGCTGAAGAGGACCACGCCGGGGCGCTGACCGAACTGATGGCGCAAATTCCCGATACGCCGa                                      >  1:489294/1‑69 (MQ=255)
                  tAACCATGCAGAAGAGGACCACGCTGGGGCGCTGACCGAACTGATGGCACAAATTCCCGATa                                            <  1:612273/62‑1 (MQ=255)
                  tAACCATGCAGAAGAGGACCACGCTGGGGCGCTGACCGAACTGATGGCACAAATTCCCGATa                                            <  1:2471268/62‑1 (MQ=255)
                  tAACCATGCAGAAGAGGACCACGCTGGGGCGCTGACCGAACTGATGGCACAAATTCCCGATa                                            >  1:2707402/1‑62 (MQ=255)
                  tAACCATGCAGAAGAGGACCACGCTGGGGCGCTGACCGAACTGATGGCACAAATTCCCGATa                                            <  1:2269968/62‑1 (MQ=255)
                  tAACCATGCAGAAGAGGACCACGCTGGGGCGCTGACCGAACTGATGGCACAAATTCCCGATa                                            <  1:3286399/62‑1 (MQ=255)
                  tAACCATGCAGAAGAGGACCACGCTGGGGCGCTGACCGAACTGATGGCACAAATTCCCGATa                                            >  1:1964582/1‑62 (MQ=255)
                  tAACCATGCAGAAGAGGACCACGCTGGGGCGCTGACCGAACTGATGGCACAAATTCCCGATa                                            <  1:3351666/62‑1 (MQ=255)
                  tAACCATGCAGAAGAGGACCACGCTGGGGCGCTGACCGAACTGATGGCACAAATTCCCGATa                                            >  1:1762628/1‑62 (MQ=255)
                  tAACCATGCAGAAGAGGACCACGCTGGGGCGCTGACCGAACTGATGGCACAAATTCCCGATa                                            >  1:72360/1‑62 (MQ=255)
                  tAACCATGCAGAAGAGGACCACGCTGGGGCGCTGACCGAACTGATGGCACAAATTCCCGATa                                            <  1:931412/62‑1 (MQ=255)
                  tAACCATGCAGAAGAGGACCACGCTGGGGCGCTGACCGAACTGATGGCACAAATTCCCGATa                                            >  1:949515/1‑62 (MQ=255)
                  tAACCATGCAGAAGAGGACCACGCTGGGGCGCTGACCGAACTGATGGCACAAATTCCCGATa                                            >  1:1382730/1‑62 (MQ=255)
                  tAACCATGCAGAAGAGGACCACGCTGGGGCGCTGACCGAACTGATGGCACAAATTCCCGATa                                            >  1:1373531/1‑62 (MQ=255)
                  tAACCATGCAGAAGAGGACCACGCTGGGGCGCTGACCGAACTGATGGCACAAATTCCCGATa                                            >  1:324613/1‑62 (MQ=255)
                  tAACCATGCAGAAGAGGACCACGCTGGGGCGCTGACCGAACTGATGGCACAAATTCCCGATa                                            <  1:3279601/62‑1 (MQ=255)
                  tAACCATGCAGAAGAGGACCACGCTGGGGCGCTGACCGAACTGATGGCACAAATTCCCGATa                                            >  1:1179745/1‑62 (MQ=255)
                          cAGAAGAGGACCACGCTGGGGCGCTGACCGAACTGATGGCACAAATTCCCGATACGCCg                                       <  1:2604154/59‑1 (MQ=255)
                          cAGAAGAGGACCACGCTGGGGCGCTGACCGAACTGATGGCACAAATTCCCGATACGCCGATCTACTg                               >  1:451224/1‑67 (MQ=255)
                          cAGAAGAGGACCACGCTGGGGCGCTGACCGAACTGATGGCACAAATTCCCGATACGCCGATCTACTg                               >  1:3302793/1‑67 (MQ=255)
                          cAGAAGAGGACCACGCTGGGGCGCTGACCGAACTGATGGCACAAATTCCCGATACGCCGATCTACTg                               >  1:46890/1‑67 (MQ=255)
                              agagGACCACGCTGGGGCGCTGACCGAACTGATGGCACAAATTCCCg                                               >  1:1136855/1‑47 (MQ=255)
                              agagGACCACGCTGGGGCGCTGACCGAACTGATGGCACAAATTCCCGATACGCCGATCTACTg                               >  1:1187129/1‑63 (MQ=255)
                              agagGACCACGCTGGGGCGCTGACCGAACTGATGGCACAAATTCCCGATACGCCGATCTACTg                               >  1:626641/1‑63 (MQ=255)
                              agagGACCACGCTGGGGCGCTGACCGAACTGATGGCACAAATTCCCGATACGCCGATCTACTg                               >  1:1114818/1‑63 (MQ=255)
                              agagGACCACGCTGGGGCGCTGACCGAACTGATGGCACAAATTCCCGATACGCCGATCTACTg                               >  1:797001/1‑63 (MQ=255)
                              agagGACCACGCTGGGGCGCTGACCGAACTGATGGCACAAATTCCCGATACGCCGATCTACTg                               >  1:2845913/1‑63 (MQ=255)
                              agagGACCACGCTGGGGCGCTGACCGAACTGATGGCACAAATTCCCGATACGCCGATCTACTg                               >  1:1785983/1‑63 (MQ=255)
                              agagGACCACGCTGGGGCGCTGACCGAACTGATGGCACAAATTCCCGATACGCCGATCTACTg                               >  1:1880233/1‑63 (MQ=255)
                              agagGACCACGCTGGGGCGCTGACCGAACTGATGGCACAAATTCCCGATACGCCGATCTACTg                               >  1:2966777/1‑63 (MQ=255)
                              agagGACCACGCTGGGGCGCTGACCGAACTGATGGCACAAATTCCAGATACGCCGATCTACt                                >  1:2262643/1‑62 (MQ=255)
                                  gACCACGCTGGGGCGCTGACCGAACTGATGGCACAAATTCCCGATACGCCGATCTACTGTACAGcc                        >  1:1016753/1‑66 (MQ=255)
                                  gACCACGCTGGGGCGCTGACCGAACTGATGGCACAAATTCCCGATACGCCGATCTACTGTACAGCCAAc                     >  1:2874267/1‑69 (MQ=255)
                                  gACCACGCTGGGGCGCTGACCGAACTGATGGCACAAATTCCCGATACGCCGATCTACTGTACAGCCAAc                     >  1:1648730/1‑69 (MQ=255)
                                  gACCACGCTGGGGCGCTGACCGAACTGATGGCACAAATTCCCGATACGCCGATCTACTGTACAGCCAAc                     >  1:3356933/1‑69 (MQ=255)
                                  gACCACGCTGGGGCGCTGACCGAACTGATGGCACAAATTCCCGATACGCCGATCTACTGTACAGCCAAc                     >  1:1923783/1‑69 (MQ=255)
                                  gACCACGCTGGGGCGCTGACCGAACTGATGGCACAAATTCCCGATACGCCGATCTACTGTACAGCCAAc                     >  1:2753869/1‑69 (MQ=255)
                                         ccggggCGCTGACCGAACTGATGGCGCAAATTCCCGATACGCCGATCTACTGTACAGCCAACGCTATCg              <  1:2900319/67‑1 (MQ=255)
                                         ccggggCGCTGACCGAACTGATGGCGCAAATTCCCGATACGCCGATCTACTGTACAGCCAACGCTATCg              <  1:1304778/67‑1 (MQ=255)
                                         ccggggCGCTGACCGAACTGATGGCGCAAATTCCCGATACGCCGATCTACTGTACAGCCAACGCTATCg              <  1:952931/67‑1 (MQ=255)
                                         ccggggCGATGACCGAACTGATGGCGCAAATTCCCGATACGCCGATCTACTGTACAGCCAACGCTATCg              <  1:2495246/67‑1 (MQ=255)
                                          tGGGGCGCTGACCGAACTGATGGCACAAATTCCCGATACGCCGATCTACTGTACAGCCAACGCTATCGa             <  1:880709/69‑1 (MQ=255)
                                          tGGGGCGCTGACCGAACTGATGGCACAAATTCCCGATACGCCGATCTACTGTACAGCCAACGCTATCGa             <  1:587673/69‑1 (MQ=255)
                                          tGGGGCGCTGACCGAACTGATGGCACAAATTCCCGATACGCCGATCTACTGTACAGCCAACGCTATCGa             <  1:1154841/69‑1 (MQ=255)
                                          tGGGGCGCTGACCGAACTGATGGCACAAATTCCCGATACGCCGATCTACTGTACAGCCAACGCTATCGa             <  1:1376270/69‑1 (MQ=255)
                                          tGGGGCGCTGACCGAACTGATGGCACAAATTCCCGATACGCCGATCTACTGTACAGCCAACGCTATCGa             <  1:3017027/69‑1 (MQ=255)
                                          tGGGGCGCTGACCGAACTGATGGCACAAATTCCCGATACGCCGATCTACTGTACAGCCAACGCTATCGa             <  1:1407157/69‑1 (MQ=255)
                                          tGGGGCGCTGACCGAACTGATGGCACAAATTCCCGATACGCCGATCTACTGTACAGCCAACGCTATCGa             <  1:1479976/69‑1 (MQ=255)
                                          tGGGGCGCTGACCGAACTGATGGCACAAATTCCCGATACGCCGATCTACTGTACAGCCAACGCTATCGa             <  1:1608049/69‑1 (MQ=255)
                                          tGGGGCGCTGACCGAACTGATGGCACAAATTCCCGATACGCCGATCTACTGTACAGCCAACGCTATCGa             <  1:2637200/69‑1 (MQ=255)
                                          cggggCGCTGACCGAACTGATGGCGCAAATTCCCGATACGCCGATCTACTGTACAGCCAACGCTATCg              <  1:2581355/67‑1 (MQ=255)
                                           ggggCGCTGACCGAACTGATGGCACAAATTCCCGATACGCCGATCTACTGTACAGCCAACGCTATCGa             <  1:258211/68‑1 (MQ=255)
                                               cgcTGACCGAACTGATGGCACAAATTCCCGATACGCCGATCTACTGTACAGCCAACGCTATCGACTc          >  1:210680/1‑67 (MQ=255)
                                               cgcTGACCGAACTGATGGCACAAATTCCCGATACGCCGATCTACTGTACAGCCAACGCTATCGACTCg         >  1:2393920/1‑68 (MQ=255)
                                               cgcTGACCGAACTGATGGCACAAATTCCCGATACGCCGATCTACTGTACAGCCAACGCTATCGACTCg         >  1:445810/1‑68 (MQ=255)
                                               cgcTGACCGAACTGATGGCACAAATTCCAGATACGCCGATCTACTGTACAGCCAACGCTATCGACTCg         >  1:3096067/1‑68 (MQ=255)
                                                 cTGACCGAACTGATGGCACAAATTCCCGATACGCCGATCTACTGTACAGCCAACGCTATCGACTCGATa      <  1:1708988/69‑1 (MQ=255)
                                                 cTGACCGAACTGATGGCACAAATTCCCGATACGCCGATCTACTGTACAGCCAACGCTATCGACTCGATa      >  1:1683861/1‑69 (MQ=255)
                                                 cTGACCGAACTGATGGCACAAATTCCCGATACGCCGATCTACTGTACAGCCAACGCTATCGACTCGATa      <  1:1549375/69‑1 (MQ=255)
                                                 cTGACCGAACTGATGGCACAAATTCCCGATACGCCGATCTACTGTACAGCCAACGCTATCGACTCGATa      >  1:948575/1‑69 (MQ=255)
                                                 cTGACCGAACTGATGGCACAAATTCCCGATACGCCGATCTACTGTACAGCCAACGCTATCGACTCGATa      <  1:3047041/69‑1 (MQ=255)
                                                 cTGACCGAACTGATGGCACAAATTCCCGATACGCCGATCTACTGTACAGCCAACGCTATCGACTCGATa      >  1:2272683/1‑69 (MQ=255)
                                                 cTGACCGAACTGATGGCACAAATTCCCGATACGCCGATCTACTGTACAGCCAACGCTATCGACTCGATa      >  1:241078/1‑69 (MQ=255)
                                                 cTGACCGAACTGATGGCACAAATTCCCGATACGCCGATCTACTGTACAGCCAACGCTATCGACTCGATa      <  1:1099639/69‑1 (MQ=255)
                                                 cTGACCGAACTGATGGCACAAATTCCCGATACGCCGATCTACTGTACAGCCAACGCTATCGACTCGATa      >  1:484597/1‑69 (MQ=255)
                                                 cTGACCGAACTGATGGCACAAATTCCCGATACGCCGATCTACTGTACAGCCAACGCTATCGACTCGATa      >  1:270272/1‑69 (MQ=255)
                                                 cTGACCGAACTGATGGCACAAATTCCCGATACGCCGATCTACTGTACAGCCAACGCTATCGACTCGATa      >  1:2796662/1‑69 (MQ=255)
                                                 cTGACCGAACTGATGGCACAAATTCCCGATACGCCGATCTACTGTACAGCCAACGCTATCGACTCGATa      <  1:412551/69‑1 (MQ=255)
                                                 cTGACCGAACTGATGGCACAAATTCCCGATACGCCGATCTACTGTACAGCCAACGCTATCGACTCGATa      <  1:392508/69‑1 (MQ=255)
                                                 cTGACCGAACTGATGGCACAAATTCCCGATACGCCGATCTACTGTACAGCCAACGCTATCGACTCGATa      >  1:3044390/1‑69 (MQ=255)
                                                  tGACCGAACTGATGGCGCAAATTCCCGATACGCCGATCTACTg                               <  1:3414057/43‑1 (MQ=255)
                                                            gATGGCACAAATTCCCGATACGCCGATCTACTGTACAGCCAACGCTATCGACTCGATAAATg  >  1:961784/1‑62 (MQ=255)
                                                                  |                                                       
TATCGATTACATCGTGATTAACCATGCAGAAGAGGACCACGCTGGGGCGCTGACCGAACTGATGGCACAAATTCCCGATACGCCGATCTACTGTACAGCCAACGCTATCGACTCGATAAATG  >  minE/1820906‑1821027

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: