Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,500,488 A→G 10.7% F38L (TTT→CTT)  nuoB ← NADH:ubiquinone oxidoreductase, chain B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,500,4880AG10.7% 249.1 / 20.1 93F38L (TTT→CTT) nuoBNADH:ubiquinone oxidoreductase, chain B
Reads supporting (aligned to +/- strand):  ref base A (34/49);  new base G (7/3);  total (41/52)
Fisher's exact test for biased strand distribution p-value = 1.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.11e-01

Only 100 of 101 total aligned reads displayed.

AGACCGAAGTTATACGGCCAAATTGAGTTTTTACGACCCCAGTTAACCATGTCATTGAGCTTGCCCATAAACACGTTTTTGTTAACTTCTTGCTCCAGAGGGTCGGTTACGATCTCCTGCTTTTGCAGGGGGTAACGG  >  minE/1500418‑1500555
                                                                      |                                                                   
aGACCGAAGTTATACGGCCAAATTGAGTTTTTACGACCCCAGTTAACCATGTCATTGAGCTTGCCCATaaa                                                                     >  1:2091535/1‑71 (MQ=255)
    cGAAGTTATACGGCCAAATTGAGTTTTTACGACCCCAGTTAACCATGTCATTGAGCTTGCCCATAAACACg                                                                 <  1:855260/71‑1 (MQ=255)
         ttATACGGCCAAATTGAGTTTTTACGACCCCAGTTAACCATGTCATTGAGCTTGCCCATAAACACGttttt                                                            <  1:1156309/71‑1 (MQ=255)
           ataCGGCCAAATTGAGTTTTTACGACCCCAGTTAACCATGTCATTGAGCTTGCCCATAAACACGTTTTTGt                                                          >  1:2273176/1‑71 (MQ=255)
           ataCGGCCAAATTGAGTTTTTACGACCCCAGTTAACCATGTCATTGAGCTTGCCCATAAACACGTTTTTGt                                                          >  1:1187106/1‑71 (MQ=255)
                gCCAAATTGAGTTTTTACGACCCCAGTTAACCATGTCATTGAGCTTGCCCATAAACACGttttt                                                            <  1:743429/64‑1 (MQ=255)
                gCCAAATTGAGTTTTTACGACCCCAGTTAACCATGTCATTGAGCTTGCCCATAAACACGTTTTTGTTAAct                                                     <  1:2174985/71‑1 (MQ=255)
                gCCAAATTGAGTTTTTACGACCCCAGTTAACCATGTCATTGAGCTTGCCCATAAACACGTTTTTGTTAAct                                                     <  1:634827/71‑1 (MQ=255)
                gCCAAATTGAGTTTTTACGACCCCAGTTAACCATGTCATTGAGCTTGCCCATAAACACGTTTTTGTTAAct                                                     <  1:625054/71‑1 (MQ=255)
                gCCAAATTGAGTTTTTACGACCCCAGTTAACCATGTCATTGAGCTTGCCCATAAACACGTTTTTGTTAAct                                                     <  1:2786200/71‑1 (MQ=255)
                 ccAAATTGAGTTTTTACGACCCCAGTTAACCATGTCATTGAGCTTGCCCATAAACACGTTTTTGTTAActt                                                    <  1:2241950/71‑1 (MQ=255)
                 ccAAATTGAGTTTTTACGACCCCAGTTAACCATGTCAATGAGCTTGCCCATAAACACGTTTTTGTTAActt                                                    <  1:1943124/71‑1 (MQ=255)
                  cAAATTGAGTTTTTACGACCCCAGTTAACCATGTCATTGAGCTTGCCCATAAACACGTTTTTGTTAActtc                                                   <  1:2051391/71‑1 (MQ=255)
                     aTTGAGTTTTTACGACCCCAGTTAACCATGTCATTGAGCTTGCCCATAAACACGTTTTTGtt                                                         <  1:2204444/62‑1 (MQ=255)
                      ttGAGTTTTTACGACCCCAGTTAACCATGTCATTGAGCTTGCCCATAAACACGTTTTTGtt                                                         <  1:3420440/61‑1 (MQ=255)
                          gTTTTTACGACCCCAGTTAACCATGTCATTGAGCTTTCCCATAAACACGTTTTTGTTAACTTCTTGCta                                             <  1:2772860/69‑2 (MQ=255)
                          gTTTTTACGACCCCAGTTAACCATGTCATTGAGCTTGCCCATAAACACGTTTTTGTTAACTTCTTGCTCCa                                           <  1:4041498/71‑1 (MQ=255)
                           tttttACGACCCCAGTTACCCATGTCATTGAGCTTGCCCATAAACACGTTTTTGTTAActtc                                                   <  1:819463/62‑1 (MQ=255)
                           tttttACGACCCCAGTTAACCATGTCATTGAGCTTGCCCATAAACACGTTTTTGTTAActtc                                                   <  1:3166936/62‑1 (MQ=255)
                             tttACGACCCCAGTTAACCATGTCATTGAGCTTGCCCATAAACACGTTTTTGTTAACTTCTTGCTCCa                                           >  1:1742555/1‑68 (MQ=255)
                             tttACGACCCCAGTTAACCATGTCATTGAGCTTGCCCATAAACACGTTTTTGTTAACTTCTTGCTCCa                                           >  1:3141378/1‑68 (MQ=255)
                              ttACGACCCCAGTTAACCATGTCATTGAGCTTGCCCATAAACACGTTTTTGTTAACTTCTTGCTCCAGAgg                                       >  1:690581/1‑71 (MQ=255)
                              ttACGACCCCAGTTAACCATGTCATTGAGCTTGCCCATAAACACGTTTTTGTTAACTTCTTGCTCCAGAgg                                       >  1:1966777/1‑71 (MQ=255)
                              ttACGACCCCAGTTAACCATGTCATTGAGCTTGCCCATAAACACGTTTTTGTTAACTTCTTGCTCCAGAgg                                       >  1:1471908/1‑71 (MQ=255)
                              ttACGACCCCAGTTAACCATGTCATTGAGCTTGCCCATAAACACGTTTTTGTTAACTTCTTGCTCCAGAgg                                       >  1:2865601/1‑71 (MQ=255)
                                 cGACCCCAGTTAACCATGTCATTGAGCTTGCCCATAAACACGTTTTTGTTAACTTCTTGCTCCAGAGGGt                                     <  1:4020683/70‑1 (MQ=255)
                                 cGACCCCAGTTAACCATGTCATTGAGCTTGCCCATAAACACGTTTTTGTTAACTTCTTGCTCCAGAGGGt                                     <  1:2219190/70‑1 (MQ=255)
                                    ccccAGTTAACCATGTCATTGAGCTTGCCCATAAACACGTTTTTGTTAACTTCTTGCTCCAGAGGGTCGGt                                 <  1:1092970/71‑1 (MQ=255)
                                    ccccAGTTAACCATGTCATTGAGCTTGCCCATAAACACGTTTTTGTTAACTTCTTGCTCCAGAGGGTCGGt                                 <  1:3096575/71‑1 (MQ=255)
                                    ccccAGTTAACCATGTCATTGAGCTTGCCCATAAACACGTTTTTGTTAACTTCTTGCTCCAGAGGGTCGGt                                 <  1:1969991/71‑1 (MQ=255)
                                    ccccAGTTAACCATGTCATTGAGCTTGCCCATAAACACGTTTTTGTTAACTTCTTGCTCCAGAGGGTCGGt                                 <  1:1844373/71‑1 (MQ=255)
                                    ccccAGTTAACCATGTCATTGAGCTTGCCCATAAACACGTTTTTGTTAACTTCTTGCTCCAGAGGGTCGGt                                 <  1:4055728/71‑1 (MQ=255)
                                    ccccAGTTAACCATGTCATTGAGCTTGCCCATAAACACGTTTTTGTTAACTTCTTGCTCCAGAGGGTCGGt                                 <  1:1177168/71‑1 (MQ=255)
                                       cAGTTAACCATGTCATTGAGCTTGCCCATAAACACGTTTTTGTTAACTTCTTGCTCCAGAGGGTCGGTTAc                              <  1:2861411/71‑1 (MQ=255)
                                       cAGTTAACCATGTCATTGAGCTTGCCCATAAACACGTTTTTGTTAACTTCTTGCTCCAGAGGGTCGGTTAc                              <  1:3345132/71‑1 (MQ=255)
                                       cAGTTAACCATGTCATTGAGCTTGCCCATAAACACGTTTTTGTTAACTTCTTGCTCCAGAGGGTCGGTTAc                              <  1:3788340/71‑1 (MQ=255)
                                          ttAACCATGTCATTGAGCTTGCCCATAAGCACGTTTTTGTTAACTTCTTGc                                               <  1:2074159/51‑1 (MQ=255)
                                          ttAACCATGTCATTGAGCTTGCCCATAAGCACGTTTTTGTTAACTTCTTGc                                               <  1:2117800/51‑1 (MQ=255)
                                          ttAACCATGTCATTGAGCTTGCCCATAAGCACGTTTTTGTTAACTTCTTGc                                               <  1:2024542/51‑1 (MQ=255)
                                          ttAACCATGTCATTGAGCTTGCCCATAAACACGtttt                                                             <  1:959995/37‑1 (MQ=255)
                                          ttAACCATGTCATTGAGCTTGCCCATAAACACGtttt                                                             <  1:2770251/37‑1 (MQ=255)
                                          ttAACCATGTCATTGAGCTTGCCCATAAACACGTTTTTGTTAActt                                                    <  1:2737422/46‑1 (MQ=255)
                                          ttAACCATGTCATTGAGCTTGCCCATAAACACGTTTTTGTTAACTTCTTGCTCCagag                                        <  1:2422948/58‑1 (MQ=255)
                                           tAACCATGTCATTGAGCTTGCCCATAAACACGTTTTTGTTAACTTCTTg                                                >  1:1486731/1‑49 (MQ=255)
                                           tAACCATGTCATTGAGCTTGCCCATAAACACGTTTTTGTTAACTTCTTg                                                >  1:2205847/1‑49 (MQ=255)
                                           tAACCATGTCATTGAGCTTGCCCATAAACACGTTTTTGTTAACTTCTTg                                                >  1:1699056/1‑49 (MQ=255)
                                           tAAACATGTCATTGAGCTTGCCCATAAACACGTTTTTGTTAActtctt                                                 >  1:2921390/1‑48 (MQ=255)
                                            aaCCATGTCATTGAGCTTGCCCATAAGCACGTTTTTGTTAACTTCTTGCt                                              >  1:1266106/1‑50 (MQ=255)
                                            aaCCATGTCATTGAGCTTGCCCATAAGCACGTTTTTGTTAACTTCTTGCTc                                             >  1:2695695/1‑51 (MQ=255)
                                            aaCCATGTCATTGAGCTTGCCCATAAGCACGTTTTTGTTAACTTCTTGCTc                                             >  1:2014739/1‑51 (MQ=255)
                                            aaCCATGTCATTGAGCTTGCCCATAAGCACGTTTTTGTTAACTTCTTGCTc                                             >  1:1870124/1‑51 (MQ=255)
                                            aaCCATGTCATTGAGCTTGCCCATAAGCACGTTTTTGTTAACTTCTTGCTc                                             >  1:2379085/1‑51 (MQ=255)
                                            aaCCATGTCATTGAGCTTGCCCATAAGCACGTTTTTGTTAACTTCTTGCTc                                             >  1:3938506/1‑51 (MQ=255)
                                            aaCCATGTCATTGAGCTTGCCCATAAGCACGTTTATGTTAACTTCTTGCTc                                             >  1:3637951/1‑51 (MQ=255)
                                                 tGTCATTGAGCTTGCCCATAAACACGTTTTTGTTAAct                                                     <  1:1029335/38‑1 (MQ=255)
                                                 tGTCATTGAGCTTGCCCATAAACACGTTTTTGTTAAct                                                     <  1:1789024/38‑1 (MQ=255)
                                                 tGTCATTGAGCTTGCCCATAAACACGTTTTTGTTAAct                                                     <  1:2033395/38‑1 (MQ=255)
                                                 tGTCATTGAGCTTGCCCATAAACACGTTTTTGTTAAct                                                     <  1:1697200/38‑1 (MQ=255)
                                                 tGTCATTGAGCTTGCCCATAAACACGTTTTTGTTAAct                                                     <  1:589492/38‑1 (MQ=255)
                                                 tGTCATTGAGCTTGCCCATAAACACGTTTTTGTTAAct                                                     <  1:25527/38‑1 (MQ=255)
                                                 tGTCATTGAGCTTGCCCATAAACACGTTTTTGTTAAct                                                     <  1:2273710/38‑1 (MQ=255)
                                                 tGTCATTGAGCTTGCCCATAAACACGTTTTTGTTAAct                                                     <  1:1429185/38‑1 (MQ=255)
                                                   tCATTGAGCTTGCCCATAAACACGTTTTTGTTAACTTCTTGCTCCAGAGGGTCGGTTACGAt                           >  1:761243/1‑62 (MQ=255)
                                                   tCATTGAGCTTGCCCATAAACACGTTTTTGTTAACTTCTTGCTCCAGAGGGTCGGTTACGAt                           >  1:2627228/1‑62 (MQ=255)
                                                   tCATTGAGCTTGCCCATAAACACGTTTTTGTTAACTTCTTGCTCCAGAGGGTCGGTTACGATCTCCTGCtt                  >  1:2912124/1‑71 (MQ=255)
                                                   tCATTGAGCTTGCCCATAAACACGTTTTTGTTAACTTCTTGCTCCAGAGGGTCGGTTACGATCTCCTGCtt                  >  1:489842/1‑71 (MQ=255)
                                                   tCATTGAGCTTGCCCATAAACACGTTTTTGTTAACTTCTTGCTCCAGAGGGTCGGTTACGATCTCCTGCtt                  >  1:3385554/1‑71 (MQ=255)
                                                   tCATTGAGCTTGCCCATAAACACGTTTTTGTTAACTTCTTGCTCCAGAGGGTCGGTTACGATCTCCTGCtt                  >  1:3202864/1‑71 (MQ=255)
                                                   tCATTGAGCTTGCCCATAAACACGTTTTTGTTAACTTCTTGCTCCAGAGGGTCGGTTACGATCTCCTGCtt                  >  1:3171944/1‑71 (MQ=255)
                                                   tCATTGAGCTTGCCATAAACACGTTTTTGTTAACTTCTTGCTCCAGAGGGTCGGTTACGAt                           >  1:2135976/1‑61 (MQ=255)
                                                      ttGAGCTTGCCCATAAACACGTTTTTGTTAACTTCTTGCTCCAGAGGGTCGGTTACGATCTCCTGCTTTTg               <  1:3294067/71‑1 (MQ=255)
                                                      ttGAGCTTGCCCATAAACACGTTTTTGTTAACTTCTTGCTCCAGAGGGTCGGTTACGATCTCCTGCTTTTg               <  1:2210931/71‑1 (MQ=255)
                                                       tGAGCTTGCCCATAAACACGTTTTTGTTAACTTCTTGCTCCAGAGGGTCGGTTACGATCTCCTGCTTTTg               >  1:2285011/1‑70 (MQ=255)
                                                       tGAGCTTGCCCATAAACACGTTTTTGTTAACTTCTTGCTCCAGAGGGTCGGTTACGATCTCCTGCTTTTg               >  1:2700962/1‑70 (MQ=255)
                                                       tGAGCTTGCCCATAAACACGTTTTTGTTAACTTCTTGCTCCAGAGGGTCGGTTACGATCTCCTGCTTTTg               >  1:1747224/1‑70 (MQ=255)
                                                       tGAGCTTGCCCATAAACACGTTTTTGTTAACTTCTTGCTCCAGAGGGTCGGTTACGATCTCCTGCTTTTg               <  1:1634181/70‑1 (MQ=255)
                                                       tGAGCTTGCCCATAAACACGTTTTTGTTAACTTCTTGCTCCAGAGGGTCGGTTACGATCTCCTGCTTTTg               >  1:3947396/1‑70 (MQ=255)
                                                        gAGCTTGCCCATAAACACGTTTTTGTTAACTTCTTGCTCCAGAGGGTCGGTTACGATCTCCTGCTTTTg               >  1:4059309/1‑69 (MQ=255)
                                                        gAGCTTGCCCATAAACACGTTTTTGTTAACTTCTTGCTCCAGAGGGTCGGTTACGATCTCCTGCTTTTg               >  1:3823685/1‑69 (MQ=255)
                                                        gAGCTTGCCCATAAACACGTTTTTGTTAACTTCTTGCTCCAGAGGGTCGGTTACGATCTCCTGCTTTTg               >  1:1835279/1‑69 (MQ=255)
                                                        gAGCTTGCCCATAAACACGTTTTTGTTAACTTCTTGCTCCAGAGGGTCGGTTACGATCTCCTGCTTTTg               >  1:3264755/1‑69 (MQ=255)
                                                        gAGCTTGCCCATAAACACGTTTTTGTTAACTTCTTGCTCCAGAGGGTCGGTTACGATCTCCTGCTTTTg               >  1:2581453/1‑69 (MQ=255)
                                                          gCTTGCCCATAAACACGTTTTTGTTAACTTCTTGCTCCAGAGGGTCGGTTACGATCTCCTGCTTTTGCAgg           <  1:2043600/71‑1 (MQ=255)
                                                              gCCCATAAACACGTTTTTGTTAACTTCTTGCTCCAGAGGGTCGGTTACGAt                           >  1:594997/1‑51 (MQ=255)
                                                              gCCCATAAACACGTTTTTGTTAACTTCTTGCTCCAGAGGGTCGGTTACGAt                           >  1:3794871/1‑51 (MQ=255)
                                                                 cATAAACACGTTTTTGTTAACTTCTTGCTCCAGAGGGTCGGTTACGATCTCCTGCTTTTGCAGGGGGt       <  1:1969784/68‑1 (MQ=255)
                                                                 cATAAACACGTTTTTGTTAACTTCTTGCTCCAGAGGGTCGGTTACGATCTCCTGCTTTTGCAGGGGGt       <  1:3923387/68‑1 (MQ=255)
                                                                   tAAACACGTTTTTGTTAACTTCTTGCTCCAGAGGGTCGGTTACGATCTCCTGCTTTTGCAgggg         <  1:719546/64‑1 (MQ=255)
                                                                   tAAACACGTTTTTGTTAACTTCTTGCTCCAGAGGGTCGGTTACGATCTCCTGCTTTTGCAgggg         <  1:1375399/64‑1 (MQ=255)
                                                                   tAAACACGTTTTTGTTAACTTCTTGCTCCAGAGGGTCGGTTACGATCTCCTGCTTTTGCAGGGGGTAACgg  >  1:906410/1‑71 (MQ=255)
                                                                   tAAACACGTTTTTGTTAACTTCTTGCTCCAGAGGGTCGGTTACGATCTCCTGCTTTTGCAGGGGGTAACgg  >  1:237164/1‑71 (MQ=255)
                                                                   tAAACACGTTTTTGTTAACTTCTTGCTCCAGAGGGTCGGTTACGATCTCCTGCTTTTGCAGGGGGTAACg   >  1:3002579/1‑70 (MQ=255)
                                                                   tAAACAAGTTTTTGTTAACTTCTTGCTCCAGAGGGTCGGTTACGATCTCCTGCTTTTGCAgggg         <  1:2124734/64‑1 (MQ=255)
                                                                    aaaCACGTTTTTGTTAACTTCTTGCTCCAGAGGGTCGGTTACGATCTCCTGCTTTTGCAg            >  1:395551/1‑60 (MQ=255)
                                                                    aaaCACGTTTTTGTTAACTTCTTGCTCCAGAGGGTCGGTTACGATCTCCTGCTTTTGCAg            >  1:847549/1‑60 (MQ=255)
                                                                    aaaCACGTTTTTGTTAACTTCTTGCTCCAGAGGGTCGGTTACGATCTCCTGCTTTTGCAg            >  1:2389714/1‑60 (MQ=255)
                                                                     aaCACGTTTTTGTTAACTTCTTGCTCCAGAGGGTCGGTTACGATCTCCTGCTTTTGCAggg          >  1:279372/1‑61 (MQ=255)
                                                                     aaCACGTTTTTGTTAACTTCTTGCTCCAGAGGGTCGGTTACGATCTCCTGCTTTTGCAggg          >  1:2696980/1‑61 (MQ=255)
                                                                     aaCACGTTTTTGTTAACTTCTTGCTCCAGAGGGTCGGTTACGATCTCCTGCTTTTGCAgg           >  1:3952341/1‑60 (MQ=255)
                                                                     aaCACGTTTTTGTTAACTTCTTGCTCCAGAGGGTCGGTTACGATCTCCTGCTATTGc              >  1:2573844/1‑57 (MQ=255)
                                                                      |                                                                   
AGACCGAAGTTATACGGCCAAATTGAGTTTTTACGACCCCAGTTAACCATGTCATTGAGCTTGCCCATAAACACGTTTTTGTTAACTTCTTGCTCCAGAGGGTCGGTTACGATCTCCTGCTTTTGCAGGGGGTAACGG  >  minE/1500418‑1500555

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