Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 718,576 G→A 46.5% intergenic (+19/‑91) mviM → / → mviN predicted oxidoreductase, NAD(P)‑binding Rossmann‑fold domain/predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE718,5760GA46.5% 8.5 / 56.3 43intergenic (+19/‑91)mviM/mviNpredicted oxidoreductase, NAD(P)‑binding Rossmann‑fold domain/predicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (6/17);  new base A (11/9);  total (17/26)
Fisher's exact test for biased strand distribution p-value = 6.77e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGC  >  minE/718512‑718642
                                                                |                                                                  
gcAAGGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCaa                                                              <  1:3883113/71‑1 (MQ=255)
gcAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCaa                                                              <  1:2487287/71‑1 (MQ=255)
 cAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAAt                                                             <  1:3948810/71‑1 (MQ=255)
   aCGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAAtt                                                            >  1:1329124/1‑70 (MQ=255)
        tCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGATAGCAATTCACCa                                                       >  1:929282/1‑70 (MQ=255)
          gTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATa                                                     <  1:1603089/70‑1 (MQ=255)
          gTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGATa                                                                  >  1:629404/1‑57 (MQ=255)
             gACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATc                                                  <  1:4144111/70‑1 (MQ=255)
                aaGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAgg                                              <  1:953578/71‑1 (MQ=255)
                 aGATCTGGCGCGATGCGATGAGTGAATAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTa                                            <  1:781902/71‑1 (MQ=255)
                  gATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTa                                            <  1:2509335/71‑1 (MQ=255)
                  gATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTa                                            <  1:2509874/71‑1 (MQ=255)
                        gcgcGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGATAGCAATTCACCATAATCCAGGTAg                                           <  1:3473962/66‑1 (MQ=255)
                        gcgcGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGATAGCAATTCACCATAATCCAGGTAGACTa                                       <  1:2831576/70‑1 (MQ=255)
                             atgcgatgAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTa                                            <  1:378997/60‑1 (MQ=255)
                                  atgAGTGAATAAACCCTGTAACATCTGGCGATAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCttt                            <  1:3692378/71‑1 (MQ=255)
                                   tgAGTGAATAAACCCTGTAACATCTGGCGATAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCttt                            >  1:3687380/1‑70 (MQ=255)
                                     aGTGAATATACCCTGTAACATCTGGCGATAGCAATTCAcc                                                        <  1:365375/40‑1 (MQ=38)
                                     aGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATc                                                  <  1:169073/46‑1 (MQ=255)
                                       tGAATAAACCCTGTAACATCTGGCGATAGCAATTCACCATAATCCAGGTAGACTa                                       <  1:3496936/55‑1 (MQ=255)
                                       tGAATAAACCCTGTAACATCTGGCGATAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAgcg                       <  1:1484121/71‑1 (MQ=255)
                                         aaTAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACt                                        >  1:3153307/1‑52 (MQ=255)
                                         aaTAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCa                          >  1:3283597/1‑66 (MQ=255)
                                         aaTAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCgcc                     >  1:3637639/1‑71 (MQ=255)
                                         aaTAAACCCTGTAACATCTGGCGATAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCgcc                     >  1:2145563/1‑71 (MQ=255)
                                           tAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAg                                           <  1:1671724/47‑1 (MQ=255)
                                           tAAACCCTGTAACATCTGGCGATAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCgcctg                   >  1:1501642/1‑71 (MQ=255)
                                             aaCCCTGTAACATCTGGCGATAGCAATTCACCATAATCCa                                                >  1:4142805/1‑40 (MQ=255)
                                             aaCCCTGTAACATCTGGCGATAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCgcctgcc                 >  1:2759564/1‑71 (MQ=255)
                                              aCCCTGTAACATCTGGCGATAGCAATTCACCATAATCCAGGTAGACTAtt                                     >  1:1196539/1‑50 (MQ=255)
                                                ccTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTg              <  1:3901170/71‑1 (MQ=255)
                                                 cTGTAACATCTGGCGATAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGc             <  1:2278125/71‑1 (MQ=255)
                                                    tAACATCTGGCGGTTGCAATTCACCATAATCCAGGTAGACTATTc                                    <  1:223754/45‑1 (MQ=255)
                                                    tAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCctc                               >  1:3907160/1‑50 (MQ=255)
                                                    tAACATCTGGCGATAGCAATTCACCATAATCCAGGTAGa                                          >  1:2562107/1‑39 (MQ=255)
                                                       cATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGt       <  1:343176/71‑1 (MQ=255)
                                                       cATCTGGCGATAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTg              <  1:1496364/64‑1 (MQ=255)
                                                        aTCTGGCGGTAGCAATTCACCATAATCCAGGTAGAc                                         <  1:1717195/36‑1 (MQ=255)
                                                        aTCTGGCGATAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGtt      <  1:3642004/71‑1 (MQ=255)
                                                          cTGGCGATAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTc                           >  1:1209423/1‑48 (MQ=255)
                                                          cTGGCGATAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGttt     >  1:2830832/1‑70 (MQ=255)
                                                           tGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTg   >  1:2799918/1‑71 (MQ=255)
                                                             gCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGc  <  1:3493422/70‑1 (MQ=255)
                                                                |                                                                  
GCAACGTATCGTTGACAAGATCTGGCGCGATGCGATGAGTGAATAAACCCTGTAACATCTGGCGGTAGCAATTCACCATAATCCAGGTAGACTATTCGCCTCTTTCAGCGCCTGCCTTGCAGGCGTTTTGC  >  minE/718512‑718642

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: