Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,906,806 G→A 43.6% A887V (GCA→GTA)  recC ← exonuclease V (RecBCD complex), gamma chain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,906,8060GA43.6% 14.3 / 39.8 55A887V (GCA→GTA) recCexonuclease V (RecBCD complex), gamma chain
Reads supporting (aligned to +/- strand):  ref base G (15/16);  new base A (15/9);  total (30/25)
Fisher's exact test for biased strand distribution p-value = 4.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.37e-01

GTAAATCCCCTGCCGCTCGGAAGCGGCGGAACAAGCGTTCGGCATCATCCTGCTCAACCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTCGGG  >  minE/1906744‑1906873
                                                              |                                                                   
gTAAATCCCCTGCCGCTCGGAAGCGGCGGAACAAGCGTTCGGCATCATCCTGCTCAACCAGTGCATTc                                                                <  1:35097/68‑1 (MQ=255)
  aaaTCCCCTGCCGCTCGGAAGCGGCGGAACAAGCGTTCGGCATCATCCTGCTCAACCAGTGCATTCAAt                                                             <  1:1502911/69‑1 (MQ=255)
  aaaTCCCCTGCCGCTCGGAAGCGGCGGAACAAGCGTTCGGCATCATCCTGCTCAACCAGTGCATTCAAt                                                             <  1:60132/69‑1 (MQ=255)
                 cGGAAGCGGCGGAACAAGCGTTCGGCATCATCCTGCTCAACCAGTGCATTCAATAACTGCtgattgat                                               >  1:1418071/1‑68 (MQ=255)
                   gAATCGGCGGAACAAGCGTTCGGCATCATCCTGCTCAACCAGTACATTTAATAACTGCTGGttgattt                                             >  1:1691023/1‑68 (MQ=255)
                   gAATCGGCGGAACAAGCGTTCGGCATCATCCTGCTCAACCAGTACATTTAATAACTGCTGGttgattt                                             >  1:2528150/1‑68 (MQ=255)
                   gAATCGGCGGAACAAGCGTTCGGCATCATCCTGCTCAACCAGTACATTTAATAACTGCTGGttgattt                                             >  1:2490447/1‑68 (MQ=255)
                   gAATCGGCGGAACAAGCGTTCGGCATCATCCTGCTCAACCAGTACATTTAATAACTGCTGGttgattt                                             >  1:3205987/1‑68 (MQ=255)
                   gAATCGGCGGAACAAGCGTTCGGCATCATCCTGCTCAACCAGTACATTTAATAACTGCTGGttgattt                                             >  1:2097579/1‑68 (MQ=255)
                   gAATCGGCGGAACAAGCGTTCGGCATCATCCTGCTCAACCAGTACATTTAATAACTGCTGGttgattt                                             >  1:1963765/1‑68 (MQ=255)
                   gAATCGGCGGAACAAGCGTTCGGCATCATCCTGCTCAACCAGTACATTTAATAACTGCTGGttgattt                                             >  1:526575/1‑68 (MQ=255)
                   gAATCGGCGGAACAAGCGTTCGGCATCATCCTGCTCAACCAGTACATTTAATAACTGCTGGttgattt                                             >  1:680342/1‑68 (MQ=255)
                   gAATCGGCGGAACAAGCGTTCGGCATCATCCTGCTCAACCAGTACATTTAATAACTGCTGGttgattt                                             >  1:1566964/1‑68 (MQ=255)
                   gAATCGGCGGAACAAGCGTTCGGCATCATCCTGCTCAACCAGTACATTTAATAACTGCTGGttgattt                                             >  1:1511111/1‑68 (MQ=255)
                   gAATCGGCGGAACAAGCGTTCGGCATCATCCTGCTCAACCAGTACATTTAATAACTGCTGGttgattt                                             >  1:812086/1‑68 (MQ=255)
                   gAATCGGCGGAACAAGCGTTCGGCATCATCCTGCTCAACCAGTACATTTAATAACTGCTGGttgattt                                             >  1:89792/1‑68 (MQ=255)
                   gAATCGGCGGAACAAGCGTTCGGCATCATCCTGCTCAACCAGTACATTTAATAACTGCTGGttgattt                                             >  1:1384200/1‑68 (MQ=255)
                   gAATCGGCGGAACAAGCGTTCGGCATCATCCTGCTCAACCAGTACATTTAATAACTGCTGGttgattt                                             >  1:943275/1‑68 (MQ=255)
                   gAATCGGCGGAACAAGCGTTCGGCATCATCCTGCTCAACCAGTACATTTAATAACTGCTGGttgattt                                             >  1:943339/1‑68 (MQ=255)
                               cAAGCGTTCGGCATCATCCTGCTCAACCAGTGCATTCAAt                                                             >  1:703378/1‑40 (MQ=255)
                                            tcatcCTGCTCAACCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATa                   >  1:1557466/1‑69 (MQ=255)
                                             catcCTGCTCAACCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATaa                  >  1:2368466/1‑69 (MQ=255)
                                              atcCTGCTCAACCAGTGCATTCAATAAGTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATaaa                 >  1:3311554/1‑69 (MQ=255)
                                              atcCTGCTCAACCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAg                                <  1:2557052/54‑1 (MQ=255)
                                              atcCTGCTCAACCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAg                                <  1:2895723/54‑1 (MQ=255)
                                              atcCTGCTCAACCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAg                                <  1:2593536/54‑1 (MQ=255)
                                              atcCTGCTCAACCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAg                                <  1:2525392/54‑1 (MQ=255)
                                              atcCTGCTCAACCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAg                                <  1:2377788/54‑1 (MQ=255)
                                              atcCTGCTCAACCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAg                                <  1:147201/54‑1 (MQ=255)
                                              atcCTGCTCAACCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAg                                <  1:2275628/54‑1 (MQ=255)
                                              atcCTGCTCAACCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATaaa                 >  1:2536018/1‑69 (MQ=255)
                                              atcCTGCTCAACCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATaaa                 >  1:1254282/1‑69 (MQ=255)
                                              atcCTGCTCAACCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATaaa                 >  1:2325553/1‑69 (MQ=255)
                                              atcCTGCTCAACCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATaaa                 >  1:1460985/1‑69 (MQ=255)
                                              atcCTGCTCAACCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATaaa                 >  1:516243/1‑69 (MQ=255)
                                              atcCTGCTCAACCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATaaa                 >  1:460544/1‑69 (MQ=255)
                                              atcCTGCTCAACCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATaaa                 >  1:1851003/1‑69 (MQ=255)
                                              atcCTGCTCAACCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATaaa                 >  1:36514/1‑69 (MQ=255)
                                              atcCTGCTCAACCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATaaa                 >  1:415870/1‑69 (MQ=255)
                                              atcCTGCTCAACCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAGTCCATCCAGAATaaa                 >  1:342458/1‑69 (MQ=255)
                                                  tGCTCAACCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAg                                <  1:865275/50‑1 (MQ=255)
                                                           aGTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAg                                <  1:3278616/41‑1 (MQ=255)
                                                           aGTACATTTAATAACTGCTGGTTGATTTGATAACGACTAAGCCCTTCAAGAATAAAAGGTTCGGTGTCg    <  1:1037745/69‑1 (MQ=255)
                                                           aGTACATTTAATAACTGCTGGTTGATTTGATAACGACTAAGCCCTTCAAGAATAAAAGGTTCGGTGTCg    <  1:881352/69‑1 (MQ=255)
                                                           aGTACATTTAATAACTGCTGGTTGATTTGATAACGACTAAGCCCTTCAAGAATAAAAGGTTCGGTGTCg    <  1:738841/69‑1 (MQ=255)
                                                           aGTACATTTAATAACTGCTGGTTGATTTGATAACGACTAAGCCCTTCAAGAATAAAAGGTTCGGTGTCg    <  1:1675536/69‑1 (MQ=255)
                                                           aGTACATTTAATAACTGCTGGTTGATTTGATAACGACTAAGCCCTTCAAGAATAAAAGGTTCGGTGTCg    <  1:673624/69‑1 (MQ=255)
                                                           aGTACATTTAATAACTGCTGGTTGATTTGATAACGACTAAGCCCTTCAAGAATAAAAGGTTCGGTGTCg    <  1:1815125/69‑1 (MQ=255)
                                                           aGTACATTTAATAACTGCTGGTTGATTTGATAACGACTAAGCCCTTCAAGAATAAAAGGTTCGGTGTCg    <  1:599713/69‑1 (MQ=255)
                                                           aGTACATTTAATAACTGCTGGTTGATTTGATAACGACTAAGCCCTTCAAGAATAAAAGGTTCGGTGTCg    <  1:1846049/69‑1 (MQ=255)
                                                           aGTACATTTAATAACTGCTGGTTGATTTGATAACGACTAAGCCCTTCAAGAATAAAAGGTTCGGTGTCg    <  1:2891039/69‑1 (MQ=255)
                                                             tGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTCggg  <  1:3363910/69‑1 (MQ=255)
                                                             tGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTCggg  <  1:3297055/69‑1 (MQ=255)
                                                             tGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTCggg  <  1:1411334/69‑1 (MQ=255)
                                                             tGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAGTCCTTACAGAATAAATGGCTCGGTGTCggg  <  1:1069896/69‑1 (MQ=255)
                                                              |                                                                   
GTAAATCCCCTGCCGCTCGGAAGCGGCGGAACAAGCGTTCGGCATCATCCTGCTCAACCAGTGCATTCAATAACTGCTGATTGATTTGATAACGGCTAAGTCCTTCCAGAATAAATGGCTCGGTGTCGGG  >  minE/1906744‑1906873

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