Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,576,174 G→C 68.6% A290A (GCC→GCG xapB ← xanthosine transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,576,1740GC68.6% 99.5 / 60.3 86A290A (GCC→GCGxapBxanthosine transporter
Reads supporting (aligned to +/- strand):  ref base G (18/9);  new base C (37/22);  total (55/31)
Fisher's exact test for biased strand distribution p-value = 8.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAGCAAAATAAATCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCAAAGCGCAGCGTCCAGGCCACCATACTCATCAGCATGACGGTTTTAATGCCAAATCGCTTTAAAAAGAA  >  minE/1576110‑1576227
                                                                |                                                     
gcagcaAAATAAATACGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCAAAGCGCAGCGTCCAGGcc                                                    >  1:936849/1‑68 (MQ=255)
 cagcaAAATAAATCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCAAAGCGCAGCGTCCAGGccac                                                  <  1:2103807/69‑1 (MQ=255)
 cagcaAAATAAATCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCAAAGCGCAGCGTCCAGGccac                                                  <  1:1258314/69‑1 (MQ=255)
 cagcaAAATAAATCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCAAAGCGCAGCGTCCAGGccac                                                  <  1:90322/69‑1 (MQ=255)
 cagcaAAATAAATCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCAAAGCGCAGCGTCCAGGccac                                                  <  1:53191/69‑1 (MQ=255)
     aaaaTAAATCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCAAAGCGCAGCGTCCAGGCCACCATa                                              <  1:1561827/69‑1 (MQ=255)
     aaaaTAAATCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCAAAGCGCAGCGTCCAGGCCACCATa                                              <  1:1341646/69‑1 (MQ=255)
     aaaaTAAATCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCAAAGCGCAGCGTCCAGGCCACCATa                                              <  1:362461/69‑1 (MQ=255)
     aaaaTAAATCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCAAAGCGCAGCGTCCAGGCCACCATa                                              <  1:1398417/69‑1 (MQ=255)
      aaataaatCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCAAAGCGCAGCGTCCAGGcc                                                    >  1:2193407/1‑62 (MQ=255)
      aaataaatCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCAAAGCGCAGCGTCCAGGcc                                                    >  1:1572764/1‑62 (MQ=255)
      aaataaatCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCAAAGCGCAGCGTCCAGGcc                                                    >  1:3182486/1‑62 (MQ=255)
      aaataaatCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCAAAGCGCAGCGTCCAGGcc                                                    >  1:3247490/1‑62 (MQ=255)
      aaataaatCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCAAAGCGCAGCGTCCAGGcc                                                    >  1:760793/1‑62 (MQ=255)
      aaataaatCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCAAAGCGCAGCGTCCAGGcc                                                    >  1:339684/1‑62 (MQ=255)
      aaataaatCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCAAAGCGCAGCGTCCAGGcc                                                    >  1:291942/1‑62 (MQ=255)
                       gatggATCGCCATAGGCAAAGAAGCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGtt                            <  1:357314/66‑1 (MQ=255)
                       gatggATCGCCATAGGCAAAGAAGCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGtt                            <  1:1800702/66‑1 (MQ=255)
                       gatggATCGCCATAGGCAAAGAAGCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGtt                            <  1:1863293/66‑1 (MQ=255)
                       gatggATCGCCATAGGCAAAGAAGCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGtt                            <  1:186537/66‑1 (MQ=255)
                       gatggATCGCCATAGGCAAAGAAGCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGtt                            <  1:1596031/66‑1 (MQ=255)
                       gatggATCGCCATAGGCAAAGAAGCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGtt                            <  1:1515360/66‑1 (MQ=255)
                       gatggATCGCCATAGGCAAAGAAGCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGtt                            <  1:3171839/66‑1 (MQ=255)
                       gatggATCGCCATAGGCAAAGAAGCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGtt                            <  1:2580671/66‑1 (MQ=255)
                       gACGGATCGCCATAGGCGAAGAAGCCAAAGCGCAGCGTCCAGGccac                                                  <  1:2309410/47‑1 (MQ=255)
                               gCCATAGGCGAAGAAGCCAAAGCGCAGCGTCCAGGCCACCATACTCATCAGCATGACGGTTTTAATGc                     >  1:1276088/1‑68 (MQ=255)
                               gCCATAGGCGAAGAAGCCAAAGCGCAGCGTCCAGGCCACCATACTCATCAGCATGACGCTTTTAATGc                     >  1:3136816/1‑68 (MQ=255)
                                    aGGCAAAGAAGCCAAAGCGCAGCGTCCACGCCACCATACacatc                                        >  1:996730/1‑44 (MQ=37)
                                    aGGCAAAGAAGCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATAccga                  >  1:1694877/1‑64 (MQ=255)
                                    aGGCAAAGAAGCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATAccga                  >  1:3321130/1‑64 (MQ=255)
                                    aGGCAAAGAAGCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATAccga                  >  1:45313/1‑64 (MQ=255)
                                    aGGCAAAGAAGCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATAccga                  >  1:2092086/1‑64 (MQ=255)
                                    aGGCAAAGAAGCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATAccga                  >  1:460814/1‑64 (MQ=255)
                                    aGGCAAAGAAGCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATAccga                  >  1:471413/1‑64 (MQ=255)
                                    aGGCAAAGAAGCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATAccga                  >  1:559691/1‑64 (MQ=255)
                                    aGGCAAAGAAGCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATAccga                  >  1:766413/1‑64 (MQ=255)
                                    aGGCAAAGAAGCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATAccga                  >  1:1369850/1‑64 (MQ=255)
                                    aGGCAAAGAAGCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATAccga                  >  1:836047/1‑64 (MQ=255)
                                    aGGCAAAGAAGCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATAccga                  >  1:3048902/1‑64 (MQ=255)
                                    aGGCAAAGAAGCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATAccga                  >  1:293713/1‑64 (MQ=255)
                                    aGGCAAAGAAGCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATAccga                  >  1:901355/1‑64 (MQ=255)
                                    aGGCAAAGAAGCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATAccga                  >  1:1123922/1‑64 (MQ=255)
                                    aGGCAAAGAAGCCAAAGCGCAGCGTCAACGCCACCATACTCATCAGCATGACGGTTTTAATAccga                  >  1:2616232/1‑64 (MQ=255)
                                     ggCGAAGAAGCCAAAGCGCAGCGTCCAGGCCACCATACTCATCAGCATGACGGTTTTAATGCCAAATCg              >  1:971656/1‑69 (MQ=255)
                                     ggCGAAGAAGCCAAAGCGCAGCGTCCAGGCCACCATACTCATCAGCATGACGGTTTTAATGCCAAATCg              >  1:447414/1‑69 (MQ=255)
                                             agCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATACCGAAACGCTTCaaaa      <  1:1292676/69‑1 (MQ=255)
                                             agCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATACCGAAACGCTTCaaaa      <  1:1417020/69‑1 (MQ=255)
                                             agCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATACCGAAACGCTTCaaaa      <  1:1474478/69‑1 (MQ=255)
                                             agCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATACCGAAACGCTTCaaaa      <  1:2457768/69‑1 (MQ=255)
                                             agCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATACCGAAACGCTTCaaaa      <  1:307447/69‑1 (MQ=255)
                                             agCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATACCGAAACGCTTCaaaa      <  1:2436243/69‑1 (MQ=255)
                                             agCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATACCGAAACGCTTCaaaa      <  1:2410078/69‑1 (MQ=255)
                                             agCCAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATACCGAAACGCTTCaaaa      <  1:3425575/69‑1 (MQ=255)
                                               ccAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATACCGAAACGCTTCaaaa      <  1:1293577/67‑1 (MQ=255)
                                                cAACGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGtt                            <  1:2079986/44‑1 (MQ=38)
                                                cAAAGCGCAGCGTCCAGGCCACCATACTCATCAGCATGACGGTTTTAATGCCGAATCGCTTTAAAAAGa   >  1:1796625/1‑69 (MQ=255)
                                                cAAAGCGCAGCGTCCAGGCCACCATACTCATCAGCATGACGGTTTTAATGCCAAATCGCTTTAAAAAGa   >  1:1130899/1‑69 (MQ=255)
                                                cAAAGCGCAGCGTCCAGGCCACCATACTCATCAGCATGACGGTTTTAATGCCAAATCGCTTTAAAAAGa   >  1:2576754/1‑69 (MQ=255)
                                                cAAAGCGCAGCGTCCAGGCCACCATACTCATCAGCATGACGGTTTTAATGCCAAATCGCTTTAAAAAGa   >  1:1642719/1‑69 (MQ=255)
                                                cAAAGCGCAGCGTCCAGGCCACCATACTCATCAGCATGACGGTTTTAATGCCAAATCGCTTTAAAAAGa   >  1:2745226/1‑69 (MQ=255)
                                                cAAAGCGCAGCGTCCAGGCCACCATACTCATCAGCATGACGGTTTTAATGCCAAATCGCTTTAAAAAGa   >  1:299675/1‑69 (MQ=255)
                                                cAAAGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGtt                            <  1:1622974/44‑1 (MQ=255)
                                                 aaaGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATACCGAAACGCTTCAAAAAGaa  >  1:585433/1‑69 (MQ=255)
                                                 aaaGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATACCGAAACGCTTCAAAAAGaa  >  1:2699936/1‑69 (MQ=255)
                                                 aaaGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATACCGAAACGCTTCAAAAAGaa  >  1:2798880/1‑69 (MQ=255)
                                                 aaaGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATACCGAAACGCTTCAAAAAGaa  >  1:2530902/1‑69 (MQ=255)
                                                 aaaGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATACCGAAACGCTTCAAAAAGaa  >  1:1313507/1‑69 (MQ=255)
                                                 aaaGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATACCGAAACGCTTCAAAAAGaa  >  1:83008/1‑69 (MQ=255)
                                                 aaaGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATACCGAAACGCTTCAAAAAGaa  >  1:807654/1‑69 (MQ=255)
                                                 aaaGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATACCGAAACGCTTCAAAAAGaa  >  1:793250/1‑69 (MQ=255)
                                                 aaaGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATACCGAAACGCTTCAAAAAGaa  >  1:709589/1‑69 (MQ=255)
                                                 aaaGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATACCGAAACGCTTCAAAAAGaa  >  1:601121/1‑69 (MQ=255)
                                                 aaaGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATACCGAAACGCTTCAAAAAGaa  >  1:1979483/1‑69 (MQ=255)
                                                 aaaGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATACCGAAACGCTTCAAAAAGaa  >  1:1558652/1‑69 (MQ=255)
                                                 aaaGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATACCGAAACGCTTCAAAAAGaa  >  1:533103/1‑69 (MQ=255)
                                                 aaaGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATACCGAAACGCTTCAAAAAGaa  >  1:1562002/1‑69 (MQ=255)
                                                 aaaGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATACCGAAACGCTTCAAAAAGaa  >  1:483524/1‑69 (MQ=255)
                                                 aaaGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATACCGAAACGCTTCAAAAAGaa  >  1:3014563/1‑69 (MQ=255)
                                                 aaaGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATACCGAAACGCTTCAAAAAGaa  >  1:3074580/1‑69 (MQ=255)
                                                 aaaGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATACCGAAACGCTTCAAAAAGaa  >  1:178617/1‑69 (MQ=255)
                                                 aaaGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATACCGAAACGCTTCAAAAAGaa  >  1:3177114/1‑69 (MQ=255)
                                                 aaaGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATACCGAAACGCTTCAAAAAGaa  >  1:2254404/1‑69 (MQ=255)
                                                 aaaGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATACCGAAACGCTTCAAAAAGaa  >  1:1878450/1‑69 (MQ=255)
                                                  aaGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGtt                            <  1:2464888/42‑1 (MQ=255)
                                                  aaGCGCAGCGTCCACGCCACCATACTCATCAGCATGACGGtt                            <  1:1164526/42‑1 (MQ=255)
                                                         gCGTCCACGCCACCATACTCATCAGCATGACGGTTTTAATACCGAAACGCTTCaaaa      <  1:164067/57‑1 (MQ=255)
                                                                |                                                     
GCAGCAAAATAAATCCGGTTGTTGACGGATCGCCATAGGCGAAGAAGCCAAAGCGCAGCGTCCAGGCCACCATACTCATCAGCATGACGGTTTTAATGCCAAATCGCTTTAAAAAGAA  >  minE/1576110‑1576227

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: