Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 809,202 T→C 5.0% K275K (AAA→AAG ychF ← predicted GTP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported809,2020TC 5.0% 235.3 / 2.7 79K275K (AAA→AAGychFpredicted GTP‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (47/28);  new base C (2/2);  total (49/30)
Fisher's exact test for biased strand distribution p-value = 6.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.17e-01

CAACCGGAATGGTCCATGCACGCACTTCTTTCACCCCAGCGGTGAAGTAAGTTTGCAGGTTCAGCAGTTTATAACCGGCACGGATCACACGGTTCAGGCCCGGCTCTTCCAGCCCAAGCTCCTGCATAAACTCG  >  Exported/809135‑809268
                                                                   |                                                                   
cAACCGGAATGGTCCATGCACGCACTTCTTTCACCCCAGCGGTGAAGTAAGTTTGCAGGTTCAGCAGtt                                                                    <  1:2236303/69‑1 (MQ=255)
cAACCGGAATGGTCCATGCACGCACTTCTTTCACCCCAGCGGTGAAGTAAGTTTGCAGGTTCAGCAGTTta                                                                  >  1:1360953/1‑71 (MQ=255)
  aCCGGAATGGTCCATGCACGCACTTCTTTCACCCCAGCGGTGAAGTAAGTTTGCAGGTTCAGCAGTTtata                                                                <  1:2105317/71‑1 (MQ=255)
  aCCGGAATGGTCCATGCACGCACTTCTTTCACCCCAGCGGTGAAGTAAGTTTGCAGGTTCAGCAGTTtata                                                                <  1:844550/71‑1 (MQ=255)
  aCCGGAATGGTCCATGCACGCACTTCTTTCACCCCAGCGGTGAAGTAAGTTTGCAGGTTCAGCAGTTtata                                                                <  1:3482861/71‑1 (MQ=255)
  aCCGGAATGGTCCATGCACGCACTTCTTTCACCCCAGCGGTGAAGTAAGTTTGCAGGTTCAGCAGTTtata                                                                <  1:3904818/71‑1 (MQ=255)
  aCCGGAATGGTCCATGCACGCACTTCTTTCACCCCAGCGGTGAAGTAAGTTTGCAGGTTCAGCAGTTtata                                                                <  1:3996607/71‑1 (MQ=255)
   ccGGAATGGTCCATGCACGCACTTCTTTCACCCCAGCGGTGAAGTAAGTTTGCAGGTTCAGCAGTTTATaa                                                               >  1:971874/1‑71 (MQ=255)
   ccGGAATGGTCCATGCACGCACTTCTTTCACCCCAGCGGTGAAGTAAGTTTGCAGGTTCAGCAGTTTATaa                                                               >  1:2873368/1‑71 (MQ=255)
   ccGGAATGGTCCATGCACGCACTTCTTTCACCCCAGCGGTGAAGTAAGTTTGCAGGTTCAGCAGTTTATaa                                                               >  1:4592320/1‑71 (MQ=255)
   ccGGAATGGTCCATGCACGCACTTCTTTCACCCCAGCGGTGAAGTAAGTTTGCAGGTTCAGCAGCTTATaa                                                               <  1:4446611/71‑1 (MQ=255)
   ccGGAATGGTCCATGCACGCACTTCTTTCACCCCAGCGGTGAAGTAAGTTTGCAGGTTCAGCAGCTTATaa                                                               <  1:1379591/71‑1 (MQ=255)
    cGGAATGGTCCATGCACGCACTTCTTTCACCCCAGCGGTGAAGTAAGTTTGCAGGTTCAGCAGTTTATAAc                                                              >  1:1478646/1‑71 (MQ=255)
            tCCATGCACGCACTTCTTTCACCCCAGCGGTGAAGTAAGTTTGCAGGTTCAGCAGTTTATAACCGGCACgg                                                      <  1:864129/71‑1 (MQ=255)
             ccATGCACGCACTTCTTTCACCCCAGCGGTGAAGTAAGTTTGCAGGTTCAGCAGTTTATAACCGGCACGGa                                                     >  1:4435342/1‑71 (MQ=255)
             ccATGCACGCACTTCTTTCACCCCAGCGGTGAAGTAAGTTTGCAGGTTCAGCAGTTTATAACCGGCACGGa                                                     >  1:4378961/1‑71 (MQ=255)
             ccATGCACGCACTTCTTTCACCCCAGCGGTGAAGTAAGTTTGCAGGTTCAGCAGTTTATAACCGGCACGGa                                                     >  1:3298644/1‑71 (MQ=255)
                      cacTTCTTTCACCCCAGCGGTGAAGTAAGTTTGCAGGTTCAGCAGTTTATAACCGGCACGGATcac                                                 >  1:1423415/1‑66 (MQ=255)
                      cacTTCTTTCACCCCAGCGGTGAAGTAAGTTTGCAGGTTCAGCAGTTTATAACCGGCACGGATcac                                                 >  1:2803651/1‑66 (MQ=255)
                      cacTTCTTTCACCCCAGCGGTGAAGTAAGTTTGCAGGTTCAGCAGTTTATAACCGGCACGGATcac                                                 >  1:2811274/1‑66 (MQ=255)
                      cacTTCTTTCACCCCAGCGGTGAAGTAAGTTTGCAGGTTCAGCAGTTTATAACCGGCACGGATcac                                                 >  1:505920/1‑66 (MQ=255)
                      cacTTCTTTCACCCCAGCGGTGAAGTAAGTTTGCAGGTTCAGCAGTTTATAACCGGCACGGATcac                                                 >  1:3933204/1‑66 (MQ=255)
                        cttcttTCACCCCAGCGGTGAAGTAAGTTTGCAGGTTCAGCAGTTTATAACCGGCACgg                                                      <  1:689897/59‑1 (MQ=255)
                        cttcttTCACCCCAGCGGTGAAGTAAGTTTGCAGGTTCAGCAGTTTATAACCGGCACgg                                                      <  1:202342/59‑1 (MQ=255)
                          tcttTCACCCCAGCGGTGAAGTAAGTTTGCAGGTTCAGCAGTTTATAACCGGCACGGATCACACGGTTCAg                                        >  1:2903481/1‑71 (MQ=255)
                          tcttTCACCCCAGCGGTGAAGTAAGTTTGCAGGTTCAGCAGTTTATAACCGGCACGGATCACACGGTTCAg                                        >  1:4955142/1‑71 (MQ=255)
                                   ccAGCGGTGAAGTAAGTTTGCAGGTTCAGCAGTTTATAAcc                                                             >  1:4222080/1‑41 (MQ=255)
                                   ccAGCGGTGAAGTAAGTTTGCAGGTTCAGCAGTTTATAACCGGCACGGATCACACGGtt                                           >  1:1109931/1‑59 (MQ=255)
                                   ccAGCGGTGAAGTAAGTTTGCAGGTTCAGCAGTTTATAACCGGCACGGATCACACGGTTCAgg                                       >  1:1877887/1‑63 (MQ=255)
                                   ccAGCGGTGAAGTAAGTTTGCAGGTTCAGCAGTTTATAACCGGCACGGATCACACGGTTCAgg                                       >  1:4202051/1‑63 (MQ=255)
                                   ccAGCGGTGAAGTAAGTTTGCAGGTTCAGCAGTTTATAACCGGCACGGATCACACGGTTCAgg                                       >  1:4267390/1‑63 (MQ=255)
                                   ccAGCGGTGAAGTAAGTTTGCAGGTTCAGCAGTTTATAACCGGCACGGATCACACGGTTCAgg                                       >  1:1888590/1‑63 (MQ=255)
                                   ccAGCGGTGAAGTAAGTTTGCAGGTTCAGCAGTTTATAACCGGCACGGATCACACGGTTCAgg                                       >  1:3212374/1‑63 (MQ=255)
                                   ccAGCGGTGAAGTAAGTTTGCAGGTTCAGCAGTTTATAACCGGCACGGATCACACGGTTCAg                                        >  1:1533818/1‑62 (MQ=255)
                                   ccAGCGGTGAAGTAAGTTTGCAGGTTCAGCAGTTTATAACCGGCACGGATCACACGGTTCAGGCCCGGctc                               >  1:2873018/1‑71 (MQ=255)
                                      gCGGTGAAGTAAGTTTGCAGGTTCAGCAGTTTATAACCGGCACGGATCACAGGGTTCAGGCCCGGctc                               >  1:3413985/1‑68 (MQ=255)
                                         gTGAAGTAAGTTTGCAGGTTCAGCAGTTTATAACCGGCACGGATCACACGGTTCAGGCCCGGCTCTTCCAg                        <  1:1904617/71‑1 (MQ=255)
                                          tGAAGTAAGTTTGCAGGTTCAGCAGTTTATAACCGGCACGGATCACACGGTTCAGGCCCGGctct                              <  1:2288095/65‑1 (MQ=255)
                                          tGAAGTAAGTTTGCAGGTTCAGCAGTTTATAACCGGCACGGATCACACGGTTCAGGCCCGGctct                              <  1:3970862/65‑1 (MQ=255)
                                          tGAAGTAAGTTTGCAGGTTCACCAGTTTATAACCGGCACGGATCACACGGTTCAGGCCCGGctct                              <  1:4774752/65‑1 (MQ=255)
                                            aagtaagtTTGAAGGTTCAGCAGTTTATAACCGGCACGGATCACACGGTTCAGGCCCGGCTCTTCCAGccc                     <  1:3595357/71‑1 (MQ=255)
                                              gtaagtTTGTAGGTTCAGCAGTTTATAACCGGCACGGATCACACGGTTCAGGCCCGGCTCTTCCAGccc                     >  1:2002028/1‑69 (MQ=255)
                                              gtaagtTTGCAGGTTCAGCAGTTTATAACCGGCAc                                                        >  1:545982/1‑35 (MQ=255)
                                              gtaagtTTGCAGGTTCAGCAGTTTATAACCGGCACGGATCACACGGTTCAGGCCCGGCTCTTCCAGccc                     >  1:4181760/1‑69 (MQ=255)
                                              gtaagtTTGCAGGTTCAGCAGTTTATAACCGGCACGGATCACACGGTTCAGGCCCGGCTCTTCCAGccc                     >  1:788391/1‑69 (MQ=255)
                                              gtaagtTTGCAGGTTCAGCAGTTTATAACCGGCACGGATCACACGGTTCAGGCCCGGCTCTTCCAGccc                     >  1:2585247/1‑69 (MQ=255)
                                                 agtTTGCAGGTTCAGCAGTTTATAACCGGCACGGATCACACGGTTCAGGCCCGGTTCTTCCAGCCCAAGc                 <  1:456227/70‑1 (MQ=255)
                                                 agtTTGCAGGTTCAGCAGTTTATAACCGGCACGGATCACACGGTTCAGGCCCGGCTCTTCCAGCCCAAGc                 <  1:4633634/70‑1 (MQ=255)
                                                 agtTTGCAGGTTCAGCAGTTTATAACCGGCACGGATCACACGGTTCAGGCCCGGCTCTTCCAGCCCAAGc                 <  1:1919890/70‑1 (MQ=255)
                                                 agtTTGCAGGTTCAGCAGTTTATAACCGGCACGGATCACACGGTTCAGGCCCGGCTCTTCCAGCCCAAGc                 <  1:3605704/70‑1 (MQ=255)
                                                 agtCTGCAGGCCCAGCAGCTTGAACGCGGCGCGGATCAGGCGATCCAGGCCCGGCTCGTCCATGCCCatat                 >  1:1524775/1‑68 (MQ=255)
                                                 agtCTGCAGGCCCAGCAGCTTGAACGCGGCGCGGATCAGGCGATCCAGGCCCGGCTCGTCCATGCCCatat                 >  1:825154/1‑68 (MQ=255)
                                                  gtTTGCAGGTTCAGCAGTTTATAACCGGCACGGATAACACGGTTCAGGCCCGGCTCTTCCAGCCCAAGc                 <  1:2866296/69‑1 (MQ=255)
                                                    ttGCAGGTTCAGCAGTTTATAACCGGCACGGATcac                                                 >  1:201131/1‑36 (MQ=255)
                                                    ttGCAGGTTCAGCAGTTTATAACCGGCACGGATCACACGGTTCAGGCCCGGCTCTTCCAGCCCAAGCt                >  1:3259237/1‑68 (MQ=255)
                                                    ttGCAGGTTCAGCAGTTTATAACCGGCACGGATCACACGGTTCAGGCCCGGCTCTTCCAGCCCAAGCt                >  1:449448/1‑68 (MQ=255)
                                                    ttGCAGGTTCAGCAGTTTATAACCGGCACGGATCACACGGTTCAGGCCCGGCTCTTCCAGCCCAAGCt                >  1:1844996/1‑68 (MQ=255)
                                                      gCAGGTTCAGCAGTTTATAACCGGCACGGATCACACGGTTCAGGCCCGGCTCTTCCAGCCCAAGCTCCTg            <  1:2460692/70‑1 (MQ=255)
                                                      gCAGGTTCAGCAGTTTATAACCGGCACGGATCACACGGTTCAGGCCCGGCTCTTCCAGCCCAAGCTCCTg            <  1:2729211/70‑1 (MQ=255)
                                                      gCAGGTTCAGCAGTTTATAACCGGCACGGATCACACGGTTCAGGCCCGGCTCTTCCAGCCCAAGCTCCTg            <  1:2947293/70‑1 (MQ=255)
                                                      gCAGGTTCAGCAGTTTATAACCGGCACGGATCACACGGTTCAGGCCCGGCTCTTCCAGCCCAAGCTCCTg            <  1:3721644/70‑1 (MQ=255)
                                                         ggTTCAGCAGTTTATAACCGGCACGGATCACACGGTTCAGGCCCGGCTCTTCCAGCCCAAGCTCCTGCATa        <  1:4746247/71‑1 (MQ=255)
                                                         ggTTCAGCAGTTTATAACCGGCACGGATCACACGGTTCAGGCCCGGCTCTTCCAGCCCAAGCTCCTGCATa        <  1:3092165/71‑1 (MQ=255)
                                                         ggTTCAGCAGTTTATAACCGGCACGGATCACACGGTTCAGGCCCGGCTCTTCCAGCCCAAGCTCCTGCATa        <  1:1721532/71‑1 (MQ=255)
                                                         ggTTCAGCAGTTTATAACCGGCACGGATCACACGGTTCAGGCCCGGCTCTTCCAGCCCAAGCTCCTGCATa        <  1:1553037/71‑1 (MQ=255)
                                                          gTTCAGCAGTTTATAACCGGCACGGATCACACGGTTCAGGCCCGGCTCTTCCAGCCCAAGCTCCTGCATa        <  1:1036713/70‑1 (MQ=255)
                                                           ttCAGCAGTTTATAACCGGCACGGATCACACGGTTCAGGCCCGGCTCTTCCAGCCCAAGCTCCTGCATa        <  1:3821927/69‑1 (MQ=255)
                                                            tCAGCAGTTTATAACCGGCACGGATCACACGGTTc                                          >  1:1173467/1‑35 (MQ=255)
                                                            tCAGCAGTTTATAACCGGCACGGATCACACGGTTc                                          >  1:3439350/1‑35 (MQ=255)
                                                            tCAGCAGTTTATAACCGGCACGGATCACACGGTTc                                          >  1:3060226/1‑35 (MQ=255)
                                                            tCAGCAGTTTATAACCGGCACGGATCACACGGTTc                                          >  1:1176371/1‑35 (MQ=255)
                                                            tCAGCAGTTTATAACCGGCACGGATCACACGGTTc                                          >  1:1858406/1‑35 (MQ=255)
                                                            tCAGCAGTTTATAACCGGCACGGATCACACGGTTc                                          >  1:1421044/1‑35 (MQ=255)
                                                            tCAGCAGTTTATAACCGGCACGGATCACACGGTTc                                          >  1:2491692/1‑35 (MQ=255)
                                                            tCAGCAGTTTATAACCGGCACGGATCACACGGTTc                                          >  1:1324625/1‑35 (MQ=255)
                                                            tCAGCAGTTTATAACCGGCACGGATCACACGGTTc                                          >  1:4179861/1‑35 (MQ=255)
                                                            tCAGCAGTTTATAACCGGCACGGATCACACGGTTc                                          >  1:644447/1‑35 (MQ=255)
                                                               gcagTTTATAACCGGCACGGATCACACGGTTCAGGCCCGGCTCTTCCAGCCCAAGCTCCTGCATAAACTCg  >  1:4116666/1‑71 (MQ=255)
                                                               gcagTTTATAACCGGCACGGATCACACGGTTCAGGCCCGGCTCTTCCAGCCCAAGCTCCTGCATAAACTCg  >  1:2003882/1‑71 (MQ=255)
                                                               gcagTTTATAACCGGCACGGATCACACGGTTCAGGCCCGGCTCTTCCAGCCCAAGCTCCTGCATAAACTCg  >  1:5051012/1‑71 (MQ=255)
                                                                   tttATAACCGGCACGGATCACACGGTTCAGGCCCGGCTCTTCCAGCCCAAGc                 <  1:1276752/52‑1 (MQ=255)
                                                                   tttATAACCGGCACGGATCACACGGTTCAGGCCCGGCTCTTCCAGCCCAAGc                 <  1:188544/52‑1 (MQ=255)
                                                                   tttATAACCGGCACGGATCACACGGTTCAGGCCCGGCTCTTCCAGCCCAAGc                 <  1:3651168/52‑1 (MQ=255)
                                                                   tttATAACCGGCACGGATCACACGGTTCAGGCCCGGCTCTTCCAGCCCAAGc                 <  1:3820867/52‑1 (MQ=255)
                                                                   |                                                                   
CAACCGGAATGGTCCATGCACGCACTTCTTTCACCCCAGCGGTGAAGTAAGTTTGCAGGTTCAGCAGTTTATAACCGGCACGGATCACACGGTTCAGGCCCGGCTCTTCCAGCCCAAGCTCCTGCATAAACTCG  >  Exported/809135‑809268

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: