Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,640,124 G→A 35.1% T41T (ACC→ACT yhhA ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,640,1240GA35.1% 50.3 / 38.7 71T41T (ACC→ACTyhhAconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (8/38);  new base A (2/23);  total (10/61)
Fisher's exact test for biased strand distribution p-value = 4.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.58e-01

TGATTCAGCATCCCTTTTTGCTGGATTTGCTGAGTCTGCATCTGGGTTTGCATCCGCTGCTGGCTGGGTATCTGATACCCCGGCTGGTTAGGGTTGTTCAGAGTATTAATG  >  minE/2640080‑2640190
                                            |                                                                  
tGATTCAGCATCCCTTTTTGCTGGATTTGCTGAGTCTGCATCTGGGTTTGCATCCGCTGCTGGCTggg                                             <  1:2566505/68‑1 (MQ=255)
tGATTCAGCATCCCTTTTTGCTGGATTTGCTGAGTCTGCATCTGGGTTTGCATCCGCTGCTGGCTggg                                             <  1:1538405/68‑1 (MQ=255)
 ggttCAGCATCCCTTTTTGCTGGATTTGCTGCGTCTGCATTTGAGTTTGCATCCGCTGCTGGCTgggg                                            <  1:2613744/66‑2 (MQ=255)
 ggttCAGCATCCCTTTTTGCTGGATTTGCTGCGTCTGCATTTGAGTTTGCATCCGCTGCTGGCTgggg                                            <  1:1929639/66‑2 (MQ=255)
 ggttCAGCATCCCTTTTTGCTGGATTTGCTGCGTCTGCATTTGAGTTTGCATCCGCTGCTGGCTgggg                                            <  1:3392958/66‑2 (MQ=255)
 ggttCAGCATCCCTTTTTGCTGGATTTGCTGCGTCTGCATTTGAGTTTGCATCCGCTGCTGGCTgggg                                            <  1:1663121/66‑2 (MQ=255)
 ggttCAGCATCCCTTTTTGCTGGATTTGCTGCGTCTGCATTTGAGTTTGCATCCGCTGCTGGCTgggg                                            <  1:142298/66‑2 (MQ=255)
 ggttCAGCATCCCTTTTTGCTGGATTTGCTGCGTCTGCATTTGAGTTTGCATCCGCTGCTGGCTgggg                                            <  1:943007/66‑2 (MQ=255)
 ggttCAGCATCCCTTTTTGCTGGATTTGCTGCGTCTGCATTTGAGTTTGCATCCGCTGCTGGCTgggg                                            <  1:1244942/66‑2 (MQ=255)
 gATTCAGCATCCCTTTTTGCTGGATTTGCTGAGTCTGCATCTGGGTTTGCATCCGCTGCTGTCTGGGTa                                           <  1:2587048/69‑1 (MQ=255)
 gATTCAGCATCCCTTTTTGCTGGATTTGCTGAGTCTGCATCTGGGTTTGCATCCGCTGCTGGCTGGGTa                                           <  1:304885/69‑1 (MQ=255)
 gATTCAGCATCCCTTTTTGCTGGATTTGCTGAGTCTGCATCTGGGTTTGCATCCGCTGCTGGCTGGGTa                                           <  1:3135859/69‑1 (MQ=255)
 gATTCAGCATCCCTTTTTGCTGGATTTGCTGAGTCTGCATCTGGGTTTGCATCCGCTGCTGGCTGGGTa                                           <  1:3398275/69‑1 (MQ=255)
 gATTCAGCATCCCTTTTTGCTGGATTTGCTGAGTCTGCATCTGGGTTTGCATCCGCTGCTGGCTGGGTa                                           <  1:3218226/69‑1 (MQ=255)
 gATTCAGCATCCCTTTTTGCTGGATTTGCTGAGTCTGCATCTGGGTTTGCATCCGCTGCTGGCTGGGTa                                           <  1:2845343/69‑1 (MQ=255)
 gATTCAGCATCCCTTTTTGCTGGATTTGCTGAGTCTGCATCTGGGTTTGCATCCGCTGCTGGCTGGGTa                                           <  1:554418/69‑1 (MQ=255)
 gATTCAGCATCCCTTTTTGCTGGATTTGCTGAGTCTGCATCTGGGTTTGCATCCGCTGCTGGCTGGGTa                                           <  1:1957502/69‑1 (MQ=255)
 gATTCAGCATCCCTTTTTGCTGGATTTGCTGAGGCTGCATCTGGGTTTGCATCCGCTGCTGGCTGGGTa                                           <  1:526671/69‑1 (MQ=255)
     cAGCATCCCTTTTTGCTGGATTTGCTGAGTCTGCATCTGGGTTTGCATCCGCTGCTGGCTGGGTATCt                                        >  1:611344/1‑68 (MQ=255)
       gCATCCCTTTTTGCTGGATTTGCTGAGTCTGCATCTGGGTTTGCATCCGCTGCTGGCTGGGTATCTGa                                      <  1:391274/68‑1 (MQ=255)
       gCATCCCTTTTTGCTGGATTTGCTGAGTCTGCATCTGGGTTTGCATCCGCTGCTGGCTGGGTATCTGa                                      <  1:577602/68‑1 (MQ=255)
       gCATCCCTTTTTGCTGGATTTGCTGAGTCTGCATCTGGGTTTGCATCCGCTGCTGGCTGGGTATCTGa                                      <  1:2209446/68‑1 (MQ=255)
       gCATCCCTTTTTGCTGGATTTGCTGAGTCTGCATCTGGGTTTGCATCCGCTGCTGGCTGGGTATCTGa                                      <  1:421790/68‑1 (MQ=255)
       gCATCCCTTTTTGCTGGATTTGCTGAGTCTGCATCTGGGTTTGCATCCGCTGCTGGCTGGGTATCTGa                                      <  1:2836901/68‑1 (MQ=255)
       gCATCCCTTTTTGCTGGATTTGCTGAGTCTGCATCTGGGTTTGCATCCGCTGCTGGCTGGGTATCTGa                                      <  1:2934609/68‑1 (MQ=255)
       gCATCCCTTTTTGCTGGATTTGCTGAGTCTGCATCTGGGTTTGCATCCGCTGCTGGCTGGGTATCTGa                                      <  1:3191934/68‑1 (MQ=255)
       gCATCCCTTTTTGCTGGATTTGCTGAGTCTGCATCTGGGTTTGCATCCGCTGCTGGCTGGGTATCTGa                                      <  1:1569732/68‑1 (MQ=255)
       gCATCCCTTTTTGCTGGATTTGCTGAGTCTGCATCTGGGTTTGCATCCGCTGCTGGCTGGGTATCTGa                                      <  1:3112868/68‑1 (MQ=255)
        cATCCCTTTTTGCTGGATTTGCTGAGTCTGCATCTGGGTTTGCATCCGCTGCTGGCTGGGTATCTGa                                      <  1:996541/67‑1 (MQ=255)
         aTCCCTTTTTGCTGGATTTGCTGCGTCTGCATTTGAGTTTGCATCCGCTGCTGGCTGGGGATCTGATat                                   >  1:3111371/1‑68 (MQ=255)
         aTCCCTTTTTGCTGGATTTGCTGCGTCTGCATTTGAGTTTGCATCCGCTGCTGGCTGGGGATCTGATat                                   >  1:469069/1‑68 (MQ=255)
            ccTTTTTGCTGGATTTGCTGAGTCTGCATCTGGGTTTGCATCCGCTGCTGGCTGGGTATCTGATAcccc                                >  1:1925114/1‑69 (MQ=255)
            ccTTTTTGCTGGATTTGCTGAGTCTGCATCTGGGTTTGCATCCGCTGCTGGCTGGGTATCTGATAccc                                 >  1:1873471/1‑68 (MQ=255)
                     tGGATTTGCTGAGTCTGCATCTGGGTTTGCATCCGCTGCTGGCTGGGTATCTGATACCCCGGCTGGtt                        <  1:2673144/68‑1 (MQ=255)
                     tGGATTTGCTGAGTCTGCATCTGGGTTTGCATCCGCTGCTGGCTGGGTATCTGATACCCCGGCTGGtt                        <  1:3236871/68‑1 (MQ=255)
                           tgctgcGTCTGCATTTGAGTTTGCATCCGCTGCTGGCTGGGGATCTGATATCCCGGCTGGTTAGGgttg                 <  1:351757/69‑1 (MQ=255)
                           tgctgcGTCTGCATTTGAGTTTGCATCCGCTGCTGGCTGGGGATCTGATATCCCGGCTGGTTAGGgttg                 <  1:628501/69‑1 (MQ=255)
                           tgctgcGTCTGCATTTGAGTTTGCATCCGCTGCTGGCTGGGGATCTGATATCCCGGCTGGTTAGGgttg                 <  1:632429/69‑1 (MQ=255)
                           tgctgcGTCTGCATTTGAGTTTGCATCCGCTGCTGGCTGGGGATCTGATATCCCGGCTGGTTAGGgttg                 <  1:774175/69‑1 (MQ=255)
                           tgctgcGTCTGCATTTGAGTTTGCATCCGCTGCTGGCTGGGGATCTGATATCCCGGCTGGTTAGGgttg                 <  1:901125/69‑1 (MQ=255)
                           tgctgcGTCTGCATTTGAGTTTGCATCCGCTGCTGGCTGGGGATCTGATATCCCGGCTGGTTAGGgttg                 <  1:966463/69‑1 (MQ=255)
                           tgctgcGTCTGCATTTGAGTTTGCATCCGCTGCTGGCTGGGGATCTGATATCCCGGCTGGTTAGGgttg                 <  1:2852326/69‑1 (MQ=255)
                           tgctgcGTCTGCATTTGAGTTTGCATCCGCTGCTGGCTGGGGATCTGATATCCCGGCTGGTTAGGgttg                 <  1:3016842/69‑1 (MQ=255)
                           tgctgcGTCTGCATTTGAGTTTGCATCCGCTGCTGGCTGGGGATCTGATATCCCGGCTGGTTAGGgttg                 <  1:1376178/69‑1 (MQ=255)
                           tgctgcGTCTGCATTTGAGTTTGCATCCGCTGCTGGCTGGGGATCTGATATCCCGGCTGGTTAGGgttg                 <  1:1567492/69‑1 (MQ=255)
                           tgctgcGTCTGCATTTGAGTTTGCATCCGCTGCTGGCTGGGGATCTGATATCCCGGCTGGTTAGGgttg                 <  1:1784800/69‑1 (MQ=255)
                           tgctgcGTCTGCATTTGAGTTTGCATCCGCTGCTGGCTGGGGATCTGATATCCCGGCTGGTTAGGgttg                 <  1:221533/69‑1 (MQ=255)
                           tgctgcGTCTGCATTTGAGTTTGCATCCGCTGCTGGCTGGGGATCTGATATCCCGGCTGGTTAGGgttg                 <  1:2078355/69‑1 (MQ=255)
                           tgctgcGTCTGCATTTGAGTTTGCATCCGCTGCTGGCTGGGGATCTGATATCCCGGCTGGTTAGGgttg                 <  1:1957197/69‑1 (MQ=255)
                            gCTGCGTCTGCATTTGAGTTTGCATCCGCTGCTGGCTGGGGATCTGATATCCCGGCTGGTTAGGgttg                 <  1:2982622/68‑1 (MQ=255)
                            gCTGCGTCTGCATTTGAGTTTGCATCCGCTGCTGGCTGGGGATCTGATATCCCGGCTGGTTAGGgttg                 <  1:356663/68‑1 (MQ=255)
                                  tCTTCATCTGGGTTTGCATCCGCTGCTGGCTGGGTATCTGATACCCCGGCTGGTTAGGGTTGTTCag            <  1:2043049/67‑1 (MQ=255)
                                  tCTGCATCTGGGTTTGCATCCGCTGCTGGCTGGGTATCTGATACCCCGGCTGGTTAGGGTTGTTCagag          <  1:1981044/69‑1 (MQ=255)
                                  tCTGCATCTGGGTTTGCATCCGCTGCTGGCTGGGTATCTGATACCCCGGCTGGTTAGGGTTGTTCagag          <  1:1135163/69‑1 (MQ=255)
                                  tCTGCATCTGGGTTTGCATCCGCTGCTGGCTGGGTATCTGATACCCCGGCTGGTTAGGGTTGTTCagag          <  1:843725/69‑1 (MQ=255)
                                  tCTGCATCTGGGTTTGCATCCGCTGCTGGCTGGGTATCTGATACCCCGGCTGGTTAGGGTTGTTCagag          <  1:722093/69‑1 (MQ=255)
                                  tCTGCATCTGGGTTTGCATCCGCTGCTGGCTGGGTATCTGATACCCCGGCTGGTTAGGGTTGTTCagag          <  1:1914880/69‑1 (MQ=255)
                                  tCTGCATCTGGGTTTGCATCCGCTGCTGGCTGGGTATCTGATACCCCGGCTGGTTAGGGTTGTTCagag          <  1:269696/69‑1 (MQ=255)
                                  tCTGCATCTGGGTTTGCATCCGCTGCTGGCTGGGTATCTGATACCCCGGCTGGTTAGGGTTGTTCagag          <  1:3241512/69‑1 (MQ=255)
                                  tCTGCATCTGGGTTTGCATCCGCTGCTGGCTGGGTATCTGATACCCCGGCTGGTTAGGGTTGTTCagag          <  1:323306/69‑1 (MQ=255)
                                  tCTGCATCTGGGTTTGCATCCGCTGCTGGCTGGGTATCTGATACCCCGGCTGGTTAGGGTTGTTCag            <  1:2265280/67‑1 (MQ=255)
                                  tCTGCATCTGGGTTTGCATCCGCTGCTGGCTGGGTATCTGATACCCCGGCTGGTTAGGGTTGTTCag            <  1:2748269/67‑1 (MQ=255)
                                    tGCATCTGGGTTTGCATCCGCTGCTGGCTGGGTATCTGATACCCCGGCTGGTTAGGGTTGTTCAGAGt         >  1:3080108/1‑68 (MQ=255)
                                        tCTGGGTTTGCATCCGCTGCTGGCTGGGTATCTGATACCCCGGCTGGTTAGGGTTGTTCAGAGTAtt      <  1:1481221/67‑1 (MQ=255)
                                        tCTGGGTTTGCATCCGCTGCTGGCTGGGTATCTGATACCCCGGCTGGTTAGGGTTGTTCAGAGTATTa     <  1:1675965/68‑1 (MQ=255)
                                        tCTGGGTTTGCATCCGCTGCTGGCTGGGTATCTGATACCCCGGCTGGTTAGGGTTGTTCAGAGTATTa     <  1:674818/68‑1 (MQ=255)
                                        tCTGGGTTTGCATCCGCTGCTGGCTGGGTATCTGATACCCCGGCTGGTTAGGGTTGTTCAGAGTATTa     <  1:872650/68‑1 (MQ=255)
                                          tGGGTTTGCATCCGCTGCTGGCTGGGTATCTGATa                                    >  1:1627916/1‑35 (MQ=255)
                                          tGGGTTTGCATCCGCTGCTGGCTGGGTATCTGATACCCCGGCTGGTTAGGGTTGTTCAGAGTATTAATg  >  1:2350791/1‑69 (MQ=255)
                                          tGGGTTTGCATCCGCTGCTGGCTGGGTATCTGATACCCCGGCTGGTTAGGGTTGTTCAGAGTATTAATg  >  1:1117413/1‑69 (MQ=255)
                                          tGGGTTTGCATCCGCTGCTGGCTGGGTATCTGATACCCCGGCTGGTTAGGGTTGTTCAGAGTATTAATg  >  1:1132880/1‑69 (MQ=255)
                                            |                                                                  
TGATTCAGCATCCCTTTTTGCTGGATTTGCTGAGTCTGCATCTGGGTTTGCATCCGCTGCTGGCTGGGTATCTGATACCCCGGCTGGTTAGGGTTGTTCAGAGTATTAATG  >  minE/2640080‑2640190

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: