Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,640,121 C→T 36.1% Q42Q (CAG→CAA yhhA ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,640,1210CT36.1% 45.4 / 39.1 69Q42Q (CAG→CAAyhhAconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (6/38);  new base T (2/23);  total (8/61)
Fisher's exact test for biased strand distribution p-value = 7.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.68e-01

GCAACTGAGTTTGCGTTTTCAGTTGCTGATTCAGCATCCCTTTTTGCTGGATTTGCTGAGTCTGCATCTGGGTTTGCATCCGCTGCTGGCTGGGTATCTGATACCCCGGCTGGTTAGGGTTGTTCAGAGTATTA  >  minE/2640054‑2640187
                                                                   |                                                                  
gCAACTGAGTTTGCGTTTTCAGTTGCTGATTCAGCATCCCTTTTTGCTGGATTTGCTGAGTCTGCATCt                                                                   >  1:1365176/1‑69 (MQ=255)
gCAACTGAGTTTGCGTTTTCAGTTGCTGATTCAGCATCCCTTTTTGCTGGATTTGCTGAGTCTGCATCt                                                                   >  1:157426/1‑69 (MQ=255)
                          tGATTCAGCATCCCTTTTTGCTGGATTTGCTGAGTCTGCATCTGGGTTTGCATCCGCTGCTGGCTggg                                          <  1:2566505/68‑1 (MQ=255)
                          tGATTCAGCATCCCTTTTTGCTGGATTTGCTGAGTCTGCATCTGGGTTTGCATCCGCTGCTGGCTggg                                          <  1:1538405/68‑1 (MQ=255)
                           ggttCAGCATCCCTTTTTGCTGGATTTGCTGCGTCTGCATTTGAGTTTGCATCCGCTGCTGGCTgggg                                         <  1:1244942/66‑2 (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:943007/66‑2 (MQ=255)
                           ggttCAGCATCCCTTTTTGCTGGATTTGCTGCGTCTGCATTTGAGTTTGCATCCGCTGCTGGCTgggg                                         <  1:142298/66‑2 (MQ=255)
                           gATTCAGCATCCCTTTTTGCTGGATTTGCTGAGTCTGCATCTGGGTTTGCATCCGCTGCTGTCTGGGTa                                        <  1:2587048/69‑1 (MQ=255)
                           gATTCAGCATCCCTTTTTGCTGGATTTGCTGAGTCTGCATCTGGGTTTGCATCCGCTGCTGGCTGGGTa                                        <  1:554418/69‑1 (MQ=255)
                           gATTCAGCATCCCTTTTTGCTGGATTTGCTGAGTCTGCATCTGGGTTTGCATCCGCTGCTGGCTGGGTa                                        <  1:3398275/69‑1 (MQ=255)
                           gATTCAGCATCCCTTTTTGCTGGATTTGCTGAGTCTGCATCTGGGTTTGCATCCGCTGCTGGCTGGGTa                                        <  1:3218226/69‑1 (MQ=255)
                           gATTCAGCATCCCTTTTTGCTGGATTTGCTGAGTCTGCATCTGGGTTTGCATCCGCTGCTGGCTGGGTa                                        <  1:1957502/69‑1 (MQ=255)
                           gATTCAGCATCCCTTTTTGCTGGATTTGCTGAGTCTGCATCTGGGTTTGCATCCGCTGCTGGCTGGGTa                                        <  1:3135859/69‑1 (MQ=255)
                           gATTCAGCATCCCTTTTTGCTGGATTTGCTGAGTCTGCATCTGGGTTTGCATCCGCTGCTGGCTGGGTa                                        <  1:2845343/69‑1 (MQ=255)
                           gATTCAGCATCCCTTTTTGCTGGATTTGCTGAGTCTGCATCTGGGTTTGCATCCGCTGCTGGCTGGGTa                                        <  1:304885/69‑1 (MQ=255)
                           gATTCAGCATCCCTTTTTGCTGGATTTGCTGAGGCTGCATCTGGGTTTGCATCCGCTGCTGGCTGGGTa                                        <  1:526671/69‑1 (MQ=255)
                               cAGCATCCCTTTTTGCTGGATTTGCTGAGTCTGCATCTGGGTTTGCATCCGCTGCTGGCTGGGTATCt                                     >  1:611344/1‑68 (MQ=255)
                                 gCATCCCTTTTTGCTGGATTTGCTGAGTCTGCATCTGGGTTTGCATCCGCTGCTGGCTGGGTATCTGa                                   <  1:2934609/68‑1 (MQ=255)
                                 gCATCCCTTTTTGCTGGATTTGCTGAGTCTGCATCTGGGTTTGCATCCGCTGCTGGCTGGGTATCTGa                                   <  1:2836901/68‑1 (MQ=255)
                                 gCATCCCTTTTTGCTGGATTTGCTGAGTCTGCATCTGGGTTTGCATCCGCTGCTGGCTGGGTATCTGa                                   <  1:3112868/68‑1 (MQ=255)
                                 gCATCCCTTTTTGCTGGATTTGCTGAGTCTGCATCTGGGTTTGCATCCGCTGCTGGCTGGGTATCTGa                                   <  1:3191934/68‑1 (MQ=255)
                                 gCATCCCTTTTTGCTGGATTTGCTGAGTCTGCATCTGGGTTTGCATCCGCTGCTGGCTGGGTATCTGa                                   <  1:391274/68‑1 (MQ=255)
                                 gCATCCCTTTTTGCTGGATTTGCTGAGTCTGCATCTGGGTTTGCATCCGCTGCTGGCTGGGTATCTGa                                   <  1:421790/68‑1 (MQ=255)
                                 gCATCCCTTTTTGCTGGATTTGCTGAGTCTGCATCTGGGTTTGCATCCGCTGCTGGCTGGGTATCTGa                                   <  1:577602/68‑1 (MQ=255)
                                 gCATCCCTTTTTGCTGGATTTGCTGAGTCTGCATCTGGGTTTGCATCCGCTGCTGGCTGGGTATCTGa                                   <  1:1569732/68‑1 (MQ=255)
                                 gCATCCCTTTTTGCTGGATTTGCTGAGTCTGCATCTGGGTTTGCATCCGCTGCTGGCTGGGTATCTGa                                   <  1:2209446/68‑1 (MQ=255)
                                  cATCCCTTTTTGCTGGATTTGCTGAGTCTGCATCTGGGTTTGCATCCGCTGCTGGCTGGGTATCTGa                                   <  1:996541/67‑1 (MQ=255)
                                   aTCCCTTTTTGCTGGATTTGCTGCGTCTGCATTTGAGTTTGCATCCGCTGCTGGCTGGGGATCTGATat                                >  1:469069/1‑68 (MQ=255)
                                   aTCCCTTTTTGCTGGATTTGCTGCGTCTGCATTTGAGTTTGCATCCGCTGCTGGCTGGGGATCTGATat                                >  1:3111371/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:1567492/69‑1 (MQ=255)
                                                     tgctgcGTCTGCATTTGAGTTTGCATCCGCTGCTGGCTGGGGATCTGATATCCCGGCTGGTTAGGgttg              <  1:628501/69‑1 (MQ=255)
                                                     tgctgcGTCTGCATTTGAGTTTGCATCCGCTGCTGGCTGGGGATCTGATATCCCGGCTGGTTAGGgttg              <  1:1784800/69‑1 (MQ=255)
                                                     tgctgcGTCTGCATTTGAGTTTGCATCCGCTGCTGGCTGGGGATCTGATATCCCGGCTGGTTAGGgttg              <  1:351757/69‑1 (MQ=255)
                                                     tgctgcGTCTGCATTTGAGTTTGCATCCGCTGCTGGCTGGGGATCTGATATCCCGGCTGGTTAGGgttg              <  1:632429/69‑1 (MQ=255)
                                                     tgctgcGTCTGCATTTGAGTTTGCATCCGCTGCTGGCTGGGGATCTGATATCCCGGCTGGTTAGGgttg              <  1:1957197/69‑1 (MQ=255)
                                                     tgctgcGTCTGCATTTGAGTTTGCATCCGCTGCTGGCTGGGGATCTGATATCCCGGCTGGTTAGGgttg              <  1:774175/69‑1 (MQ=255)
                                                     tgctgcGTCTGCATTTGAGTTTGCATCCGCTGCTGGCTGGGGATCTGATATCCCGGCTGGTTAGGgttg              <  1:2078355/69‑1 (MQ=255)
                                                     tgctgcGTCTGCATTTGAGTTTGCATCCGCTGCTGGCTGGGGATCTGATATCCCGGCTGGTTAGGgttg              <  1:3016842/69‑1 (MQ=255)
                                                     tgctgcGTCTGCATTTGAGTTTGCATCCGCTGCTGGCTGGGGATCTGATATCCCGGCTGGTTAGGgttg              <  1:901125/69‑1 (MQ=255)
                                                     tgctgcGTCTGCATTTGAGTTTGCATCCGCTGCTGGCTGGGGATCTGATATCCCGGCTGGTTAGGgttg              <  1:2852326/69‑1 (MQ=255)
                                                     tgctgcGTCTGCATTTGAGTTTGCATCCGCTGCTGGCTGGGGATCTGATATCCCGGCTGGTTAGGgttg              <  1:221533/69‑1 (MQ=255)
                                                     tgctgcGTCTGCATTTGAGTTTGCATCCGCTGCTGGCTGGGGATCTGATATCCCGGCTGGTTAGGgttg              <  1:1376178/69‑1 (MQ=255)
                                                     tgctgcGTCTGCATTTGAGTTTGCATCCGCTGCTGGCTGGGGATCTGATATCCCGGCTGGTTAGGgttg              <  1:966463/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:722093/69‑1 (MQ=255)
                                                            tCTGCATCTGGGTTTGCATCCGCTGCTGGCTGGGTATCTGATACCCCGGCTGGTTAGGGTTGTTCagag       <  1:843725/69‑1 (MQ=255)
                                                            tCTGCATCTGGGTTTGCATCCGCTGCTGGCTGGGTATCTGATACCCCGGCTGGTTAGGGTTGTTCagag       <  1:1914880/69‑1 (MQ=255)
                                                            tCTGCATCTGGGTTTGCATCCGCTGCTGGCTGGGTATCTGATACCCCGGCTGGTTAGGGTTGTTCagag       <  1:3241512/69‑1 (MQ=255)
                                                            tCTGCATCTGGGTTTGCATCCGCTGCTGGCTGGGTATCTGATACCCCGGCTGGTTAGGGTTGTTCagag       <  1:323306/69‑1 (MQ=255)
                                                            tCTGCATCTGGGTTTGCATCCGCTGCTGGCTGGGTATCTGATACCCCGGCTGGTTAGGGTTGTTCagag       <  1:1981044/69‑1 (MQ=255)
                                                            tCTGCATCTGGGTTTGCATCCGCTGCTGGCTGGGTATCTGATACCCCGGCTGGTTAGGGTTGTTCagag       <  1:1135163/69‑1 (MQ=255)
                                                            tCTGCATCTGGGTTTGCATCCGCTGCTGGCTGGGTATCTGATACCCCGGCTGGTTAGGGTTGTTCagag       <  1:269696/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)
                                                                   |                                                                  
GCAACTGAGTTTGCGTTTTCAGTTGCTGATTCAGCATCCCTTTTTGCTGGATTTGCTGAGTCTGCATCTGGGTTTGCATCCGCTGCTGGCTGGGTATCTGATACCCCGGCTGGTTAGGGTTGTTCAGAGTATTA  >  minE/2640054‑2640187

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: