Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,260,264 G→C 6.3% T220T (ACC→ACG torZ ← trimethylamine N‑oxide reductase system III, catalytic subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,260,2640GC 6.3% 225.9 / 3.8 79T220T (ACC→ACGtorZtrimethylamine N‑oxide reductase system III, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base G (45/29);  new base C (2/3);  total (47/32)
Fisher's exact test for biased strand distribution p-value = 3.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

TGCCAGATTTTTTCAGCTGATGGAAGTATTCCAGCCCTTGCTCATCGGTACTGCTCCAGGCAATTTTTAAGGTATTCAGCGGGTTCATTCCCCACAGCACCACAACCTGGCTGTTTTCCAGAATCAGCGGCCAGGAGGTCT  >  minE/1260194‑1260334
                                                                      |                                                                      
tGCCAGATTTTTTCAGCTGATGGAAGTATTCCAGCCCTTGCTCATCGGTACTGCTCCAGGCAATTTTTAAg                                                                        >  1:3825630/1‑71 (MQ=255)
tGCCAGATTTTTTCAGCTGATGGAAGTATTCCAGCCCTTGCTCATCGGTACTGCTCCAGGCAATTTTTAAg                                                                        >  1:2003855/1‑71 (MQ=255)
 gCCAGATTTTTTCAGCTGATGGAAGTATTCCAGCCCTTGCTCATCGGTACTGCTCCAGGCAATTTTTAAg                                                                        <  1:972653/70‑1 (MQ=255)
   cAGATTTTTTCAGCTGATGGAAGTATTCCAGCCCTTGCTCATCGGTACTGCTCCAGGCAATTTTTAAGGTa                                                                     <  1:2561546/71‑1 (MQ=255)
    aGATTTTTTCAGCTGATGGAAGTATTCCAGCCCTTGCTCATCGGTACTGCTCCAGGCAATTTTTAAGGTAt                                                                    >  1:349523/1‑71 (MQ=255)
    aGATTTTTTCAGCTGATGGAAGTATTCCAGCCCTTGCTCATCGGTACTGCTCCAGGCAATTTTTAAGGTAt                                                                    >  1:656751/1‑71 (MQ=255)
     gATTTTTTCAGCTGATGGAAGTATTCCAGCCCTTGCTCATCGGTACTTCTCCAGGCAATTTTTAAGGTAt                                                                    >  1:3996536/1‑70 (MQ=255)
         ttttCAGCTGATGGAAGTATTCCAGCCCTTGCTCATCGGTACTGCTCCAGGCAATTTTTAAGGTATTCAg                                                                <  1:3159600/70‑1 (MQ=255)
          tttCAGCTGATGGAAGTATTCCAGCCCTTGCTCATCGGTACTGCTCCAGGCAATTTTTAAGGTATTCAGCg                                                              >  1:2553243/1‑71 (MQ=255)
          tttCAGCTGATGGAAGTATTCCAGCCCTTGCTCATCGGTACTGCTCCAGGCAATTTTTAAGGTATTCAGCg                                                              >  1:1779133/1‑71 (MQ=255)
              aGCTGATGGAAGTATTCCAGCCCTTGCTCATCGGTACTGCTCCAGGCAATTTTTAAGGTATTCAGCGGGtt                                                          <  1:1694908/71‑1 (MQ=255)
              aGCTGATGGAAGTATTCCAGCCCTTGCTCATCGGTACTGCTCCAGGCAATTTTTAAGGTATTCAGCGGGtt                                                          <  1:2379281/71‑1 (MQ=255)
              aGCTGATGGAAGTATTCCAGCCCTTGCTCATCGGTACTGCTCCAGGCAATTTTTAAGGTATTCAGCGGGt                                                           <  1:3855275/70‑1 (MQ=255)
                   aTGGAAGTATTCCAGCCCTTGCTCATCGGTACTGCTCCAGGCAATTTTTAAGGTATTCAGCGGGTTCAtt                                                      >  1:1681416/1‑70 (MQ=255)
                    tGGAAGTATTCCAGCCCTTGCTCATCGGTACTGCTCCAGGCAATTTTTAAGGTATTCAGCGGGTTCATTcc                                                    >  1:4041160/1‑71 (MQ=255)
                     ggAAGTATTCCAGCCCTTGCTCATCGGTACTGCTCCAGGCAATTTTGAGCGTGTTCAGCGGGTTCAtcc                                                     >  1:1953086/1‑67 (MQ=255)
                      gAAGTATTCCAGCCCTTGCTCATCGGTACTGCTCCAGGCAATTTTTAAGGTATTCAGCGGGTTCATTccc                                                   <  1:3975682/70‑1 (MQ=255)
                       aaGTATTCCAGCCCTTGCTCATCGGTACTGCTCCAGGCAATTTTTAAGGTATTCAGCGGGTTCATTCCCca                                                 <  1:1761453/71‑1 (MQ=255)
                        aGTATTCCAGCCCTTGCTCATCGGTACTGCTCCAGGCAATTTTTAAGGTATTCAGCGGGTTCATTCCCca                                                 >  1:517849/1‑70 (MQ=255)
                           aTTCCAGCCCTTGCTCATAGGTACTGCTCCAGGCAATTTTTAAg                                                                        >  1:1856159/1‑44 (MQ=255)
                            ttCCAGCCCTTGCTCATCGGTACTGCTCCAGGCAATTTTTAAGGTa                                                                     <  1:794059/46‑1 (MQ=255)
                            ttCCAGCCCTTGCTCATCGGTACTGCTCCAGGCAATTTTTAAGGTATTCAGc                                                               >  1:2844393/1‑52 (MQ=255)
                            ttCCAGCCCTTGCTCATCGGTACTGCTCCAGGCAATTTTTAAGGTATTCAGc                                                               >  1:1242440/1‑52 (MQ=255)
                            ttCCAGCCCTTGCTCATCGGTACTGCTCCAGGCAATTTTTAAGGTATTCAGc                                                               >  1:668058/1‑52 (MQ=255)
                            ttCCAGCCCTTGCTCATCGGTACTGCTCCAGGCAATTTTTAAGGTATTCAGCGGGTTCAt                                                       <  1:3309137/60‑1 (MQ=255)
                            ttCCAGCCCTTGCTCATCGGTACTGCTCCAGGCAATTTTTAAGGTATTCAGCGGGTTCATTcc                                                    >  1:2033046/1‑63 (MQ=255)
                            ttCCAGCCCTTGCTCATCGGTACTGCTCCAGGCAATTTTTAAGGTATTCAGCGGGTTCATTCCCCACAGca                                            >  1:2859108/1‑71 (MQ=255)
                              ccAGCCCTTGCTCATCGGTACTGCTCCAGGCAATTTTTAAGGTATTCAGCgg                                                             >  1:196136/1‑52 (MQ=255)
                              ccAGCCCTTGCTCATCGGTACTGCTCCAGGCAATTTTTAAGGTATTCAGCGGGTTCa                                                        >  1:3009869/1‑57 (MQ=255)
                              ccAGCCCTTGCTCATCGGTACTGCTCCAGGCAATTTTTAAGGTATTCAGCGGGTTCATTCCCCACAGcac                                           >  1:869309/1‑70 (MQ=255)
                              ccAGCCCTTGCTCATCGGTACTGCTCCAGGCAATTTTTAAGGTATTCAGCGGGTTCATTCCCCACAGcac                                           >  1:993297/1‑70 (MQ=255)
                              ccAGCCCTTGCTCATCGGTACTGCTCCAGGCAATTTTTAAGGTATTCAGCGGGTTCATTCCCCACAGcac                                           >  1:1332250/1‑70 (MQ=255)
                              ccAGCCCTTGCTCATCGGTACTGCTCCAGGCAATTTTTAAGGTATTCAGCGGGTTCATTCCCCACAGca                                            >  1:726931/1‑69 (MQ=255)
                               cAGCCCTTGCTCATCGGTACTGCTCCAGGCAATTTTTAAGGTATTCAGCGGGTTCATTccc                                                   <  1:2669945/61‑1 (MQ=255)
                               cAGCCCTTGCTCATCGGTACTGCTCCAGGCAATTTTTAAGGTATTCAGCGGGTTCATTCCCCACAGcacc                                          <  1:1289954/70‑1 (MQ=255)
                               cAGCCCTTGCTCATCGGTACTGCTCCAGGCAATTTTTAAGGTATTCAGCGGGTTCATTCCCCACAGcacc                                          >  1:2149879/1‑70 (MQ=255)
                                   ccTTGCTCATCGGTACTGCTCCAGGCAATTTTTAAGGTATTCAGCGGGTTCATTCCCCACAGCACCACAAc                                     >  1:3968348/1‑71 (MQ=255)
                                    cTTGCTCATCGGTACTGCTCCAGGCAATTTTTAAGGTATTCAGCGGGtt                                                          <  1:3531775/49‑1 (MQ=255)
                                    cTTGCTCATCGGTACTGCTCCAGGCAATTTTTAAGGTATTCAGCGGGTTCAt                                                       >  1:2605433/1‑52 (MQ=255)
                                     ttGCTCATCGGTACTGCTCCAGGCAATTTTTAAGGt                                                                      <  1:285115/36‑1 (MQ=255)
                                      tGCTCATCGGTACTGCTCCAGGCAATTTTTAAGGTATTCAGCGGGTTCATTCCCCACAGCACCACAACCTg                                  >  1:1314854/1‑71 (MQ=255)
                                      tGCTCATCGGTACTGCTCCAGGCAATTTTTAAGGTATTCAGCGGGTTCATTCCCCACAGCACCACAACCTg                                  >  1:2656355/1‑71 (MQ=255)
                                       gCTCATTGGTACTGCTACAGGCAATTTTGAGCGTGTTCAGCGGGTTCATCCCCCACAGCACCACCACCTgg                                 >  1:2089805/1‑71 (MQ=255)
                                       gCTCATCGGTACTGCTCCAGGCAATTTTTAAGGTATTCa                                                                 >  1:81478/1‑39 (MQ=255)
                                       gCTCATCGGTACTGCTCCAGGCAATTTTTAAGGTATTCAGCGGGTTCATTCCCCACAGCACCACAACCTgg                                 >  1:302596/1‑71 (MQ=255)
                                         tCATCGGTACTGCTCCAGGCAATTTTTAAGGTATTCAGCGGGtt                                                          <  1:686505/44‑1 (MQ=255)
                                           aTCGGTACTGCTCCAGGCAATTTTTAAGGTATTCAGCGGGTTCCTTCCCca                                                 >  1:3453387/1‑51 (MQ=255)
                                           aTCGGTACTGCTCCAGGCAATTTTTAAGGTATTCAGCGGGTTCATTCCCca                                                 >  1:3431545/1‑51 (MQ=255)
                                             cGGTACTGCTCCAGGCAATTTTTAAGGTATTCAGCgg                                                             <  1:4459859/37‑1 (MQ=255)
                                             cGGTACTGCTCCAGGCAATTTTTAAGGTATTCAGCGGGTTCATTCCCCACAGCACCACAACCTGGCTGttt                           >  1:926221/1‑71 (MQ=255)
                                                tACTGCTCCAGGCAATTTTTAAGGTATTCAGCGGGtt                                                          <  1:564853/37‑1 (MQ=255)
                                                 aCTGCTCCAGGCAATTTTTAAGGTATTCAGCGGGTTCATTCCCcaca                                               >  1:979210/1‑47 (MQ=255)
                                                   tGCTCCAGGCAATTTTTAAGGTATTCAGCGGGTTCATTccc                                                   <  1:423902/41‑1 (MQ=255)
                                                    gCTCCAGGCAATTTTTAAGGTATTCAGCGGGTTCATTCCCCACAGCACCACAACCTGGCTGTTTTCCAg                      >  1:2495951/1‑69 (MQ=255)
                                                    gCTCCAGGCAATTTTTAAGGTATTCAGCGGGTTCATTCCCCACAGCACCACAACCTGGCTGTTTTCCAGaa                    >  1:4697953/1‑71 (MQ=255)
                                                    gCTCCAGGCAATTTTTAAGGTATTCAGCGGGTTCATTCCCCACAGCACCACAACCTGGCTGTTTTCCAGaa                    >  1:4658312/1‑71 (MQ=255)
                                                    gCTCCAGGCAATTTTTAAGGTATTCAGCGGGTTCATTCCCCACAGCACCACAACCTGGCTGTTTTCCAGa                     >  1:652890/1‑70 (MQ=255)
                                                    attCCACGCAATTTTCAGCGTATTGAGTGGGTTAGCACTCCACAGCACCACGACATCGCTATGTTCCAg                      <  1:2550735/67‑1 (MQ=255)
                                                    attCCACGCAATTTTCAGCGTATTGAGTGGGTTAGCACTCCACAGCACCACGACATCGCTATGTTCCAg                      <  1:1832019/67‑1 (MQ=255)
                                                     ttCCACGCAATTTTCAGCGTATTGAGTGGGTTAGCACTCCACAGCACCACGACATCGCTATGTTCCAg                      <  1:3821885/67‑1 (MQ=255)
                                                     cTCCAGGCAATTTTTAAGGTATTCAGCGGGTTCATTCCCCACAg                                              <  1:2248890/44‑1 (MQ=255)
                                                     cTCCAGGCAATTTTTAAGGTATTCAGCGGGTTCATTCCCCACAg                                              <  1:1053278/44‑1 (MQ=255)
                                                     cTCCAGGCAATTTTTAAGGTATTCAGCGGGTTCATTCCCCACAGcaccac                                        >  1:2320586/1‑50 (MQ=255)
                                                     cTCCAGGCAATTTTTAAGGTATTCAGCGGGTTCATTCCCCACAGCACCACAACCTGGCTGTTTTCCAGa                     <  1:2045695/69‑1 (MQ=255)
                                                     cTCCAGGCAATTTTTAAGGTATTCAGCGGGTACATTCCCCACAGCACCACAAc                                     >  1:2576371/1‑53 (MQ=255)
                                                     cTCCAGGAAATTTTTAAGGTATTCAGCGGGTTCATTCCCCACAGCACCACaa                                      >  1:4329194/1‑52 (MQ=255)
                                                      tCCAGGCAATTTTTAAGGTATTCAGCGGGTTCATTCCCCACAGCACCACAACCTGGCTGTTTTCCAGAAt                   >  1:917040/1‑70 (MQ=255)
                                                        cAGGCAATTTTTAAGGTATTCAGCGGGTTCATTCCCCACAGCACCACAACCTGGCTGTTTTCCAGAATc                  <  1:911493/69‑1 (MQ=255)
                                                        cAGGCAATTTTTAAGGTATTCAGCGGGTTCATTCCCCACAGCACCACAACCTGGCTGTTTTCCAGAATCAg                >  1:2762794/1‑71 (MQ=255)
                                                           gCAATTTTTAAGGTATTCAGCGGGTTCATTCCCCACAGCACCACAACCTGGCTg                              <  1:3645029/54‑1 (MQ=255)
                                                            cAATTTTTAAGGTATTCAGCGGGTTCATTCCCCACAGCACCACAACCTGGCTGTTTTCCAGAATCAGCGGc            >  1:3199586/1‑71 (MQ=255)
                                                             aaTTTTTAAGGTATTCAGCGGGTTCATTCCCCACAGCACCACAACCTGGCTGTTTTCCAGAATCAGCGGcc           <  1:1271933/71‑1 (MQ=255)
                                                               tttttAAGGTATTCAGCGGGTTCATTCCCCACAGCACCACAACCTGGCTGtt                            <  1:4597774/52‑1 (MQ=255)
                                                               tttttAAGGTATTCAGCGGGTTCATTCCCCACAGCACCACAACCTGGCTGtt                            <  1:2998249/52‑1 (MQ=255)
                                                               tttttAAGGTATTCAGCGGGTTCATTCCCCACAGCACCACAACCTGGCTGTTTTCCAGAATCAGCGGCCag         >  1:2586623/1‑71 (MQ=255)
                                                               tttttAAGGTATTCAGCGGGTTCATTCCCCACAGCACCACAACCTGGCTGTTTTCCAGAATCAGCGGCCag         <  1:1732167/71‑1 (MQ=255)
                                                                ttttAAGGTATTCAGCGGGTTCATTCCCCACAGCACCACAACCTGGCTGTTTTCCAGaa                    >  1:1924374/1‑59 (MQ=255)
                                                                 tttAAGGTATTCAGCGGGTTCATTCCCCACAGcacca                                         >  1:1992032/1‑37 (MQ=255)
                                                                 tttAAGGTATTCAGCGGGTTCATTCCCCACAGcacca                                         >  1:2004359/1‑37 (MQ=255)
                                                                 tttAAGGTATTCAGCGGGTTCATTCCCCACAGCACCACAACCTGGCTGTTTTCCAGAATCAGCg              >  1:2334487/1‑64 (MQ=255)
                                                                   tAAGGTATTCAGCGGGTTCATTCCCCACAGCACCACAACCTGGCTGTTTTCCAGAATCAGCGGCCaggag      <  1:4109742/70‑1 (MQ=255)
                                                                   tAAGGTATTCAGCGGGTTCATTCCCCACAGCACCACAACCTGGCTGTTTTCCAGAATCAGCGGCCaggag      <  1:91989/70‑1 (MQ=255)
                                                                     aGGTATTCAGCGGGTTCATTCCCCACAGCACCACAACCTGGCTGTTTTCCAGAATCAGCGGCCAGGAGGt    <  1:2645756/70‑1 (MQ=255)
                                                                      ggTATTCAGCGGGTTCATTCCCCACAGCACCACAACCTGGCTGTTTTCCAGAATCAGCGGCCAGGAGGTct  >  1:264036/1‑71 (MQ=255)
                                                                      |                                                                      
TGCCAGATTTTTTCAGCTGATGGAAGTATTCCAGCCCTTGCTCATCGGTACTGCTCCAGGCAATTTTTAAGGTATTCAGCGGGTTCATTCCCCACAGCACCACAACCTGGCTGTTTTCCAGAATCAGCGGCCAGGAGGTCT  >  minE/1260194‑1260334

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: