Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,854,662 T→C 64.6% I356M (ATA→ATG) ‡ yjeS ← predicted Fe‑S electron transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,854,6620TC64.6% 43.7 / 33.1 62I356M (ATA→ATG) ‡yjeSpredicted Fe‑S electron transport protein
Reads supporting (aligned to +/- strand):  ref base T (5/17);  new base C (9/31);  total (14/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.68e-01

AACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCGACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTG  >  minE/2854606‑2854730
                                                        |                                                                    
aaCCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCGACTATGCACGCAttt                                                          >  1:3325733/1‑69 (MQ=255)
           aaTCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCGACTATGCACGCATTTCGTCTCTCTa                                                <  1:794428/68‑1 (MQ=255)
           aaTCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCGACTATGCACGCATTTCGTCTCTCTa                                                <  1:1158844/68‑1 (MQ=255)
           aaTCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCGACTATGCACGCATTTCGTCTCTCTa                                                <  1:123378/68‑1 (MQ=255)
           aaTCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCGACTATGCACGCATTTCGTCTCTCTa                                                <  1:2552945/68‑1 (MQ=255)
           aaTCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCGACTATGCACGCATTTCGTCTCTCTa                                                <  1:2418921/68‑1 (MQ=255)
           aaTCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCGACTATGCACGCATTTCGTCTCTCTa                                                <  1:2116875/68‑1 (MQ=255)
           aaTCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCGACTATGCACGCATTTCGTCTCTCTa                                                <  1:1891284/68‑1 (MQ=255)
           aaTCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCGACTATGCACGCATTTCGTCTCTCTa                                                <  1:1865837/68‑1 (MQ=255)
           aaTCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCGACTATGCACGCATTTCGTCTCTCTa                                                <  1:1603856/68‑1 (MQ=255)
            aTCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCGACTATGCACGCATTTCGTCTCTCTa                                                <  1:1127586/67‑1 (MQ=255)
            aTCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCGACTATGCACGCATTTCGTCTCTCTa                                                <  1:835764/67‑1 (MQ=255)
            aTCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCGACTATGCACGCATTTCGTCTCTCTa                                                <  1:1586240/67‑1 (MQ=255)
             tCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCGACTATGCACGCATTTCGTCTCTCTAttt                                             >  1:777116/1‑69 (MQ=255)
               aCCCGAACCAGACGCTGTTTTTTCGGTAATTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTtgc                                           >  1:2766498/1‑69 (MQ=255)
               aCCCGAACCAGACGCTGTTTTTTCGGTAATTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTtgc                                           >  1:296791/1‑69 (MQ=255)
               aCCCGAACCAGACGCTGTTTTTTCGGTAATTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTtgc                                           >  1:2653249/1‑69 (MQ=255)
               aCCCGAACCAGACGCTGTTTTTTCGGTAATTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTtgc                                           >  1:2648339/1‑69 (MQ=255)
               aCCCGAACCAGACGCTGTTTTTTCGGTAATTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTtgc                                           >  1:2642364/1‑69 (MQ=255)
               aCCCGAACCAGACGCTGTTTTTTCGGTAATTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTtgc                                           >  1:2568991/1‑69 (MQ=255)
               aCCCGAACCAGACGCTGTTTTTTCGGTAATTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTtgc                                           >  1:131832/1‑69 (MQ=255)
               aCCCGAACCAGACGCTGTTTTTTCGGTAATTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTtgc                                           >  1:2472206/1‑69 (MQ=255)
               aCCCGAACCAGACGCTGTTTTTTCGGTAATTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTtgc                                           >  1:2258613/1‑69 (MQ=255)
                    aaCCAGACGCTGTTTTTTCGGCAATTGCACTTCGACTATGCACGCATTTCGTCTCTCTa                                                <  1:2068060/59‑1 (MQ=255)
                       cAGACGCTGTTTTTTCGGCAATTGCACTTCGACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAAt                                   <  1:1848239/69‑1 (MQ=255)
                            gCTGTTTTTTCGGCAATTGCACTTCGACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGcc                               >  1:239483/1‑68 (MQ=255)
                                        gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                   <  1:2751448/66‑1 (MQ=255)
                                        gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                   <  1:3333228/66‑1 (MQ=255)
                                        gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                   <  1:1145575/66‑1 (MQ=255)
                                        gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                   <  1:3420124/66‑1 (MQ=255)
                                        gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                   <  1:3172502/66‑1 (MQ=255)
                                        gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                   <  1:760022/66‑1 (MQ=255)
                                        gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                   <  1:777811/66‑1 (MQ=255)
                                        gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                   <  1:3150039/66‑1 (MQ=255)
                                        gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                   <  1:1074309/66‑1 (MQ=255)
                                        gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                   <  1:866135/66‑1 (MQ=255)
                                        gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                   <  1:978742/66‑1 (MQ=255)
                                        gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                   <  1:994408/66‑1 (MQ=255)
                                        gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                   <  1:2832643/66‑1 (MQ=255)
                                        gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                   <  1:2808574/66‑1 (MQ=255)
                                        gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                   <  1:998628/66‑1 (MQ=255)
                                        gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                   <  1:2755457/66‑1 (MQ=255)
                                        gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                   <  1:1639766/66‑1 (MQ=255)
                                        gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                   <  1:1873320/66‑1 (MQ=255)
                                        gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                   <  1:1874225/66‑1 (MQ=255)
                                        gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                   <  1:1579033/66‑1 (MQ=255)
                                        gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                   <  1:1919569/66‑1 (MQ=255)
                                        gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                   <  1:194372/66‑1 (MQ=255)
                                        gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                   <  1:2095978/66‑1 (MQ=255)
                                        gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                   <  1:2133725/66‑1 (MQ=255)
                                        gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                   <  1:2399941/66‑1 (MQ=255)
                                        gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                   <  1:1429939/66‑1 (MQ=255)
                                        gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                   <  1:2516528/66‑1 (MQ=255)
                                        gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                   <  1:2598956/66‑1 (MQ=255)
                                        gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                   <  1:1612705/66‑1 (MQ=255)
                                        gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                   <  1:2751734/66‑1 (MQ=255)
                                        gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                   <  1:2801920/66‑1 (MQ=255)
                                         cAATTGCACTTCGACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTc                 <  1:1523506/69‑1 (MQ=255)
                                             tGCACTTCGACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGc                          >  1:779006/1‑56 (MQ=255)
                                             tGCACTTCGACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGc                          >  1:1139063/1‑56 (MQ=255)
                                                 cTTCGACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCgaga          <  1:773294/68‑1 (MQ=255)
                                                 cTTCGACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCgaga          <  1:2974026/68‑1 (MQ=255)
                                                        tATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTg  <  1:846820/69‑1 (MQ=255)
                                                        |                                                                    
AACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCGACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTG  >  minE/2854606‑2854730

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: