Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,260,317 T→C 13.7% I203V (ATT→GTT) ‡ 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,3170TC13.7% 228.4 / 23.7 95I203V (ATT→GTT) ‡torZtrimethylamine N‑oxide reductase system III, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base T (40/42);  new base C (9/4);  total (49/46)
Fisher's exact test for biased strand distribution p-value = 2.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

CAGGCAATTTTTAAGGTATTCAGCGGGTTCATTCCCCACAGCACCACAACCTGGCTGTTTTCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCACCTGTG  >  minE/1260250‑1260385
                                                                   |                                                                    
cAGGCAATTTTTAAGGTATTCAGCGGGTTCATTCCCCACAGCACCACAACCTGGCTGTTTTCGAGAATc                                                                     <  1:692/69‑1 (MQ=255)
cAGGCAATTTTTAAGGTATTCAGCGGGTTCATTCCCCACAGCACCACAACCTGGCTGTTTTCCAGAATCa                                                                    <  1:4729602/70‑1 (MQ=255)
    cAATTTTTAAGGTATTCAGCGGGTTCATTCCCCACAGCACCACAACCTGGCTGTTTTCCAGAATCAGCGGc                                                               <  1:3450325/71‑1 (MQ=255)
    cAATTTTCAGCGTATTGAGTGGGTTAGCACTCCACAGCACCACGACATCGCTATGTTCCAGCACCAGCGGc                                                               >  1:4677745/1‑71 (MQ=255)
     aaTTTTTAAGGTATTCAGCGGGTTCATTCCCCACAGCACCACAACCTGGCTGTTTTCCAGAATCAGCGGcc                                                              <  1:5603068/71‑1 (MQ=255)
     aaTTTTTAAGGTATTCAGCGGGTTCATTCCCCACAGCACCACAACCTGGCTGTTTTCCAGAATCAGCGGcc                                                              <  1:3640974/71‑1 (MQ=255)
     aaTTTTTAAGGTATTCAGCGGGTTCATTCCCCACAGCACCACAACCTGGCTGTTTTCCAGAATCAGCGGcc                                                              <  1:1686427/71‑1 (MQ=255)
       tttttAAGGTATTCAGCGGGTTCATTCCCCACAGCACCACAACCTGGGTGTTTTCCAGAATCAGCGGCCa                                                             <  1:2295390/70‑1 (MQ=255)
       tttttAAGGTATTCAGCGGGTTCATTCCCCACAGCACCACAACCTGGCTGTTTTCCAGAATCAGCGGc                                                               <  1:145547/68‑1 (MQ=255)
         tttAAGGTATTCAGCGGGTTCATTCCCCACAGCACCACAACCTGGCTGTTTTCCAGAATCAGCGGCCagg                                                           <  1:4674443/70‑1 (MQ=255)
         tttAAGGTATTCAGCGGGTTCATTCCCCACAGCACCACAACCTGGCTGTTTTCCAGAATCAGCGGCCag                                                            >  1:5569303/1‑69 (MQ=255)
            aaGGTATTCAGCGGGTTCATTCCCCACAGCACCACAACCTGGCTGTTTTCCAGAATCAGCGGCCaggag                                                         >  1:1400898/1‑69 (MQ=255)
             aGGTATTCAGCGGGTTCATTCCCCACAGCACCACAACCTGGCTGTTTTCCATAATCAGCGGCCAGGAGGTc                                                      >  1:3443078/1‑71 (MQ=255)
             aGGTATTCAGCGGGTTCATTCCCCACAGCACCACAACCTGGCTGTTTTCCAGAATCAGCGGCCAGGAGGTc                                                      >  1:711866/1‑71 (MQ=255)
             aGGTATTCAGCGGGTTCATTCCCCACAGCACCAAAACCTGGCTGTTTTCCAGAATCAGCGGCCAGGAGGt                                                       <  1:2277717/70‑1 (MQ=255)
                tATTCAGCGGGTTCATTCCCCACAGCACCACAACCTGGCTGTTTTCCAGAATCAGCGGCCAGGAGGTctgc                                                   <  1:4918965/71‑1 (MQ=255)
                 aTTCAGCGGGTTCATTCCCCACAGCACCACAACCTGGCTGTTTTCCAGAATCAGCGGCCAGGAGGTctgc                                                   >  1:4336443/1‑70 (MQ=255)
                  ttCAGCGGGTTCATTCCCCACAGCACCACAACCTGGCTGTTTTCCAGAATCAGCGGCCAGGAGGTctgctg                                                 <  1:5566174/71‑1 (MQ=255)
                    cAGCGGGTTCATTCCCCACAGCACCACAACCTGGCTGTTTTCCAGAATCAGCGGCCag                                                            >  1:3571836/1‑58 (MQ=255)
                    cAGCGGGTTCATTCCCCACAGCACCACAACCTGGCTGTTTTCCAGAATCAGCGGCCAGGAGGTCTGCTGt                                                >  1:4141498/1‑70 (MQ=255)
                     aGCGGGTTCATTCCCCACAGCACCACAACCTGGCTGTTTTCCAGAAt                                                                      >  1:375176/1‑47 (MQ=255)
                      gCGGGTTCATTCCCCACAGCACCACAACCTGGCTGTTTTCCAGAATCAGCGGCCagg                                                           <  1:3114892/57‑1 (MQ=255)
                      gCGGGTTCATTCCCCACAGCACCACAACCTGGCTGTTTTCCAGAATCAGCGGCCagg                                                           <  1:1378839/57‑1 (MQ=255)
                      gCGGGTTCATTCCCCACAGCACCACAACCTGGCTGTTTTCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCa                                             <  1:5631559/71‑1 (MQ=255)
                       cGGGTTCATTCCCCACAGCACCACAACCTGGCTGTTTTCCAGAATCAGCGGCCAGGAGGt                                                       >  1:3354627/1‑60 (MQ=255)
                       cGGGTTCATTCCCCACAGCACCACAACCTGGCTGTTTTCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCa                                             >  1:1336383/1‑70 (MQ=255)
                         ggTTCATTCCCCACAGCACCACAACCTGGCTGTTTTCCAGAATCAGCGGCCaggag                                                         >  1:3823476/1‑56 (MQ=255)
                         ggTTCATTCCCCACAGCACCACAACCTGGCTGTTTTCCAGAATCAGCGGCCaggag                                                         >  1:3723437/1‑56 (MQ=255)
                         ggTTCATTCCCCACAGCACCACAACCTGGCTGTTTTCCAGAATCAGCGGCCaggag                                                         >  1:1223074/1‑56 (MQ=255)
                         ggTTCATTCCCCACAGCACCACAACCTGGCTGTTTTCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATac                                          <  1:4160766/71‑1 (MQ=255)
                         ggTTCATTCCCCACAGCACCACAACCTGGCTGTTTTCCAGAATAAGCGGCCaggag                                                         >  1:3482302/1‑56 (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:5071132/1‑67 (MQ=255)
                              aTTCCCCACAGCACCACAACCTGGCTGTTTTCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATaca                                         >  1:2900489/1‑67 (MQ=255)
                              aTTCCCCACAGCACCACAACCTGGCTGTTTTCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATa                                           >  1:4038360/1‑65 (MQ=255)
                               ttCCCCACAGCACCACAACCTGGCTGTTTTCCAGAATCAGCGGc                                                               <  1:4256170/44‑1 (MQ=255)
                                tCCCCACAGCACCACAACCTGGCTGTTTTCCAGAATCAg                                                                   <  1:677708/39‑1 (MQ=255)
                                tCCCCACAGCACCACAACCTGGCTGTTTTCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTa                                   >  1:4327800/1‑71 (MQ=255)
                                tCCCCACAGCACCACAACCTGGCTGTTTTCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTa                                   >  1:4819010/1‑71 (MQ=255)
                                  cccACAGCACCACAACCTGGCTGTTTTCCAGAATCAGCg                                                                 >  1:5255032/1‑39 (MQ=255)
                                  cccACAGCACCACAACCTGGCTGTTTTCCAGAATCAGCGGCCa                                                             >  1:2272889/1‑43 (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)
                                           ccacAACCTGGCTGTTTTCCAGAATCAGCGGCCAGGAGGTCTGCTGt                                                <  1:5568185/47‑1 (MQ=255)
                                           ccacAACCTGGCTGTTTTCCAGAATCAGCGGCCAGGAGGTCTGCTGt                                                <  1:1189137/47‑1 (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:3638016/53‑1 (MQ=255)
                                                         ttttCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAAcc                            <  1:380373/53‑1 (MQ=255)
                                                         ttttCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACcaca                     >  1:2118921/1‑60 (MQ=255)
                                                         ttttCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACcaca                     >  1:2125503/1‑60 (MQ=255)
                                                         ttttCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCa              <  1:3013841/67‑1 (MQ=255)
                                                         ttttCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCa              <  1:4982701/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:4888777/1‑70 (MQ=255)
                                                           ttCCAGCACCAGCGGCCAACTGGTCTGCTGTTGATAAACTTCACTACCACCCACGACATACGGCATGATcg        <  1:3771724/71‑2 (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:1905987/68‑1 (MQ=255)
                                                             ccAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGATc         <  1:4622432/68‑1 (MQ=255)
                                                             ccAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCAcc      >  1:2463294/1‑71 (MQ=255)
                                                             ccAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCAc       <  1:3754421/70‑1 (MQ=255)
                                                             ccAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCAc       <  1:4802583/70‑1 (MQ=255)
                                                             ccAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCAc       >  1:1843578/1‑70 (MQ=255)
                                                             ccAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCAc       <  1:65811/70‑1 (MQ=255)
                                                               aGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGa                               <  1:2130681/44‑1 (MQ=255)
                                                                gAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACa                                 <  1:4680426/41‑1 (MQ=255)
                                                                 aaTCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCACCtgtg  >  1:636680/1‑71 (MQ=255)
                                                                 aaTCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCACCtgtg  >  1:2739886/1‑71 (MQ=255)
                                                                   |                                                                    
CAGGCAATTTTTAAGGTATTCAGCGGGTTCATTCCCCACAGCACCACAACCTGGCTGTTTTCCAGAATCAGCGGCCAGGAGGTCTGCTGTTCATACACCTCTACAGAACCGACCACATGCGGCATGATCACCTGTG  >  minE/1260250‑1260385

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