Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,985,612 T→C 28.6% Q316Q (CAA→CAG yggM ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,985,6120TC28.6% 102.6 / 44.1 101Q316Q (CAA→CAGyggMconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (53/19);  new base C (13/16);  total (66/35)
Fisher's exact test for biased strand distribution p-value = 1.03e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 110 total aligned reads displayed.

TATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACTTTTTTACTGTTTATTTCTGAAACTTTCAT  >  minE/1985544‑1985677
                                                                    |                                                                 
tatCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCt                                                                   >  1:1559974/1‑69 (MQ=255)
tatCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCt                                                                   >  1:793101/1‑69 (MQ=255)
tatCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCt                                                                   >  1:603674/1‑69 (MQ=255)
tatCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCt                                                                   >  1:50300/1‑69 (MQ=255)
tatCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCt                                                                   >  1:1999718/1‑69 (MQ=255)
tatCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCt                                                                   >  1:2790491/1‑69 (MQ=255)
tatCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCt                                                                   >  1:2806314/1‑69 (MQ=255)
tatCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCt                                                                   >  1:2816816/1‑69 (MQ=255)
  tCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTg                                                                 >  1:1049070/1‑69 (MQ=255)
  tCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTg                                                                 >  1:1654197/1‑69 (MQ=255)
  tCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTg                                                                 >  1:275807/1‑69 (MQ=255)
  tCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTg                                                                 >  1:42931/1‑69 (MQ=255)
    aaaCTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTgg                                                                >  1:1170712/1‑68 (MQ=255)
    aaaCTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTgg                                                                >  1:2756532/1‑68 (MQ=255)
          acttaATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTgg                                                                >  1:1665282/1‑62 (MQ=255)
          acttaATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTg                                                             >  1:1214318/1‑65 (MQ=255)
          acttaATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTg                                                             >  1:921089/1‑65 (MQ=255)
            ttaATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATc                                                        >  1:3337103/1‑68 (MQ=255)
            ttaATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATc                                                        >  1:2398623/1‑68 (MQ=255)
            ttaATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATc                                                        >  1:467352/1‑68 (MQ=255)
            ttaATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATc                                                        >  1:2629873/1‑68 (MQ=255)
            ttaATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATc                                                        >  1:1503854/1‑68 (MQ=255)
            ttaATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATc                                                        >  1:134277/1‑68 (MQ=255)
                   ccccATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAAtt                                                 >  1:1135267/1‑68 (MQ=255)
                    cccATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCCTGGTTGAATCGTCAATTTtc                                              <  1:3365198/69‑1 (MQ=255)
                     ccATTGCGCCGCAGTTGTTCTGTATCTTTTTTCCAGCGTTCATTAGCCTGGTTGATATCGTCAATTTtc                                              <  1:40479/69‑1 (MQ=255)
                     ccATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCCTGGTTGATATCGTCAATTTtc                                              <  1:2992399/69‑1 (MQ=255)
                     ccATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCCTGGTTGATATCGTCAATTTtc                                              <  1:529131/69‑1 (MQ=255)
                     ccATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCCTGGTTGATATCGTCAATTTtc                                              <  1:2913647/69‑1 (MQ=255)
                     ccATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCCTGGTTGATATCGTCAATTTtc                                              <  1:2688098/69‑1 (MQ=255)
                     ccATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCCTGGTTGATATCGTCAATTTtc                                              <  1:213984/69‑1 (MQ=255)
                     ccATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCCTGGTTGATATCGTCAATTTtc                                              <  1:824262/69‑1 (MQ=255)
                     ccATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCCTGGTTGATATCGTCAATTTtc                                              <  1:375321/69‑1 (MQ=255)
                     ccATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCCTGGTTGATATCGTCAATTTtc                                              <  1:1922632/69‑1 (MQ=255)
                     ccATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCCTGGTTGATATCGTCAATTTtc                                              <  1:1810931/69‑1 (MQ=255)
                     ccATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCCTGGTTGATATCGTCAATTTtc                                              <  1:3077403/69‑1 (MQ=255)
                     ccATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCCTGGTTGATATCGTCAATTTtc                                              <  1:1476329/69‑1 (MQ=255)
                      cATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCCTGGTTGATATCGTCAATTTtc                                              <  1:1850614/68‑1 (MQ=255)
                      cATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCCTGGTTGATATCGTCAATTTtc                                              <  1:1369242/68‑1 (MQ=255)
                         tGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTtcatc                                           >  1:1914037/1‑68 (MQ=255)
                         tGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTtcatc                                           >  1:1910183/1‑68 (MQ=255)
                          gcgcCGCAGTTTTTCGGTATCTTTTTTCCAGCGTTCATTAGCCTGGTTGATATCGTCAATTTTCATCaa                                         >  1:930925/1‑69 (MQ=255)
                          gcgcCGCAGTTGTTCGGTATCTTTTTTCCTGCGTTCATTAGCCTGGTTGATATCGTCAATTTTCATCaa                                         >  1:250745/1‑69 (MQ=255)
                          gcgcCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCaa                                         >  1:2832101/1‑69 (MQ=255)
                          gcgcCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCaa                                         >  1:368381/1‑69 (MQ=255)
                          gcgcCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCaa                                         >  1:116163/1‑69 (MQ=255)
                          gcgcCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCaa                                         >  1:1529158/1‑69 (MQ=255)
                          gcgcCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCaa                                         >  1:3398061/1‑69 (MQ=255)
                          gcgcCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCaa                                         >  1:2949848/1‑69 (MQ=255)
                          gcgcCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCaa                                         >  1:886951/1‑69 (MQ=255)
                          gcgcCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCaa                                         >  1:1685525/1‑69 (MQ=255)
                          gcgcCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCCTGGTTGATATCGTCAATTTtcataa                                          >  1:13994/1‑66 (MQ=255)
                          gcgcCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCCTGGTTGATATCGTCAATTTTCATCaa                                         >  1:599114/1‑69 (MQ=255)
                          gcgcCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCCTGGTTGATATCGTCAATTTTCATCaa                                         >  1:1127965/1‑69 (MQ=255)
                          gcgcCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCCTGGTTGATATCGTCAATTTTCATCaa                                         >  1:2258232/1‑69 (MQ=255)
                          gcgcCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCCTGGTTGATATCGTCAATTTTCATCaa                                         >  1:2245561/1‑69 (MQ=255)
                          gcgcCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCCTGGTTGATATCGTCAATTTTCATCaa                                         >  1:610098/1‑69 (MQ=255)
                          gcgcCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCCTGGTTGATATCGTCAATTTTCATCaa                                         >  1:82640/1‑69 (MQ=255)
                          gcgcCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCCTGGTTGATATCGTCAATTTTCATCaa                                         >  1:2741983/1‑69 (MQ=255)
                          gcgcCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCCTGGTTGATATCGTCAATTTTCATCaa                                         >  1:1809747/1‑69 (MQ=255)
                          gcgcCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCCTGGTTGATATCGTCAATTTTCATCaa                                         >  1:1363201/1‑69 (MQ=255)
                                    tgttCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACt                               >  1:553014/1‑69 (MQ=255)
                                    tgttCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACt                               >  1:2689075/1‑69 (MQ=255)
                                    tgttCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACt                               >  1:142267/1‑69 (MQ=255)
                                    tgttCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACt                               >  1:1962847/1‑69 (MQ=255)
                                    tgttCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACt                               >  1:978102/1‑69 (MQ=255)
                                          gTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACtttttt                          <  1:1642559/68‑1 (MQ=255)
                                          gTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACtttttt                          <  1:3167997/68‑1 (MQ=255)
                                          gTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACtttttt                          <  1:983642/68‑1 (MQ=255)
                                          gTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACtttttt                          <  1:2904555/68‑1 (MQ=255)
                                          gTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACtttttt                          <  1:1573848/68‑1 (MQ=255)
                                          gTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACtttttt                          <  1:2491016/68‑1 (MQ=255)
                                          gTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACtttttt                          <  1:27227/68‑1 (MQ=255)
                                           tATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACtttttt                          <  1:2053203/67‑1 (MQ=255)
                                            aTCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACtttttt                          <  1:2512241/66‑1 (MQ=255)
                                               ttttttCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACTTTTTTACTGtt                    >  1:2870637/1‑69 (MQ=255)
                                               ttttttCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACTTTTTTACTGtt                    >  1:1304936/1‑69 (MQ=255)
                                               ttttttCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACTTTTTTACTGtt                    >  1:1455401/1‑69 (MQ=255)
                                               ttttttCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACTTTTTTACTGtt                    >  1:378299/1‑69 (MQ=255)
                                               ttttttCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACTTTTTTACTGtt                    >  1:1435955/1‑69 (MQ=255)
                                               ttttttCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACTTTTTTACTGtt                    >  1:485248/1‑69 (MQ=255)
                                               ttttttCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACTTTTTTACTGtt                    >  1:3186579/1‑69 (MQ=255)
                                               ttttttCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACTTTTTTACTGtt                    >  1:2053749/1‑69 (MQ=255)
                                                                 aGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACTTTTTTACTGtt                    <  1:2954764/51‑1 (MQ=255)
                                                                 aGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACTTTTTTACTGtt                    <  1:1718569/51‑1 (MQ=255)
                                                                 aGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACTTTTTTACTGtt                    <  1:1119207/51‑1 (MQ=255)
                                                                 aGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACTTTTTTACTGtt                    <  1:737923/51‑1 (MQ=255)
                                                                 aGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACTTTTTTACTGtt                    <  1:731207/51‑1 (MQ=255)
                                                                 aGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACTTTTTTACTGtt                    <  1:2391190/51‑1 (MQ=255)
                                                                 aGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACTTTTTTACTGtt                    <  1:2769293/51‑1 (MQ=255)
                                                                 aGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACTTTTTTACTGTTTATTTCTGAAACTTTCAt  >  1:673911/1‑69 (MQ=255)
                                                                 aGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACTTTTTTACTGTTTATTTCTGAAACTTTCAt  >  1:2407731/1‑69 (MQ=255)
                                                                 aGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACTTTTTTACTGTTTATTTCTGAAACTTTCAt  >  1:2511015/1‑69 (MQ=255)
                                                                 aGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACTTTTTTACTGTTTATTTCTGAAACTTTCAt  >  1:1020472/1‑69 (MQ=255)
                                                                 aGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACTTTTTTACTGTTTATTTCTGAAACTTTCAt  >  1:1223726/1‑69 (MQ=255)
                                                                 aGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACTTTTTTACTGTTTATTTCTGAAACTTTCAt  >  1:2708296/1‑69 (MQ=255)
                                                                 aGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACTTTTTTACTGTTTATTTCTGAAACTTTCAt  >  1:1723569/1‑69 (MQ=255)
                                                                 aGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACTTTTTTACTGTTTATTTCTGAAACTTTCAt  >  1:1102280/1‑69 (MQ=255)
                                                                 aGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACTTTTTTACTGTTTATTTCTGAAACTTTCAt  >  1:1280580/1‑69 (MQ=255)
                                                                 aGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACTTTTTTACTGTTTATTTCTGAAACTTTCAt  >  1:3355869/1‑69 (MQ=255)
                                                                    |                                                                 
TATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACTTTTTTACTGTTTATTTCTGAAACTTTCAT  >  minE/1985544‑1985677

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: