Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 854,891 A→G 9.6% F33F (TTT→TTC)  yciG ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE854,8910AG 9.6% 140.4 / 6.0 52F33F (TTT→TTC) yciGhypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (19/28);  new base G (2/3);  total (21/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAA  >  minE/854823‑854959
                                                                    |                                                                    
gATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTaa                                                                     >  1:3311015/1‑70 (MQ=255)
 aTTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTaaaa                                                                   >  1:2608063/1‑71 (MQ=255)
 aTTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTaaaa                                                                   <  1:1728286/71‑1 (MQ=255)
 aTTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTaa                                                                     <  1:3055749/69‑1 (MQ=255)
 aTTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTaa                                                                     <  1:1100104/69‑1 (MQ=255)
  tttATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAAt                                                                  >  1:880963/1‑71 (MQ=255)
  tttATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAAt                                                                  <  1:775819/71‑1 (MQ=255)
  tttATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAAt                                                                  <  1:1195799/71‑1 (MQ=255)
  tttATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAAt                                                                  >  1:2845390/1‑71 (MQ=255)
  tttATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAAt                                                                  >  1:2621739/1‑71 (MQ=255)
   ttattaCCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAAt                                                                  >  1:3186597/1‑70 (MQ=255)
   ttattaCCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAAt                                                                  >  1:2959527/1‑70 (MQ=255)
         ccaccGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTAccgcc                                                           >  1:3602794/1‑71 (MQ=255)
           accGCTTTTTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACcgccgc                                                         <  1:2905950/71‑1 (MQ=255)
            ccGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACcgccgc                                                         >  1:1423511/1‑70 (MQ=255)
              gCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTAc                                                               <  1:1423181/62‑1 (MQ=255)
               cTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTAccg                                                             >  1:3628490/1‑63 (MQ=255)
               cTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTAccg                                                             >  1:976829/1‑63 (MQ=255)
                  tgttgACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTaa                                                                     <  1:3636562/52‑1 (MQ=255)
                     ttaccaccTTTTTTACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAAt                                                                  >  1:3810581/3‑52 (MQ=255)
                     ttaccaccTTTTTTACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAAt                                                                  >  1:1305664/3‑52 (MQ=255)
                          gCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCtt                                          <  1:1599748/71‑1 (MQ=255)
                           ccTTTTTTACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAATTACCCCCGCTGTGCTGTCCACCttt                                         <  1:4034793/71‑1 (MQ=255)
                           ccTTTTTTACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAATTACCCCCGCTGTGCTGTCCACCttt                                         <  1:3075387/71‑1 (MQ=255)
                            cTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAAtt                                                                 <  1:3144155/46‑1 (MQ=255)
                             ttttttACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAAt                                                                  <  1:2300730/44‑1 (MQ=255)
                             ttttttACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTAt                                                      <  1:2961839/56‑1 (MQ=255)
                             ttttttACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTAt                                                      <  1:335368/56‑1 (MQ=255)
                               ttttACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAATTACCCCCGCTGTGCTGTCCACCtttttt                                      <  1:1365589/70‑4 (MQ=255)
                               ttttACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTAc                                       <  1:2023585/69‑1 (MQ=255)
                               ttttACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACgg                                     >  1:3378161/1‑71 (MQ=255)
                                        cTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCg                            <  1:64729/71‑1 (MQ=255)
                                        cTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCg                            <  1:3970837/71‑1 (MQ=255)
                                          tCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGa                          >  1:854799/1‑71 (MQ=255)
                                            aGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATg                        <  1:1524835/71‑1 (MQ=255)
                                             gATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTATGGCCTGCGTCGGATGc                       <  1:3409029/71‑1 (MQ=255)
                                              aTGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGcc                      >  1:834722/1‑71 (MQ=255)
                                                gcgcGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCg                               <  1:2461211/60‑1 (MQ=255)
                                                gcgcGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCtt                    >  1:3489900/1‑71 (MQ=255)
                                                 cgcgTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCtt                    <  1:3442763/70‑1 (MQ=255)
                                                  gcgTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCtt                    <  1:3753653/69‑1 (MQ=255)
                                                  gcgTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTtct                  <  1:2298393/71‑1 (MQ=255)
                                                   cgTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTtctc                 <  1:3162706/71‑1 (MQ=255)
                                                     ttGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCAc               <  1:1642648/71‑1 (MQ=255)
                                                     ttGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCAc               <  1:1705479/71‑1 (MQ=255)
                                                         ggATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCtt                                          <  1:27010/40‑1 (MQ=255)
                                                          gATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCg                            <  1:319862/53‑1 (MQ=255)
                                                            tCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTtct                  <  1:1585511/61‑1 (MQ=255)
                                                            tCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCtt         >  1:885141/1‑70 (MQ=255)
                                                            tCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTc        >  1:467385/1‑71 (MQ=255)
                                                             cATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCg       >  1:1068802/1‑71 (MQ=255)
                                                              aTTTTTAAAATTACCGCCGCTATGCTGACCGCCtt                                          <  1:3278000/35‑1 (MQ=255)
                                                              aTTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCgg      >  1:2516019/1‑71 (MQ=255)
                                                               tttttAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGc     <  1:2007114/71‑1 (MQ=255)
                                                                ttttAAATTTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGc     <  1:3925395/70‑1 (MQ=255)
                                                                ttttAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGc     >  1:541787/1‑70 (MQ=255)
                                                                atttAAAATAACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTtc                   >  1:3664662/2‑56 (MQ=255)
                                                                   tAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGaa  <  1:2890543/70‑1 (MQ=255)
                                                                    aaaaTTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGa                          >  1:3120658/1‑45 (MQ=255)
                                                                    aaaaTTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGa                          >  1:1759880/1‑45 (MQ=255)
                                                                    |                                                                    
GATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAA  >  minE/854823‑854959

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: ‑