Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,488,647 A→C 9.5% V127V (GTT→GTG 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,6470AC 9.5% 232.8 / 20.4 95V127V (GTT→GTGnuoMNADH:ubiquinone oxidoreductase, membrane subunit M
Reads supporting (aligned to +/- strand):  ref base A (55/31);  new base C (4/5);  total (59/36)
Fisher's exact test for biased strand distribution p-value = 2.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 104 total aligned reads displayed.

ATCGGCACCAGCATCATTTCCCAGAAGAAGAAGAACAGGAACATGTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTCCACGAA  >  minE/1488579‑1488716
                                                                    |                                                                     
aTCGGCACCAGCATCATTTCCCAGAAGAAGAAGAACAGGAACATGTCGATGGCAAGGAACACGCCGATaa                                                                      >  1:915908/1‑70 (MQ=255)
aTCGGCACCAGCATCATTTCCCAGAAGAAGAAGAACAGGAACATGTCGATGGCAAGGAACACGCCGATAAc                                                                     >  1:4331577/1‑71 (MQ=255)
         aGCATCATTTCCCAGAAGAAGAAGAACAGGAACATGTCGATGGCAAGGAACACGCCGATaa                                                                      >  1:3902023/1‑61 (MQ=255)
         aGCATCATTTCCCAGAAGAAGAAGAACAGGAACATGTCGATGGCAAGGAACACGCCGATaa                                                                      >  1:2596952/1‑61 (MQ=255)
         aGCATCATTTCCCAGAAGAAGAAGAACAGGAACATGTCGATGGCAAGGAACACGCCGATaa                                                                      >  1:2610378/1‑61 (MQ=255)
         aGCATCATTTCCCAGAAGAAGAAGAACAGGAACATGTCGATGGCAAGGAACACGCCGATaa                                                                      >  1:1655390/1‑61 (MQ=255)
         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:2768686/62‑1 (MQ=255)
                        aagaagaagaaCAGGAACAGGTCGATGGCGAGGAACACGCCGATCACACCGCCCAGGATCca                                                      <  1:3713450/62‑1 (MQ=255)
                        aagaagaagaaCAGGAACAGGTCGATGGCGAGGAACACGCCGATCACACCGCCCAGGATCca                                                      <  1:4091278/62‑1 (MQ=255)
                        aagaagaagaaCAGGAACAGGTCGATGGCGAGGAACACGCCGATCACACCGCCCAGGATCca                                                      <  1:2188976/62‑1 (MQ=255)
                          gaagaagaaCAGGAACATGTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGa                                                          >  1:5109751/1‑56 (MQ=255)
                             gaagaaCAGGAACATGTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATTCACATCAgg                                               >  1:2747887/1‑64 (MQ=255)
                             gaagaaCAGGAACATGTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATCAgg                                               >  1:2887115/1‑64 (MQ=255)
                             gaagaaCAGGAACATGTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATCAgg                                               >  1:612069/1‑64 (MQ=255)
                             gaagaaCAGGAACATGTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATCAgg                                               >  1:2136153/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:2050930/1‑50 (MQ=255)
                                        aCATGTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATc                                                  >  1:1544512/1‑50 (MQ=255)
                                        aCATGTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATCa                                                 >  1:4016406/1‑51 (MQ=255)
                                        aCATGTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATCa                                                 >  1:495916/1‑51 (MQ=255)
                                        aCATGTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATCa                                                 >  1:2153989/1‑51 (MQ=255)
                                        aCATGTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATCa                                                 >  1:4399548/1‑51 (MQ=255)
                                        aCATGTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATCa                                                 >  1:1643898/1‑51 (MQ=255)
                                          aTGTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAGGTgg                                      <  1:1334506/60‑1 (MQ=255)
                                            gTCGATGGCAAGGAACACGCCGATAACGCCGCCCAgg                                                           >  1:1296419/1‑37 (MQ=255)
                                            gTCGATGGCAAGGAACACGCCGATAACGCCGCCCAgg                                                           >  1:4944848/1‑37 (MQ=255)
                                            gTCGATGGCAAGGAACACGCCGATAACGCCGCCCAgg                                                           >  1:1745620/1‑37 (MQ=255)
                                            gTCGATGGCAAGGAACACGCCGATAACGCCGCCCAgg                                                           >  1:4543982/1‑37 (MQ=255)
                                            gTCGATGGCAAGGAACACGCCGATAACGCCGCCCAgg                                                           >  1:293323/1‑37 (MQ=255)
                                            gTCGATGGCAAGGAACACGCCGATAACGCCGCCCAgg                                                           >  1:4587768/1‑37 (MQ=255)
                                            gTCGATGGCAAGGAACACGCCGATAACGCCGCCCAgg                                                           >  1:1687650/1‑37 (MQ=255)
                                            gTCGATGGCAAGGAACACGCCGATAACGCCGCCCAgg                                                           >  1:166177/1‑37 (MQ=255)
                                            gTCGATGGCAAGGAACACGCCGATAACGCCGCCCAg                                                            >  1:3094202/1‑36 (MQ=255)
                                            gTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATCAgg                                               <  1:3997151/49‑1 (MQ=255)
                                            gTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATCAgg                                               <  1:1340821/49‑1 (MQ=255)
                                            gTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATCAgg                                               <  1:4022782/49‑1 (MQ=255)
                                            gTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATCAgg                                               <  1:5114574/49‑1 (MQ=255)
                                            gTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATCAgg                                               <  1:1720408/49‑1 (MQ=255)
                                            gTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAg                                          >  1:4506119/1‑54 (MQ=255)
                                            gTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAg                                          >  1:1699321/1‑54 (MQ=255)
                                            gTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAg                                          >  1:1526321/1‑54 (MQ=255)
                                            gTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAg                                          >  1:1523374/1‑54 (MQ=255)
                                             tCGATGGCGAGGAACACGCCGATCACACCGCCCAGGATCCACAGCAGGTTCAGATGga                                     >  1:3359346/1‑58 (MQ=255)
                                             tCGATGGCGAGGAACACGCCGATCACACCGCCCAGGATCCACAGCAGGTTCAGATGga                                     >  1:4972047/1‑58 (MQ=255)
                                             tCGATGGCGAGGAACACGCCGATCACACCGCCCAGGATCCACAGCAGGTTCAGATGga                                     >  1:4902287/1‑58 (MQ=255)
                                             tCGATGGCGAGGAACACGCCGATCACACCGCCCAGGATCCACAGCAGGTTCAGATGga                                     >  1:1095425/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)
                                                        gAACACGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAGGt                                        <  1:2876952/44‑1 (MQ=255)
                                                        gAACACGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAGGt                                        <  1:1343273/44‑1 (MQ=255)
                                                        gAACACGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAGGt                                        <  1:4559179/44‑1 (MQ=255)
                                                         aaCACGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGAtctc            >  1:1430015/1‑71 (MQ=255)
                                                         aaCACGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGAtctc            >  1:4056653/1‑71 (MQ=255)
                                                         aaCACGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGAt               >  1:435225/1‑68 (MQ=255)
                                                         aaCACGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGAt               >  1:4000879/1‑68 (MQ=255)
                                                         aaCACGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGAt               >  1:154742/1‑68 (MQ=255)
                                                         aaCACGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGAt               >  1:891052/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:3915376/71‑1 (MQ=255)
                                                            acGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCttt         <  1:3657517/71‑1 (MQ=255)
                                                            acGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCttt         <  1:4598837/71‑1 (MQ=255)
                                                             cGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTc        >  1:4636644/1‑71 (MQ=255)
                                                             cGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTc        >  1:2929432/1‑71 (MQ=255)
                                                             cGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTc        >  1:2551157/1‑71 (MQ=255)
                                                             cGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTc        >  1:2454803/1‑71 (MQ=255)
                                                             cGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTc        >  1:2228706/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:75664/61‑1 (MQ=255)
                                                                  aTAACGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGAtct             <  1:2669801/61‑1 (MQ=255)
                                                                  aTAACGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGAtct             <  1:2660206/61‑1 (MQ=255)
                                                                  aTAACGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGAtct             <  1:201660/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)
                                                                   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:168476/1‑70 (MQ=255)
                                                                    aaCGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTCCACGaa  >  1:2610114/1‑70 (MQ=255)
                                                                    aaCGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTCCACGaa  >  1:3538139/1‑70 (MQ=255)
                                                                    aaCGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTCCACGaa  >  1:4127212/1‑70 (MQ=255)
                                                                    |                                                                     
ATCGGCACCAGCATCATTTCCCAGAAGAAGAAGAACAGGAACATGTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTCCACGAA  >  minE/1488579‑1488716

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: