Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,854,664 T→C 63.6% I356V (ATA→GTA) ‡ 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,6640TC63.6% 39.7 / 35.6 63I356V (ATA→GTA) ‡yjeSpredicted Fe‑S electron transport protein
Reads supporting (aligned to +/- strand):  ref base T (5/18);  new base C (9/31);  total (14/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.47e-01

AACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCGACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCT  >  minE/2854606‑2854732
                                                          |                                                                    
aaCCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCGACTATGCACGCAttt                                                            >  1:3325733/1‑69 (MQ=255)
           aaTCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCGACTATGCACGCATTTCGTCTCTCTa                                                  <  1:2116875/68‑1 (MQ=255)
           aaTCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCGACTATGCACGCATTTCGTCTCTCTa                                                  <  1:1891284/68‑1 (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:2418921/68‑1 (MQ=255)
           aaTCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCGACTATGCACGCATTTCGTCTCTCTa                                                  <  1:1865837/68‑1 (MQ=255)
           aaTCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCGACTATGCACGCATTTCGTCTCTCTa                                                  <  1:2552945/68‑1 (MQ=255)
           aaTCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCGACTATGCACGCATTTCGTCTCTCTa                                                  <  1:1603856/68‑1 (MQ=255)
            aTCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCGACTATGCACGCATTTCGTCTCTCTa                                                  <  1:1586240/67‑1 (MQ=255)
            aTCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCGACTATGCACGCATTTCGTCTCTCTa                                                  <  1:835764/67‑1 (MQ=255)
            aTCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCGACTATGCACGCATTTCGTCTCTCTa                                                  <  1:1127586/67‑1 (MQ=255)
             tCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCGACTATGCACGCATTTCGTCTCTCTAttt                                               >  1:777116/1‑69 (MQ=255)
               aCCCGAACCAGACGCTGTTTTTTCGGTAATTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTtgc                                             >  1:2568991/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:2766498/1‑69 (MQ=255)
               aCCCGAACCAGACGCTGTTTTTTCGGTAATTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTtgc                                             >  1:296791/1‑69 (MQ=255)
               aCCCGAACCAGACGCTGTTTTTTCGGTAATTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTtgc                                             >  1:2472206/1‑69 (MQ=255)
               aCCCGAACCAGACGCTGTTTTTTCGGTAATTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTtgc                                             >  1:2258613/1‑69 (MQ=255)
               aCCCGAACCAGACGCTGTTTTTTCGGTAATTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTtgc                                             >  1:131832/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:760022/66‑1 (MQ=255)
                                        gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                     <  1:2801920/66‑1 (MQ=255)
                                        gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                     <  1:2808574/66‑1 (MQ=255)
                                        gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                     <  1:1429939/66‑1 (MQ=255)
                                        gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                     <  1:3150039/66‑1 (MQ=255)
                                        gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                     <  1:3172502/66‑1 (MQ=255)
                                        gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                     <  1:3333228/66‑1 (MQ=255)
                                        gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                     <  1:3420124/66‑1 (MQ=255)
                                        gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                     <  1:2832643/66‑1 (MQ=255)
                                        gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                     <  1:777811/66‑1 (MQ=255)
                                        gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                     <  1:1145575/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:998628/66‑1 (MQ=255)
                                        gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                     <  1:2516528/66‑1 (MQ=255)
                                        gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                     <  1:2095978/66‑1 (MQ=255)
                                        gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                     <  1:194372/66‑1 (MQ=255)
                                        gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                     <  1:2133725/66‑1 (MQ=255)
                                        gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                     <  1:1919569/66‑1 (MQ=255)
                                        gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                     <  1:2399941/66‑1 (MQ=255)
                                        gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                     <  1:1874225/66‑1 (MQ=255)
                                        gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                     <  1:1873320/66‑1 (MQ=255)
                                        gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                     <  1:1074309/66‑1 (MQ=255)
                                        gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                     <  1:2598956/66‑1 (MQ=255)
                                        gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                     <  1:1639766/66‑1 (MQ=255)
                                        gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                     <  1:1612705/66‑1 (MQ=255)
                                        gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                     <  1:1579033/66‑1 (MQ=255)
                                        gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                     <  1:2751448/66‑1 (MQ=255)
                                        gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                     <  1:2751734/66‑1 (MQ=255)
                                        gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                     <  1:2755457/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)
                                                          tGCACGCATTTCTTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCt  >  1:941051/1‑69 (MQ=255)
                                                          tGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGc   >  1:1402848/1‑68 (MQ=255)
                                                          tGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCt  >  1:629692/1‑69 (MQ=255)
                                                          |                                                                    
AACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCGACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCT  >  minE/2854606‑2854732

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: