Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,022,941 A→G 38.3% C322C (TGT→TGC ygiQ ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,022,9410AG38.3% 35.7 / 47.4 73C322C (TGT→TGCygiQconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (10/35);  new base G (9/19);  total (19/54)
Fisher's exact test for biased strand distribution p-value = 4.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.34e-01

GGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCAC  >  minE/2022891‑2023009
                                                  |                                                                    
gtgTTGATCCATACATAGCGGTCGCCATGTTTTTGCATCAATGCGCGGGCGCAGCCTGGGTTAGTTTCg                                                    <  1:2842892/67‑1 (MQ=255)
gggTTGATCCATACATAGCGGTCGCCATGTTTTTGCATCAATGCGCGGGCGCAGCCTGGGTTAGTTTCg                                                    <  1:657783/69‑1 (MQ=255)
gggTTGATCCATACATAGCGGTCGCCATGTTTTTGCATCAATGCGCGGGCGCAGCCTGGGTTAGTTTCg                                                    <  1:1812754/69‑1 (MQ=255)
gggTTGATCCATACATAGCGGTCGCCATGTTTTTGCATCAATGCGCGGGCGCAGCCTGGGTTAGTTTCg                                                    <  1:2840562/69‑1 (MQ=255)
gggTTGATCCATACATAGCGGTCGCCATGTTTTTGCATCAATGCGCGGGCGCAGCCTGGGTTAGTTTCg                                                    <  1:2816681/69‑1 (MQ=255)
gggTTGATCCATACATAGCGGTCGCCATGTTTTTGCATCAATGCGCGGGCGCAGCCTGGGTTAGTTTCg                                                    <  1:1969538/69‑1 (MQ=255)
gggTTGATCCATACATAGCGGTCGCCATGTTTTTGCATCAATGCGCGGGCGCAGCCTGGGTTAGTTTCg                                                    <  1:161787/69‑1 (MQ=255)
gggTTGATCCATACATAGCGGTCGCCATGTTTTTGCATCAATGCGCGGGCGCAGCCTGGGTTAGTTTCg                                                    <  1:337131/69‑1 (MQ=255)
gggTTGATCCATACATAGCGGTCGCCATGTTTTTGCATCAATGCGCGGGCGCAGCCTGGGTTAGTTTCg                                                    <  1:2687296/69‑1 (MQ=255)
gggTTGATCCATACATAGCGGTCGCCATGTTTTTGCATCAATGCGCGGGCGCAGCCTGGGTTAGTTTCg                                                    <  1:1591747/69‑1 (MQ=255)
gggTTGATCCATACATAGCGGTCGCCATGTTTTTGCATCAATGCGCGGGCGCAGCCTGGGTTAGTTTCg                                                    <  1:398853/69‑1 (MQ=255)
gggTTGATCCATACATAGCGGTCGCCATGTTTTTGCATCAATGCGCGGGCGCAGCCTGGGTTAGTTTCg                                                    <  1:661782/69‑1 (MQ=255)
gggTTGATCCATACATAGCGGTCGCCATGTTTTTGCATCAATGCGCGGGCGCAGCCTGGGTTAGTTTCg                                                    <  1:688120/69‑1 (MQ=255)
gggTTGATCCATACATAGCGGTCGCCATGTTTTTGCATCAATGCGCGGGCGCAGCCTGGGTTAGTTTCg                                                    <  1:1388037/69‑1 (MQ=255)
gggTTGATCCATACATAGCGGTCGCCATGTTTTTGCATCAATGCGCGGGCGCAGCCTGGGTTAGTTTCg                                                    <  1:769479/69‑1 (MQ=255)
gggTTGATCCATACATAGCGGTCGCCATGTTTTTGCATCAATGCGCGGGCGCAGCCTGGGTTAGTTTCg                                                    <  1:2363335/69‑1 (MQ=255)
gggTTGATCCATACATAGCGGTCGCCATGTTTTTGCATCAATGCGCGGGCGCAGCCTGGGTTAGTTTCg                                                    <  1:2086360/69‑1 (MQ=255)
gggTTGATCCATACATAGCGGTCGCCATGTTTTTGCATCAATGCGCGGGCGCAGCCTGGGTTAGTTTCg                                                    <  1:866687/69‑1 (MQ=255)
gggTTGATCCATACATAGCGGTCGCCATGTTTTTGCATCAATGCGCGGGCGCAGCCTGGGTTAGTTTCg                                                    <  1:20975/69‑1 (MQ=255)
gggTTGATCCACACATAGCTGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGggttggt                                                        <  1:1585849/65‑1 (MQ=255)
gggTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGggttggt                                                        <  1:3134171/65‑1 (MQ=255)
gggTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGggttggt                                                        <  1:2558357/65‑1 (MQ=255)
gggTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGggttggt                                                        <  1:2392297/65‑1 (MQ=255)
gggTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGggttggt                                                        <  1:2924466/65‑1 (MQ=255)
gggTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGggttggt                                                        <  1:1602754/65‑1 (MQ=255)
gggTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGggttggt                                                        <  1:1430781/65‑1 (MQ=255)
gggTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGggttggt                                                        <  1:820201/65‑1 (MQ=255)
 ggTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGggttggt                                                        <  1:794522/64‑1 (MQ=255)
        ccACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAg                                            <  1:895340/69‑1 (MQ=255)
            acaTAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATa                                        <  1:1769760/69‑1 (MQ=255)
            acaTAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATa                                        <  1:1755383/69‑1 (MQ=255)
            acaTAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATa                                        <  1:2734385/69‑1 (MQ=255)
            acaTAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATa                                        <  1:1451155/69‑1 (MQ=255)
            acaTAGCGGTCGCCATGTTTTTGCATCAATGCGCGGGCGCAGCCTGggtta                                                          >  1:286992/1‑50 (MQ=38)
            acaTAGCGGTCGCCATGTTTTTGCATCAATGCGCGGGCGCAGCCTGggtta                                                          >  1:1086223/1‑50 (MQ=38)
            acaTAGCGGTCGCCATGTTTTTGCATCAATGCGCGGGCGCAGCCTGggtta                                                          >  1:2257181/1‑50 (MQ=38)
            acaTAGCGGTCGCCATGTTTTTGCATCAATGCGCGGGCGCAGCCTGggtta                                                          >  1:1379510/1‑50 (MQ=38)
            acaTAGCGGTCGCCATGTTTTTGCATCAATGCGCGGGCGCAGCCTGggtta                                                          >  1:1466352/1‑50 (MQ=38)
            acaTAGCGGTCGCCATGTTTTTGCATCAATGCGCGGGCGCAGCCTGggtta                                                          >  1:147416/1‑50 (MQ=38)
            acaTAGCGGTCGCCATGTTTTTGCATCAATGCGCGGGCGCAGCCTGggtta                                                          >  1:1563738/1‑50 (MQ=38)
            acaTAGCGGTCGCCATGTTTTTGCATCAATGCGCGGGCGCAGCCTGggtta                                                          >  1:1714260/1‑50 (MQ=38)
            acaTAGCGGTCGCCATGTTTTTGCATCAATGCGCGGGCGCAGCCTGggtta                                                          >  1:1716454/1‑50 (MQ=38)
                 gCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAg                                            >  1:2069125/1‑60 (MQ=255)
                 gCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAg                                            >  1:960381/1‑60 (MQ=255)
                 gCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAg                                            >  1:3299418/1‑60 (MQ=255)
                 gCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAg                                            >  1:310085/1‑60 (MQ=255)
                 gCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAg                                            >  1:1903925/1‑60 (MQ=255)
                 gCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAg                                            >  1:2471142/1‑60 (MQ=255)
                   ggTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAg                                 <  1:7740/69‑1 (MQ=255)
                   ggTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAg                                 <  1:453399/69‑1 (MQ=255)
                   ggTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAg                                 <  1:1613834/69‑1 (MQ=255)
                   ggTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAg                                 <  1:3217155/69‑1 (MQ=255)
                    gTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGc                                <  1:1971212/69‑1 (MQ=255)
                    gTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGc                                <  1:3235737/69‑1 (MQ=255)
                    gTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGc                                <  1:3196709/69‑1 (MQ=255)
                     tCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGc                                <  1:2974537/68‑1 (MQ=255)
                     tCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGc                                <  1:1950271/68‑1 (MQ=255)
                       gCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTg                                                               >  1:3242510/1‑35 (MQ=255)
                       gCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATa                                        >  1:75015/1‑58 (MQ=255)
                       gCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATg                             >  1:1159902/1‑69 (MQ=255)
                       gCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATg                             >  1:2361747/1‑69 (MQ=255)
                           tgtTTTTGCATCAATGCGCTGGCACAGCCTGGGTTGGTTTCGTTGTGCAGAATACGCGAAGCATGGGCg                         <  1:1972485/69‑1 (MQ=255)
                           tgtTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCATAATACGCGAAGCATGGGCg                         <  1:528051/69‑1 (MQ=255)
                           tgtTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCg                         <  1:2249340/69‑1 (MQ=255)
                           tgtTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCg                         <  1:2882664/69‑1 (MQ=255)
                           tgtTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCg                         <  1:1228851/69‑1 (MQ=255)
                           tgtTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCg                         <  1:759979/69‑1 (MQ=255)
                           tgtTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCg                         <  1:737919/69‑1 (MQ=255)
                             tttttGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTa                       <  1:3068084/69‑1 (MQ=255)
                             tttttGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTa                       <  1:635350/69‑1 (MQ=255)
                             tttttGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTa                       <  1:1574798/69‑1 (MQ=255)
                             tttttGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTa                       <  1:3278566/69‑1 (MQ=255)
                             tttttGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTa                       <  1:2613201/69‑1 (MQ=255)
                                                  acaGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGc         <  1:579663/62‑1 (MQ=255)
                                                  acaGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGc         <  1:604175/62‑1 (MQ=255)
                                                  acaGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGc         <  1:2064350/62‑1 (MQ=255)
                                                  acaGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGc         <  1:1834098/62‑1 (MQ=255)
                                                  acaGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGc         <  1:1829965/62‑1 (MQ=255)
                                                  acaGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGc         <  1:1322065/62‑1 (MQ=255)
                                                  acaGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGc         <  1:259401/62‑1 (MQ=255)
                                                  acaGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGc         <  1:1021764/62‑1 (MQ=255)
                                                  acaGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGc         <  1:1701119/62‑1 (MQ=255)
                                                  acaGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCAc  <  1:2820006/69‑1 (MQ=255)
                                                  acaGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCAc  <  1:2419129/69‑1 (MQ=255)
                                                  acaGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCAc  <  1:2012668/69‑1 (MQ=255)
                                                  acaGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCAc  <  1:3233790/69‑1 (MQ=255)
                                                  |                                                                    
GGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCAC  >  minE/2022891‑2023009

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: