Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 123,550 A→G 12.6% intergenic (+177/+25) cueO → / ← gcd multicopper oxidase (laccase)/glucose dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE123,5500AG12.6% 269.4 / 25.8 103intergenic (+177/+25)cueO/gcdmulticopper oxidase (laccase)/glucose dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (48/42);  new base G (9/4);  total (57/46)
Fisher's exact test for biased strand distribution p-value = 3.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.04e-01

Only 100 of 104 total aligned reads displayed.

ATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAG  >  minE/123480‑123614
                                                                      |                                                                
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)
           ttGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGc                                                        >  1:2512564/1‑70 (MQ=255)
                  aTTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAg                                                <  1:2960430/71‑1 (MQ=255)
                   ttttGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAg                                                <  1:3380507/70‑1 (MQ=255)
                   ttttGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGc                                               <  1:4778690/71‑1 (MQ=255)
                   ttttGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGc                                               <  1:4246953/71‑1 (MQ=255)
                              ggATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCt                                                   <  1:1582692/56‑1 (MQ=255)
                              ggATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCt                                                   <  1:4003144/56‑1 (MQ=255)
                              ggATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCt                                                   <  1:3419434/56‑1 (MQ=255)
                                  aaGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTcaca                                >  1:4399778/1‑71 (MQ=255)
                                  aaGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTcaca                                >  1:5044388/1‑71 (MQ=255)
                                  aaGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTcaca                                >  1:3201378/1‑71 (MQ=255)
                                  aaGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTcaca                                >  1:2159653/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:1671422/63‑1 (MQ=255)
                                        ttCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTca                                  <  1:3379220/63‑1 (MQ=255)
                                        ttCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTca                                  <  1:2995317/63‑1 (MQ=255)
                                        ttCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTca                                  <  1:4823033/63‑1 (MQ=255)
                                        ttCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTca                                  <  1:2615997/63‑1 (MQ=255)
                                           aCGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGc                       <  1:4112905/71‑1 (MQ=255)
                                                   tCCGGCATGAACAACTCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCAcatca                             <  1:492172/57‑1 (MQ=255)
                                                   tCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCAcatca                             <  1:2445230/57‑1 (MQ=255)
                                                   tCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCAcatca                             <  1:2618008/57‑1 (MQ=255)
                                                   tCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCAcatca                             <  1:1245504/57‑1 (MQ=255)
                                                   tCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCAcatca                             <  1:3083657/57‑1 (MQ=255)
                                                   tCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCAcatca                             <  1:124561/57‑1 (MQ=255)
                                                    ccGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGc                       <  1:1362091/62‑1 (MQ=255)
                                                    ccGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGc                       <  1:180119/62‑1 (MQ=255)
                                                    ccGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGc                       <  1:2661195/62‑1 (MQ=255)
                                                    ccGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGc                       <  1:1959027/62‑1 (MQ=255)
                                                     cGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACatcatc                           >  1:2228695/1‑57 (MQ=255)
                                                     cGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACatcatc                           >  1:2059409/1‑57 (MQ=255)
                                                     cGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACatcatc                           >  1:1080901/1‑57 (MQ=255)
                                                     cGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACatcatc                           >  1:4274750/1‑57 (MQ=255)
                                                     cGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGc                       <  1:3257209/61‑1 (MQ=255)
                                                      ggCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACtt                                    <  1:1834637/47‑1 (MQ=255)
                                                      ggCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACtt                                    <  1:2885687/47‑1 (MQ=255)
                                                      ggCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAg                     >  1:3908905/1‑62 (MQ=255)
                                                      ggCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAg                     >  1:2126313/1‑62 (MQ=255)
                                                      ggCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAg                     >  1:4033520/1‑62 (MQ=255)
                                                      ggCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAg                     >  1:3121950/1‑62 (MQ=255)
                                                      ggCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAg                     >  1:3308964/1‑62 (MQ=255)
                                                      ggCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAg                     >  1:4152990/1‑62 (MQ=255)
                                                      ggCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAg                     >  1:3444632/1‑62 (MQ=255)
                                                      ggCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAg                     >  1:3515051/1‑62 (MQ=255)
                                                         ataaacaacGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCca          >  1:2222117/4‑70 (MQ=255)
                                                         aTGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCca          >  1:5037588/1‑70 (MQ=255)
                                                         aTGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCca          >  1:1705526/1‑70 (MQ=255)
                                                          tGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCt                                         <  1:2383851/38‑1 (MQ=255)
                                                          tGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCt                                         <  1:1799818/38‑1 (MQ=255)
                                                          tGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCt                                         <  1:4722535/38‑1 (MQ=255)
                                                          tGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCcaca        <  1:4251523/71‑1 (MQ=255)
                                                          tGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCcaca        <  1:2796833/71‑1 (MQ=255)
                                                          tGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCcac         >  1:4176877/1‑70 (MQ=255)
                                                          tGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCcac         >  1:4145266/1‑70 (MQ=255)
                                                           tgacaacGTGCGCGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAgc                    <  1:5194924/56‑1 (MQ=38)
                                                           tgacaacGTGCGCGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAgc                    <  1:805606/56‑1 (MQ=38)
                                                           tgacaacGTGCGCGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAgc                    <  1:2698849/56‑1 (MQ=38)
                                                           tgacaacGTGCGCGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAgc                    <  1:2108003/56‑1 (MQ=38)
                                                            aacaacGCGCACGTTGTCAGCAATCTGAGCAAGTc                                          >  1:1728994/1‑35 (MQ=255)
                                                            aacaacGCGCACGTTGTCAGCAATCTGAGCAAGTc                                          >  1:4790261/1‑35 (MQ=255)
                                                            aacaacGCGCACGTTGTCAGCAATCTGAGCAAGTc                                          >  1:2636692/1‑35 (MQ=255)
                                                            aacaacGCGCACGTTGTCAGCAATCTGAGCAAGTc                                          >  1:1633322/1‑35 (MQ=255)
                                                            aacaacGCGCACGTTGTCAGCAATCTGAGCAAGTc                                          >  1:2218350/1‑35 (MQ=255)
                                                            aacaacGCGCACGTTGTCAGCAATCTGAGCAAGTc                                          >  1:4666991/1‑35 (MQ=255)
                                                            aacaacGCGCACGTTGTCAGCAATCTGAGCAAGTc                                          >  1:4969919/1‑35 (MQ=255)
                                                            aacaacGCGCACGTTGTCAGCAATCTGAGCAAGTc                                          >  1:5021884/1‑35 (MQ=255)
                                                            aacaacGCGCACGTTGTCAGCAATCTGAGCAAGTc                                          >  1:5031380/1‑35 (MQ=255)
                                                            aacaacGCGCACGTTGTCAGCAATCTGAGCAAGTc                                          >  1:752411/1‑35 (MQ=255)
                                                            aacaacGCGCACGTTGTCAGCAATCTGAGCAAGTc                                          >  1:418076/1‑35 (MQ=255)
                                                            aacaacGCGCACGTTGTCAGCAATCTGAGCAAGTc                                          >  1:3663052/1‑35 (MQ=255)
                                                            aacaacGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACaa       <  1:4267466/70‑1 (MQ=255)
                                                            aacaacGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACaa       <  1:1919328/70‑1 (MQ=255)
                                                             acaacGTGCGCGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAgcgc                  >  1:4748103/1‑58 (MQ=255)
                                                             acaacGTGCGCGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAgcgc                  >  1:1194893/1‑58 (MQ=255)
                                                             acaacGTGCGCGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAgcgc                  >  1:750273/1‑58 (MQ=255)
                                                             acaacGTGCGCGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAgcgc                  >  1:2619464/1‑58 (MQ=255)
                                                             acaacGTGCGCGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAgcgc                  >  1:4718154/1‑58 (MQ=255)
                                                             acaacGTGCGCGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAgcgc                  >  1:1794515/1‑58 (MQ=255)
                                                             acaacGTGCGCGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAgcgc                  >  1:3461024/1‑58 (MQ=255)
                                                             acaacGTGCGCGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAgcgc                  >  1:1514934/1‑58 (MQ=255)
                                                             acaacGTGCGCGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAgcgc                  >  1:4183298/1‑58 (MQ=255)
                                                             acaacGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAat      <  1:605641/70‑1 (MQ=255)
                                                                acGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCa                 >  1:2188052/1‑56 (MQ=255)
                                                                acGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCa                 >  1:1240637/1‑56 (MQ=255)
                                                                acGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGc            <  1:2961292/61‑1 (MQ=255)
                                                                acGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGc            <  1:4381169/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)
                                                                      |                                                                
ATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAG  >  minE/123480‑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: