Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,759,854 C→T 21.9% L98L (CTC→CTT smg → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,759,8540CT21.9% 106.7 / 20.6 77L98L (CTC→CTTsmgconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (50/10);  new base T (13/4);  total (63/14)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-02

TGAACGACTGGATGCCAGCTGCCGTGGGTTTCTGCTTTTCCTTGAGCAGATTCAGGTGCTCAACCTTGAAACTCGTGAAATGGTGATAGAGCGAGTGCTCGCGCTGGATAACGCAGAGTTCGAAC  >  minE/2759794‑2759918
                                                            |                                                                
tgAACGACTGGATGCCAGCTGCCGTGGGTTTCTGCTTTTCCTTGAGCAGATTCAGGTGCTCAACCTTGa                                                          >  1:626551/1‑69 (MQ=255)
  aaagACTGGATGCCAGCTGCCGTGGGTTTCTGCTTTTCCTTGAGCAGATTCAGGTGCTCAACCTTGaaa                                                        >  1:3293107/4‑69 (MQ=255)
  aaCGACTGGATGCCAGCTGCCGTGGGTTTCTGCTTTTCCTTGAGCAGATTCAGGTGCTCAACCTTGaaa                                                        >  1:972833/1‑69 (MQ=255)
  aaCGACTGGATGCCAGCTGCCGTGGGTTTCTGCTTTTCCTTGAGCAGATTCAGGTGCTCAACCTTGaaa                                                        >  1:903879/1‑69 (MQ=255)
  aaCGACTGGATGCCAGCTGCCGTGGGTTTCTGCTTTTCCTTGAGCAGATTCAGGTGCTCAACCTTGaaa                                                        >  1:2166710/1‑69 (MQ=255)
  aaCGACTGGATGCCAGCTGCCGTGGGTTTCTGCTTTTCCTTGAGCAGATTCAGGTGCTCAACCTTGaaa                                                        >  1:337257/1‑69 (MQ=255)
  aaCGACTGGATGCCAGCTGCCGTGGGTTTCTGCTTTTCCTTGAGCAGATTCAGGTGCTCAACCTTGaa                                                         >  1:2700973/1‑68 (MQ=255)
  aaCGACTGGATGCCAGCTGCCGTGGGTTTCTGCTTTTCCTTGAGCAGATTCAGGTGCTCAACCTTGaa                                                         >  1:176939/1‑68 (MQ=255)
      aCTGGATGCCAGCTGCCGTGGGTTTCTGCTTTTCCTTGAGCAGATTCAGGTGCTCAACCTTGAAACt                                                      >  1:257288/1‑67 (MQ=255)
      aCTGGATGCCAGCTGCCGTGGGTTTCTGCTTTTCCTTGAGCAGATTCAGGTGCTCAACCTTGAAACTc                                                     >  1:1813808/1‑68 (MQ=255)
      aCTGGATGCCAGCTGCCGTGGGTTTCTGCTTTTCCTTGAGCAGATTCAGGTGCTCAACCTTGAAACTc                                                     >  1:2463365/1‑68 (MQ=255)
      aCTGGATGCCAGCTGCCGTGGGTTTCTGCTTTTCCTTGAGCAGATTCAGGTGCTCAACCTTGAAACTc                                                     >  1:2263112/1‑68 (MQ=255)
      aCTGGATGCCAGCTGCCGTGGGTTTCTGCTTTTCCTTGAGCAGATTCAGGTGCTCAACCTTGAAACTc                                                     >  1:1279677/1‑68 (MQ=255)
      aCTGGATGCCAGCTGCCGTGGGTTTCTGCTTTTCCTTGAGCAGATTCAGGTGCTCAACCTTGAAACTc                                                     >  1:1453759/1‑68 (MQ=255)
      aCTGGATGCCAGCTGCCGTGGGTTTCTGCTTTTCCTTGAGCAGATTCACGTGCTCAACCTTGAAACTc                                                     >  1:2654051/1‑68 (MQ=255)
          gATGCCAGCTGCCGTGTGTTTCTGCTTTTCCTTGAGCAGATTCAGGTGCTCAACCTTGAAACTCGTGaa                                                >  1:3221758/1‑69 (MQ=255)
          gATGCCAGCTGCCGTGGGTTTCTGCTTTTCCTTGAGCAGATTCAGGTGCTCAACCTTGaaa                                                        >  1:2534467/1‑61 (MQ=255)
          gATGCCAGCTGCCGTGGGTTTCTGCTTTTCCTTGAGCAGATTCAGGTGCTCAACCTTGAAACTCGTGaa                                                >  1:1599572/1‑69 (MQ=255)
          gATGCCAGCTGCCGTGGGTTTCTGCTTTTCCTTGAGCAGATTCAGGTGCTCAACCTTGAAACTCGTGaa                                                >  1:3061628/1‑69 (MQ=255)
          gATGCCAGCTGCCGTGGGTTTCTGCTTTTCCTTGAGCAGATTCAGGTGCTCAACCTTGAAACTCGTGaa                                                >  1:1746743/1‑69 (MQ=255)
                aGCTGTCGTGGTTTTCTGCTTTTCCTTGAGCAAATTCAGGTGCTTAACCTTGAAACCCGTGAAATgg                                            >  1:1959298/1‑67 (MQ=255)
                aGCTGTCGTGGTTTTCTGCTTTTCCTTGAGCAAATTCAGGTGCTTAACCTTGAAACCCGTGAAATTGt                                           >  1:3378190/1‑68 (MQ=255)
                aGCTGTCGTGGTTTTCTGCTTTTCCTTGAGCAAATTCAGGTGCTTAACCTTGAAACCCGTGAAATGGt                                           >  1:3134509/1‑68 (MQ=255)
                aGCTGTCGTGGTTTTCTGCTTTTCCTTGAGCAAATTCAGGTGCTTAACCTTGAAACCCGTGAAATGGt                                           >  1:2103864/1‑68 (MQ=255)
                aGCTGTCGTGGTTTTCTGCTTTTCCTTGAGCAAATTCAGGTGCTTAACCTTGAAACCCGTGAAATGGt                                           >  1:1846463/1‑68 (MQ=255)
                aGCTGTCGTGGTTTTCTGCTTTTCCTTGAGCAAATTCAGGTGCTTAACCTTGAAACCCGTGAAATGGt                                           >  1:2235242/1‑68 (MQ=255)
                aGCTGTCGTGGTTTTCTGCTTTTCCTTGAGCAAATTCAGGTGCTTAACCTTGAAACCCGTGAAATGGt                                           >  1:3393276/1‑68 (MQ=255)
                aGCTGTCGTGGTTTTCTGCTTTTCCTTGAGCAAATTCAGGTGCTTAACCTTGAAACCCGTGAAATGGt                                           >  1:3143663/1‑68 (MQ=255)
                aGCTGTCGTGGTTTTCTGCTTTTCCTTGAGCAAATTCAGGTGCTTAACCTTGAAACCCGTGAAATGGt                                           >  1:289756/1‑68 (MQ=255)
                aGCTGTCGTGGTTTTCTGCTTTTCCTTGAGCAAATTCAGGTGCTTAACCTTGAAACCCGTGAAATGGt                                           >  1:2759781/1‑68 (MQ=255)
                aGCTGTCGTGGTTTTCTGCTTTTCCTTGAGCAAATTCAGGTGCTTAACCTTGAAACCCGTGAAATGGt                                           >  1:1273290/1‑68 (MQ=255)
                aGCTGTCGTGGTTTTCTGCTTTTCCTTGAGCAAATTCAGGTGCTTAACCTTGAAACCCGTGAAATGGt                                           >  1:888069/1‑68 (MQ=255)
                aGCTGTCGTGGTTTTCTGCTTTTCCTTGAGCAAATTCAGGTGCTTAACCTTGAAACCCGTGAAATGGt                                           >  1:1494657/1‑68 (MQ=255)
                    gCCGTGGGTTTCTGCTTTTCCTTGAGCAGATTCAGGTGCTCAACCTTGAAACTCGTGAAATGGTGATag                                      >  1:2887856/1‑69 (MQ=255)
                       gTGGGTTTCTGCTTTTCCTTGAGCAGATTCAGGTGCTCAACCTTGAAACTCGTGAAATGGTGATAGAGc                                   <  1:677303/69‑1 (MQ=255)
                       gTGGGTTTCTGCTTTTCCTTGAGCAGATTAAGGTGCTCAACCTTGAAACTCGTGAAATGGTGATAGAGc                                   <  1:3051129/69‑1 (MQ=255)
                             ttCTGCTTTTCCTTGAGCAAATTCAGGTGCTTAACCTTGAAACCCGTGAAATGGTGATAGAGCGAGTGc                             <  1:2356514/69‑1 (MQ=255)
                             ttCTGCTTTTCCTTGAGCAAATTCAGGTGCTTAACCTTGAAACCCGTGAAATGGTGATAGAGCGAGTGc                             <  1:411909/69‑1 (MQ=255)
                             ttCTGCTTTTCCTTGAGCAAATTCAGGTGCTTAACCTTGAAACCCGTGAAATGGTGATAGAGCGAGTGc                             <  1:1588190/69‑1 (MQ=255)
                             ttCTGCTTTTCCTTGAGCAAATTCAGGTGCTTAACCTTGAAACCCGTGAAATGGTGATAGAGCGAGTGc                             <  1:2018727/69‑1 (MQ=255)
                              tCTGCTTTTCCTTGAGCAGATTCAGGTGCTCAACCTTGAAACTCGTGaaa                                               >  1:3111575/1‑50 (MQ=255)
                              tCTGCTTTTCCTTGAGCAGATTCAGGTGCTCAACCTTGAAACTCGTGAAATGGTGa                                         >  1:2417749/1‑56 (MQ=255)
                              tCTGCTTTTCCTTGAGCAGATTCAGGTGCTCAACCTTGAAACTCGTGAAATGGTGa                                         >  1:709709/1‑56 (MQ=255)
                              tCTGCTTTTCCTTGAGCAGATTCAGGTGCTCAACCTTGAAACTCGTGAAATGGTGa                                         >  1:419356/1‑56 (MQ=255)
                              tCTGCTTTTCCTTGAGCAGATTCAGGTGCTCAACCTTGAAACTCGTGAAATGGTGATAGAGc                                   <  1:917910/62‑1 (MQ=255)
                              tCTGCTTTTCCTTGAGCAGATTCAGGTGCTCAACCTTGAAACTCGTGAAATGGTGATAGAGCGAGTGc                             >  1:3361434/1‑68 (MQ=255)
                              tCTGCTTTTCCTTGAGCAGATTCAGGTGCTCAACCTTGAAACTCGTGAAATGGTGATAGAGCGAGTGCt                            >  1:1165704/1‑69 (MQ=255)
                              tCTGCTTTTCCTTGAGCAGATTCAGGTGCTCAACCTTGAAACTCGTGAAATGGTGATAGAGCGAGTGCt                            >  1:95377/1‑69 (MQ=255)
                              tCTGCTTTTCCTTGAGCAGATTCAGGTGCTCAACCTTGAAACTCGTGAAATGGTGATAGAGCGAGTGCt                            >  1:1780765/1‑69 (MQ=255)
                              tCTGCTTTTCCTTGAGCAGATTCAGGTGCTCAACCTTGAAACTCGTGAAATGGTGATAGAGCGAGTGCt                            >  1:3188863/1‑69 (MQ=255)
                              tCTGCTTTTCCTTGAGCAGATTCAGGTGCTCAACCTTGAAACTCGTGAAATGGTGATAGAGCGAGTGCt                            >  1:1825799/1‑69 (MQ=255)
                              tCTGCTTTTCCTTGAGCAGATTCAGGTGCTCAACCTTGAAACTCGTGAAATGGTGATAGAGCGAGTGCt                            >  1:1161432/1‑69 (MQ=255)
                              tCTGCTTTTCCTTGAGCAGATTCAGGTGCTCAACCTTGAAACTCGTGAAATGGTGATAGAGCGAGTGCt                            >  1:2359166/1‑69 (MQ=255)
                              tCTGCTTTTCCTTGAGCAGATTCAGGTGCTCAACCTTGAAACTCGTGAAATGGTGATAGAGCGAGTGCt                            >  1:3123933/1‑69 (MQ=255)
                                   ttttCCTTGAGCAGATTCAGGTGCTCAACCTTGAAACTCGTGAAATGGTGATAGAGCGAGTGCTCgcgc                       >  1:310072/1‑69 (MQ=255)
                                                   tCAGGTGCTCAACCTTGAAACTCGTGAAATGGTGATAGAGCGAGTGCTCGCGCTGGATAACGCAGAGt        >  1:517752/1‑68 (MQ=255)
                                                   tCAGGTGCTCAACCTTGAAACTCGTGAAATGGTGATAGAGCGAGTGCTCGCGCTGGATAACGCAGAGt        >  1:1249710/1‑68 (MQ=255)
                                                   tCAGGTGCTCAACCTTGAAACTCGTGAAATGGTGATAGAGCGAGTGCTCGCGCTGGATAACGCAGAGt        >  1:659377/1‑68 (MQ=255)
                                                   tCAGGTGCTCAACCTTGAAACTCGTGAAATGGTGATAGAGCGAGTGCTCGCGCTGGATAACGCAGAGt        >  1:3160170/1‑68 (MQ=255)
                                                   tCAGGTGCTCAACCTTGAAACTCGTGAAATGGTGATAGAGCGAGTGCTCGCGCTGCATAACGCagag         >  1:2854326/1‑67 (MQ=255)
                                                    cAGGTGCTCAACCTTGAAACTCGTGAAATGGTGATAGAGCGAGTGCTCGCGCTGGATAACGCAGAGtt       >  1:3112725/1‑68 (MQ=255)
                                                    cAGGTGCTCAACCTTGAAACTCGTGAAATGGTGATAGAGCGAGTGCTCGCGCTGGATAACGCAGAGtt       >  1:548476/1‑68 (MQ=255)
                                                    cAGGTGCTCAACCTTGAAACTCGTGAAATGGTGATAGAGCGAGTGCTCGCGCTGGATAACGCAGAGt        >  1:3124735/1‑67 (MQ=255)
                                                      ggTGCTCAACCTTGAAACTCGTGAAATGGTGATAGAGCGAGTGCTCGCGCTGGATAACGCAGAGTTCg     <  1:1674412/68‑1 (MQ=255)
                                                      ggTGCTCAACCTTGAAACTCGTGAAATGGTGATAGAGCGAGTGCTCGCGCTGGATAACGCAGAGTTCg     <  1:478591/68‑1 (MQ=255)
                                                      ggTGCTCAACCTTGAAACTCGTGAAATGGTGATAGAGCGAGTGCTCGCGCTGGATAACGCAGAGTTCg     <  1:334555/68‑1 (MQ=255)
                                                      ggTGCTCAACCTTGAAACTCGTGAAATGGTGATAGAGCGAGTGCTCGCGCTGGATAACGCAGAGTTCg     <  1:706346/68‑1 (MQ=255)
                                                      ggTGCTCAACCTTGAAACTCGTGAAATGGTGATAGAGCGAGTGCTCGCGCTGGATAACGCAGAGTTCg     <  1:1189832/68‑1 (MQ=255)
                                                      ggTGCTCAACCTTGAAACTCGTGAAATGGTGATAGAGCGAGTGCTCGCGCTGGATAACGCAGAGTTCg     <  1:2597154/68‑1 (MQ=255)
                                                       gTGCTCAACCTTGAAACTCGTGAAATGGTGATAGAGCGAGTGCTCGCGCTGGATAACGCAGAGTTCg     <  1:1729673/67‑1 (MQ=255)
                                                        tGCTCAACCTTGAAACTCGTGAAATGGTGATAGAGCGAGTGCTCGCGCTGGATAACGCAGAGTTCGAAc  >  1:1360535/1‑69 (MQ=255)
                                                        tGCTCAACCTTGAAACTCGTGAAATGGTGATAGAGCGAGTGCTCGCGCTGGATAACGCAGAGTTCGAAc  >  1:650410/1‑69 (MQ=255)
                                                        tGCTCAACCTTGAAACTCGTGAAATGGTGATAGAGCGAGTGCTCGCGCTGGATAACGCAGAGTTCGAAc  >  1:1389383/1‑69 (MQ=255)
                                                        tGCTCAACCTTGAAACTCGTGAAATGGTGATAGAGCGAGTGCTCGCGCTGGATAACGCAGAGTTCGAAc  >  1:2481414/1‑69 (MQ=255)
                                                         gCTCAACCTTGAAACTCGTGAAATGGTGATAGAGCGAGTGCTCGCGCTGGATAACGCAGAGTTCGAAc  >  1:524823/1‑68 (MQ=255)
                                                         gCTCAACCTTGAAACTCGTGAAATGGTGATAGAGCGAGTGCTCGCGCTGGATAACGCAGAGTTCGAAc  >  1:368506/1‑68 (MQ=255)
                                                         gCTCAACCTTGAAACTCGTGAAATGGTGATAGAGCGAGTGCTCGCGCTGGATAACGCAGAGTTAGAAc  >  1:3314883/1‑68 (MQ=255)
                                                            |                                                                
TGAACGACTGGATGCCAGCTGCCGTGGGTTTCTGCTTTTCCTTGAGCAGATTCAGGTGCTCAACCTTGAAACTCGTGAAATGGTGATAGAGCGAGTGCTCGCGCTGGATAACGCAGAGTTCGAAC  >  minE/2759794‑2759918

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: