Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 518,289 G→A 18.6% D105N (GAT→AAT)  bioC → predicted methltransferase, enzyme of biotin synthesis

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE518,2890GA18.6% 159.9 / 28.3 70D105N (GAT→AAT) bioCpredicted methltransferase, enzyme of biotin synthesis
Reads supporting (aligned to +/- strand):  ref base G (34/23);  new base A (6/7);  total (40/30)
Fisher's exact test for biased strand distribution p-value = 5.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

CAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAGCTGTATCG  >  minE/518220‑518359
                                                                     |                                                                      
cAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGa                                                                       <  1:3230205/71‑1 (MQ=255)
  gAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGAt                                                                      >  1:2408037/1‑70 (MQ=255)
  gAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGAt                                                                      >  1:1423549/1‑70 (MQ=255)
    aGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCtt                                                                   >  1:4138256/1‑71 (MQ=255)
      gATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCAATCTTg                                                                  <  1:2389395/70‑1 (MQ=255)
           cGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTA‑‑GACTGCGACGTTCGATCTTGCATGGAg                                                           >  1:1944329/1‑70 (MQ=255)
             cAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGc                                                          <  1:561804/71‑1 (MQ=255)
             cAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGc                                                          <  1:4014692/71‑1 (MQ=255)
             cAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCAATCTTGc                                                                 <  1:2662252/64‑1 (MQ=255)
              aggccATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCa                                                         <  1:451400/68‑1 (MQ=255)
              aGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCa                                                         <  1:1449823/71‑1 (MQ=255)
                aCCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCCACTGCGACGTTCg                                                                        <  1:1699943/54‑1 (MQ=255)
                aCCATTATCTGGCGGGAGATATAGAATCCCTGCCGTTAGCGACTGCGACGTTCGATc                                                                     >  1:815493/1‑57 (MQ=255)
                 ccATAATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAAtc                                                      >  1:2811345/1‑71 (MQ=255)
                  cATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCAATCTTGCa                                                                >  1:3533218/1‑60 (MQ=255)
                    ttATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATg                                                              <  1:1141432/60‑1 (MQ=255)
                    ttATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAAtct                                                     >  1:1536406/1‑69 (MQ=255)
                    ttATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCg                                                   <  1:2634982/71‑1 (MQ=255)
                    ttATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCg                                                   >  1:2007936/1‑71 (MQ=255)
                    ttATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCg                                                   >  1:345076/1‑71 (MQ=255)
                    ttATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCAATCTTGCATGGAGCAATCTCg                                                   <  1:69013/71‑1 (MQ=255)
                          ggCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCg                                                                        >  1:367088/1‑44 (MQ=255)
                           gCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGc                                                                 >  1:1968469/1‑50 (MQ=255)
                           gCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCgcagtg                                              <  1:3434001/69‑1 (MQ=255)
                             gggAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCgcagtgcagt                                          <  1:455413/71‑1 (MQ=255)
                              ggAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCgca                                                 >  1:492385/1‑63 (MQ=255)
                                 gaTATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGAt                                                                      >  1:2225484/1‑39 (MQ=255)
                                  atatCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCgcagtgcag                                           <  1:2153417/65‑1 (MQ=255)
                                  atatCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCgcagtgcag                                           <  1:2269026/65‑1 (MQ=255)
                                    atCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCa                                                                >  1:2889542/1‑42 (MQ=255)
                                    atCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCa                                                                >  1:2743168/1‑42 (MQ=255)
                                    atCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCa                                                                >  1:1313528/1‑42 (MQ=255)
                                    atCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCgg                                   >  1:1612664/1‑71 (MQ=255)
                                     tCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAAt                                                       >  1:1079925/1‑50 (MQ=255)
                                     tCGAATCCCTGCCGTTAGCGACTGCGACGTTCAATCTTGCATGGAGCAATCTCgcagtgcag                                           <  1:2022379/62‑1 (MQ=255)
                                       gAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGa                                                            >  1:742500/1‑43 (MQ=255)
                                       gAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGa                                                            >  1:4130244/1‑43 (MQ=255)
                                       gAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCg                                                   <  1:3956721/52‑1 (MQ=255)
                                       gAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCg                                                   <  1:676633/52‑1 (MQ=255)
                                       gAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTaa                                >  1:2192894/1‑71 (MQ=255)
                                       gAATCCCTGCCGTTAGCGACTGCGACGTTCAATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTaa                                >  1:1573534/1‑71 (MQ=255)
                                         aTCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAAtt                              <  1:3709868/71‑1 (MQ=255)
                                         aTCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAAtt                              >  1:2878295/1‑71 (MQ=255)
                                          tCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAAttt                             <  1:3253522/71‑1 (MQ=255)
                                          tCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAAtt                              >  1:3385211/1‑70 (MQ=255)
                                          tCCCTGCCGTTAGCGACTGCGACGTTCGATCTTCATGGa                                                            >  1:1576807/1‑39 (MQ=39)
                                              tGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAAttt                             <  1:1907932/67‑1 (MQ=255)
                                               gCCGTTAGCGACTGCGACGTTCGATCTTGCATGGGGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCa                        >  1:3357977/1‑71 (MQ=255)
                                               gCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCa                        >  1:3271016/1‑71 (MQ=255)
                                               gCCGTTAGCGACTGCGACGTTCAATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCa                        >  1:1263228/1‑71 (MQ=255)
                                                ccGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAAtt                              >  1:1128314/1‑64 (MQ=255)
                                                 cGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTa                                 >  1:3620925/1‑60 (MQ=255)
                                                   ttAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCgcagtg                                              <  1:1594129/45‑1 (MQ=255)
                                                   ttAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACgg                     >  1:338412/1‑70 (MQ=255)
                                                   ttAGCGACTGCGACGTTCAATCTTGCATTGAGCAATCTCgca                                                 >  1:2055005/1‑42 (MQ=255)
                                                   ttAGCGACTGCGACGTTCAATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATc                          >  1:122175/1‑65 (MQ=255)
                                                   ttAGCGACTGCGACGTTCAATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACgg                     >  1:1589864/1‑70 (MQ=255)
                                                     aGCGACTGCGACGTTCGATCTTGCATGGAGCAAtct                                                     >  1:2385726/1‑36 (MQ=255)
                                                     aGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCAc                  >  1:389640/1‑71 (MQ=255)
                                                     aGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCAc                  >  1:3125618/1‑71 (MQ=255)
                                                       cGACTGCGACGTTCAATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTc                <  1:248502/71‑1 (MQ=255)
                                                       cGACTGCGACGTTCAATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTc                <  1:2560008/71‑1 (MQ=255)
                                                         aCTGCGACGTTCGATCTTGCATGGAGCAATCTCgcag                                                >  1:273611/1‑37 (MQ=255)
                                                           tGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATc                          <  1:3699364/57‑1 (MQ=255)
                                                            gCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTaa                                >  1:86602/1‑50 (MQ=255)
                                                              gACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAAttt                             <  1:2738943/51‑1 (MQ=255)
                                                              gACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATc                          <  1:3555310/54‑1 (MQ=255)
                                                                cGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTAt                           >  1:847082/1‑51 (MQ=255)
                                                                 gTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAGCTGt      <  1:634729/71‑1 (MQ=255)
                                                                  ttCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTa                            <  1:3837582/48‑1 (MQ=255)
                                                                  ttCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCTGCGAGCTGTa     <  1:1089957/71‑1 (MQ=255)
                                                                  ttCAATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATc                          <  1:2770073/50‑1 (MQ=255)
                                                                     gATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAGCTGTATCg  >  1:546575/1‑71 (MQ=255)
                                                                     |                                                                      
CAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGGTAATTTATCCACGGCACTCCGCGAGCTGTATCG  >  minE/518220‑518359

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: