Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,874,661 A→G 39.9% G80G (GGT→GGC ygcF ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,874,6610AG39.9% 17.6 / 32.4 51G80G (GGT→GGCygcFconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (8/22);  new base G (14/6);  total (23/28)
Fisher's exact test for biased strand distribution p-value = 3.72e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAGCAAATCATGAATGCAAGGCTCACCACCCGTAATCACCACATGCCGCGCGGTGTATCCCTGGCGACCAATGACAGCCAGCAAATCTTCACTGCTCGCAGCCCCCCACTTATCACTCTCTTTGGTCTTCGCCAG  >  minE/1874594‑1874729
                                                                   |                                                                    
gcagcaAATCATGAATGCAAGGCTCACCACCCGTAATCACCACATGCCGCGCGGTGTATCCCTGGCGa                                                                      >  1:1643269/1‑68 (MQ=255)
gcagcaAATCATGAATGCAAGGCTCACCACCCGTAATCACCACATGCCGCGCGGTGTATCCCTGGCGa                                                                      >  1:561879/1‑68 (MQ=255)
gcagcaAATCATGAATGCAAGGCTCACCACCCGTAATCACCACATGCCGCGCGGTGTATCCCTGGCGa                                                                      >  1:396138/1‑68 (MQ=255)
    caAATCATGAATGCAAGGCTCACCACCCGTAATCACCACATGCCGCGCGGTGTATCCCTGGCGACCAAt                                                                 <  1:2138262/69‑1 (MQ=255)
      aaTCATGAATGCAAGGCTCACCACCCGTAATCACCACATGACGCGCGGTGTATCCCTGACGGCCAATGa                                                               <  1:895232/69‑1 (MQ=255)
      aaTCATGAATGCAAGGCTCACCACCCGTAATCACCACATGACGCGCGGTGTATCCCTGACGGCCAATGa                                                               <  1:2028702/69‑1 (MQ=255)
      aaTCATGAATGCAAGGCTCACCACCCGTAATCACCACATGACGCGCGGTGTATCCCTGACGGCCAATGa                                                               <  1:1606259/69‑1 (MQ=255)
      aaTCATGAATGCAAGGCTCACCACCCGTAATCACCACATGACGCGCGGTGTATCCCTGACGGCCAATGa                                                               <  1:1272269/69‑1 (MQ=255)
      aaTCATGAATGCAAGGCTCACCACCCGTAATCACCACATGACGCGCGGTGTATCCCTGACGGCCAATGa                                                               <  1:1236407/69‑1 (MQ=255)
       aTCATGAATGCAAGGCTCACCACCCGTAATCACCACATGACGCGCGGTGTATCCCTGACGGCCAATGa                                                               <  1:1931490/68‑1 (MQ=255)
          aTGAATGCAAGGCTCACCACCCGTAATCACCACATGCCGCGCGGTGTATCCCTGGCGACCAATGaaa                                                             >  1:611790/1‑65 (MQ=255)
          aTGAATGCAAGGCTCACCACCCGTAATCACCACATGCCGCGCGGTGTATCCCTGGCGACCAATGAca                                                             >  1:1216268/1‑67 (MQ=255)
          aTGAATGCAAGGCTCACCACCCGTAATCACCACATGCCGCGCGGTGTATCCCTGGCGACCAATGAca                                                             >  1:1615310/1‑67 (MQ=255)
          aTGAATGCAAGGCTCACCACCCGTAATCACCACATGCCGCGCGGTGTATCCCTGGCGACCAATGAca                                                             >  1:2993280/1‑67 (MQ=255)
               tGCAAGGCTCACCACCCGTAATCACCACATGCCGCGCGGTGTATCCCTGGCGACCAATGACAGCCAGCa                                                      <  1:2152122/69‑1 (MQ=255)
               tGCAAGGCTCACCACCCGTAATCACCACATGCCGCGCGGTGTATCCCTGGCGACCAATGACAGCCAGCa                                                      <  1:2529125/69‑1 (MQ=255)
               tGCAAGGCTCACCACCCGTAATCACCACATGCCGCGCGGTGTATCCCTGGCGACCAATGACAGCCAGCa                                                      <  1:1459656/69‑1 (MQ=255)
               tGCAAGGCTCACCACCCGTAATCACCACATGCCGCGCGGTGTATCCCTGGCGACCAATGACAGCCAGCa                                                      <  1:1772166/69‑1 (MQ=255)
               tGCAAGGCTCACCACCCGTAATCACCACATGCCGCGCGGTGTATCCCTGGCGACCAATGACAGCCAGCa                                                      <  1:605000/69‑1 (MQ=255)
               tGCAAGGCTCACCACCCGTAATCACCACATGCCGCGCGGTGTATCCCTGGCGACCAATGACAGCCAGCa                                                      <  1:1704249/69‑1 (MQ=255)
                gCAAGGCTCACCACCCGTAATCACCACATGCCGCGCGGTGTATCCCTGGCGACCAATGACAGCCAGCa                                                      <  1:2570068/68‑1 (MQ=255)
                   aGGCTCACCACCCGTAATCACCACATGCCGCGCGGTGTATCCCTGGCGACCAATGACAGCCAGCAAATc                                                  <  1:877618/69‑1 (MQ=255)
                      cTCACCACCCGTAATCACCACATGCCGCGCGGTGTATCCCTGGCGACCAATGACAGCCAGCAAATCTTc                                               <  1:2946219/69‑1 (MQ=255)
                      cTCACCACCCGTAATCACCACATGCCGCGCGGTGTATCCCTGGCGACCAATGACAGCCAGCAAATCTTc                                               <  1:338510/69‑1 (MQ=255)
                      cTCACCACCCGTAATCACCACATGCCGCGCGGTGTATCCCTGGCGACCAATGACAGCCAGCAAATCTTc                                               <  1:2623337/69‑1 (MQ=255)
                      cTCACCACCCGTAATCACCACATGCCGCGCGGTGTATCCCTGGCGACCAATGACAGCCAGCAAATCTTc                                               <  1:1844686/69‑1 (MQ=255)
                      cTCACCACCCGTAATCACCACATGCCGCGCGGTGTATCCCTGGCGACCAATGACAGCCAGCAAATCTTc                                               <  1:1728858/69‑1 (MQ=255)
                      cTCACCACCCGTAATCACCACATGCCGCGCGGTGTATCCCTGGCGACCAATGACAGCCAGCAAATCTTc                                               <  1:3194171/69‑1 (MQ=255)
                      cTCACCACCCGTAATCACCACATGCCGCGCGGTGTATCCCTGGCGACCAATGACAGCCAGCAAATCTTc                                               <  1:363965/69‑1 (MQ=255)
                                      accacATGCCGCGCGGTGTATCCCTGGCGACCAATGACAGCCAGCAAATCTTCACTGCTCGCAGccccc                               >  1:1514705/1‑69 (MQ=255)
                                      accacATGCCGCGCGGTGTATCCCTGGCGACCAATGACAGCCAGCAAATCTTCACTGCTCGCAGccccc                               >  1:1380270/1‑69 (MQ=255)
                                       ccacATGACGCGCGGTGTATCCCTGACGGCCAATGACAGCCAGCAAATCtt                                                >  1:1575527/1‑51 (MQ=255)
                                       ccacATGACGCGCGGTGTATCCCTGACGGCCAATGACAGCCAGCAAATCTTCACTGc                                          >  1:3055221/1‑57 (MQ=255)
                                       ccacATGACGCGCGGTGTATCCCTGACGGCCAATGACAGCCAGCAAATCTTCACTGc                                          >  1:255814/1‑57 (MQ=255)
                                       ccacATGACGCGCGGTGTATCCCTGACGGCCAATGACAGCCAGCAAATCTTCACTGc                                          >  1:3237963/1‑57 (MQ=255)
                                       ccacATGACGCGCGGTGTATCCCTGACGGCCAATGACAGCCAGCAAATCTTCACTGc                                          >  1:3247514/1‑57 (MQ=255)
                                       ccacATGACGCGCGGTGTATCCCTGACGGCCAATGACAGCCAGCAAATCTTCACTGc                                          >  1:2400786/1‑57 (MQ=255)
                                       ccacATGACGCGCGGTGTATCCCTGACGGCCAATGACAGCCAGCAAATCTTCACTGc                                          >  1:1195847/1‑57 (MQ=255)
                                       ccacATGACGCGCGGTGTATCCCTGACGGCCAATGACAGCCAGCAAATCTTCACTGc                                          >  1:210558/1‑57 (MQ=255)
                                       ccacATGACGCGCGGTGTATCCCTGACGGCCAATGACAGCCAGCAAATCTTCACTGc                                          >  1:1067383/1‑57 (MQ=255)
                                       ccacATGACGCGCGGTGTATCCCTGACGGCCAATGACAGCCAGCAAATCTTCACTGc                                          >  1:1486934/1‑57 (MQ=255)
                                       ccacATGACGCGCGGTGTATCCCTGACGGCCAATGACAGCCAGCAAATCTTCACTGc                                          >  1:1853926/1‑57 (MQ=255)
                                       ccacATGACGCGCGGTGTATCCCTGACGGCCAATGACAGCCAGCAAATCTTCACTGc                                          >  1:1765361/1‑57 (MQ=255)
                                       ccacATGACGCGCGGTGTATCCCTGACGGCCAATGACAGCCAGCAAATCTTCACTGc                                          >  1:1679283/1‑57 (MQ=255)
                                       ccacATGACGCGCGGTGTATCCCTGACGGCCAATGACAGCCAGCAAATCTTCACTGc                                          >  1:1757782/1‑57 (MQ=255)
                                            tGCCGCGCGGTGTATCCCTGGCGACCAATGACAGCCAGCAAATCTTCACTGCTCGCAGCCCCCCACTTa                         <  1:2093395/69‑1 (MQ=255)
                                            tGCCGCGCGGTGTATCCCTGGCGACCAATGACAGCCAGCAAATCTTCACTGCTCGCAGCCCCCCACTTa                         <  1:3389335/69‑1 (MQ=255)
                                            tGCCGCGCGGTGTATCCCTGGCGACCAATGACAGCCAGCAAATCTTCACTGCTCGCAGCCCCCCACTTa                         <  1:318513/69‑1 (MQ=255)
                                            tGCCGCGCGGTGTATCCCTGGCGACCAATGACAGCCAGCAAATCTTCACTGCTCGCAGCCCCCCACTTa                         <  1:601445/69‑1 (MQ=255)
                                                 cgcgGTGTATCCCTGGCGACCAATGACAGCCAGCAAATCTTCACTGCTCGCAGCCCCCCACTTATCAct                    <  1:433231/69‑1 (MQ=255)
                                                 cgcgGTGTATCCCTGGCGACCAATGACAGCCAGCAAATCTTCACTGCTCGCAGCCCCCCACTTATCAct                    <  1:360758/69‑1 (MQ=255)
                                                          tCCCTGGCGCCCAATGACAGCCAGCAAATCTTCACTGCTCGCAGCCGCCCACTTATCACTCTCTTTGGt           >  1:3318832/1‑69 (MQ=255)
                                                          tCCCTGGCGACCAATGACAGCCAGCAAATCTTCACTGCTCGCAGCCCCCCACTTATCACTCTCTTTGGt           >  1:1152158/1‑69 (MQ=255)
                                                          tCCCTGGCGACCAATGACAGCCAGCAAATCTTCACTGCTCGCAGCCCCCCACTTATCACTCTCTTTGGt           >  1:3325508/1‑69 (MQ=255)
                                                                   aCCAATGACAGCCAGCTAATCTTCACTGCTCGCAGCCCCCCACTTATCACTCTCTTTGGTCTTCGCCAg  <  1:441234/69‑1 (MQ=255)
                                                                   aCCAATGACAGCCAGCAAATCTTCACTGCTCGCAGCCCCCCACTTATCACTCTCTTTGGTCTTCGCCAg  <  1:1109688/69‑1 (MQ=255)
                                                                   aCCAATGACAGCCAGCAAATCTTCACTGCTCGCAGCCCCCCACTTATCACTCTCTTTGGTCTTCGCCAg  <  1:1152554/69‑1 (MQ=255)
                                                                   aCCAATGACAGCCAGCAAATCTTCACTGCTCGCAGCCCCCCACTTATCACTCTCTTTGGTCTTCGCCAg  <  1:38758/69‑1 (MQ=255)
                                                                   aCCAATGACAGCCAGCAAATCTTCACTGCTCGCAGCCCCCCACTTATCACTCTCTTTGGTCTTCGCCAg  <  1:1080366/69‑1 (MQ=255)
                                                                   aCCAATGACAGCCAGCAAATCTTCACTGCTCGCAGCCCCCCACTTATCACTCTCTTTGGTCTTCGCCAg  <  1:1165182/69‑1 (MQ=255)
                                                                   aCCAATGACAGCCAGCAAATCTTCACTGCTCGCAGCCCCCCACTTATCACTCTCTTTGGTCTTCGCCAg  <  1:863866/69‑1 (MQ=255)
                                                                   aCCAATGACAGCCAGCAAATCTTCACTGCTCGCAGCCCCCCACTTATCACTCTCTTTGGTCTTCGCCAg  <  1:350964/69‑1 (MQ=255)
                                                                   aCCAATGACAGCCAGCAAATCTTCACTGCTCGCAGCCCCCCACTTATCACTCTCTTTGGTCTTCGCCAg  <  1:887718/69‑1 (MQ=255)
                                                                   aCCAATGACAGCCAGCAAATCTTCACTGCTCGCAGCCCCCCACTTATCACTCTCTTTGGTCTTCGCCAg  <  1:1038945/69‑1 (MQ=255)
                                                                   aCCAATGACAGCCAGCAAATCTTCACTGCTCGCAGCCCCCCACTTATCACTCTCTTTGGTCTTCGCCAg  <  1:3304644/69‑1 (MQ=255)
                                                                   aCCAATGACAGCCAGCAAATCTTCACTGCTCGCAGCCCCCCACTTATCACTCTCTTTGGTCTTCGCCAg  <  1:3265453/69‑1 (MQ=255)
                                                                   aCCAATGACAGCCAGCAAATCTTCACTGCTCGCAGCCCCCCACTTATCACTCTCTTTGGTCTTCGCCAg  <  1:1321226/69‑1 (MQ=255)
                                                                   aCCAATGACAGCCAGCAAATCTTCACTGCTCGCAGCCCCCCACTTATCACTCTCTTTGGTCTTCGCCAg  <  1:1497251/69‑1 (MQ=255)
                                                                   aCCAATGACAGCCAGCAAATCTTCACTGCTCGCAGCCCCCCACTTATCACTCTCTTTGGTCTTCGCCAg  <  1:3108656/69‑1 (MQ=255)
                                                                   aCCAATGACAGCCAGCAAATCTTCACTGCTCGCAGCCCCCCACTTATCACTCTCTTTGGTCTTCGCCAg  <  1:1622345/69‑1 (MQ=255)
                                                                   aCCAATGACAGCCAGCAAATCTTCACTGCTCGCAGCCCCCCACTTATCACTCTCTTTGGTCTTCGCCAg  <  1:3031147/69‑1 (MQ=255)
                                                                   aCCAATGACAGCCAGCAAATCTTCACTGCTCGCAGCCCCCCACTTATCACTCTCTTTGGTCTTCGCCAg  <  1:3026422/69‑1 (MQ=255)
                                                                   aCCAATGACAGCCAGCAAATCTTCACTGCTCGCAGCCCCCCACTTATCACTCTCTTTGGTCTTCGCCAg  <  1:2995664/69‑1 (MQ=255)
                                                                   aCCAATGACAGCCAGCAAATCTTCACTGCTCGCAGCCCCCCACTTATCACTCTCTTTGGTCTTCGCCAg  <  1:2969954/69‑1 (MQ=255)
                                                                   aCCAATGACAGCCAGCAAATCTTCACTGCTCGCAGCCCCCCACTTATCACTCTCTTTGGTCTTCGCCAg  <  1:2842364/69‑1 (MQ=255)
                                                                   aCCAATGACAGCCAGCAAATCTTCACTGCTCGCAGCCCCCCACTTATCACTCTCTTTGGTCTTCGCCAg  <  1:1905189/69‑1 (MQ=255)
                                                                   aCCAATGACAGCCAGCAAATCTTCACTGCTCGCAGCCCCCCACTTATCACTCTCTTTGGTCTTCGCCAg  <  1:2531027/69‑1 (MQ=255)
                                                                   aCCAATGACAGCCAGCAAATCTTCACTGCTCGCAGCCCCCCACTTATCACTCTCTTTGGTCTTCGCCAg  <  1:2040729/69‑1 (MQ=255)
                                                                   aCCAATGACAGCCAGCAAATCTTCACTGCTCGCAGCCCCCCACTTATCACTCTCTTTGGTCTTCGCCAg  <  1:1029075/69‑1 (MQ=255)
                                                                   aCCAATGACAGCCAGCAAATCTTCACTGCTCGCAGCCCCCCACTTATCACTCTCTTTGGTCTTCGCCAg  <  1:2160485/69‑1 (MQ=255)
                                                                   |                                                                    
GCAGCAAATCATGAATGCAAGGCTCACCACCCGTAATCACCACATGCCGCGCGGTGTATCCCTGGCGACCAATGACAGCCAGCAAATCTTCACTGCTCGCAGCCCCCCACTTATCACTCTCTTTGGTCTTCGCCAG  >  minE/1874594‑1874729

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: