Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 224,636 G→A 20.4% F153F (TTC→TTT pepD ← aminoacyl‑histidine dipeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE224,6360GA20.4% 162.4 / 27.7 107F153F (TTC→TTTpepDaminoacyl‑histidine dipeptidase
Reads supporting (aligned to +/- strand):  ref base G (59/26);  new base A (13/9);  total (72/35)
Fisher's exact test for biased strand distribution p-value = 4.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.30e-01

Only 100 of 110 total aligned reads displayed.

TTTCACCTTCTTCTTCGGAGTCGGTGTTAATCAGAATATCAGCCTGCAACCAGTTGCCCTGTAAGCCGAACGCACCGTCCATACCGGCTTCTTCGGTCATGGTCAGCAGCACTTCCAGCGGGCCGTGAACCACGT  >  minE/224569‑224703
                                                                   |                                                                   
tttCACCTTCTTCTTCGGAGTCGGTGTTAATCAGAATATCAGCCTGCAACCAGTTGCCCTGTAAGCCg                                                                     <  1:1295454/68‑1 (MQ=255)
tttCACCTTCTTCTTCGGAGTCGGTGTTAATCAGAATATCAGCCTGCAACCAGTTGCCCTGTAAGCCg                                                                     <  1:2912125/68‑1 (MQ=255)
     ccTTCTTCTTCGGAGTCGGTGTTAATCAGAATATCAGCCTGCAACCAGTTGCCCTGTAAGCCGAACGc                                                                <  1:661871/68‑1 (MQ=255)
      cttcttcttCGGAGTCGGTGTTAATCAGAATATCAGCCTGCAACCAGTTGCCCTGTAAGCCGAACGc                                                                <  1:1539024/67‑1 (MQ=255)
       ttcttcttcGGAGTCGGTGTTAATCAGAATATCAGCCTGCAACCAGTTGCCCTGTAAGCCGAACGc                                                                >  1:38698/1‑66 (MQ=255)
       ttcttcttcGGAGTCGGTGTTAATCAGAATATCAGCCTGCAACCAGTTGCCCTGTAAGCCGAACGCa                                                               >  1:1004836/1‑67 (MQ=255)
       ttcttcttcGGAGTCGGTGTTAATCAGAATATCAGCCTGCAACCAGTTGCCCTGTAAGCCGAACGCa                                                               >  1:871337/1‑67 (MQ=255)
       ttcttcttcGGAGTCGGTGTTAATCAGAATATCAGCCTGCAACCAGTTGCCCTGTAAGCCGAACGCa                                                               >  1:1526417/1‑67 (MQ=255)
       ttcttcttcGGAGTCGGTGTTAATCAGAATATCAGCCTGCAACCAGTTGCCCTGTAAGCCGAACGCa                                                               >  1:1391482/1‑67 (MQ=255)
       ttcttcttcGGAGTCGGTGTTAATCAGAATATCAGCCTGCAACCAGTTGCCCTGTAAGCCGAACGCa                                                               >  1:2294718/1‑67 (MQ=255)
       ttcttcttcGGAGTCGGTGTTAATCAGAATATCAGCCTGCAACCAGTTGCCCTGTAAGCCGAACGCa                                                               >  1:2688715/1‑67 (MQ=255)
       ttcttcttcGGAGTCGGTGTTAATCAGAATATCAGCCTGCAACCAGTTGCCCTGTAAGCCGAACGCa                                                               >  1:1348386/1‑67 (MQ=255)
           tcttcGGAGTCGGTGTTAATCAGAATATCAGCCTGCAACCAGTTGCCCTGTAAGCCGAACGCACCGTc                                                          <  1:1654193/68‑1 (MQ=255)
           tcttcGGAGTCGGTGTTAATCAGAATATCAGCCTGCAACCAGTTGCCCTGTAAGCCGAACGCACCGTc                                                          <  1:150064/68‑1 (MQ=255)
           tcttcGGAGTCGGTGTTAATCAGAATATCAGCCTGCAACCAGTTGCCCTGTAAGCCGAACGCACCGTc                                                          <  1:2124441/68‑1 (MQ=255)
           tcttcGGAGTCGGTGTTAATCAGAATATCAGCCTGCAACCAGTTGCCCTGTAAGCCGAACGCACCGTc                                                          <  1:2377038/68‑1 (MQ=255)
           tcttcGGAGTCGGTGTTAATCAGAATATCAGCCTGCAACCAGTTGCCCTGTAAGCCGAACGCACCGTc                                                          <  1:3326092/68‑1 (MQ=255)
           tcttcGGAGTCGGTGTTAATCAGAATATCAGCCTGCAACCAGTTGCCCTGTAAGCCGAACGCACCGTc                                                          <  1:1137689/68‑1 (MQ=255)
           tcttcGGAGTCGGTGTTAATCAGAATATCAGCCTGCAACCAGTTGCCCTGTAAGCCGAACGCACCGTc                                                          <  1:739599/68‑1 (MQ=255)
           tcttcGGAGTCGGTGTTAATCAGAATATCAGCCTGCAACCAGTTGCCCTGTAAGCCGAACGCACCGTc                                                          <  1:1094765/68‑1 (MQ=255)
                   gTCGGTGTTAATCAGAATATCAGCCTGCAACCAGTTGCTCTGTAAGCCAAAAGCACCGTCCATACCGGc                                                 >  1:859298/1‑69 (MQ=255)
                   gTCGGTGTTAATCAGAATATCAGCCTGCAACCAGTTGCTCTGTAAGCCAAAAGCACCGTCCATACCGGc                                                 >  1:1807969/1‑69 (MQ=255)
                   gTCGGTGTTAATCAGAATATCAGCCTGCAACCAGTTCTCTGTAAGCCAAAAGCACCGTCCATACCgg                                                  >  1:682618/1‑67 (MQ=255)
                          ttAATCAGGATATCAGCTTGCAACCAGTTGCCCTGTAAGCCGAACGCACCGTCCATACCGGCTTCTTCg                                          <  1:3219628/69‑1 (MQ=255)
                          ttAATCAGGATATCAGCCTGCAACCAGTTGCCCTGTAAGCCGAACGCACCGTCCATACCGGCTTCTTCg                                          <  1:2744463/69‑1 (MQ=255)
                          ttAATCAGGATATCAGCCTGCAACCAGTTGCCCTGTAAGCCGAACGCACCGTCCATACCGGCTTCTTCg                                          <  1:3045671/69‑1 (MQ=255)
                          ttAATCAGGATATCAGCCTGCAACCAGTTGCCCTGTAAGCCGAACGCACCGTCCATACCGGCTTCTTCg                                          <  1:573628/69‑1 (MQ=255)
                          ttAATCAGGATATCAGCCTGCAACCAGTTGCCCTGTAAGCCGAACGCACCGTCCATACCGGCTTCTTCg                                          <  1:1815160/69‑1 (MQ=255)
                          ttAATCAGGATATCAGCCTGCAACCAGTTGCCCTGTAAGCCGAACGCACCGTCCATACCGGCTTCTTCg                                          <  1:938018/69‑1 (MQ=255)
                          ttAATCAGGATATCAGCCTGCAACCAGTTGCCCTGTAAGCCGAACGCACCGTCCATACCGGCTTCTTCg                                          <  1:1303504/69‑1 (MQ=255)
                          ttAATCAGGATATCAGCCTGCAACCAGTTGCCCTGTAAGCCGAACGCACCGTCCATACCGGCTTCTTCg                                          <  1:1166439/69‑1 (MQ=255)
                          ttAATCAGGATATCAGCCTGCAACCAGTTGCCCTGTAAGCCGAACGCACCGTCCATACCGGCTTCTTCg                                          <  1:1543365/69‑1 (MQ=255)
                          ttAATCAGGATATCAGCCTGCAACCAGTTGCCCTGTAAGCCGAACGCACCGTCCATACCGGCTTCTACg                                          <  1:1548927/69‑1 (MQ=255)
                          ttAATCAGAATATCAGCCTGCAACCAGTTGCCCTGTAAGCCGAACGCACCGTCCATACCGGCTTCTTCg                                          >  1:1856362/1‑69 (MQ=255)
                          ttAATCAGAATATCAGCCTGCAACCAGTTGCCCTGTAAGCCGAACGCACCGTCCATACCGGCTTCTTCg                                          >  1:1270680/1‑69 (MQ=255)
                            aaTCAGGATATCAGCCTGCAACCAGTTGCCCTGTAAGCCGAACGCACCGTCCa                                                        >  1:4259/1‑53 (MQ=255)
                            aaTCAGAATATCAGCCTGCAACGAGTTGCCCTGTAAGCCGAACGCACCGTCCa                                                        >  1:2165520/1‑53 (MQ=255)
                            aaTCAGAATATCAGCCTGCAACCAGTTGCCCTGTAAGCCGACCGCACCGTCCa                                                        >  1:788314/1‑53 (MQ=255)
                            aaTCAGAATATCAGCCTGCAACCAGTTGCCCTGTAAGCCGAACGCAcc                                                             >  1:1098183/1‑48 (MQ=255)
                            aaTCAGAATATCAGCCTGCAACCAGTTGCCCTGTAAGCCGAACGCACCg                                                            >  1:908817/1‑49 (MQ=255)
                            aaTCAGAATATCAGCCTGCAACCAGTTGCCCTGTAAGCCGAACGCACCGt                                                           >  1:2358023/1‑50 (MQ=255)
                            aaTCAGAATATCAGCCTGCAACCAGTTGCCCTGTAAGCCGAACGCACCGt                                                           >  1:2024897/1‑50 (MQ=255)
                            aaTCAGAATATCAGCCTGCAACCAGTTGCCCTGTAAGCCGAACGCACCGTc                                                          >  1:3289568/1‑51 (MQ=255)
                            aaTCAGAATATCAGCCTGCAACCAGTTGCCCTGTAAGCCGAACGCACCGTc                                                          >  1:3276936/1‑51 (MQ=255)
                            aaTCAGAATATCAGCCTGCAACCAGTTGCCCTGTAAGCCGAACGCACCGTc                                                          >  1:2565480/1‑51 (MQ=255)
                            aaTCAGAATATCAGCCTGCAACCAGTTGCCCTGTAAGCCGAACGCACCGTc                                                          >  1:713385/1‑51 (MQ=255)
                            aaTCAGAATATCAGCCTGCAACCAGTTGCCCTGTAAGCCGAACGCACCGTc                                                          >  1:2931017/1‑51 (MQ=255)
                            aaTCAGAATATCAGCCTGCAACCAGTTGCCCTGTAAGCCGAACGCACCGTc                                                          >  1:1504254/1‑51 (MQ=255)
                            aaTCAGAATATCAGCCTGCAACCAGTTGCCCTGTAAGCCGAACGCACCGTc                                                          >  1:2061079/1‑51 (MQ=255)
                            aaTCAGAATATCAGCCTGCAACCAGTTGCCCTGTAAGCCGAACGCACCGTc                                                          >  1:3152705/1‑51 (MQ=255)
                            aaTCAGAATATCAGCCTGCAACCAGTTGCCCTGTAAGCCGAACGCACCGTc                                                          >  1:1970493/1‑51 (MQ=255)
                            aaTCAGAATATCAGCCTGCAACCAGTTGCCCTGTAAGCCGAACGCACCGTCCa                                                        >  1:1322524/1‑53 (MQ=255)
                            aaTCAGAATATCAGCCTGCAACCAGTTGCCCTGTAAGCCGAACGCACCGTCCa                                                        >  1:3118988/1‑53 (MQ=255)
                            aaTCAGAATATCAGCCTGCAACCAGTTGCCCTGTAAGCCGAACGCACCGTCCa                                                        >  1:2515598/1‑53 (MQ=255)
                            aaTCAGAATATCAGCCTGCAACCAGTTGCCCTGTAAGCCGAACGCACCGTCCa                                                        >  1:2679248/1‑53 (MQ=255)
                            aaTCAGAATATCAGCCTGCAACCAGTTGCCCTGTAAGCCGAACGCACCGTCCa                                                        >  1:1358281/1‑53 (MQ=255)
                            aaTCAGAATATCAGCCTGCAACCAGTTGCCCTGTAAGCCGAACGCACCGTCCATACCGGCTTCTTCGGt                                        >  1:2694453/1‑69 (MQ=255)
                            aaTCAGAATATCAGCCTGCAACCAGTTGCCCTGTAAGCAGAACGCACCGTCCa                                                        >  1:1665854/1‑53 (MQ=255)
                              tCAGGATATCAGCCTGCAACCAGTTGCCCTGTAAGCCGAACGCACCGTCCATACCGGCTTCTTCg                                          <  1:1830333/65‑1 (MQ=255)
                                        aGCCTGCAACCAGTTGCCCTGTAAGCCGAACGCACCGTCCATACCGGCTTCTTCGGTCATGGTcagcag                            >  1:449136/1‑69 (MQ=255)
                                        aGCCTGCAACCAGTTGCCCTGTAAGCCGAACGCACCGTCCATACCGGCTTCTTCGGTCATGGTcagcag                            >  1:2642801/1‑69 (MQ=255)
                                        aGCCTGCAACCAGTTGCCCTGTAAGCCGAACGCACCGTCCATACCGGCTTCTTCGGTCATGGTcagca                             >  1:1249129/1‑68 (MQ=255)
                                          ccTGCAACCAGTTGCTCTGTAAGCCAAAAGCACCg                                                            >  1:1800126/1‑35 (MQ=25)
                                          ccTGCAACCAGTTGCTCTGTAAGCCAAAAGCACCg                                                            >  1:2727189/1‑35 (MQ=25)
                                          ccTGCAACCAGTTGCTCTGTAAGCCAAAAGCACCg                                                            >  1:2572212/1‑35 (MQ=25)
                                          ccTGCAACCAGTTGCTCTGTAAGCCAAAAGCACCg                                                            >  1:3384165/1‑35 (MQ=25)
                                          ccTGCAACCAGTTGCTCTGTAAGCCAAAAGCACCg                                                            >  1:1378884/1‑35 (MQ=25)
                                          ccTGCAACCAGTTGCTCTGTAAGCCAAAAGCACCg                                                            >  1:1445314/1‑35 (MQ=25)
                                          ccTGCAACCAGTTGCTCTGTAAGCCAAAAGCACCg                                                            >  1:1926898/1‑35 (MQ=25)
                                          ccTGCAACCAGTTGCTCTGTAAGCCAAAAGCACCg                                                            >  1:10409/1‑35 (MQ=25)
                                          ccTGCAACCAGTTGCTCTGTAAGCCAAAAGCACCg                                                            >  1:88005/1‑35 (MQ=25)
                                                          cTGTAAGCCAAAAGCACCGTCCATACCGGCTTCTTCGGTCATGGTCAGCAGTACTTCCAGCGGGCCGTg          <  1:2711436/69‑1 (MQ=255)
                                                          cTGTAAGCCAAAAGCACCGTCCATACCGGCTTCTTCGGTCATGGTCAGCAGCACTTCCAGCGGGCCGTg          <  1:1115941/69‑1 (MQ=255)
                                                          cTGTAAGCCAAAAGCACCGTCCATACCGGCTTCTTCGGTCATGGTCAGCAGCACTTCCAGCGGGCCGTg          <  1:508695/69‑1 (MQ=255)
                                                          cTGTAAGCCAAAAGCACCGTCCATACCGGCTTCTTCGGTCATGGTCAGCAGCACTTCCAGCGGGCCGTg          <  1:451160/69‑1 (MQ=255)
                                                          cTGTAAGCCAAAAGCACCGTCCATACCGGCTTCTTCGGTCATGGTCAGCAGCACTTCCAGCGGGCCGTg          <  1:197814/69‑1 (MQ=255)
                                                          cTGTAAGCCAAAAGCACCGTCCATACCGGCTTCTTCGGTCATGGTCAGCAGCACTTCCAGCGGGCCGTg          <  1:1073706/69‑1 (MQ=255)
                                                          cTGTAAGCCAAAAGCACCGTCCATACCGGCTTCTTCGGTCATGGTCAGCAGCACTTCCAGCGGGCCGTg          <  1:1360808/69‑1 (MQ=255)
                                                           tGTAAGCCAAAAGCACCGTCCATACCGGCTTCTTCGGTCATGGTCAGCAGCACTTCCAGCGGGCCGTg          <  1:2986785/68‑1 (MQ=255)
                                                           tGTAAGCCAAAAGCACCGTCCATACCGGCTTCTTCGGTCATGGTCAGCAGCACTTCCAGCGGGCCGTg          <  1:184149/68‑1 (MQ=255)
                                                             tAAGCCGAACGCACCGTCCATACCGGCTTCTTCGGTCATGGTCAGCAGCACTTCCAGCGGGCCGTGAAc       >  1:1018518/1‑69 (MQ=255)
                                                             tAAGCCGAACGCACCGTCCATACCGGCTTCTTCGGTCATGGTCAGCAGCACTTCCAGCGGGCCGTGAAc       >  1:2176721/1‑69 (MQ=255)
                                                             tAAGCCGAACGCACCGTCCATACCGGCTTCTTCGGTCATGGTCAGCAGCACTTCCAGCGGGCCGTGAAc       >  1:124576/1‑69 (MQ=255)
                                                             tAAGCCGAACGCACCGTCCATACCGGCTTCTTCGGTCATGGTCAGCAGCACTTCCAGCGGGCCGTGAAc       >  1:2476439/1‑69 (MQ=255)
                                                             tAAGCCGAACGCACCGTCCATACCGGCTTCTTCGGTCATGGTCAGCAGCACTTCCAGCGGGCCGTGAAc       >  1:2815213/1‑69 (MQ=255)
                                                             tAAGCCGAACGCACCGTCCATACCGGCTTCTTCGGTCATGGTCAGCAGCACTTCCAGCGGGCCGTGAAc       >  1:3193324/1‑69 (MQ=255)
                                                             tAAGCCGAACGCACCGTCCATACCGGCTTCTTCGGTCATGGTCAGCAGCACTTCCAGCGGGCCGTGAAc       >  1:1317149/1‑69 (MQ=255)
                                                                  cGAACGCACCGTCCATACCGGCTTCTTCGGTCATGGTcag                               >  1:2701149/1‑40 (MQ=255)
                                                                  cGAACGCACCGTCCATACCGGCTTCTTCGGTCATGGTCAgcagca                          >  1:1706019/1‑45 (MQ=255)
                                                                  cGAACGCACCGTCCATACCGGCTTCTTCGGTCATGGTCAgcagca                          >  1:1699401/1‑45 (MQ=255)
                                                                  cGAACGCACCGTCCATACCGGCTTCTTCGGTCATGGTCAgcagca                          >  1:2319922/1‑45 (MQ=255)
                                                                  cGAACGCACCGTCCATACCGGCTTCTTCGGTCATGGTCAgcagca                          <  1:1098095/45‑1 (MQ=255)
                                                                  cGAACGCACCGTCCATACCGGCTTCTTCGGTCATGGTCAgcagca                          >  1:1988907/1‑45 (MQ=255)
                                                                  cGAACGCACCGTCCATACCGGCTTCTTCGGTCATGGTCAgcagca                          >  1:2110755/1‑45 (MQ=255)
                                                                  cGAACGCACCGTCCATACCGGCTTCTTCGGTCATGGTCAgcagca                          >  1:2360534/1‑45 (MQ=255)
                                                                  cGAACGCACCGTCCATACCGGCTTCTTCGGTCATGGTCAgcagca                          <  1:3343151/45‑1 (MQ=255)
                                                                  cGAACGCACCGTCCATACCGGCTTCTTCGGTCATGGTCAgcagca                          >  1:2901023/1‑45 (MQ=255)
                                                                  cGAACGCACCGTCCATACCGGCTTCTTCGGTCATGGTCAgcagca                          >  1:3115906/1‑45 (MQ=255)
                                                                  cGAACGCACCGTCCATACCGGCTTCTTCGGTCATGGTCAGCAGCACTTCCAGCGGGCCGTGAACCACGt  >  1:1577625/1‑69 (MQ=255)
                                                                  cGAACCCACCGTCCATACCGGCTTCTTCGGTCATGGTaagc                              >  1:386396/1‑41 (MQ=255)
                                                                   |                                                                   
TTTCACCTTCTTCTTCGGAGTCGGTGTTAATCAGAATATCAGCCTGCAACCAGTTGCCCTGTAAGCCGAACGCACCGTCCATACCGGCTTCTTCGGTCATGGTCAGCAGCACTTCCAGCGGGCCGTGAACCACGT  >  minE/224569‑224703

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: