Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,488,467 A→G 73.5% I187I (ATT→ATC 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,4670AG73.5% 108.6 / 35.5 94I187I (ATT→ATCnuoMNADH:ubiquinone oxidoreductase, membrane subunit M
Reads supporting (aligned to +/- strand):  ref base A (18/7);  new base G (24/45);  total (42/52)
Fisher's exact test for biased strand distribution p-value = 2.02e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.36e-01

Only 100 of 105 total aligned reads displayed.

AGCTCTTCATAGTTGAAGGTCCAGACGCCGGTCGCATTGTAGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCG  >  minE/1488399‑1488531
                                                                    |                                                                
agcTCTTCATAGTTGAAGGGCCAGACGCCGGTCGCATTGTAGTGAACAAAAACCAGCGCCAGGATGGCa                                                                  <  1:473727/69‑1 (MQ=255)
  ctctTCATAGTTGAAGGTCCAGACGCCGGTCGCATTGTAGTGAACAAAAACCAGCGCCAGGATTGCaa                                                                 >  1:3066831/1‑68 (MQ=255)
  ctctTCATAGTTGAAGGTCCAGACGCCGGTCGCATTGTAGTGAACAAAAACCAGCGCCAGGATGGCaa                                                                 >  1:886075/1‑68 (MQ=255)
  ctctTCATAGTTGAAGGTCCAGACGCCGGTCGCATTGTAGTGAACAAAAACCAGCGCCAGGATGGCaa                                                                 >  1:157863/1‑68 (MQ=255)
  ctctTCATAGTTGAAGGTCCAGACGCCGGTCGCATTGTAGTGAACAAAAACCAGCGCCAGGATGGCaa                                                                 >  1:2386903/1‑68 (MQ=255)
  ctctTCATAGTTGAAGGTCCAGACGCCGGTCGCATTGTAGTGAACAAAAACCAGCGCCAGGATGGCaa                                                                 >  1:105241/1‑68 (MQ=255)
  ctctTCATAGTTGAAGGTACAGACGCCGGTCGCATTGTAGTGAACAAAAACCAGCGCCAGGATGGCa                                                                  >  1:3248532/1‑67 (MQ=255)
       cATAGTTGAAGGTCCAGACGCCGGTCGCATTGTAGTGAACAAAAACCAGCGCCAGGATGGCAATCAAc                                                            >  1:865269/1‑68 (MQ=255)
         tAGTTGAAGGTCCAGACGCCGGTCGCATTGTAGTGAACAAAAACCAGCGCCAGGATGGCAATCAAc                                                            <  1:1620597/66‑1 (MQ=255)
             tGAAGGTCCAGACGCCGGTCGCATTGTAGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCAc                                                       <  1:1621049/67‑1 (MQ=255)
                    ccAGACGCCGGTCGCATTGTAGTGAACAAAAACCAGCGCCAGGATGGCa                                                                  <  1:669929/49‑1 (MQ=255)
                    ccAGACGCCGGTCGCATTGTAGTGAACAAAAACCAGCGCCAGGATGGCa                                                                  <  1:2504428/49‑1 (MQ=255)
                    ccAGACGCCGGTCGCATTGTAGTGAACAAAAACCAGCGCCAGGATGGCa                                                                  <  1:1852558/49‑1 (MQ=255)
                    ccAGACGCCGGTCGCATTGTAGTGAACAAAAACCAGCGCCAGGATGGCa                                                                  <  1:1661692/49‑1 (MQ=255)
                                      gTAGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATga                             >  1:666265/1‑68 (MQ=255)
                                      gTAGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATga                             >  1:1931296/1‑68 (MQ=255)
                                      gTAGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATga                             >  1:1831518/1‑68 (MQ=255)
                                      gTAGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATga                             >  1:1021568/1‑68 (MQ=255)
                                        aGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATgaag                           <  1:2225037/68‑1 (MQ=255)
                                         gTGAACGAAAACCAGCGCCAGGATGGCGATCAACATCACCAGACCACT‑‑CCTGGGTGTAAATGAAGAAc                        <  1:405860/68‑1 (MQ=255)
                                           gAACGAAAACCAGCGCCAGGATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAAc                        <  1:3365258/68‑1 (MQ=39)
                                           gAACGAAAACCAGCGCCAGGATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAAc                        <  1:2858643/68‑1 (MQ=39)
                                           gAACGAAAACCAGCGCCAGGATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAAc                        <  1:2705890/68‑1 (MQ=39)
                                           gAACGAAAACCAGCGCCAGGATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAAc                        <  1:2564492/68‑1 (MQ=39)
                                           gAACGAAAACCAGCGCCAGGATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAAc                        <  1:2981725/68‑1 (MQ=39)
                                           gAACGAAAACCAGCGCCAGGATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAAc                        <  1:386196/68‑1 (MQ=39)
                                           gAACGAAAACCAGCGCCAGGATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAAc                        <  1:592725/68‑1 (MQ=39)
                                           gAACGAAAACCAGCGCCAGGATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAAc                        <  1:1894943/68‑1 (MQ=39)
                                           gAACGAAAACCAGCGCCAGGATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAAc                        <  1:13136/68‑1 (MQ=39)
                                           gAACGAAAACCAGCGCCAGGATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAAc                        <  1:1410425/68‑1 (MQ=39)
                                           gAACGAAAACCAGCGCCAGGATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAAc                        <  1:167634/68‑1 (MQ=39)
                                           gAACGAAAACCAGCGCCAGGATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAAc                        <  1:147032/68‑1 (MQ=39)
                                           gAACGAAAACCAGCGCCAGGATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAAc                        <  1:1409243/68‑1 (MQ=39)
                                           gAACGAAAACCAGCGCCAGGATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAAc                        <  1:1588851/68‑1 (MQ=39)
                                            aaCAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACtt                      >  1:738888/1‑69 (MQ=255)
                                            aaCAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACtt                      >  1:2726750/1‑69 (MQ=255)
                                                 aaaCCAGCGCCAGGATGGCGATCAACATCACCAGACCACTCGc                                           >  1:3412963/1‑43 (MQ=37)
                                                 aaaCCAGCGCCAGGATGGCGATCAACATCACCAGACCACTCGCCTgg                                       >  1:463226/1‑47 (MQ=37)
                                                 aaaCCAGCGCCAGGATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGTTg                 >  1:461440/1‑69 (MQ=255)
                                                 aaaCCAGCGCCAGGATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGTTg                 >  1:1156237/1‑69 (MQ=255)
                                                 aaaCCAGCGCCAGGATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGTTg                 >  1:2379688/1‑69 (MQ=255)
                                                 aaaCCAGCGCCAGGATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGTTg                 >  1:3013666/1‑69 (MQ=255)
                                                 aaaCCAGCGCCAGGATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGTTg                 >  1:869816/1‑69 (MQ=255)
                                                 aaaCCAGCGCCAGGATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGTTg                 >  1:1257403/1‑69 (MQ=255)
                                                 aaaCCAGCGCCAGGATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGTTg                 >  1:1790918/1‑69 (MQ=255)
                                                 aaaCCAGCGCCAGGATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGTTg                 >  1:2649396/1‑69 (MQ=255)
                                                 aaaCCAGCGCCAGGATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGTTg                 >  1:2888628/1‑69 (MQ=255)
                                                 aaaCCAGCGCCAGGATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGTTg                 >  1:2280432/1‑69 (MQ=255)
                                                 aaaCCAGCGCCAGGATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGTTg                 >  1:1746668/1‑69 (MQ=255)
                                                 aaaCCAGCGCCAGGATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGTTg                 >  1:1851817/1‑69 (MQ=255)
                                                 aaaCCAGCGCCAGGATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGTTg                 >  1:1948416/1‑69 (MQ=255)
                                                 aaaCCAGCGCCAGGATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGTTg                 >  1:381033/1‑69 (MQ=255)
                                                 aaaCCAGCGCCAGGATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGTTg                 >  1:1084670/1‑69 (MQ=255)
                                                 aaaCCAGCGCCAGGATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGTTg                 >  1:1118372/1‑69 (MQ=255)
                                                 aaaCCAGCGCCAGGATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGTTg                 >  1:6562/1‑69 (MQ=255)
                                                 aaaCCAGCGCCAGGATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGTTg                 >  1:638602/1‑69 (MQ=255)
                                                 aaaCCAGCGCCAGGATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGTTg                 >  1:2273504/1‑69 (MQ=255)
                                                 aaaCCAGCGCCAGGATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGTTg                 >  1:2304703/1‑69 (MQ=255)
                                                 aaaCCAGCGCCAGGATGGCGATCAACATCACAAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGTTg                 >  1:730549/1‑69 (MQ=255)
                                                    ccAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTgccg              >  1:1000563/1‑69 (MQ=255)
                                                    ccAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTgccg              >  1:1287423/1‑69 (MQ=255)
                                                    ccAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTgccg              >  1:428037/1‑69 (MQ=255)
                                                    ccAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTgccg              >  1:2835546/1‑69 (MQ=255)
                                                    ccAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTgccg              >  1:2839325/1‑69 (MQ=255)
                                                    ccAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTgcc               >  1:936533/1‑68 (MQ=255)
                                                              gATGGCGATCAACATCCCCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGtt                  <  1:3104970/55‑1 (MQ=38)
                                                              gATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGtt                  <  1:1126520/55‑1 (MQ=38)
                                                              gATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGtt                  <  1:561947/55‑1 (MQ=38)
                                                              gATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGtt                  <  1:583025/55‑1 (MQ=38)
                                                              gATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGtt                  <  1:682596/55‑1 (MQ=38)
                                                              gATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGtt                  <  1:754199/55‑1 (MQ=38)
                                                              gATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGtt                  <  1:764556/55‑1 (MQ=38)
                                                              gATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGtt                  <  1:1538960/55‑1 (MQ=38)
                                                              gATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGtt                  <  1:3389345/55‑1 (MQ=38)
                                                              gATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGtt                  <  1:3378263/55‑1 (MQ=38)
                                                              gATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGtt                  <  1:1506893/55‑1 (MQ=38)
                                                              gATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGtt                  <  1:1601037/55‑1 (MQ=38)
                                                              gATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGtt                  <  1:1427115/55‑1 (MQ=38)
                                                              gATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGtt                  <  1:1631176/55‑1 (MQ=38)
                                                              gATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGtt                  <  1:1914240/55‑1 (MQ=38)
                                                              gATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGtt                  <  1:215467/55‑1 (MQ=38)
                                                              gATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGtt                  <  1:2255036/55‑1 (MQ=38)
                                                              gATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGtt                  <  1:2330103/55‑1 (MQ=38)
                                                              gATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGtt                  <  1:1284125/55‑1 (MQ=38)
                                                              gATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGtt                  <  1:2479723/55‑1 (MQ=38)
                                                              gATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGtt                  <  1:2569311/55‑1 (MQ=38)
                                                              gATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGtt                  <  1:2936488/55‑1 (MQ=38)
                                                              gATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGtt                  <  1:3031068/55‑1 (MQ=38)
                                                              gATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGtt                  <  1:3049248/55‑1 (MQ=38)
                                                              gATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGtt                  <  1:1186067/55‑1 (MQ=38)
                                                              gATGGCGATCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGGtt                  <  1:3387425/55‑1 (MQ=38)
                                                              gATGGCGAGCAACATCACCAGACCACTCGCCTGGGTGTAAATGAAGAACTTGAtt                  <  1:2309748/55‑1 (MQ=37)
                                                                tGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATAcgcg  >  1:2250844/1‑69 (MQ=255)
                                                                tGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATAcgcg  >  1:2749390/1‑69 (MQ=255)
                                                                tGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATAcgcg  >  1:1972794/1‑69 (MQ=255)
                                                                tGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGCGATAcgcg  >  1:1612076/1‑69 (MQ=255)
                                                                  gCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGa        <  1:24515/61‑1 (MQ=255)
                                                                  gCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGa        <  1:129519/61‑1 (MQ=255)
                                                                  gCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGa        <  1:1072069/61‑1 (MQ=255)
                                                                   cAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGa        <  1:2264489/60‑1 (MQ=255)
                                                                    |                                                                
AGCTCTTCATAGTTGAAGGTCCAGACGCCGGTCGCATTGTAGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCG  >  minE/1488399‑1488531

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: