Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,483,575 T→C 69.3% T307T (ACA→ACG trmE ← GTPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,483,5750TC69.3% 20.0 / 8.3 26T307T (ACA→ACGtrmEGTPase
Reads supporting (aligned to +/- strand):  ref base T (2/6);  new base C (13/5);  total (15/11)
Fisher's exact test for biased strand distribution p-value = 3.84e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

TTCGCTGGCAGACGGGCAATAAATTCCGGCCAGATCTCTGCCGGATCCACGGCGTCTGTTGTGGTGCCATCGACCATAAACAGCACGCGGTCGGCCTGTTCAATTTCCTGCCACGCGCGCTC  >  minE/2483516‑2483637
                                                           |                                                              
ttCGCTGGCAGACGGGCAATAAATTCCGGCCAGATCTCTGCCGGATCCACGGCGTCtgttgt                                                              <  1:661935/62‑1 (MQ=255)
ttCGCTGGCAGACGGGCAATAAATTCCGGCCAGATCTCTGCCGGATCCACGGCGTCtgttgt                                                              <  1:2177102/62‑1 (MQ=255)
      ggCAGACGGGCAATAAATTCCGGCCAGATCTCTGCCGGATCCACGGCGTCTGTTgtggt                                                           >  1:2094733/1‑59 (MQ=255)
      ggCAGACGGGCAATAAATTCCGGCCAGATCTCTGCCGGATCCACGGCGTCTGTTGTGGTGCCATCGAc                                                  >  1:3170401/1‑68 (MQ=255)
         aGACGGGTGATAAATTCCGGCCAGATCTCTGCCGGATCAACGTCGTCTGTCGTGGTGCCATCGACCATa                                              <  1:334976/69‑1 (MQ=255)
         aGACGGGCGATAAATTCCGGCCAGATCTCTGCCGGATCAACGTCGTCTGTCGTGGTGCCATCGACCATa                                              <  1:257295/69‑1 (MQ=255)
         aGACGGGCGATAAATTCCGGCCAGATCTCTGCCGGATCAACGTCGTCTGTCGTGGTGCCATCGACCATa                                              <  1:92771/69‑1 (MQ=255)
         aGACGGGCGATAAATTCCGGCCAGATCTCTGCCGGATCAACGTCGTCTGTCGTGGTGCCATCGACCATa                                              <  1:2281587/69‑1 (MQ=255)
         aGACGGGCGATAAATTCCGGCCAGATCTCTGCCGGATCAACGTCGTCTGTCGTGGTGCCATCGACCATa                                              <  1:1835681/69‑1 (MQ=255)
         aGACGGGCAATAAATTCCGGCCAGATCTCTGCCGGATCCACGGCGTCTGTTGTGGTGCCATCGACCATa                                              <  1:2426295/69‑1 (MQ=255)
            cGGGCAATAAATTCCGGCCAGATCTCTGCCGGATCCACGGCGTCTGTTGTGGTGCCATCGACCATAAAc                                           <  1:3102064/69‑1 (MQ=255)
            cGGGCAATAAATTCCGGCCAGATCTCTGCCGGATCCACGGCGTCTGTTGTGGTGCCATCGACCATAAAc                                           <  1:2187614/69‑1 (MQ=255)
                        tCCGGCCAGATCTCTGCCGGATCAACGGCGTCTGTCGTGGTGCCATCGACCATaa                                             >  1:3102699/1‑55 (MQ=255)
                        tCCGGCCAGATCTCTGCCGGATCAACGGCGTCTGTCGTGGTGCCATCGACCATaa                                             >  1:1482692/1‑55 (MQ=255)
                        tCCGGCCAGATCTCTGCCGGATCAACGGCGTCTGTCGTGGTGCCATCGACCATaa                                             >  1:915690/1‑55 (MQ=255)
                        tCCGGCCAGATCTCTGCCGGATCAACGGCGTCTGTCGTGGTGCCATCGACCATaa                                             >  1:745260/1‑55 (MQ=255)
                        tCCGGCCAGATCTCTGCCGGATCAACGGCGTCTGTCGTGGTGCCATCGACCATaa                                             >  1:678624/1‑55 (MQ=255)
                        tCCGGCCAGATCTCTGCCGGATCAACGGCGTCTGTCGTGGTGCCATCGACCATaa                                             >  1:1484087/1‑55 (MQ=255)
                        tCCGGCCAGATCTCTGCCGGATCAACGGCGTCTGTCGTGGTGCCATCGACCATaa                                             >  1:514860/1‑55 (MQ=255)
                        tCCGGCCAGATCTCTGCCGGATCAACGGCGTCTGTCGTGGTGCCATCGACCATaa                                             >  1:280552/1‑55 (MQ=255)
                        tCCGGCCAGATCTCTGCCGGATCAACGGCGTCTGTCGTGGTGCCATCGACCATaa                                             >  1:1086402/1‑55 (MQ=255)
                        tCCGGCCAGATCTCTGCCGGATCAACGGCGTCTGTCGTGGTGCCATCGACCATaa                                             >  1:2558059/1‑55 (MQ=255)
                        tCCGGCCAGATCTCTGCCGGATCAACGGCGTCTGTCGTGGTGCCATCGACCATaa                                             >  1:2379900/1‑55 (MQ=255)
                        tCCGGCCAGATCTCTGCCGGATCAACGGCGTCTGTCGTGGTGCCATCGACCATaa                                             >  1:2374309/1‑55 (MQ=255)
                        tCCGGCCAGATCTCTGCCGGATCAACGGCGTCTGTCGTGGTGCCATCGACCATa                                              >  1:2739243/1‑54 (MQ=39)
                                 aTCTCTGCCGGATCCACGGCGTCTGTTGTGGTGCCATCGACCATAAACAGCACGCGGTCGGCCTGTTCa                      <  1:155531/69‑1 (MQ=255)
                                                     gTCTGTTGTGGTGCCATCGACCATAAACAGCACGCGGTCGGCCTGTTCAATTTCCTGCCACGCGCGCTc  <  1:3327922/69‑1 (MQ=255)
                                                     gTCTGTTGTGGTGCCATCGACCATAAACAGCACGCGGTCGGCCTGTTCAATTTCCTGCCACGCGCGCTc  <  1:218048/69‑1 (MQ=255)
                                                           |                                                              
TTCGCTGGCAGACGGGCAATAAATTCCGGCCAGATCTCTGCCGGATCCACGGCGTCTGTTGTGGTGCCATCGACCATAAACAGCACGCGGTCGGCCTGTTCAATTTCCTGCCACGCGCGCTC  >  minE/2483516‑2483637

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: