Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 20,146 C→G 49.2% P26P (CCC→CCG ileS → isoleucyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE20,1460CG49.2% ‑1.2 / 82.2 67P26P (CCC→CCGileSisoleucyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base C (9/25);  new base G (20/13);  total (29/38)
Fisher's exact test for biased strand distribution p-value = 6.75e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCCGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCG  >  minE/20082‑20197
                                                                |                                                   
actCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCTAACCCgg                                                   <  1:2603158/65‑1 (MQ=255)
aaTCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCCgg                                                   <  1:3311931/67‑1 (MQ=255)
aaTCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCCgg                                                   <  1:1109148/67‑1 (MQ=255)
aaTCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCCgg                                                   <  1:1986150/67‑1 (MQ=255)
aaTCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCCgg                                                   <  1:2521511/67‑1 (MQ=255)
aaTCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCCgg                                                   <  1:3186027/67‑1 (MQ=255)
aaTCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCCgg                                                   <  1:918004/67‑1 (MQ=255)
aaTCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCCgg                                                   <  1:2739172/67‑1 (MQ=255)
aaTCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCCgg                                                   <  1:1601079/67‑1 (MQ=255)
aaTCAACACTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCCgg                                                   <  1:2292560/67‑1 (MQ=255)
 aTCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAAccc                                                     >  1:3047416/1‑64 (MQ=255)
 aTCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAAccc                                                     >  1:2592568/1‑64 (MQ=255)
 aTCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCCg                                                    >  1:2060799/1‑65 (MQ=255)
 aTCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCCg                                                    >  1:700492/1‑65 (MQ=255)
 aTCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCCg                                                    >  1:1584538/1‑65 (MQ=255)
 aTCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCCg                                                    >  1:2844030/1‑65 (MQ=255)
 aTCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCCg                                                    >  1:577636/1‑65 (MQ=255)
 aTCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCCg                                                    >  1:968405/1‑65 (MQ=255)
 aTCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCCg                                                    >  1:3288858/1‑65 (MQ=255)
 aTCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCCg                                                    >  1:982970/1‑65 (MQ=255)
  tgaaCCCTGAATTTGCCGGAACCAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCCg                                                    <  1:2226482/62‑1 (MQ=255)
  tCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCCg                                                    <  1:355673/64‑1 (MQ=255)
  tCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCCg                                                    <  1:3376138/64‑1 (MQ=255)
  tCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCCg                                                    <  1:982313/64‑1 (MQ=255)
             tttGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCTCGAACCCGGAATGCTGGCGCGTTg                                    <  1:2304709/69‑1 (MQ=255)
             tttGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCCGGAATGCTGGCGCGTTg                                    <  1:2691925/69‑1 (MQ=255)
             tttGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCCGGAATGCTGGCGCGTTg                                    <  1:2791547/69‑1 (MQ=255)
              ttGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCCGGAATGCTGGCGCGTTgg                                   <  1:888183/69‑1 (MQ=255)
              ttGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCCGGAATGCTGGCGCGTTg                                    <  1:549314/68‑1 (MQ=255)
                  cGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCCGGAATGCTGGCGCGTTGGACt                                <  1:2266069/68‑1 (MQ=255)
                  cGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCCGGAATGCTGGCGCGTTGGACt                                <  1:1283375/68‑1 (MQ=255)
                  cGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCCGGAATGCTGGCGCGTTGGACt                                <  1:1648750/68‑1 (MQ=255)
                  cGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCCGGAATGCTGGCGCGTTGGACt                                <  1:2778855/68‑1 (MQ=255)
                     aaaCAGGGTTCCCGATGCGTGGCGATCTCACCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACtgatg                            <  1:2084455/69‑1 (MQ=255)
                     aaaCAGGGTTCCCGATGCGTGGCGATCTCACCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACtgatg                            <  1:388390/69‑1 (MQ=255)
                     aaaCAGGGTTCCCGATGCGTGGCGATCTCACCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACtgatg                            <  1:992164/69‑1 (MQ=255)
                     aaaCAGGGTTCCCGATGCGTGGCGATCTCACCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACtgatg                            <  1:1542276/69‑1 (MQ=255)
                     aaaCAGGGTTCCCGATGCGTGGCGATCTCACCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACtgatg                            <  1:1047663/69‑1 (MQ=255)
                     aaaCAGGGTTCCCGATGCGTGGCGATCTCACCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACtgatg                            <  1:1834420/69‑1 (MQ=255)
                     aaaCAGGGTTCCCGATGCGTGGCGATCTCACCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACtgatg                            <  1:2040546/69‑1 (MQ=255)
                      aaCAGGGTTCCCGATGCGTGGCGATCTCACCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACtgatg                            <  1:815649/68‑1 (MQ=255)
                               cccGATGCGTGGCGATCTCGCCAAGCGCGAACCCgg                                                   <  1:1905725/36‑1 (MQ=255)
                                   aTGCGTGGCGATCTCGCCAAGCGCGCACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTa                  >  1:1225877/1‑65 (MQ=255)
                                   aTGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACgg               >  1:16794/1‑68 (MQ=255)
                                   aTGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACgg               >  1:169608/1‑68 (MQ=255)
                                   aTGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACgg               >  1:2953645/1‑68 (MQ=255)
                                   aTGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACgg               >  1:1995898/1‑68 (MQ=255)
                                   aTGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGc              >  1:993796/1‑69 (MQ=255)
                                   aTGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGc              >  1:625094/1‑69 (MQ=255)
                                   aTGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGc              >  1:698547/1‑69 (MQ=255)
                                   aTGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGc              >  1:884010/1‑69 (MQ=255)
                                   aTGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGc              >  1:329527/1‑69 (MQ=255)
                                   aTGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGc              >  1:1121240/1‑69 (MQ=255)
                                   aTGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGc              >  1:1199779/1‑69 (MQ=255)
                                   aTGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGc              >  1:1240134/1‑69 (MQ=255)
                                   aTGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGc              >  1:1421508/1‑69 (MQ=255)
                                   aTGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGc              >  1:2207895/1‑69 (MQ=255)
                                   aTGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGc              >  1:2717864/1‑69 (MQ=255)
                                   aTGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGc              >  1:2798419/1‑69 (MQ=255)
                                   aTGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGc              >  1:3016026/1‑69 (MQ=255)
                                   aTGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGc              >  1:3074654/1‑69 (MQ=255)
                                   aTGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGc              >  1:3363649/1‑69 (MQ=255)
                                            gATCTCGCCAAGCGCGAACCCGGAATGCTGGCGCGTTGGACt                                <  1:1035736/42‑1 (MQ=255)
                                               ctcGCCAAGCGCGAACCCGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCg  >  1:2404270/1‑69 (MQ=255)
                                                    cAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATTCGTGc   <  1:2949751/63‑1 (MQ=255)
                                                    cAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATTCGTGc   <  1:1315643/63‑1 (MQ=255)
                                                    cAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATTCGTGc   <  1:3155561/63‑1 (MQ=255)
                                                     aaGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATTCGTGc   <  1:923708/62‑1 (MQ=255)
                                                      aGCGGGAACCGGGAATGCTGGCGCGTTGGACtgatg                            <  1:1517713/36‑1 (MQ=37)
                                                                |                                                   
AATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCCGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCG  >  minE/20082‑20197

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: