Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 936,833 G→A 34.0% G131G (GGG→GGA fdnI → formate dehydrogenase‑N, cytochrome B556 (gamma) subunit, nitrate‑inducible

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE936,8330GA34.0% 69.2 / 47.1 88G131G (GGG→GGAfdnIformate dehydrogenase‑N, cytochrome B556 (gamma) subunit, nitrate‑inducible
Reads supporting (aligned to +/- strand):  ref base G (21/37);  new base A (18/12);  total (39/49)
Fisher's exact test for biased strand distribution p-value = 4.26e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.44e-01

Only 100 of 105 total aligned reads displayed.

TACAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGC  >  minE/936765‑936887
                                                                    |                                                      
tACAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTTTCGTGCTGCTGGTGACCggg                                                        >  1:1972140/1‑69 (MQ=255)
tACAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTTTCGTGCTGCTGGTGACCggg                                                        >  1:2457304/1‑69 (MQ=255)
tACAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTTTCGTGCTGCTGGTGACCggg                                                        >  1:1752189/1‑69 (MQ=255)
tACAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTTTCGTGCTGCTGGTGACCggg                                                        >  1:2627600/1‑69 (MQ=255)
tACAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTTTCGTGCTGCTGGTGACCggg                                                        >  1:2722918/1‑69 (MQ=255)
tACAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTTTCGTGCTGCTGGTGACCggg                                                        >  1:1439101/1‑69 (MQ=255)
tACAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTTTCGTGCTGCTGGTGACCggg                                                        >  1:3035871/1‑69 (MQ=255)
 aCAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTTTCGTGCTGCTGGTGACCggg                                                        >  1:2177211/1‑68 (MQ=255)
 aCAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTTTCGTGCTGCTGGTGACCggg                                                        >  1:2038690/1‑68 (MQ=255)
 aCAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTTTCGTGCTGCTGGTGACCggg                                                        >  1:255210/1‑68 (MQ=255)
 aCAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTTTCGTGCTGCTGGTGACCggg                                                        >  1:172453/1‑68 (MQ=255)
 aCAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTTTCGTGCTGCTGGTGACCggg                                                        >  1:2617401/1‑68 (MQ=255)
 aCAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTTTCGTGCTGCTGGTGACCggg                                                        >  1:2661393/1‑68 (MQ=255)
 aCAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTTTCGTGCTGCTGGTGACCggg                                                        >  1:2913962/1‑68 (MQ=255)
 aCAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTTTCGTGCTGCTGGTGACCggg                                                        >  1:1036188/1‑68 (MQ=255)
 aCAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTTTCGTGCTGCTGGTGACCggg                                                        >  1:79657/1‑68 (MQ=255)
 aCAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTTTCGTGCTGCTGGTGACCggg                                                        >  1:657683/1‑68 (MQ=255)
  cAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTTTCGTGCTTCTGGTGACCGGGGt                                                      >  1:829673/1‑69 (MQ=255)
  cAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGGGt                                                      >  1:2645247/1‑69 (MQ=255)
  cAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTTTCGTGCTGCTGGTGACAGGGGt                                                      >  1:176011/1‑69 (MQ=255)
     cGCCGGGCAAAAGATGATGTTCTGGTCAATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGAGTGAt                                                   >  1:1444764/1‑69 (MQ=255)
     cGCCGGGCAAAAGATGATGTTCTGGTCAATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGAGTGAt                                                   >  1:2321300/1‑69 (MQ=255)
     cGCCGGGCAAAAGATGATGTTCTGGTCAATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGAGTGAt                                                   >  1:243821/1‑69 (MQ=255)
     cGCCGGGCAAAAGATGATGTTCTGGTCAATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGAGTGAt                                                   >  1:743311/1‑69 (MQ=255)
     cGCCGGGCAAAAGATGATGTTCTGGTCAATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGAGTGAt                                                   >  1:1104724/1‑69 (MQ=255)
     cGCCGGGCAAAAGATGATGTTCTGGTCAATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGAGTGAt                                                   >  1:666398/1‑69 (MQ=255)
     cGCCGGGCAAAAGATGATGTTCTGGTCAATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGAGTGAt                                                   >  1:2562863/1‑69 (MQ=255)
     cGCCGGGCAAAAGATGATGTTCTGGTCAATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGAGTGAt                                                   >  1:1168608/1‑69 (MQ=255)
     cGCCGGGCAAAAGATGATGTTCTGGTCAATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGAGTGAt                                                   >  1:171127/1‑69 (MQ=255)
     cGCCGGGCAAAAGATGATGTTCTGGTCAATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGAGTGAt                                                   >  1:2630346/1‑69 (MQ=255)
     cGCCGGGCAAAAGATGATGTTCTGGTCAATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGAGTGAt                                                   >  1:3263424/1‑69 (MQ=255)
     cGCCGGGCAAAAGATGATGTTCTGGTCAATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGAGTGAt                                                   >  1:1515590/1‑69 (MQ=255)
     cGCCGGGCAAAAGATGATGTTCTGGTCAATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGAGTGAt                                                   >  1:1369050/1‑69 (MQ=255)
     cGCCGGGCAAAAGATGATGTTCTGGTCAATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGAGTGAt                                                   >  1:3256345/1‑69 (MQ=255)
     cGCCGGGCAAAAGATGATGTTCTGGTCAATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGAGTGAt                                                   >  1:2981727/1‑69 (MQ=255)
     cGCCGGGCAAAAGATGATGTTCTGGTCAATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGAGTGAt                                                   >  1:31871/1‑69 (MQ=255)
     cGCCGGGCAAAAGATGATATTCTGGTCAATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGAGTGAt                                                   >  1:862254/1‑69 (MQ=255)
         gggCAAAAGATGATGTTCTGGTCAATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGAGTGATTATc                                               <  1:3423384/69‑1 (MQ=255)
         gggCAAAAGATGATGTTCTGGTCAATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGAGTGATTATc                                               <  1:532838/69‑1 (MQ=255)
         gggCAAAAGATGATGTTCTGGTCAATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGAGTGATTATc                                               <  1:3010799/69‑1 (MQ=255)
         gggCAAAAGATGATGTTCTGGTCAATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGAGTGATTATc                                               <  1:569020/69‑1 (MQ=255)
         gggCAAAAGATGATGTTCTGGTCAATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGAGTGATTATc                                               <  1:2976192/69‑1 (MQ=255)
         gggCAAAAGATGATGTTCTGGTCAATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGAGTGATTATc                                               <  1:2972919/69‑1 (MQ=255)
         gggCAAAAGATGATGTTCTGGTCAATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGAGTGATTATc                                               <  1:621179/69‑1 (MQ=255)
         gggCAAAAGATGATGTTCTGGTCAATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGAGTGATTATc                                               <  1:2774501/69‑1 (MQ=255)
         gggCAAAAGATGATGTTCTGGTCAATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGAGTGATTATc                                               <  1:1666512/69‑1 (MQ=255)
         gggCAAAAGATGATGTTCTGGTCAATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGAGTGATTATc                                               <  1:2527304/69‑1 (MQ=255)
         gggCAAAAGATGATGTTCTGGTCAATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGAGTGATTATc                                               <  1:1206162/69‑1 (MQ=255)
          ggCAAAAGATGATGTTCTGGTCAATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGATTGATTATc                                               <  1:2169753/68‑1 (MQ=255)
           gCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTTTCGTGCTTCTGGTGACCGGGGTGATTATCTg                                             <  1:1141148/69‑1 (MQ=255)
           gCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTg                                             <  1:1236106/69‑1 (MQ=255)
           gCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTg                                             <  1:308241/69‑1 (MQ=255)
           gCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTg                                             <  1:2968661/69‑1 (MQ=255)
           gCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTg                                             <  1:1455318/69‑1 (MQ=255)
           gCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTg                                             <  1:2439790/69‑1 (MQ=255)
           gCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTg                                             <  1:604910/69‑1 (MQ=255)
            cAAAAGATGATGTTCTGGTCGATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTgg                                            >  1:1713775/1‑69 (MQ=255)
            cAAAAGATGATGTTCTGGTCGATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTgg                                            >  1:401671/1‑69 (MQ=255)
            cAAAAGATGATGTTCTGGTCGATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTgg                                            >  1:3316805/1‑69 (MQ=255)
            cAAAAGATGATGTTCTGGTCGATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTgg                                            >  1:2056914/1‑69 (MQ=255)
            cAAAAGATGATGTTCTGGTCGATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTgg                                            >  1:55730/1‑69 (MQ=255)
            cAAAAGATGATGTTCTGGTCGATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTgg                                            >  1:2079369/1‑69 (MQ=255)
            cAAAAGATGATGTTCTGGTCGATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTgg                                            >  1:2562800/1‑69 (MQ=255)
            cAAAAGATGATGTTCTGGTCGATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTgg                                            >  1:1705270/1‑69 (MQ=255)
            cAAAAGATGATGTTCTGGTCGATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTgg                                            >  1:2743105/1‑69 (MQ=255)
            cAAAAGATGATGTTCTGGTCGATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTgg                                            >  1:1097634/1‑69 (MQ=255)
            cAAAAGATGATGTTCTGGTCGATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTgg                                            >  1:773253/1‑69 (MQ=255)
            cAAAAGATGATGTTCTGGTCGATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTgg                                            >  1:1485734/1‑69 (MQ=255)
            cAAAAGATGATGTTCTGGTCGATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTgg                                            >  1:2652830/1‑69 (MQ=255)
            cAAAAGATGATGTTCTGGTCGATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTg                                             >  1:3126897/1‑68 (MQ=255)
            cAAAAGATGATGTTCTGGTCGATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTg                                             <  1:509704/68‑1 (MQ=255)
                aGATGATGTTCTGGTCGATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTc                                        <  1:2023336/69‑1 (MQ=255)
                aGATGATGTTCTGGTCGATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTc                                        <  1:1712804/69‑1 (MQ=255)
                aGATGATGTTCTGGTCGATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTc                                        <  1:287819/69‑1 (MQ=255)
                  atgatgTTCTGGTCGATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCg                                      <  1:3137687/69‑1 (MQ=255)
                  atgatgTTCTGGTCGATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCg                                      <  1:3191467/69‑1 (MQ=255)
                  atgatgTTCTGGTCGATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCg                                      <  1:4748/69‑1 (MQ=255)
                      tgTTCTTGTCGATGATGAGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTg                                             <  1:1381294/58‑1 (MQ=38)
                               cGATCATGAGCATGATTTTGGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAg                          <  1:2750882/68‑1 (MQ=255)
                               cGATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAg                          <  1:634768/68‑1 (MQ=255)
                               cGATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAg                          <  1:1081168/68‑1 (MQ=255)
                               cGATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAg                          <  1:829083/68‑1 (MQ=255)
                               cGATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAg                          <  1:1065336/68‑1 (MQ=255)
                               cGATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAg                          <  1:3191152/68‑1 (MQ=255)
                               cGATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAg                          <  1:2281167/68‑1 (MQ=255)
                               cGATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAg                          <  1:2312204/68‑1 (MQ=255)
                               cGATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAg                          <  1:2408709/68‑1 (MQ=255)
                               cGATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAg                          <  1:1883471/68‑1 (MQ=255)
                               cGATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAg                          <  1:1873666/68‑1 (MQ=255)
                               cGATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAg                          <  1:1829736/68‑1 (MQ=255)
                               cGATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAg                          <  1:2788064/68‑1 (MQ=255)
                               cGATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAg                          <  1:1447619/68‑1 (MQ=255)
                               cGATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAg                          <  1:1402148/68‑1 (MQ=255)
                               cGATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAg                          <  1:138261/68‑1 (MQ=255)
                                gATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTa                        <  1:1690821/69‑1 (MQ=255)
                                       aGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCg                 >  1:1085605/1‑69 (MQ=255)
                                           tGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGc             >  1:3212541/1‑69 (MQ=255)
                                             aTTTTCGTGCTGCTGGTGACCGGGGTGATTATCTg                                             <  1:2856570/35‑1 (MQ=255)
                                                      ctgctgGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGc  <  1:1648568/69‑1 (MQ=255)
                                                      ctgctgGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGc  <  1:703294/69‑1 (MQ=255)
                                                                    |                                                      
TACAACGCCGGGCAAAAGATGATGTTCTGGTCGATCATGAGCATGATTTTCGTGCTGCTGGTGACCGGGGTGATTATCTGGCGTCCGTACTTTGCGCAGTACTTCCCGATGCAGGTTGTTCGC  >  minE/936765‑936887

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: