Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 684,711 G→A 18.7% D23N (GAC→AAC)  yccU → predicted CoA‑binding protein with NAD(P)‑binding Rossmann‑fold domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE684,7110GA18.7% 66.6 / 9.8 32D23N (GAC→AAC) yccUpredicted CoA‑binding protein with NAD(P)‑binding Rossmann‑fold domain
Reads supporting (aligned to +/- strand):  ref base G (14/12);  new base A (3/3);  total (17/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

ATAATGAAAGAAACCGATATTGCAGGCATTTTAACGTCCACTCACACCATTGCGCTGGTGGGGGCAAGCGACAAACCCGATCGTCCCAGCTATCGCGTGATGAAATATCTACTCGACCAGGGGTATCACGTCATTCCG  >  minE/684642‑684779
                                                                     |                                                                    
ataatGAAAGAAACCGATATTGCAGGCATTTTAACGTCCACTCACACCATTGCGCTGGTGGGGGCAAGCGa                                                                     >  1:397991/1‑71 (MQ=255)
            aCCGATATTGCAGGCATTTTAACGTCCACTCACACCATTGCGCTGGTGGGGGCAAGCGa                                                                     >  1:425407/1‑59 (MQ=255)
               gATATTGCAGGCATTTTAACGTCCACTCACACCATTGCGCTGGTGGGGGCAAGCGACAAACCCGATCGTcc                                                      <  1:88174/71‑1 (MQ=255)
               gATATTGCAGGCATTTTAACGTCCACTCACACCATTGCGCTGGTGGGGGCAAGCAACAAACCCGATCGTcc                                                      <  1:3034960/71‑1 (MQ=255)
                    tGCAGGCATTTTAACGTCCACTCACACCATTGCGCTGGTGGGGGCAAGCGa                                                                     >  1:454807/1‑51 (MQ=255)
                     gCAGGCATTTTAACGTCCACTCACACCATTGCGCTGGTGGGGGCAAGCGACa                                                                   <  1:3495377/52‑1 (MQ=255)
                         gCATTTTAACGTCCACTCACACCATTGCGCTGGTGGGGGCAAGCAACAAACCCGATCGTccc                                                     >  1:864848/1‑62 (MQ=255)
                             tttAACGTCCACTCACACCATTGCGCTGGTGGGGGCAAGCAACaaa                                                                 >  1:420787/1‑46 (MQ=255)
                                 aCGTCCACTCACACCATTGCGCTGGTGGGGGCAAGCGACAAACCCGATCGTCCCAGCTATCGCGTGATGaa                                    <  1:224133/71‑1 (MQ=255)
                                   gTCCACTCACACCATTGCGCTGGTGGGGGCAAGCGACAAACCCGATCGTCCCAGCTATCGCGTGATGAAat                                  <  1:2899892/71‑1 (MQ=255)
                                    tCCACTCACACCATTGCGCTGGTGGGGGCAAGCGACAAACCCGATCGTCCCAGCTATCGCGTGATGAAata                                 <  1:3774997/71‑1 (MQ=255)
                                          cacaCCATTGCGCTGGTGGGGGCAAGCGACAAACCCGATCGTCCCAGCTATcgcg                                           >  1:2397122/1‑55 (MQ=255)
                                          cacaCCATTGCGCTGGTGGGGGCAAGCAACAAACCCGATCGTCCCAGCTATCGCGTGATGAAATATCTACt                           >  1:945453/1‑71 (MQ=255)
                                           acacCATTGCGCTGGTGGGGGCAAGCGACAAACCCGAt                                                           >  1:2508580/1‑38 (MQ=255)
                                              ccATTGCGCTGGTGGGGGCAAGCAACAAACCCGATCGTCCCAGCTATCGCGTGATGAAATATCTACTCGAc                       <  1:581241/71‑1 (MQ=255)
                                                aTTGCGCTGGTGGGGGCAAGCGACAAACCCGATCGTCCCAGCTATCGCGTGATGAAatat                                <  1:1925019/60‑1 (MQ=255)
                                                aTTGCGCTGGTGGGGGCAAGCAACAAACCCGATCGTCCCAGCTATCGCGTGATGAAATATCTACTCGACCa                     <  1:985529/71‑1 (MQ=255)
                                                  tGCGCTGGTGGGGGCAAGCGACAAACCCGATCGTCCCAGCTATCGCGTGATGAAATATCt                              >  1:1699369/1‑60 (MQ=255)
                                                  tGCGCTGGTGGGGGCAAGCGACAAACCCGATCGTCCCAGCTATCGCGTGATGAAATATCt                              >  1:1058337/1‑60 (MQ=255)
                                                      cTGGTGGGGGCAAGCGACAAACCCGATCGTCCCAGCTATCGCGTGATGAAATATCTACTCGACCa                     <  1:1824242/65‑1 (MQ=255)
                                                      cTGGTGGGGGCAAGCGACAAACCCGATCGTCCCAGCTATCGCGTGATGAAATATCTACTCGACCAGGGGTa               <  1:324572/71‑1 (MQ=255)
                                                             gggCAAGCGACAAACCCGATCGTCCCAGCTATCGCGTGATGAAATATCTACTCGACCAGGGGTATCACg          >  1:2476250/1‑69 (MQ=255)
                                                             gggCAAGCGACAAACCCGATCGTCCCAGCTATCGCGTGATGAAATATCTACTCGACCAGGGGTATCACGt         <  1:595595/70‑1 (MQ=255)
                                                               gcaagcaaCAAACCCGATCGTCCCAGCTATCGCGtgatg                                      <  1:1050319/39‑1 (MQ=255)
                                                                cAAGCGACAAACCCGATCGTCCCAGCTATCGCGTGATGAAATATCTACTCGACCAGGGGTATCACg          <  1:1101572/66‑1 (MQ=255)
                                                                cAAGCGACAAACCCGATCGTCCCAGCTATCGCGTGATGAAATATCTACTCGACCAGGGGTATCACGTCAtt     <  1:263143/71‑1 (MQ=255)
                                                                 aaGCGACAAACCCGATCGTCCCAGCTATCGCGtgatga                                     <  1:572550/38‑1 (MQ=255)
                                                                   gCGACAAACCCGATCGTCCCAGCTATCGCGTGATGaaa                                   >  1:2087527/1‑38 (MQ=255)
                                                                   gCGACAAACCCGATCGTCCCAGCTATCGCGTGATGaaa                                   >  1:567741/1‑38 (MQ=255)
                                                                   gCGACAAACCCGATCGTCCCAGCTATCGCGTGATGaaa                                   >  1:585505/1‑38 (MQ=255)
                                                                   gCGACAAACCCGATCGTCCCAGCTATCGCGTGATGAAATATCTACTCGACCAGGGGTATCACGTCATTCCg  >  1:3608129/1‑71 (MQ=255)
                                                                   gCGACAAACCCGATCGTCCCAGCTATCGCGTGATGAAATATCTACTCGACCAGGGGTATCACGTCATTCCg  >  1:839735/1‑71 (MQ=255)
                                                                   gCGACAAACCCGATCGTCCCAGCTATCGCGTGATGAAATATCTACTCGACCAGGGGTATCACGTCATTCCg  >  1:1418023/1‑71 (MQ=255)
                                                                     |                                                                    
ATAATGAAAGAAACCGATATTGCAGGCATTTTAACGTCCACTCACACCATTGCGCTGGTGGGGGCAAGCGACAAACCCGATCGTCCCAGCTATCGCGTGATGAAATATCTACTCGACCAGGGGTATCACGTCATTCCG  >  minE/684642‑684779

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: