Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,969,927 G→A 17.5% R161C (CGC→TGC)  lplA ← lipoate‑protein ligase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,969,9270GA17.5% 90.0 / 11.3 40R161C (CGC→TGC) lplAlipoate‑protein ligase A
Reads supporting (aligned to +/- strand):  ref base G (9/24);  new base A (3/4);  total (12/28)
Fisher's exact test for biased strand distribution p-value = 4.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

CGGGAACGTACCGACGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCGGCATTGAGTAGCAAGGTGCCGTGGTGGAAGCCGCGATCTTTGGTTTCGCGATAGGCCGAG  >  minE/2969857‑2969997
                                                                      |                                                                      
cgGGAACGTACCGACGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAggcg                                                                        >  1:175149/1‑71 (MQ=255)
cgGGAACGTACCGACGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAggcg                                                                        >  1:2407291/1‑71 (MQ=255)
        tACCGACGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGt                                                                >  1:3130741/1‑71 (MQ=255)
              cGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCGGCAt                                                          <  1:995092/71‑1 (MQ=255)
              cGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCGGCAt                                                          <  1:598449/71‑1 (MQ=255)
              cGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCGGCAt                                                          <  1:3162554/71‑1 (MQ=255)
                      cTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCGGCATTGAGTAGc                                                  <  1:4222307/71‑1 (MQ=255)
                         tCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCGGCATTGAGTAGCaa                                                <  1:4568464/70‑1 (MQ=255)
                         tCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCGGCATTGAGTAGCAAg                                               >  1:3019297/1‑71 (MQ=255)
                         tCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCGGCATTGAGTAGCAAg                                               >  1:3749224/1‑71 (MQ=255)
                         tCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCGGCATTGAGTAGCAAg                                               >  1:2408622/1‑71 (MQ=255)
                         tCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCGGCATTGAGTAGCAAg                                               >  1:1628663/1‑71 (MQ=255)
                         tCGCCGCCAATTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCGGCATTGAGTAGCaa                                                <  1:5024958/70‑1 (MQ=255)
                                   tttCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCGGCATTGAGTAGc                                                  <  1:4617311/58‑1 (MQ=255)
                                   tttCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCGGCATTGAGTAGc                                                  <  1:362834/58‑1 (MQ=255)
                                   tttCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCGGCATTGAGTAGc                                                  <  1:2829803/58‑1 (MQ=255)
                                   tttCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCGGCATTGAGTAGc                                                  <  1:3058542/58‑1 (MQ=255)
                                     tCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCGGCATTGAGTAGCAAGGTGCCgtgg                                      >  1:1390307/1‑68 (MQ=255)
                                       tttttATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCGGc                                                            <  1:2837458/44‑1 (MQ=255)
                                       tttttATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCGGc                                                            <  1:1617240/44‑1 (MQ=255)
                                       tttttATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCGGc                                                            <  1:2026354/44‑1 (MQ=255)
                                       tttttATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCGGc                                                            <  1:2597386/44‑1 (MQ=255)
                                       tttttATCCGGATTGAGATAGTTTGCCAGGCAGCTGAGGTCGGCATTGAGTAGCAAgg                                              <  1:128278/58‑1 (MQ=255)
                                       tttttATCCGGATTGAGATAGTTTGCCAGGCAGCTGAGGTCGGCATTGAGTAGCAAgg                                              <  1:930088/58‑1 (MQ=255)
                                       tttttATCCGGATTGAGATAGTTTGCCAGGCAGCTGAGGTCGGCATTGAGTAGCAAgg                                              <  1:2644564/58‑1 (MQ=255)
                                       tttttATCCGGATTGAGATAGTTTGCCAGGCAGCTGAGGTCGGCATTGAGTAGCAAgg                                              <  1:3166303/58‑1 (MQ=255)
                                         tttATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCGGCATTGAGTAGCAAGGTGCCGTGGTGGa                                  <  1:491241/68‑1 (MQ=255)
                                         tttATCCGGATTGAGATAGTTTGCCAGGCAGCTGAGGTCGGCATTGAGTAGCAAGGTg                                            >  1:749839/1‑58 (MQ=255)
                                         tttATCCGGATTGAGATAGTTTGCCAGGCAGCTGAGGTCGGCATTGAGTAGCAAGGTg                                            >  1:4039071/1‑58 (MQ=255)
                                         tttATCCGGATTGAGATAGTTTGCCAGGCAGCTGAGGTCGGCATTGAGTAGCAAGGTg                                            >  1:2250038/1‑58 (MQ=255)
                                                     gagaTAGTTTGCCAGGCGGCTGAGGTCGGCATTGAGTAGCAAGGTGCCGTGGTGGaa                                 >  1:3925752/1‑57 (MQ=255)
                                                     gagaTAGTTTGCCAGGCGGCTGAGGTCGGCATTGAGTAGCAAGGTGCCGTGGTGGaa                                 >  1:4410173/1‑57 (MQ=255)
                                                           gTTTGCCAGGCGGCTGAGGTCGGCATTGAGTAGCAAGGTGCCgtg                                       >  1:3276183/1‑45 (MQ=255)
                                                            tttGCCAGGCGGCTGAGGTCGGCATTGAGTAGCAAGGTGCCGTGGTGGaa                                 <  1:2850491/50‑1 (MQ=255)
                                                            tttGCCAGGCGGCTGAGGTCGGCATTGAGTAGCAAGGTGCCGTGGTGGaa                                 <  1:4612019/50‑1 (MQ=255)
                                                            tttGCCAGGCGGCTGAGGTCGGCATTGAGTAGCAAGGTGCCGTGGTGGaa                                 <  1:3695728/50‑1 (MQ=255)
                                                             ttGCCAGGCGGCTGAGGTCGGCATTGAGTAGCAAGGTGCCGTGGTGGAAGCCGCGATCTTTGGTTTc               <  1:3179365/67‑1 (MQ=255)
                                                             ttGCCAGGCGGCTGAGGTCGGCATTGAGTAGCAAGGTGCCGTGGTGGAAGCCGCGATCTTTGGTTTc               <  1:47077/67‑1 (MQ=255)
                                                             ttGCCAGGCGGCTGAGGTCGGCATTGAGTAGCAAGGTGCCGTGGTGGAAGCCGCGATCTTTGGTTTc               <  1:2532567/67‑1 (MQ=255)
                                                              tGCCAGGCGGCTGAGGTCGGCATTGAGTAGCAAGGTGCCGTGGTGGAAGCCGCGATCTTTGGTTTCGCGa           <  1:5190545/70‑1 (MQ=255)
                                                              tGCCAGGCGGCTGAGGTCGGCATTGAGTAGCAAGGTGCCGTGGTGGAAGCCGCGATCTTTGGTTTCGCGa           <  1:1458459/70‑1 (MQ=255)
                                                              tGCCAGGCGGCTGAGGTCGGCATTGAGTAGCAAGGTGCCGTGGTGGAAGCCGCGATCTTTGGTTTCGCGa           <  1:3308629/70‑1 (MQ=255)
                                                                   ggcggcTGAGGTCGGCATTGAGTAGCAAGGTGCCGTGGTGGAAGCCGCGATCTTTGGTTTCGCGATAGGcc     >  1:2521496/1‑71 (MQ=255)
                                                                    gcggcTGAGGTCGGCATTGAGTAGCAAGGTGCCGTGGTGGAAGCCGCGATCTTTGGTTTCGCGATAGGcc     >  1:4827366/1‑70 (MQ=255)
                                                                    gcggcTGAGGTCGGCATTGAGTAGCAAGGTGCCGTGGTGGAAGCCGCGATCTTTGGTTTCGCGATAGGcc     >  1:4953448/1‑70 (MQ=255)
                                                                    gcggcTGAGGTCGGCATTGAGTAGCAAGGTGCCGTGGTGGAAGCCGCGATCTTTGGTTTCGCGATAGGcc     >  1:5191612/1‑70 (MQ=255)
                                                                     cggcTGAGGTCGGCATTGAGTAGCAAGGTGCCGTGGTGGAAGCCGCGATCTTTGGTTTCGCGATAGGc      >  1:4214044/1‑68 (MQ=255)
                                                                      ggcTGAGGTCGGCATTGAGTAGCAAGGTGCCGTGGTGGAAGCCGCGATCTTTGGTTTCGCGATAGGCCGAg  >  1:3959687/1‑71 (MQ=255)
                                                                      |                                                                      
CGGGAACGTACCGACGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCGGCATTGAGTAGCAAGGTGCCGTGGTGGAAGCCGCGATCTTTGGTTTCGCGATAGGCCGAG  >  minE/2969857‑2969997

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: