Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,488,719 T→G 10.0% L103L (CTA→CTC nuoM ← NADH:ubiquinone oxidoreductase, membrane subunit M

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,488,7190TG10.0% 248.2 / 21.6 99L103L (CTA→CTCnuoMNADH:ubiquinone oxidoreductase, membrane subunit M
Reads supporting (aligned to +/- strand):  ref base T (70/19);  new base G (4/6);  total (74/25)
Fisher's exact test for biased strand distribution p-value = 1.53e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 103 total aligned reads displayed.

AGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTCCACGAACATAGTACCGCCAGCACACCGAGCAGACCGGTCAGCACGACCATCAGCAGCGACAGCCCGTCAATGGCGAGAT  >  minE/1488657‑1488789
                                                              |                                                                      
agcaTCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTCCACGAACATACTACCGc                                                                <  1:2455536/68‑1 (MQ=255)
aGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTCCACGAACATAGTAcc                                                                  <  1:1105239/69‑1 (MQ=255)
aGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTCCACGAACATAGTAcc                                                                  <  1:324182/69‑1 (MQ=255)
aGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTCCACGAACATAGTAcc                                                                  <  1:3581336/69‑1 (MQ=255)
aGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTCCACGAACATAGTAcc                                                                  <  1:696209/69‑1 (MQ=255)
  gATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTCCACGAACATAGt                                                                     >  1:2451859/1‑64 (MQ=255)
  gATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTCCACGAACATAGTACCGc                                                                >  1:2416326/1‑69 (MQ=255)
  gATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTCCACGAACATAGTACCGc                                                                >  1:4696703/1‑69 (MQ=255)
  gATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTCCACGAACATAGTACCGc                                                                >  1:4737967/1‑69 (MQ=255)
  gATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTCCACGAACATAGTACCGc                                                                >  1:1545618/1‑69 (MQ=255)
  gATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTCCACGAACATAGTACCGc                                                                >  1:3848426/1‑69 (MQ=255)
    tCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTCCACGAACATAGTACCGCCAg                                                             >  1:1652934/1‑70 (MQ=255)
    tCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTCCACGAACATAGTACCGCCAg                                                             <  1:1194506/70‑1 (MQ=255)
    tCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTCCACGAACATAGTACCGCCAg                                                             >  1:4560682/1‑70 (MQ=255)
    tCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTCCACGAACATAGTACCGCCAg                                                             >  1:4768977/1‑70 (MQ=255)
    tCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTCCACGAACATAGTACCGCCAg                                                             >  1:4719017/1‑70 (MQ=255)
    tCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTCAACGAACATAGTACCGCCa                                                              >  1:87844/1‑69 (MQ=255)
       acaTCAGGTTGAGGTGGAAGAAGGCCTGATATTTTTCGATCTCTTTCCACGAACAGAGTACCGccgg                                                             >  1:119239/1‑65 (MQ=255)
       acaTCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTCCACGAACAGAGTACCGCCAg                                                             >  1:504440/1‑67 (MQ=255)
       acaTCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTCCACGAACAGAGTACCGCCAg                                                             >  1:4808637/1‑67 (MQ=255)
       acaTCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTCCACGAACAGAGTACCGCCAg                                                             >  1:4725458/1‑67 (MQ=255)
        caTCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTCCACGAACATAGTACCGCCAGCacac                                                        >  1:34939/1‑71 (MQ=255)
        caTCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTCCACGAACATAGTACCGCCAGCacac                                                        >  1:869044/1‑71 (MQ=255)
            aGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTCCACGAACATAGTACCGCCAGCACACCGAg                                                    >  1:4325456/1‑71 (MQ=255)
             ggTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTCCACGAACATAGTACCGCCACCACACCGAGc                                                   <  1:464702/71‑1 (MQ=255)
                tGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTCCACGAACATAGTACCGCCAGCACACCGAGCa                                                  >  1:2169873/1‑69 (MQ=255)
                tGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTCCACGAACATAGTACCGCCAGCACACCGAGCa                                                  >  1:1054262/1‑69 (MQ=255)
                tGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTCCACGAACATAGTACCGCCAGCACACCGAGCa                                                  >  1:3443263/1‑69 (MQ=255)
                tGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTCCACGAACATAGTACCGCCAGCACACCGAGCa                                                  >  1:5073500/1‑69 (MQ=255)
                 gAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTCCACGAACATa                                                                       >  1:135274/1‑47 (MQ=255)
                 gAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTCCACGAACATa                                                                       >  1:2224119/1‑47 (MQ=255)
                 gAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTCCACGAACATAGTACCGCCAGCACACCGAGCAGAc                                               >  1:2821955/1‑71 (MQ=255)
                       gaagaaGCCCTGATATTTTTCGATCTCTTTCCACGAACATAGTACCGCCAGCACACCGAGCAGACCGGt                                           >  1:4071143/1‑69 (MQ=255)
                         agaagCCCTGATATTTTTCGATCTCTTTCCACGAACATAGTACCGCCAGCACACGGAGCAGACCGGTCAg                                        >  1:3792873/1‑70 (MQ=255)
                         agaagCCCTGATATTTTTCGATCTCTTTCCACGAACATAGTACCGCCAGCACACCGAg                                                    >  1:1417368/1‑58 (MQ=255)
                         agaagCCCTGATATTTTTCGATCTCTTTCCACGAACATAGTACCGCCAGCACACCGAg                                                    >  1:3389323/1‑58 (MQ=255)
                         agaagCCCTGATATTTTTCGATCTCTTTCCACGAACATAGTACCGCCAGCACACCGAg                                                    >  1:2180938/1‑58 (MQ=255)
                         agaagCCCTGATATTTTTCGATCTCTTTCCACGAACATAGTACCGCCAGCACACCGAg                                                    >  1:702188/1‑58 (MQ=255)
                             gCCCTGATATTTTTCGATCTCTTTCCACGAACATAGTACCGCCAGCACACCGAGCAGACCGGTCAGCACGa                                   <  1:1106812/71‑1 (MQ=255)
                               ccTGATATTTTTCGATCTCTTTCCACGAACATAGTACCGCCAGCACACCGAGCAGACCGGTCAGCACGAc                                  >  1:4641240/1‑70 (MQ=255)
                                 tGATATTTTTCGATCTCTTTCCACGAACATAGTACCGCCAGCACACCGAGCAGACCGGTCAg                                        <  1:1832852/62‑1 (MQ=255)
                                  gATATTTTTCTATCTCTTTCCACGAACAGAGTACCGCCAGCACACCGAGCAGACCGGTCAGCACGAc                                  <  1:1611153/67‑1 (MQ=255)
                                  gATATTTTTCTATCTCTTTCCACGAACAGAGTACCGCCAGCACACCGAGCAGACCGGTCAGCACGAc                                  <  1:2968461/67‑1 (MQ=255)
                                  gATATTTTTCGATCTCTTTCCACGAACAGAGTACCGCCAGCACACCGAGCAGACCGGTCAGCACGAc                                  <  1:2289720/67‑1 (MQ=255)
                                  gATATTTTTCGATCTCTTTCCACGAACAGAGTACCGCCAGCACACCGAGCAGACCGGTCAGCACGAc                                  <  1:4478611/67‑1 (MQ=255)
                                  gATATTTTTCGATCTCTTTCCACGAACAGAGTACCGCCAGCACACCGAGCAGACCGGTCAGCACGAc                                  <  1:1680038/67‑1 (MQ=255)
                                  gATATTTTTCGATCTCTTTCCACGAACAGAGTACCGCCAGCACACCGAGCAGACCGGTCAGCACGAc                                  <  1:4588041/67‑1 (MQ=255)
                                    tatTTTTCGATCTCTTTCCACGAACATAGTACCGCCAGCACACCGAGCAGACCGGTCAGCACGACCATcag                            >  1:3432250/1‑71 (MQ=255)
                                     atTTTTCGATCTCTTTCCACGAACATAGTACCGCCAGCACACCGAGCAGACCGGTCAGCACGACCATcag                            >  1:2681099/1‑70 (MQ=255)
                                      tttttCGATCTCTTTCCACGAACATAGTACCGCCAGCACACCGAGCAGACCGGTCAGCa                                      >  1:4672537/1‑59 (MQ=255)
                                      tttttCGATCTCTTTCCACGAACATAGTACCGCCAGCACACCGAGCAGACCGGTCAGCACGACCATcag                            >  1:13005/1‑69 (MQ=255)
                                      tttttCGATCTCTTTCCACGAACATAGTACCGCCAGCACACCGAGCAGACCGGTCAGCACGACCATcag                            >  1:2153947/1‑69 (MQ=255)
                                      tttttCGATCTCTTTCCACGAACATAGTACCGCCAGCACACCGAGCAGACCGGTCAGCACGACCATcag                            >  1:1954598/1‑69 (MQ=255)
                                      tttttCGATCTCTTTCCACGAACATAGTACCGCCAGCACACCGAGCAGACCGGTCAGCACGACCATcag                            >  1:4177999/1‑69 (MQ=255)
                                      tttttCGATCTCTTTCCACGAACATAGTACCGCCAGCACACCGAGCAGACCGGTCAGCACGACCATcag                            >  1:4918578/1‑69 (MQ=255)
                                      tttttCGATCTCTTTCCACGAACATAGTACAGCCAGCACACCGAGCAGACCGGTCa                                         >  1:1558065/1‑56 (MQ=255)
                                        tttCGATCTCTTTCCACGAACATAGTACCGCCAGCACACCGAGCAGACCGGTCAGCACGACCATCagcagc                        >  1:815723/1‑71 (MQ=255)
                                         ttCGATCTCTTTCCACGAACATAGTACCGCCAGCACACCGAGCAGACCGGTCAGCACGACCATCagcagc                        >  1:4401873/1‑70 (MQ=255)
                                         ttCGATCTCTTTCCACGAACATAGTACCGCCAGCACACCGAGCAGACCGGTCAGCACGACCATCagcagc                        >  1:2719397/1‑70 (MQ=255)
                                          tCGATCTCTTTCGACGAACATAGTACCGCCAGCACACCGAGCAGACCGGTCAGCACGACCATCAGCAGCg                       <  1:4474548/70‑1 (MQ=255)
                                          tCGATCTCTTTCCACGAACATAGTACCGCCAGCACACCGAGCAGa                                                >  1:4982629/1‑45 (MQ=255)
                                          tCGATCTCTTTCCACGAACATAGTACCGCCAGCACACCGAGCAGa                                                >  1:592949/1‑45 (MQ=255)
                                          tCGATCTCTTTCCACGAACATAGTACCGCCAGCACACCGAGCAGa                                                >  1:842601/1‑45 (MQ=255)
                                          tCGATCTCTTTCCACGAACATAGTACCGCCAGCACACCGAGCAGa                                                >  1:1120763/1‑45 (MQ=255)
                                          tCGATCTCTTTCCACGAACATAGTACCGCCAGCACACCGAGCAGa                                                >  1:879114/1‑45 (MQ=255)
                                          tCGATCTCTTTCCACGAACATAGTACCGCCAGCACACCGAGCAGACCGGTCAGCACGACCATCAGCAGCg                       <  1:2731391/70‑1 (MQ=255)
                                          tCGATCTCTTTCCACGAACATAGTACCGCCAGCACACCGAGCAGACCGGTCAGCACGACCATCAGCAGCg                       <  1:2920931/70‑1 (MQ=255)
                                          tCGATCTCTTTCCACGAACATAGTACCGCCAGCACACCGAGCAGACCGGTCAGCACGACCATCAGCAGCg                       <  1:37006/70‑1 (MQ=255)
                                          tCGATCTCTTTCCACGAACATAGTACCGCCAGCACACCGAGCAGACCGGTCAGCACGACCATCAGCAGCg                       <  1:1356107/70‑1 (MQ=255)
                                          tCGATCTATTTCCACGAACATAGTACCGCCAGCACACCGAGCAGa                                                >  1:3911443/1‑45 (MQ=255)
                                                tctTTCCACGAACATAGTACCGCCAGCACACCGAGCAGACCGGTCAGCACGACCATCAGCAGCGACAGccc                >  1:1323555/1‑71 (MQ=255)
                                                tctTTCCACGAACATAGTACCGCCAGCACACCGAGCAGACCGGTCAGCACGACCATCAGCAGCGACAGccc                >  1:4589345/1‑71 (MQ=255)
                                                 ctTTCCACGAACATAGTACCGCCAGCACACCGAGCAGACCGGTCAGCACGACCATCAGCAGCGACAGCCCg               >  1:315333/1‑71 (MQ=255)
                                                   ttCCACGAACATAGTACTGCCAGCACACCGAGCAGACCGGTCAGCACGACCATcagcag                         >  1:4729546/1‑59 (MQ=255)
                                                   ttCCACGAACATAGTACCGCCAGCACACCGAGCAGACCGGTCAGCa                                      >  1:969395/1‑46 (MQ=255)
                                                   ttCCACGAACATAGTACCGCCAGCACACCGAGCAGACCGGTCAGCACGACCATcagcag                         >  1:2134388/1‑59 (MQ=255)
                                                   ttCCACGAACATAGTACCGCCAGCACACCGAGCAGACCGGTCAGCACGACCATcagcag                         >  1:2058543/1‑59 (MQ=255)
                                                   ttCCACGAACATAGTACCGCCAGCACACCGAGCAGACCGGTCAGCACGACCATcagcag                         >  1:2070248/1‑59 (MQ=255)
                                                   ttCCACGAACATAGTACCGCCAGCACACCGAGCAGACCGGTCAGCACGACCATcagcag                         >  1:336271/1‑59 (MQ=255)
                                                   ttCCACGAACATAGTACCGCCAGCACACCGAGCAGACCGGTCAGCACGACCATcagcag                         >  1:3169019/1‑59 (MQ=255)
                                                   ttCCACGAACATAGTACCGCCAGCACACCGAGCAGACCGGTCAGCACGACCATcagcag                         >  1:3483060/1‑59 (MQ=255)
                                                   ttCCACGAACATAGTACCGCCAGCACACCGAGCAGACCGGTCAGCACGACCATcagcag                         >  1:2351797/1‑59 (MQ=255)
                                                   ttCCACGAACATAGTACCGCCAGCACACCGAGCAGACCGGTCAGCACGACCATcagcag                         >  1:3693314/1‑59 (MQ=255)
                                                   ttCCACGAACATAGTACCGCCAGCACACCGAGCAGACCGGTCAGCACGACCATcagcag                         >  1:4130358/1‑59 (MQ=255)
                                                   ttCCACGAACATAGTACCGCCAGCACACCGAGCAGACCGGTCAGCACGACCATcagcag                         >  1:2065212/1‑59 (MQ=255)
                                                    tCCACGAACATAGTACCGCCAGCACACCGAGCAGACCGGTCAGCACGACCATcagca                          >  1:3917408/1‑57 (MQ=255)
                                                    tCCACGAACATAGTACCGCCAGCACACCGAGCAGACCGGTCAGCACGACCATcagca                          >  1:843104/1‑57 (MQ=255)
                                                    tCCACGAACATAGTACCGCCAGCACACCGAGCAGACCGGTCAGCACGACCATcagca                          >  1:5063782/1‑57 (MQ=255)
                                                    tCCACGAACATAGTACCGCCAGCACACCGAGCAGACCGGTCAGCACGACCATcagca                          >  1:2442525/1‑57 (MQ=255)
                                                    tCCACGAACATAGTACCGCCAGCACACCGAGCAGACCGGTCAGCACGACCATCAGCAGCGACAGCCCGTc             <  1:1981334/70‑1 (MQ=255)
                                                    tCCACGAAAATAGTACCGCCAGCACACCGAGCAGACCGGTCAGCACGACCATcagca                          >  1:3134759/1‑57 (MQ=255)
                                                      cACGAACATAGTACCGCCAGCACACCGAGCAGACCGGTCAGCACGACCATCAGCAGCGACAGCCCGTCAAt          >  1:3365362/1‑71 (MQ=255)
                                                        cGAACATAGTACCGCCAGCACACCGAGCAGACCGGTCAGCACGACCATCAGCAGCGACAGCCCGTCAATgg        <  1:2516882/71‑1 (MQ=255)
                                                        cGAACATAGTACCGCCAGCACACCGAGCAGACCGGTCAGCACGACCATCAGCAGCGACAGCCCGTCAATgg        <  1:3524152/71‑1 (MQ=255)
                                                        cGAACATAGTACCGCCAGCACACCGAGCAGACCGGTCAGCACGACCATCAGCAGCGACAGCCCGTCAATgg        <  1:909864/71‑1 (MQ=255)
                                                        cGAACATAGTACCGCCAGCACACCGAGCAGACCGGTCAGCACGACCATCAGCAGCGACAGCCCGTCAATgg        <  1:3364593/71‑1 (MQ=255)
                                                              tAGTACCGCCAGCACACCGAGCAGACCGGTCAGCACGACCATCAGCAGCGACAGCCCGTCAATGGCGAGAt  >  1:143480/1‑71 (MQ=255)
                                                              tAGTACCGCCAGCACACCGAGCAGACCGGTCAGCACGACCATCAGCAGCGACAGCCCGTCAATGGCGAGAt  >  1:714730/1‑71 (MQ=255)
                                                              tAGTACCGCCAGCACACCGAGCAGACCGGTCAGCACGACCATCAGCAGCGACAGCCCGTCAATGGCGAGAt  >  1:3516667/1‑71 (MQ=255)
                                                              tAGTACCGCCAGCACACCGAGCAGACCGGTCAGCACGACCATCAGCAGCGACAGCCCGTCAATGGCGAGAt  >  1:2055226/1‑71 (MQ=255)
                                                              |                                                                      
AGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTCCACGAACATAGTACCGCCAGCACACCGAGCAGACCGGTCAGCACGACCATCAGCAGCGACAGCCCGTCAATGGCGAGAT  >  minE/1488657‑1488789

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