Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,260,263 A→G 6.1% L221L (TTA→CTA)  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,2630AG 6.1% 221.2 / 5.2 82L221L (TTA→CTA) torZtrimethylamine N‑oxide reductase system III, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base A (48/29);  new base G (2/3);  total (50/32)
Fisher's exact test for biased strand distribution p-value = 3.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

TGCCAGATTTTTTCAGCTGATGGAAGTATTCCAGCCCTTGCTCATCGGTACTGCTCCAGGCAATTTTTAAGGTATTCAGCGGGTTCATTCCCCACAGCACCACAACCTGGCTGTTTTCCAGAATCAGCGGCCAGGAGGT  >  minE/1260194‑1260332
                                                                     |                                                                     
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:656751/1‑71 (MQ=255)
    aGATTTTTTCAGCTGATGGAAGTATTCCAGCCCTTGCTCATCGGTACTGCTCCAGGCAATTTTTAAGGTAt                                                                  >  1:349523/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)
                 tGATGGAAGTATTCCAGCCCTTGCTCATCGGTACTGCTCCAGGCAATTTTTaa                                                                       <  1:4154378/53‑1 (MQ=255)
                  gATGGAAGTATTCCAGCCCTTGCTCATCGGTACTGCTCCAGGCAATTTTTaa                                                                       >  1:615133/1‑52 (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:668058/1‑52 (MQ=255)
                            ttCCAGCCCTTGCTCATCGGTACTGCTCCAGGCAATTTTTAAGGTATTCAGc                                                             >  1:1242440/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)
                              ccAGCCCTTGCTCATCGGTACTGCTCCAGGCAATTTTTaa                                                                       >  1:549112/1‑40 (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:2149879/1‑70 (MQ=255)
                               cAGCCCTTGCTCATCGGTACTGCTCCAGGCAATTTTTAAGGTATTCAGCGGGTTCATTCCCCACAGcacc                                        <  1:1289954/70‑1 (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:2656355/1‑71 (MQ=255)
                                      tGCTCATCGGTACTGCTCCAGGCAATTTTTAAGGTATTCAGCGGGTTCATTCCCCACAGCACCACAACCTg                                >  1:1314854/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:1832019/67‑1 (MQ=255)
                                                    attCCACGCAATTTTCAGCGTATTGAGTGGGTTAGCACTCCACAGCACCACGACATCGCTATGTTCCAg                    <  1:2550735/67‑1 (MQ=255)
                                                     ttCCACGCAATTTTCAGCGTATTGAGTGGGTTAGCACTCCACAGCACCACGACATCGCTATGTTCCAg                    <  1:3821885/67‑1 (MQ=255)
                                                     cTCCAGGCAATTTTTAAGGTATTCAGCGGGTTCATTCCCCACAg                                            <  1:1053278/44‑1 (MQ=255)
                                                     cTCCAGGCAATTTTTAAGGTATTCAGCGGGTTCATTCCCCACAg                                            <  1:2248890/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:1732167/71‑1 (MQ=255)
                                                               tttttAAGGTATTCAGCGGGTTCATTCCCCACAGCACCACAACCTGGCTGTTTTCCAGAATCAGCGGCCag       >  1:2586623/1‑71 (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:91989/70‑1 (MQ=255)
                                                                   tAAGGTATTCAGCGGGTTCATTCCCCACAGCACCACAACCTGGCTGTTTTCCAGAATCAGCGGCCaggag    <  1:4109742/70‑1 (MQ=255)
                                                                     aGGTATTCAGCGGGTTCATTCCCCACAGCACCACAACCTGGCTGTTTTCCAGAATCAGCGGCCAGGAGGt  <  1:2645756/70‑1 (MQ=255)
                                                                     |                                                                     
TGCCAGATTTTTTCAGCTGATGGAAGTATTCCAGCCCTTGCTCATCGGTACTGCTCCAGGCAATTTTTAAGGTATTCAGCGGGTTCATTCCCCACAGCACCACAACCTGGCTGTTTTCCAGAATCAGCGGCCAGGAGGT  >  minE/1260194‑1260332

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: