Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,260,341 C→G 15.1% E195Q (GAA→CAA)  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,3410CG15.1% 208.2 / 26.4 86E195Q (GAA→CAA) torZtrimethylamine N‑oxide reductase system III, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base C (39/34);  new base G (9/4);  total (48/38)
Fisher's exact test for biased strand distribution p-value = 3.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

GCGGGTTCATTCCCCACAGCACCACAACCTGGCTGTTTTCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGAATAATCGCCGC  >  minE/1260272‑1260409
                                                                     |                                                                    
gCGGGTTCATTCCCCACAGCACCACAACCTGGCTGTTTTCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCa                                                                     <  1:5631559/71‑1 (MQ=255)
 cGGGTTCATTCCCCACAGCACCACAACCTGGCTGTTTTCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCa                                                                     >  1:1336383/1‑70 (MQ=255)
   ggTTCATTCCCCACAGCACCACAACCTGGCTGTTTTCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATac                                                                  <  1:4160766/71‑1 (MQ=255)
    gTTCATTCCCCACAGCACCACAACCTGGCTGTTTTCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATaca                                                                 <  1:352442/71‑1 (MQ=255)
      tCATTCCCCACAGCACCACAACCTGGCTGTTTTCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACAcc                                                               >  1:1785152/1‑71 (MQ=255)
        aTTCCCCACAGCACCACAACCTGGCTGTTTTCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCa                                                                     <  1:3573935/63‑1 (MQ=255)
        aTTCCCCACAGCACCACAACCTGGCTGTTTTCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATaca                                                                 >  1:2900489/1‑67 (MQ=255)
        aTTCCCCACAGCACCACAACCTGGCTGTTTTCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATaca                                                                 >  1:5071132/1‑67 (MQ=255)
        aTTCCCCACAGCACCACAACCTGGCTGTTTTCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATa                                                                   >  1:4038360/1‑65 (MQ=255)
          tCCCCACAGCACCACAACCTGGCTGTTTTCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTa                                                           >  1:4819010/1‑71 (MQ=255)
          tCCCCACAGCACCACAACCTGGCTGTTTTCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTa                                                           >  1:4327800/1‑71 (MQ=255)
             ccACAGCACCACAACCTGGCTGTTTTCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAg                                                        <  1:4607476/71‑1 (MQ=255)
             ccACAGCACCACAACCTGGCTGTTTTCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAg                                                        <  1:5527836/71‑1 (MQ=255)
              cacaGCACCACGACATCGCTATGTTCCAGCACCAGCGGCCAACTGGTCTGCTGTTGATaaac                                                                <  1:1399712/62‑1 (MQ=255)
              cacaGCACCACGACATCGCTATGTTCCAGCACCAGCGGCCAACTGGTCTGCTGTTGATAAACttcac                                                           >  1:2175164/1‑65 (MQ=255)
                caGCACCACGACATCGCTATGTTCCAGCACCAGCGGCCAACTGGTCTGCTGTTGATaaact                                                               >  1:4406122/1‑60 (MQ=255)
                 aGCACCACGACATCGCTATGTTCCAGCACCAGCGGTCAACTGGTCTGCTGTTGATaaa                                                                 >  1:4046215/1‑56 (MQ=255)
                 aGCACCACGACATCGCTATGTTCCAGCACCAGCGGCCAACTGGTCTGCTGTTGATaaa                                                                 >  1:5234272/1‑56 (MQ=255)
                 aGCACCAAGACATCGCTATGTTCCAGCACCAGCGGCCAACTGGTCTGCTGTTGATaaa                                                                 >  1:4937699/1‑56 (MQ=255)
                  gCACCACAACCTGGCTGTTTTCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAAcc                                                    >  1:4129664/1‑70 (MQ=255)
                   caccacGACATCGCTATGTTCCAGCACCAGCGGCCAACTGGTCTGCTGTTGATaaa                                                                 >  1:3883834/1‑54 (MQ=255)
                      cacaACCTGGCTGTTTTCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACa                                                         <  1:3424540/61‑1 (MQ=255)
                      cacaACCTGGCTGTTTTCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACca                                               <  1:657866/71‑1 (MQ=255)
                       acaACCTGGCTGTTTTCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACcac                                              <  1:1199900/71‑1 (MQ=255)
                        caACCTGGCTGTTTTCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACctc                                                             <  1:2356329/55‑1 (MQ=255)
                        caACCTGGCTGTTTTCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACcaca                                             >  1:2621935/1‑71 (MQ=255)
                         aaCCTGGCTGTTTTCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACAt                                            >  1:745218/1‑71 (MQ=255)
                          aCCTGGCTGTTTTCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACcaca                                             <  1:5497130/69‑1 (MQ=255)
                          aCCTGGCTGTTTTCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATg                                           >  1:1607392/1‑71 (MQ=255)
                          aCCTGGCTGTTTTCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATg                                           <  1:1657138/71‑1 (MQ=255)
                            cTGGCTGTTTTCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCg                                         <  1:1819726/71‑1 (MQ=255)
                              ggCTGTTTTCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCgg                                        >  1:5348665/1‑70 (MQ=255)
                                   ttttCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAAcc                                                    <  1:380373/53‑1 (MQ=255)
                                   ttttCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAAcc                                                    <  1:3638016/53‑1 (MQ=255)
                                   ttttCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACcaca                                             >  1:2118921/1‑60 (MQ=255)
                                   ttttCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACcaca                                             >  1:2125503/1‑60 (MQ=255)
                                   ttttCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCa                                      <  1:4982701/67‑1 (MQ=255)
                                   ttttCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCa                                      <  1:3013841/67‑1 (MQ=255)
                                   ttttCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGa                                   >  1:2053805/1‑70 (MQ=255)
                                   tgttCCAGCACCAGCGGCCAACTGGTCTGCTGTTGATAAACTTCACTACCACCCACGACATACGGCATGa                                   <  1:3837518/68‑1 (MQ=255)
                                    tttCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATa                                                                   >  1:1862993/1‑37 (MQ=255)
                                    tttCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATa                                                                   >  1:4256130/1‑37 (MQ=255)
                                     ttCCAGCACCAGCGGCCAACTGGTCTGCTGTTGATAAACTTCACTACCACCCACGACATACGGCATGAt                                  >  1:3448164/1‑69 (MQ=255)
                                     ttCCAGCACCAGCGGCCAACTGGTCTGCTGTTGATAAACTTCACTACCACCCACGACATACGGCATGATcg                                <  1:5635217/71‑2 (MQ=255)
                                     ttCCAGCACCAGCGGCCAACTGGTCTGCTGTTGATAAACTTCACTACCACCCACGACATACGGCATGATcg                                <  1:3771724/71‑2 (MQ=255)
                                     ttCCAGCACCAGCGGCCAACTGGTCTGCTGTTGATAAACTTCACTACCACCCACGACATACGGCATGATcg                                >  1:4888777/1‑70 (MQ=255)
                                     ttCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACcaca                                             >  1:5386442/1‑58 (MQ=255)
                                     ttCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACcac                                              >  1:2792368/1‑57 (MQ=255)
                                     ttCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATg                                    >  1:5605029/1‑67 (MQ=255)
                                     ttCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATg                                    >  1:3667326/1‑67 (MQ=255)
                                     ttCCAGAATCAGCGGCCAGGAGGTCTGCAGTTCATACACCTCTACAGAACCGACCACATGCGGCATg                                    >  1:1734160/1‑67 (MQ=255)
                                       ccAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGATc                                 <  1:4622432/68‑1 (MQ=255)
                                       ccAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGATc                                 <  1:1905987/68‑1 (MQ=255)
                                       ccAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCAcc                              >  1:2463294/1‑71 (MQ=255)
                                       ccAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCAc                               <  1:3754421/70‑1 (MQ=255)
                                       ccAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCAc                               <  1:65811/70‑1 (MQ=255)
                                       ccAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCAc                               >  1:1843578/1‑70 (MQ=255)
                                       ccAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCAc                               <  1:4802583/70‑1 (MQ=255)
                                         aGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGa                                                       <  1:2130681/44‑1 (MQ=255)
                                          gAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACa                                                         <  1:4680426/41‑1 (MQ=255)
                                           aaTCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCACCtgtg                          >  1:2739886/1‑71 (MQ=255)
                                           aaTCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCACCtgtg                          >  1:636680/1‑71 (MQ=255)
                                             ccAGCGGCCAGGAGGTTTGCTGTTCATACACTTCCATCGAACCGACTACGCGCGGCAGGATCAc                               >  1:2469112/2‑64 (MQ=255)
                                             ccAGCGGCCAGGAGGTTTGCTGTTCATACACTTCCATCGAACCGACTACGCGCGGCAGGATCAc                               >  1:4428977/2‑64 (MQ=255)
                                             ccAGCGGCCAGGAGGTTTGCTGTTCATACACTTCCATCGAACCGACTACGCGCGGCAGGATCACCtgcg                          >  1:2004066/2‑67 (MQ=255)
                                                 cGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCACCtgt                           <  1:1127941/64‑1 (MQ=255)
                                                 cGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCACCtgt                           <  1:2812304/64‑1 (MQ=255)
                                                 cGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCACCTGTGCTGccc                    <  1:3558249/71‑1 (MQ=255)
                                                   gCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCg                  >  1:3156104/1‑71 (MQ=255)
                                                   gCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCg                  >  1:2095323/1‑71 (MQ=255)
                                                     cAGGAGGTTTGCTGTTCATACACTTCCATCGAACCGACTACGCGCGGCAGGATCACCAGCGCAGCACCgg                 >  1:1902837/1‑70 (MQ=255)
                                                      aggaggTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGtt               >  1:2762715/1‑71 (MQ=255)
                                                        gaggTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCAc                               >  1:2985231/1‑53 (MQ=255)
                                                         tggTCTGCTGTTGATAAACTTCACTACCACCCACGACATACGGCATGAt                                  >  1:4584112/2‑49 (MQ=255)
                                                          ggTTTGCTGTTCATACACTTCCATCGAACCGACTAAGCGCGGCAGGATCACCTGCGCAGCACAGGTAGag            >  1:866146/1‑69 (MQ=255)
                                                            tCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGaataa         <  1:3517141/71‑1 (MQ=255)
                                                             ctgctgTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGaataat        >  1:5492888/1‑71 (MQ=255)
                                                             ctgctgTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGaataat        >  1:2371063/1‑71 (MQ=255)
                                                              tgctgTTCATACACTTCCATCGAACCGACTACGCGCGGCAGGATCACCTGCGCAGCACCGGTAGagtaat        <  1:5009371/70‑1 (MQ=255)
                                                              tgctgTTCATACACCTCTACAGAACCGACCACATGCGGCATGATc                                 <  1:2666496/45‑1 (MQ=255)
                                                              tgctgTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGAATAATc       >  1:3645477/1‑71 (MQ=255)
                                                                 tgTTCATACACCTCTACAGAACCGACCACATGCgg                                        <  1:702544/35‑1 (MQ=255)
                                                                 tgTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGAATAATcgcc    >  1:181128/1‑71 (MQ=255)
                                                                   ttCATACACCTCTACAGCACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGaat           <  1:4559352/62‑1 (MQ=255)
                                                                   ttCATACACCTCTACAGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGaat           <  1:3231834/62‑1 (MQ=255)
                                                                   ttCATACACCTCTACAGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGAATAATcgccgc  <  1:1894208/71‑1 (MQ=255)
                                                                     cATACACCTCTACAGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGaataa         >  1:4765835/1‑62 (MQ=255)
                                                                     cATACACCTCTACAGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGaataa         >  1:440420/1‑62 (MQ=255)
                                                                     |                                                                    
GCGGGTTCATTCCCCACAGCACCACAACCTGGCTGTTTTCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCACCTGTGCTGCCCCGGTTGAATAATCGCCGC  >  minE/1260272‑1260409

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: