Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,854,683 T→G 52.8% I349I (ATA→ATC 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,6830TG52.8% 24.1 / 74.9 76I349I (ATA→ATCyjeSpredicted Fe‑S electron transport protein
Reads supporting (aligned to +/- strand):  ref base T (9/27);  new base G (9/31);  total (18/58)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCGACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGA  >  minE/2854617‑2854743
                                                                  |                                                            
aaTCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCGACTATGCACGCATTTCGTCTCTCTa                                                             <  1:1865837/68‑1 (MQ=255)
aaTCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCGACTATGCACGCATTTCGTCTCTCTa                                                             <  1:1891284/68‑1 (MQ=255)
aaTCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCGACTATGCACGCATTTCGTCTCTCTa                                                             <  1:794428/68‑1 (MQ=255)
aaTCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCGACTATGCACGCATTTCGTCTCTCTa                                                             <  1:2116875/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:1603856/68‑1 (MQ=255)
aaTCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCGACTATGCACGCATTTCGTCTCTCTa                                                             <  1:2552945/68‑1 (MQ=255)
 aTCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCGACTATGCACGCATTTCGTCTCTCTa                                                             <  1:1127586/67‑1 (MQ=255)
 aTCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCGACTATGCACGCATTTCGTCTCTCTa                                                             <  1:1586240/67‑1 (MQ=255)
 aTCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCGACTATGCACGCATTTCGTCTCTCTa                                                             <  1:835764/67‑1 (MQ=255)
  tCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCGACTATGCACGCATTTCGTCTCTCTAttt                                                          >  1:777116/1‑69 (MQ=255)
    aCCCGAACCAGACGCTGTTTTTTCGGTAATTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTtgc                                                        >  1:2642364/1‑69 (MQ=255)
    aCCCGAACCAGACGCTGTTTTTTCGGTAATTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTtgc                                                        >  1:2648339/1‑69 (MQ=255)
    aCCCGAACCAGACGCTGTTTTTTCGGTAATTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTtgc                                                        >  1:2472206/1‑69 (MQ=255)
    aCCCGAACCAGACGCTGTTTTTTCGGTAATTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTtgc                                                        >  1:2653249/1‑69 (MQ=255)
    aCCCGAACCAGACGCTGTTTTTTCGGTAATTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTtgc                                                        >  1:2766498/1‑69 (MQ=255)
    aCCCGAACCAGACGCTGTTTTTTCGGTAATTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTtgc                                                        >  1:296791/1‑69 (MQ=255)
    aCCCGAACCAGACGCTGTTTTTTCGGTAATTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTtgc                                                        >  1:2258613/1‑69 (MQ=255)
    aCCCGAACCAGACGCTGTTTTTTCGGTAATTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTtgc                                                        >  1:2568991/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:2751448/66‑1 (MQ=255)
                             gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                                <  1:1874225/66‑1 (MQ=255)
                             gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                                <  1:2751734/66‑1 (MQ=255)
                             gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                                <  1:2755457/66‑1 (MQ=255)
                             gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                                <  1:2801920/66‑1 (MQ=255)
                             gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                                <  1:2808574/66‑1 (MQ=255)
                             gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                                <  1:2832643/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:760022/66‑1 (MQ=255)
                             gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                                <  1:1145575/66‑1 (MQ=255)
                             gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                                <  1:777811/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:1873320/66‑1 (MQ=255)
                             gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                                <  1:1639766/66‑1 (MQ=255)
                             gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                                <  1:2133725/66‑1 (MQ=255)
                             gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                                <  1:1429939/66‑1 (MQ=255)
                             gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                                <  1:2399941/66‑1 (MQ=255)
                             gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                                <  1:1612705/66‑1 (MQ=255)
                             gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                                <  1:1579033/66‑1 (MQ=255)
                             gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                                <  1:2516528/66‑1 (MQ=255)
                             gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                                <  1:194372/66‑1 (MQ=255)
                             gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                                <  1:1074309/66‑1 (MQ=255)
                             gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                                <  1:2598956/66‑1 (MQ=255)
                             gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                                <  1:2095978/66‑1 (MQ=255)
                             gtaaTTGCACTTCGACCACGCACGCATTTCGCCTCTCGATTTGCTGCGCAATCGCCCACGCTATGTGc                                <  1:1919569/66‑1 (MQ=255)
                              cAATTGCACTTCGACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTc                              <  1:1523506/69‑1 (MQ=255)
                                  tGCACTTCGACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGc                                       >  1:1139063/1‑56 (MQ=255)
                                  tGCACTTCGACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGc                                       >  1:779006/1‑56 (MQ=255)
                                      cTTCGACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCgaga                       <  1:2974026/68‑1 (MQ=255)
                                      cTTCGACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCgaga                       <  1:773294/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)
                                                 cACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTc            <  1:659666/68‑1 (MQ=255)
                                                 cACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTc            <  1:652520/68‑1 (MQ=255)
                                                 cACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTc            <  1:2272015/68‑1 (MQ=255)
                                                 cACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTc            <  1:3038407/68‑1 (MQ=255)
                                                 cACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTc            <  1:1129183/68‑1 (MQ=255)
                                                 cACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTc            <  1:2646612/68‑1 (MQ=255)
                                                        ttCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCgaga                       <  1:1721678/50‑1 (MQ=255)
                                                         tCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACg   >  1:359246/1‑69 (MQ=255)
                                                         tCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACg   >  1:1479686/1‑69 (MQ=255)
                                                           gTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCgaga                       <  1:1538250/47‑1 (MQ=255)
                                                           gTCTCTATATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGa  <  1:271564/68‑1 (MQ=255)
                                                                  |                                                            
AATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCGACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGA  >  minE/2854617‑2854743

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: