Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 854,888 T→C 9.4% K34K (AAA→AAG)  yciG ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE854,8880TC 9.4% 143.0 / 6.8 53K34K (AAA→AAG) yciGhypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (19/29);  new base C (2/3);  total (21/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.25e-01

GATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGC  >  minE/854823‑854956
                                                                 |                                                                    
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:1305664/3‑52 (MQ=255)
                     ttaccaccTTTTTTACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAAt                                                               >  1:3810581/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:335368/56‑1 (MQ=255)
                             ttttttACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTAt                                                   <  1:2961839/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)
                                                                 |                                                                    
GATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGC  >  minE/854823‑854956

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