Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 123,547 C→T 12.5% intergenic (+174/+28) cueO → / ← gcd multicopper oxidase (laccase)/glucose dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE123,5470CT12.5% 280.7 / 25.1 104intergenic (+174/+28)cueO/gcdmulticopper oxidase (laccase)/glucose dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (49/42);  new base T (9/4);  total (58/46)
Fisher's exact test for biased strand distribution p-value = 3.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

Only 100 of 106 total aligned reads displayed.

CAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAG  >  minE/123478‑123614
                                                                     |                                                                   
cAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAcgcg                                                                    <  1:4716372/71‑1 (MQ=255)
 aaTCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAcgcg                                                                    <  1:528352/70‑1 (MQ=255)
  aTCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCa                                                                  >  1:2988647/1‑71 (MQ=255)
        tGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGt                                                            >  1:1171462/1‑71 (MQ=255)
        tGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGt                                                            >  1:4209872/1‑71 (MQ=255)
        tGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGt                                                            >  1:1287221/1‑71 (MQ=255)
        tGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGt                                                            >  1:535671/1‑71 (MQ=255)
             ttGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGc                                                        >  1:1495503/1‑70 (MQ=255)
             ttGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGc                                                        >  1:279495/1‑70 (MQ=255)
                    aTTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAg                                                <  1:2960430/71‑1 (MQ=255)
                     ttttGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAg                                                <  1:3380507/70‑1 (MQ=255)
                     ttttGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGc                                               <  1:4246953/71‑1 (MQ=255)
                                ggATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCt                                                   <  1:1582692/56‑1 (MQ=255)
                                ggATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCt                                                   <  1:3419434/56‑1 (MQ=255)
                                ggATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCt                                                   <  1:4003144/56‑1 (MQ=255)
                                    aaGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTcaca                                >  1:5044388/1‑71 (MQ=255)
                                    aaGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTcaca                                >  1:2159653/1‑71 (MQ=255)
                                    aaGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTcaca                                >  1:4399778/1‑71 (MQ=255)
                                    aaGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTcaca                                >  1:1871547/1‑71 (MQ=255)
                                    aaGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGAAAGTCTTACTTcaca                                >  1:4471958/1‑71 (MQ=255)
                                      ggCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCAcat                               >  1:3863334/1‑70 (MQ=255)
                                      ggCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCAcat                               >  1:2462121/1‑70 (MQ=255)
                                       gCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCAcatca                             <  1:2066781/71‑1 (MQ=255)
                                          ttCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTca                                  <  1:4266970/63‑1 (MQ=255)
                                          ttCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTca                                  <  1:2615997/63‑1 (MQ=255)
                                          ttCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTca                                  <  1:2995317/63‑1 (MQ=255)
                                          ttCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTca                                  <  1:4823033/63‑1 (MQ=255)
                                          ttCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTca                                  <  1:1671422/63‑1 (MQ=255)
                                          ttCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTca                                  <  1:3379220/63‑1 (MQ=255)
                                             aCGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGc                       <  1:4112905/71‑1 (MQ=255)
                                                     tCCGGCATGAACAACTCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCAcatca                             <  1:492172/57‑1 (MQ=255)
                                                     tCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCAcatca                             <  1:124561/57‑1 (MQ=255)
                                                     tCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCAcatca                             <  1:2618008/57‑1 (MQ=255)
                                                     tCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCAcatca                             <  1:2445230/57‑1 (MQ=255)
                                                     tCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCAcatca                             <  1:3083657/57‑1 (MQ=255)
                                                     tCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCAcatca                             <  1:1245504/57‑1 (MQ=255)
                                                     tCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCAt                >  1:3172786/1‑70 (MQ=255)
                                                      ccGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGc                       <  1:1362091/62‑1 (MQ=255)
                                                      ccGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGc                       <  1:2661195/62‑1 (MQ=255)
                                                      ccGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGc                       <  1:1959027/62‑1 (MQ=255)
                                                       cGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACatcatc                           >  1:1080901/1‑57 (MQ=255)
                                                       cGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACatcatc                           >  1:2228695/1‑57 (MQ=255)
                                                       cGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACatcatc                           >  1:2059409/1‑57 (MQ=255)
                                                       cGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACatcatc                           >  1:4274750/1‑57 (MQ=255)
                                                       cGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGc                       <  1:3257209/61‑1 (MQ=255)
                                                        ggCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACtt                                    <  1:2885687/47‑1 (MQ=255)
                                                        ggCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACtt                                    <  1:1834637/47‑1 (MQ=255)
                                                        ggCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAg                     >  1:3908905/1‑62 (MQ=255)
                                                        ggCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAg                     >  1:4033520/1‑62 (MQ=255)
                                                        ggCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAg                     >  1:2126313/1‑62 (MQ=255)
                                                        ggCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAg                     >  1:3515051/1‑62 (MQ=255)
                                                        ggCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAg                     >  1:3444632/1‑62 (MQ=255)
                                                        ggCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAg                     >  1:3308964/1‑62 (MQ=255)
                                                        ggCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAg                     >  1:3121950/1‑62 (MQ=255)
                                                           ataaacaacGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCca          >  1:2222117/4‑70 (MQ=255)
                                                           aTGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCca          >  1:1705526/1‑70 (MQ=255)
                                                           aTGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCca          >  1:5037588/1‑70 (MQ=255)
                                                            tGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCt                                         <  1:1799818/38‑1 (MQ=255)
                                                            tGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCt                                         <  1:2383851/38‑1 (MQ=255)
                                                            tGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCt                                         <  1:4722535/38‑1 (MQ=255)
                                                            tGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCcaca        <  1:2119942/71‑1 (MQ=255)
                                                            tGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCcaca        <  1:4251523/71‑1 (MQ=255)
                                                            tGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCcaca        <  1:2796833/71‑1 (MQ=255)
                                                            tGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCcac         >  1:4145266/1‑70 (MQ=255)
                                                            tGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCcac         >  1:4176877/1‑70 (MQ=255)
                                                             tgacaacGTGCGCGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAgc                    <  1:805606/56‑1 (MQ=38)
                                                             tgacaacGTGCGCGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAgc                    <  1:2698849/56‑1 (MQ=38)
                                                             tgacaacGTGCGCGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAgc                    <  1:2108003/56‑1 (MQ=38)
                                                             tgacaacGTGCGCGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAgc                    <  1:5194924/56‑1 (MQ=38)
                                                              aacaacGCGCACGTTGTCAGCAATCTGAGCAAGTc                                          >  1:1633322/1‑35 (MQ=255)
                                                              aacaacGCGCACGTTGTCAGCAATCTGAGCAAGTc                                          >  1:4790261/1‑35 (MQ=255)
                                                              aacaacGCGCACGTTGTCAGCAATCTGAGCAAGTc                                          >  1:4969919/1‑35 (MQ=255)
                                                              aacaacGCGCACGTTGTCAGCAATCTGAGCAAGTc                                          >  1:2636692/1‑35 (MQ=255)
                                                              aacaacGCGCACGTTGTCAGCAATCTGAGCAAGTc                                          >  1:2218350/1‑35 (MQ=255)
                                                              aacaacGCGCACGTTGTCAGCAATCTGAGCAAGTc                                          >  1:1728994/1‑35 (MQ=255)
                                                              aacaacGCGCACGTTGTCAGCAATCTGAGCAAGTc                                          >  1:4666991/1‑35 (MQ=255)
                                                              aacaacGCGCACGTTGTCAGCAATCTGAGCAAGTc                                          >  1:5021884/1‑35 (MQ=255)
                                                              aacaacGCGCACGTTGTCAGCAATCTGAGCAAGTc                                          >  1:5031380/1‑35 (MQ=255)
                                                              aacaacGCGCACGTTGTCAGCAATCTGAGCAAGTc                                          >  1:418076/1‑35 (MQ=255)
                                                              aacaacGCGCACGTTGTCAGCAATCTGAGCAAGTc                                          >  1:3663052/1‑35 (MQ=255)
                                                              aacaacGCGCACGTTGTCAGCAATCTGAGCAAGTc                                          >  1:752411/1‑35 (MQ=255)
                                                              aacaacGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACaa       <  1:4267466/70‑1 (MQ=255)
                                                              aacaacGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACaa       <  1:1919328/70‑1 (MQ=255)
                                                               acaacGTGCGCGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAgcgc                  >  1:1514934/1‑58 (MQ=255)
                                                               acaacGTGCGCGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAgcgc                  >  1:1194893/1‑58 (MQ=255)
                                                               acaacGTGCGCGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAgcgc                  >  1:750273/1‑58 (MQ=255)
                                                               acaacGTGCGCGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAgcgc                  >  1:4718154/1‑58 (MQ=255)
                                                               acaacGTGCGCGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAgcgc                  >  1:4748103/1‑58 (MQ=255)
                                                               acaacGTGCGCGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAgcgc                  >  1:1794515/1‑58 (MQ=255)
                                                               acaacGTGCGCGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAgcgc                  >  1:2619464/1‑58 (MQ=255)
                                                               acaacGTGCGCGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAgcgc                  >  1:3461024/1‑58 (MQ=255)
                                                               acaacGTGCGCGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAgcgc                  >  1:4183298/1‑58 (MQ=255)
                                                               acaacGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAat      <  1:605641/70‑1 (MQ=255)
                                                                  acGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCa                 >  1:1240637/1‑56 (MQ=255)
                                                                  acGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCa                 >  1:2188052/1‑56 (MQ=255)
                                                                  acGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGc            <  1:4381169/61‑1 (MQ=255)
                                                                  acGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGc            <  1:2961292/61‑1 (MQ=255)
                                                                  acGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGc            <  1:3486199/61‑1 (MQ=255)
                                                                  acGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGc            <  1:566306/61‑1 (MQ=255)
                                                                  acGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAg  >  1:3744951/1‑71 (MQ=255)
                                                                     |                                                                   
CAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAG  >  minE/123478‑123614

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: