Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,290,390 C→T 37.2% G74G (GGC→GGT yeeN → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,290,3900CT37.2% 69.6 / 68.8 110G74G (GGC→GGTyeeNconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (26/43);  new base T (21/20);  total (47/63)
Fisher's exact test for biased strand distribution p-value = 2.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.82e-01

Only 100 of 110 total aligned reads displayed.

CGTGCAAAGCAGGCACAAGTTCCAAAGCACGTTATTGATAAAGCAATTGATAAAGCCAAAGGCGGCGGAGATGAAACGTTCGTGCAGGGACGTTATGAAGGCTTTGGTCCTAATGGCTCAATG  >  minE/1290325‑1290447
                                                                 |                                                         
cGTGCAAAGCAGGCACAAGTTCCAAAGCACGTTATTGATAAAGCAATTGATAAAGCCAAAGGCGGCGga                                                        <  1:722324/69‑1 (MQ=255)
cGTGCAAAGCAGGCACAAGTTCCAAAGCACGTTATTGATAAAGCAATTGATAAAGCCAAAGGCGGCGga                                                        <  1:1990854/69‑1 (MQ=255)
         cAGGCACAAGTTCCAAAGCACGTTATTGATAAAGCAATTGATAAAGCCAAAGGCGGCGGAGATGAAACg                                               <  1:3436209/69‑1 (MQ=255)
         cAGGCACAAGTTCCAAAGCACGTTATTGATAAAGCAATTGATAAAGCCAAAGGCGGCGGAGATGAAACg                                               <  1:2891444/69‑1 (MQ=255)
         cAGGCACAAGTTCCAAAGCACGTTATTGATAAAGCAATTGATAAAGCCAAAGGCGGCGGAGATGAAACg                                               <  1:2578023/69‑1 (MQ=255)
         cAGGCACAAGTTCCAAAGCACGTTATTGATAAAGCAATTGATAAAGCCAAAGGCGGCGGAGATGAAACg                                               <  1:1066938/69‑1 (MQ=255)
         cAGGCACAAGTTCCAAAGCACGTTATTGATAAAGCAATTGATAAAGCCAAAGGCGGCGGAGATGAAACg                                               <  1:2373412/69‑1 (MQ=255)
         cAGGCACAAGTTCCAAAGCACGTTATTGATAAAGCAATTGATAAAGCCAAAGGCGGCGGAGATGAAACg                                               <  1:841144/69‑1 (MQ=255)
         cAGGCACAAGTTCCAAAGCACGTTATTGATAAAGCAATTGATAAAGCCAAAGGCGGCGGAGATGAAACg                                               <  1:2092818/69‑1 (MQ=255)
         cAGGCACAAGTTCCAAAGCACGTTATTGATAAAGCAATTGATAAAGCCAAAGGCGGCGGAGATGAAACg                                               <  1:1937309/69‑1 (MQ=255)
         cAGGCACAAGTTCCAAAGCACGTTATTGATAAAGCAATTGATAAAGCCAAAGGCGGCGGAGATGAAACg                                               <  1:1898713/69‑1 (MQ=255)
           tgCACAAGTTCCAAAGCACGTTATTGATAAAGCAATTGATAAAGCCAAAGGCGGTGGAGATGAAACGtt                                             >  1:1039562/2‑69 (MQ=255)
           ggCACAAGTTCCAAAGCACGTTATTGATAAAGCAATTGATAAAGCCAAAGGCGGTGGAGATGAAACGtt                                             >  1:1765582/1‑69 (MQ=255)
           ggCACAAGTTCCAAAGCACGTTATTGATAAAGCAATTGATAAAGCCAAAGGCGGTGGAGATGAAACGtt                                             >  1:1358913/1‑69 (MQ=255)
           ggCACAAGTTCCAAAGCACGTTATTGATAAAGCAATTGATAAAGCCAAAGGCGGTGGAGATGAAACGtt                                             >  1:2358453/1‑69 (MQ=255)
           ggCACAAGTTCCAAAGCACGTTATTGATAAAGCAATTGATAAAGCCAAAGGCGGTGGAGATGAAACGtt                                             >  1:812994/1‑69 (MQ=255)
           ggCACAAGTTCCAAAGCACGTTATTGATAAAGCAATTGATAAAGCCAAAGGCGGTGGAGATGAAACGtt                                             >  1:2151763/1‑69 (MQ=255)
           ggCACAAGTTCCAAAGCACGTTATTGATAAAGCAATTGATAAAGCCAAAGGCGGTGGAGATGAAACGtt                                             >  1:1188097/1‑69 (MQ=255)
           ggCACAAGTTCCAAAGCACGTTATTGATAAAGCAATTGATAAAGCCAAAGGCGGTGGAGATGAAACGtt                                             >  1:1937871/1‑69 (MQ=255)
           ggCACAAGTTCCAAAGCACGTTATTGATAAAGCAATTGATAAAGCCAAAGGCGGTGGAGATGAAACGtt                                             >  1:1239282/1‑69 (MQ=255)
           ggCACAAGTTCCAAAGCACGTTATTGATAAAGCAATTGATAAAGCCAAAGGCGGTGGAGATGAAACGtt                                             >  1:1922675/1‑69 (MQ=255)
           ggCACAAGTTCCAAAGCACGTTATTGATAAAGCAATTGATAAAGCCAAAGGCGGTGGAGATGAAACGtt                                             >  1:2593636/1‑69 (MQ=255)
           ggCACAAGTTCCAAAGCACGTTATTGATAAAGCAATTGATAAAGCCAAAGGCGGTGGAGATGAAACGtt                                             >  1:189802/1‑69 (MQ=255)
           ggCACAAGTTCCAAAGCACGTTATTGATAAAGCAATTGATAAAGCCAAAGGCGGTGGAGATGAAACGtt                                             >  1:277081/1‑69 (MQ=255)
           ggCACAAGTTCCAAAGCACGTTATTGATAAAGCAATTGATAAAGCCAAAGGCGGTGGAGATGAAACGtt                                             >  1:1668516/1‑69 (MQ=255)
           ggCACAAGTTCCAAAGCACGTTATTGATAAAGCAATTGATAAAGCCAAAGGCGGTGGAGATGAAACGtt                                             >  1:1667707/1‑69 (MQ=255)
           ggCACAAGTTCCAAAGCACGTTATTGATAAAGCAATTGATAAAGCCAAAGGCGGTGGAGATGAAACGtt                                             >  1:1393352/1‑69 (MQ=255)
           ggCACAAGTTCCAAAGCACGTTATTGATAAAGCAATTGATAAAGCCAAAGGCGGTGGAGATGAAACGtt                                             >  1:1434187/1‑69 (MQ=255)
           ggCACAAGTTCCAAAGCACGTTATTGATAAAGCAATTGATAAAGCCAAAGGCGGTGGAGATGAAACGtt                                             >  1:1466020/1‑69 (MQ=255)
           ggCACAAGTTCCAAAGCACGTTATTGATAAAGCAATTGATAAAGCCAAAGGCGGTGGAGATGAAACGtt                                             >  1:3352620/1‑69 (MQ=255)
           ggCACAAGTTCCAAAGCACGTTATTGATAAAGCAATTGATAAAGCCAAAGGCGGTGGAGATGAAACGtt                                             >  1:3290707/1‑69 (MQ=255)
           ggCACAAGTTCCAAAGCACGTTATTGATAAAGCAATTGATAAAGCCAAAGGCGGTGGAGATGAAACGtt                                             >  1:1576377/1‑69 (MQ=255)
           ggCACAAGTTCCAAAGCACGTTATTGATAAAGCAATTGATAAAGCCAAAGGCGGCGGAGATGAAACGtt                                             <  1:340144/69‑1 (MQ=255)
            gCACAAGTTCCAAAGCACGTTATTGATAAAGCAATTGATAAAGCCAAAGGCGGCGGAGATGAAACGtt                                             <  1:373278/68‑1 (MQ=255)
            gCACAAGTTCCAAAGCACGTTATTGATAAAGCAATTGATAAAGCCAAAGGCGGCGGAGATGAAACGtt                                             <  1:1058159/68‑1 (MQ=255)
            gCACAAGTTCCAAAGCACGTTATTGATAAAGCAATTGATAAAGCCAAAGGCGGCGGAGATGAAACGtt                                             <  1:260931/68‑1 (MQ=255)
            gCACAAGTTCCAAAGCACGTTATTGATAAAGCAATTGATAAAGCCAAAGGCGGCGGAGATGAAACGtt                                             <  1:3062925/68‑1 (MQ=255)
                        aaacACGTTATTGATAAAGCAATTGATAAAGCCAAAGGCGGCGGAGATGAAACGTTCGTGCAGGGACGt                                <  1:1840990/66‑1 (MQ=255)
                        aaGCACGTTATTGATAAAGCAATTGATAAAGCCAAAGGCGGTGGAGATGAAACGTTCGTACAggg                                    <  1:1871545/65‑1 (MQ=255)
                        aaGCACGTTATTGATAAAGCAATTGATAAAGCCAAAGGCGGTGGAGATGAAACGTTCGTACAggg                                    <  1:1859383/65‑1 (MQ=255)
                        aaGCACGTTATTGATAAAGCAATTGATAAAGCCAAAGGCGGTGGAGATGAAACGTTCGTACAggg                                    <  1:2031395/65‑1 (MQ=255)
                        aaGCACGTTATTGATAAAGCAATTGATAAAGCCAAAGGCGGTGGAGATGAAACGTTCGTACAggg                                    <  1:3352842/65‑1 (MQ=255)
                        aaGCACGTTATTGATAAAGCAATTGATAAAGCCAAAGGCGGTGGAGATGAAACGTTCGTACAggg                                    <  1:1045625/65‑1 (MQ=255)
                        aaGCACGTTATTGATAAAGCAATTGATAAAGCCAAAGGCGGCGGAGATGAAACGTTCGTGCAGGGACGt                                <  1:1374823/69‑1 (MQ=255)
                        aaGCACGTTATTGATAAAGCAATTGATAAAGCCAAAGGCGGCGGAGATGAAACGTTCGTGCAGGGACGt                                <  1:2593215/69‑1 (MQ=255)
                        aaGCACGTTATTGATAAAGCAATTGATAAAGCCAAAGGCGGCGGAGATGAAACGTTCGTGCAGGGACGt                                <  1:3279854/69‑1 (MQ=255)
                        aaGCACGTTATTGATAAAGCAATTGATAAAGCCAAAGGCGGCGGAGATGAAACGTTCGTGCAGGGACGt                                <  1:243755/69‑1 (MQ=255)
                        aaGCACGTTATTGATAAAGCAATTGATAAAGCCAAAGGCGGCGGAGATGAAACGTTCGTGCAGGGACGt                                <  1:3226902/69‑1 (MQ=255)
                        aaGCACGTTATTGATAAAGCAATTGATAAAGCCAAAGGCGGCGGAGATGAAACGTTCGTGCAGGGACGt                                <  1:3181142/69‑1 (MQ=255)
                        aaGCACGTTATTGATAAAGCAATTGATAAAGCCAAAGGCGGCGGAGATGAAACGTTCGTGCAGGGACGt                                <  1:1721782/69‑1 (MQ=255)
                        aaGCACGTTATTGATAAAGCAATTGATAAAGCCAAAGGCGGCGGAGATGAAACGTTCGTGCAGGGACGt                                <  1:186533/69‑1 (MQ=255)
                         aGCACGTTATTGATAAAGCAATTGATAAAGCCAAAGGCGGCGGAGATGAAACGTTCGTGCAGGGACGt                                <  1:1195186/68‑1 (MQ=255)
                         aGCACGTTATTGATAAAGAAATTGATAAAGCCAAAGGCGGTGGAGATGAAACGTTCGTACAggg                                    <  1:157588/64‑1 (MQ=255)
                          gCACGTTATTGATAAAGCAATTGATAAAGCCAAAGGCGGCGGAGATGAAACGTTCGTGCAGGGACGtt                               >  1:2645067/1‑68 (MQ=255)
                                 aTTGATAAAGCAATTGATAAAGCCAAAGGCGGCGGAGATGAAACGTTCGt                                          >  1:3137568/1‑50 (MQ=255)
                                 aTTGATAAAGCAATTGATAAAGCCAAAGGCGGCGGAGATGAAACGTTCGTGCAGGGACGTTATGAAGGc                       >  1:1626617/1‑69 (MQ=255)
                                 aTTGATAAAGCAATTGATAAAGCCAAAGGCGGCGGAGATGAAACGTTCGTGCAGGGACGTTATGAAGGc                       >  1:1088587/1‑69 (MQ=255)
                                 aTTGATAAAGCAATTGATAAAGCCAAAGGCGGCGGAGATGAAACGTTCGTGCAGGGACGTTATGAAGGc                       >  1:442981/1‑69 (MQ=255)
                                 aTTGATAAAGCAATTGATAAAGCCAAAGGCGGCGGAGATGAAACGTTCGTGCAGGGACGTTATGAAGGc                       >  1:3095640/1‑69 (MQ=255)
                                 aTTGATAAAGCAATTGATAAAGCCAAAGGCGGCGGAGATGAAACGTTCGTGCAGGGACGTTATGAAGGc                       >  1:3051646/1‑69 (MQ=255)
                                 aTTGATAAAGCAATTGATAAAGCCAAAGGCGGCGGAGATGAAACGTTCGTGCAGGGACGTTATGAAGGc                       >  1:2513255/1‑69 (MQ=255)
                                      tAAAGCAATTGATAAAGCCAAAGGCGGCGGAGATGAAACGTTCGTGCAGGGACGTTATGAAGGCTTTg                   <  1:878058/68‑1 (MQ=255)
                                      tAAAGCAATTGATAAAGCCAAAGGCGGCGGAGATGAAACGTTCGTGCAGGGACGTTATGAAGGCTTTg                   <  1:1591634/68‑1 (MQ=255)
                                       aaaGCAATTGATAAAGCCAAAGGCGGCGGAGATGAAACGTTCGTGCAGGGACGTTATGAAGGCTTTGGt                 <  1:3186945/69‑1 (MQ=255)
                                        aaGCAATTGATAAAGCCAAAGGCGGCGGAGATGAAACGTTCGTGCAGGGACGTTATGAAGGCTTTGGt                 <  1:1388768/68‑1 (MQ=255)
                                        aaGCAATTGATAAAGCCAAAGGCGGCGGAGATGAAACGTTCGTGCAGGGACGTTATGAAGGCTTTGGt                 <  1:2009953/68‑1 (MQ=255)
                                          ggaaTTGATAAAGCCAAAGGCGGCGGAGATGAAACGTTCGTGCAGGGACGTTATGAAGGCTTTGGTCCt              <  1:2586542/67‑1 (MQ=255)
                                          gCAATTGATAAAGCCAAAGGCGGCGGAGATGAAACGTTCGTGCAGGGACGTTATGAAGGCTTTGGTCCt              <  1:1605034/69‑1 (MQ=255)
                                           cAATTGATAAAGCCAAAGGCGGCGGAGATGAAACGTTCGTGCAGGGACGTTATGAAGGCTTTGGTCCTa             <  1:155493/69‑1 (MQ=255)
                                           cAATTGATAAAGCCAAAGGCGGCGGAGATGAAACGTTCGTGCAGGGACGTTATGAAGGCTTTGGTCCTa             <  1:908904/69‑1 (MQ=255)
                                           cAATTGATAAAGCCAAAGGCGGCGGAGATGAAACGTTCGTGCAGGGACGTTATGAAGGCTTTGGTCCTa             <  1:3267842/69‑1 (MQ=255)
                                           cAATTGATAAAGCCAAAGGCGGCGGAGATGAAACGTTCGTGCAGGGACGTTATGAAGGCTTTGGTCCTa             <  1:141774/69‑1 (MQ=255)
                                           cAATTGATAAAGCCAAAGGCGGCGGAGATGAAACGTTCGTGCAGGGACGTTATGAAGGCTTTGGTCCTa             <  1:1113084/69‑1 (MQ=255)
                                            aaTTGATAAGCCAAAGGCGGTGGAGATGAAACGTTCGTACAGGGACGTTATGAAGGCTTTGGTCCTAAt           <  1:2420134/69‑1 (MQ=255)
                                             aTTGATAAAGCCAAAGGCGGTGGAGATGAAACGTTCGTACAGGGACGTTATGAAGGCTTTGGTCCTAAt           <  1:3227947/69‑1 (MQ=255)
                                             aTTGATAAAGCCAAAGGCGGTGGAGATGAAACGTTCGTACAGGGACGTTATGAAGGCTTTGGTCCTAAt           <  1:1407275/69‑1 (MQ=255)
                                             aTTGATAAAGCCAAAGGCGGTGGAGATGAAACGTTCGTACAGGGACGTTATGAAGGCTTTGGTCCTAAt           <  1:1011351/69‑1 (MQ=255)
                                             aTTGATAAAGCCAAAGGCGGTGGAGATGAAACGTTCGTACAGGGACGTTATGAAGGCTTTGGTCCTAAt           <  1:3206027/69‑1 (MQ=255)
                                             aTTGATAAAGCCAAAGGCGGTGGAGATGAAACGTTCGTACAGGGACGTTATGAAGGCTTTGGTCCTAAt           <  1:909737/69‑1 (MQ=255)
                                             aTTGATAAAGCCAAAGGCGGTGGAGATGAAACGTTCGTACAGGGACGTTATGAAGGCTTTGGTCCTAAt           <  1:2612501/69‑1 (MQ=255)
                                             aTTGATAAAGCCAAAGGCGGTGGAGATGAAACGTTCGTACAGGGACGTTATGAAGGCTTTCGTCCTAAt           <  1:3313324/69‑1 (MQ=255)
                                             aTTGATAAAGCCAAAGGCGGCGGAGATGAAACGTTCGTGCAGGGACGTTATGAAGGCTTTGGTCCTaa            >  1:740675/1‑68 (MQ=255)
                                             aTTGATAAAGCCAAAGGCGGCGGAGATGAAACGTTCGTGCAGGGACGTTATGAAGGCTTTGGTCCTaa            >  1:2326289/1‑68 (MQ=255)
                                             aTTGATAAAGCCAAAGGCGGCGGAGATGAAACGTTCGTGCAGGGACGTTATGAAGGCTTTGGTCCTaa            >  1:81662/1‑68 (MQ=255)
                                             aTTGATAAAGCCAAAGGCGGCGGAGATGAAACGTTCGTGCAGGGACGTTATGAAGGCTTTGGTCCTAAt           >  1:1689315/1‑69 (MQ=255)
                                             aTTGATAAAGCCAAAGGCGGCGGAGATGAAACGTTCGTGCAGGGACGTTATGAAGGCTTTGGTCCTAAt           >  1:2752144/1‑69 (MQ=255)
                                             aTTGATAAAGCCAAAGGCGGCGGAGATGAAACGTTCGTGCAGGGACGTTATGAAGGCTTTGGTCCTAAt           >  1:1296584/1‑69 (MQ=255)
                                             aTTGATAAAGCCAAAGGCGGCGGAGATGAAACGTTCGTGCAGGGACGTTATGAAGGCTTTGGTCCTAAt           >  1:771280/1‑69 (MQ=255)
                                             aTTGATAAAGACAAAGGCGGTGGAGATGAAACGTTCGTACAGGGACGTTATGAAGGCTTTGGTCCTAAt           <  1:2296673/69‑1 (MQ=255)
                                              ttGATAAAGCCAAAGGCGGTGGAGATGAAACGTTCGTACAGGGACGTTATGAAGGCTTTGGTCCTAAt           <  1:2154290/68‑1 (MQ=255)
                                              ttGATAAAGCCAAAGGCGGCGGAGATGAAACGTTCGTGCAGGGACGTTATGAAGGCTTTGGTCCt              <  1:2268628/65‑1 (MQ=255)
                                                gATAAAGCCAAAGGCGGTGGAGATGAAACGTTCGTACAggg                                    <  1:2401105/41‑1 (MQ=255)
                                                gATAAAGCCAAAGGCGGCGGAGATGAAACGTTCGTGCAGGGACGTTATGAAGGCTTTGGTCCt              >  1:2651123/1‑63 (MQ=255)
                                                gATAAAGCCAAAGGCGGCGGAGATGAAACGTTCGTGCAGGGACGTTATGAAGGCTTTGGTCCt              >  1:2526131/1‑63 (MQ=255)
                                                gATAAAGCCAAAGGCGGCGGAGATGAAACGTTCGTGCAGGGACGTTATGAAGGCTTTGGTCCt              >  1:2230240/1‑63 (MQ=255)
                                                gATAAAGCCAAAGGCGGCGGAGATGAAACGTTCGTGCAGGGACGTTATGAAGGCTTTGGTCCt              >  1:971776/1‑63 (MQ=255)
                                                    aaGCCAAAGGCGGCGGAGATGAAACGTTCGTGCAGGGACGTTATGAAGGCTTTGGTCCTAATGGCTc      >  1:2225553/1‑67 (MQ=255)
                                                    aaGCCAAAGGCGGCGGAGATGAAACGTTCGTGCAGGGACGTTATGAAGGCTTTGGTCCTAATGGCTCa     >  1:2221647/1‑68 (MQ=255)
                                                      gCCAAAGGCGGCGGAGATGAAACGTTCGTGCAGGGACGTTATGAAGGCTTTGGTCCTAATGGCTCAATg  >  1:1053791/1‑69 (MQ=255)
                                                       ccAAAGGCGGCGGAGATGAAACGTTCGTGCAGGGACGTTATGAAGGCTTTGGTCCTAATGGCTAAATg  >  1:1971527/1‑68 (MQ=255)
                                                                 |                                                         
CGTGCAAAGCAGGCACAAGTTCCAAAGCACGTTATTGATAAAGCAATTGATAAAGCCAAAGGCGGCGGAGATGAAACGTTCGTGCAGGGACGTTATGAAGGCTTTGGTCCTAATGGCTCAATG  >  minE/1290325‑1290447

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: