Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,488,650 G→A 9.3% G126G (GGC→GGT 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,6500GA 9.3% 275.9 / 13.7 96G126G (GGC→GGTnuoMNADH:ubiquinone oxidoreductase, membrane subunit M
Reads supporting (aligned to +/- strand):  ref base G (57/30);  new base A (4/5);  total (61/35)
Fisher's exact test for biased strand distribution p-value = 2.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCATCATTTCCCAGAAGAAGAAGAACAGGAACATGTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTCCACGAACA  >  minE/1488588‑1488718
                                                              |                                                                    
aGCATCATTTCCCAGAAGAAGAAGAACAGGAACATGTCGATGGCAAGGAACACGCCGATAACGCCGCCCAg                                                              >  1:4533502/1‑71 (MQ=255)
    tcatTTCCCAGAAGAAGAAGAACAGGAACATGTCGATGGCAAGGAACACGCCCATAACGCCGCCCAGGATc                                                          <  1:1945249/71‑1 (MQ=255)
      atTTCCCAGAAGAAGAAGAACAGGAACATGTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATcc                                                         <  1:1659856/70‑1 (MQ=255)
       tttCCCAGAAGAAGAAGAACAGGAACAGGTCGATGGCGAGGAACACGCCGATCACACCGCCCAGGATCcac                                                       <  1:3030360/71‑1 (MQ=255)
             agaagaagaagaACAGGAACATGTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATCAgg                                                 >  1:3551469/1‑71 (MQ=255)
               aagaagaagaaCAGGAACAGGTCGATGGCGAGGAACACGCCGATCACACCGCCCAGGATCca                                                        <  1:2188976/62‑1 (MQ=255)
               aagaagaagaaCAGGAACAGGTCGATGGCGAGGAACACGCCGATCACACCGCCCAGGATCca                                                        <  1:2768686/62‑1 (MQ=255)
               aagaagaagaaCAGGAACAGGTCGATGGCGAGGAACACGCCGATCACACCGCCCAGGATCca                                                        <  1:4091278/62‑1 (MQ=255)
               aagaagaagaaCAGGAACAGGTCGATGGCGAGGAACACGCCGATCACACCGCCCAGGATCca                                                        <  1:3713450/62‑1 (MQ=255)
                 gaagaagaaCAGGAACATGTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGa                                                            >  1:5109751/1‑56 (MQ=255)
                    gaagaaCAGGAACATGTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATTCACATCAgg                                                 >  1:2747887/1‑64 (MQ=255)
                    gaagaaCAGGAACATGTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATCAgg                                                 >  1:2136153/1‑64 (MQ=255)
                    gaagaaCAGGAACATGTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATCAgg                                                 >  1:612069/1‑64 (MQ=255)
                    gaagaaCAGGAACATGTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATCAgg                                                 >  1:2887115/1‑64 (MQ=255)
                    gaagaaCAGGAACATGTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATCAgg                                                 >  1:2982546/1‑64 (MQ=255)
                    gaagaaCAGGAACATGTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATCAgg                                                 >  1:1752023/1‑64 (MQ=255)
                    gaagaaCAGGAACATGTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAGGt                                          >  1:2288601/1‑71 (MQ=255)
                     aagaaCAGGAACATGTCGATGGCAAGGAACACGCCGATAAcg                                                                      <  1:1366484/42‑1 (MQ=255)
                      agaaCAGGAACATGTCGATGGCAAGGAACACGCCGATAAcg                                                                      <  1:2370373/41‑1 (MQ=255)
                      agaaCAGGAACATGTCGATGGCAAGGAACACGCCGATAAcg                                                                      <  1:2654791/41‑1 (MQ=255)
                               aCATGTCGATGGCAAGGAACACGCCGATAAcgccgc                                                                  >  1:332863/1‑36 (MQ=255)
                               aCATGTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCca                                                        >  1:3236265/1‑46 (MQ=255)
                               aCATGTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATc                                                    >  1:1544512/1‑50 (MQ=255)
                               aCATGTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATc                                                    >  1:2050930/1‑50 (MQ=255)
                               aCATGTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATCa                                                   >  1:4016406/1‑51 (MQ=255)
                               aCATGTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATCa                                                   >  1:1643898/1‑51 (MQ=255)
                               aCATGTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATCa                                                   >  1:4399548/1‑51 (MQ=255)
                               aCATGTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATCa                                                   >  1:2153989/1‑51 (MQ=255)
                               aCATGTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATCa                                                   >  1:495916/1‑51 (MQ=255)
                                 aTGTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAGGTgg                                        <  1:1334506/60‑1 (MQ=255)
                                 aTGTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAGGTgg                                        <  1:952947/60‑1 (MQ=255)
                                   gTCGATGGCAAGGAACACGCCGATAACGCCGCCCAgg                                                             >  1:4543982/1‑37 (MQ=255)
                                   gTCGATGGCAAGGAACACGCCGATAACGCCGCCCAgg                                                             >  1:4587768/1‑37 (MQ=255)
                                   gTCGATGGCAAGGAACACGCCGATAACGCCGCCCAgg                                                             >  1:4944848/1‑37 (MQ=255)
                                   gTCGATGGCAAGGAACACGCCGATAACGCCGCCCAgg                                                             >  1:293323/1‑37 (MQ=255)
                                   gTCGATGGCAAGGAACACGCCGATAACGCCGCCCAgg                                                             >  1:1745620/1‑37 (MQ=255)
                                   gTCGATGGCAAGGAACACGCCGATAACGCCGCCCAgg                                                             >  1:1687650/1‑37 (MQ=255)
                                   gTCGATGGCAAGGAACACGCCGATAACGCCGCCCAgg                                                             >  1:1296419/1‑37 (MQ=255)
                                   gTCGATGGCAAGGAACACGCCGATAACGCCGCCCAgg                                                             >  1:166177/1‑37 (MQ=255)
                                   gTCGATGGCAAGGAACACGCCGATAACGCCGCCCAg                                                              >  1:3094202/1‑36 (MQ=255)
                                   gTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATCAgg                                                 <  1:5114574/49‑1 (MQ=255)
                                   gTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATCAgg                                                 <  1:1720408/49‑1 (MQ=255)
                                   gTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATCAgg                                                 <  1:3997151/49‑1 (MQ=255)
                                   gTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATCAgg                                                 <  1:4022782/49‑1 (MQ=255)
                                   gTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATCAgg                                                 <  1:1340821/49‑1 (MQ=255)
                                   gTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAg                                            >  1:4506119/1‑54 (MQ=255)
                                   gTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAg                                            >  1:1526321/1‑54 (MQ=255)
                                   gTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAg                                            >  1:1523374/1‑54 (MQ=255)
                                   gTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAg                                            >  1:29197/1‑54 (MQ=255)
                                   gTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAg                                            >  1:1699321/1‑54 (MQ=255)
                                    tCGATGGCGAGGAACACGCCGATCACACCGCCCAGGATCCACAGCAGGTTCAGATGga                                       >  1:4902287/1‑58 (MQ=255)
                                    tCGATGGCGAGGAACACGCCGATCACACCGCCCAGGATCCACAGCAGGTTCAGATGga                                       >  1:1095425/1‑58 (MQ=255)
                                    tCGATGGCGAGGAACACGCCGATCACACCGCCCAGGATCCACAGCAGGTTCAGATGga                                       >  1:4972047/1‑58 (MQ=255)
                                    tCGATGGCGAGGAACACGCCGATCACACCGCCCAGGATCCACAGCAGGTTCAGATGga                                       >  1:3359346/1‑58 (MQ=255)
                                        tGGCAAGGAACACGCCGATAACGCCGCCCAGGATcc                                                         <  1:2641569/36‑1 (MQ=255)
                                        tGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATCAgg                                                 <  1:1258148/44‑1 (MQ=255)
                                            aaGGAACACGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGa                  >  1:2522150/1‑71 (MQ=255)
                                            aaGGAACACGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGa                  >  1:1898301/1‑71 (MQ=255)
                                               gAACACGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAGGt                                          <  1:4559179/44‑1 (MQ=255)
                                               gAACACGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAGGt                                          <  1:1343273/44‑1 (MQ=255)
                                               gAACACGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAGGt                                          <  1:2876952/44‑1 (MQ=255)
                                                aaCACGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGAtctc              >  1:1430015/1‑71 (MQ=255)
                                                aaCACGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGAtctc              >  1:4056653/1‑71 (MQ=255)
                                                aaCACGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGAt                 >  1:4000879/1‑68 (MQ=255)
                                                aaCACGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGAt                 >  1:891052/1‑68 (MQ=255)
                                                aaCACGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGAt                 >  1:435225/1‑68 (MQ=255)
                                                aaCACGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGAt                 >  1:154742/1‑68 (MQ=255)
                                                aaCACGCCGATAACGCCGCCCAGGATCCACATCAGGTTAGGTGGAAGAAGCCCTGATATTTTTCGATctct             >  1:3489106/1‑71 (MQ=255)
                                                 acacGCCGATAACGCCGCCCAGGATCCACATCAgg                                                 <  1:4350580/35‑1 (MQ=255)
                                                 acacGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAGGGGGAAGAAGCCCTGATATTTTTCg                   <  1:3751305/65‑1 (MQ=255)
                                                   acGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCttt           <  1:3657517/71‑1 (MQ=255)
                                                   acGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCttt           <  1:4598837/71‑1 (MQ=255)
                                                   acGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCttt           <  1:3915376/71‑1 (MQ=255)
                                                    cGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTc          >  1:2551157/1‑71 (MQ=255)
                                                    cGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTc          >  1:2228706/1‑71 (MQ=255)
                                                    cGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTc          >  1:2929432/1‑71 (MQ=255)
                                                    cGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTc          >  1:4636644/1‑71 (MQ=255)
                                                    cGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTc          >  1:2454803/1‑71 (MQ=255)
                                                    cGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCACttac          >  1:181104/1‑69 (MQ=255)
                                                      ccGATAACGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATAtttt                      >  1:1223194/1‑57 (MQ=255)
                                                       cGATAACGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTCCAc       <  1:3501276/71‑1 (MQ=255)
                                                         aTAACGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGAtct               <  1:201660/61‑1 (MQ=255)
                                                         aTAACGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGAtct               <  1:75664/61‑1 (MQ=255)
                                                         aTAACGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGAtct               <  1:2669801/61‑1 (MQ=255)
                                                         aTAACGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGAtct               <  1:2660206/61‑1 (MQ=255)
                                                         aTAACGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTCCACg      >  1:3128120/1‑70 (MQ=255)
                                                         aTAACGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTCCACg      >  1:700406/1‑70 (MQ=255)
                                                         aTAACGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTCCACGa     <  1:2762213/71‑1 (MQ=255)
                                                         aTAACGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTCCACGa     <  1:5189226/71‑1 (MQ=255)
                                                          tAACGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTCTACGaa    <  1:4261188/71‑1 (MQ=255)
                                                          tAACGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTCCACGaa    >  1:2494152/1‑71 (MQ=255)
                                                          tAACGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTCCACGaa    <  1:4813241/71‑1 (MQ=255)
                                                          tAACGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTCCACGa     >  1:1235916/1‑70 (MQ=255)
                                                          tAACGCCGACCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTCCACGaa    >  1:4397262/1‑71 (MQ=255)
                                                           aaCGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTCCACGaa    >  1:2610114/1‑70 (MQ=255)
                                                           aaCGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTCCACGaa    >  1:4127212/1‑70 (MQ=255)
                                                           aaCGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTCCACGaa    >  1:168476/1‑70 (MQ=255)
                                                           aaCGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTCCACGaa    >  1:3538139/1‑70 (MQ=255)
                                                            aCGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTCCAc       <  1:3478351/66‑1 (MQ=255)
                                                            aCGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTCCACGAACa  <  1:4347952/71‑1 (MQ=255)
                                                              |                                                                    
AGCATCATTTCCCAGAAGAAGAAGAACAGGAACATGTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTCCACGAACA  >  minE/1488588‑1488718

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: