Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,259,817 C→G 19.2% V369V (GTG→GTC 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,259,8170CG19.2% 169.8 / 32.4 78V369V (GTG→GTCtorZtrimethylamine N‑oxide reductase system III, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base C (29/34);  new base G (8/7);  total (37/41)
Fisher's exact test for biased strand distribution p-value = 7.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

GGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGC  >  minE/1259748‑1259887
                                                                     |                                                                      
ggAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTc                                                                        <  1:1265038/70‑1 (MQ=255)
 gAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCa                                                                       <  1:1543201/70‑1 (MQ=255)
 gAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCAc                                                                      <  1:423356/71‑1 (MQ=255)
        gATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATc                                                               >  1:311156/1‑71 (MQ=255)
         aTAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATcc                                                              <  1:3902408/71‑1 (MQ=255)
             gagaAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGc                                                                  <  1:1773659/63‑1 (MQ=255)
                  aCCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAgt                                                           >  1:1623018/1‑65 (MQ=255)
                   ccAAAACCGCCGCCCGGTGTACCAATCTGCCCTATCATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTg                                                    <  1:1866248/71‑1 (MQ=255)
                    cAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGt                                                   >  1:237533/1‑71 (MQ=255)
                    cAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGt                                                   >  1:1587592/1‑71 (MQ=255)
                      aaaCCGCCGCCGGGTGTGCCGATTTGCCCCAACATTGCTGCCAGCGTGACGATCATCCAGTGTTTTTgct                                                  >  1:1019307/1‑68 (MQ=255)
                          cgccgcCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCAc                                                                      >  1:2249374/1‑46 (MQ=255)
                          cgccgcCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCa                                             >  1:956166/1‑71 (MQ=255)
                          cgccgcCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCa                                             <  1:2071491/71‑1 (MQ=255)
                            ccgccCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCAc                                                                      >  1:1219358/1‑44 (MQ=255)
                               cccGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATcc                                                              <  1:814624/49‑1 (MQ=255)
                               cccGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATAtt                                         <  1:507695/70‑1 (MQ=255)
                                 cGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGttt                                                       <  1:1074327/54‑1 (MQ=255)
                                 cGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGc                                       >  1:589207/1‑70 (MQ=255)
                                 cGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCt                                      >  1:1012164/1‑71 (MQ=255)
                                 cGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCt                                      >  1:1824658/1‑71 (MQ=255)
                                  ggTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTg                                     <  1:1504740/71‑1 (MQ=255)
                                  ggTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTg                                     >  1:1156444/1‑71 (MQ=255)
                                   gtgtACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGtt                                                  >  1:2807651/1‑57 (MQ=255)
                                     gtACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGgcg                                  >  1:3950268/1‑71 (MQ=255)
                                      tACCAATCTGCCCTAACATCGCCTCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGgcg                                  <  1:1955581/70‑1 (MQ=255)
                                        ccAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGtt                                                  <  1:332172/52‑1 (MQ=255)
                                        ccAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGt                                                   <  1:1372389/51‑1 (MQ=255)
                                           atttGCCCCAACATTGCTGCCAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCAt                            >  1:1223111/4‑71 (MQ=255)
                                           aTCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAgtgt                                                         >  1:233247/1‑42 (MQ=255)
                                           aTCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGtt                                                  <  1:2919854/49‑1 (MQ=255)
                                           aTCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGtt                                                  <  1:3436580/49‑1 (MQ=255)
                                            tCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCa                                                >  1:2617609/1‑50 (MQ=255)
                                              tGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTg                               >  1:3121549/1‑65 (MQ=255)
                                              tGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTg                               >  1:3416274/1‑65 (MQ=255)
                                               gCCCCAACATTGCTGCCAGCGTGACGATCATCCAGTGTTTTTgc                                                   <  1:3303665/44‑2 (MQ=12)
                                               gCCCCAACATTGCTGCCAGCGTGACGATCATCCAGTGTTTTTgc                                                   <  1:3576895/44‑2 (MQ=12)
                                                cccTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATAtt                                         <  1:3882319/53‑1 (MQ=255)
                                                 ccTAACATCGCCGCCAGCGTCACCAGCATCCAgtgt                                                         <  1:2709702/36‑1 (MQ=255)
                                                 ccTAACATCGCCGCCAGCGTCACCAGCATCCAgtgt                                                         <  1:1739990/36‑1 (MQ=255)
                                                 ccTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGt                                                   <  1:3317346/42‑1 (MQ=255)
                                                 ccTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGacgc                                 >  1:2492214/1‑60 (MQ=255)
                                                 ccTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCa                      <  1:843594/71‑1 (MQ=255)
                                                 ccTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCa                      <  1:3514365/71‑1 (MQ=255)
                                                   tAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTg                                        <  1:1012511/51‑1 (MQ=255)
                                                   caaCATTGCTGCCAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTg                               <  1:2524134/59‑1 (MQ=255)
                                                   caaCATTGCTGCCAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGc                    >  1:3033525/2‑71 (MQ=255)
                                                    aaCATTGCTGCCAGCGTGACGATCATCCAGTGTTTTTGCTCACCaaa                                           >  1:2259757/1‑45 (MQ=12)
                                                    aaCATCGCCGCCAGCGTCACCAGCATCCAGTGttt                                                       <  1:1163994/35‑1 (MQ=255)
                                                    aaCATCGCCGCCAGCGTCACCAGCATCCAGTGttt                                                       <  1:1159233/35‑1 (MQ=255)
                                                    aaCATCGCCGCCAGCGTCACCAGCATCCAGTGttt                                                       <  1:3115781/35‑1 (MQ=255)
                                                    aaCATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGtt                                                  >  1:3326263/1‑40 (MQ=255)
                                                     aCATTGCTGCCAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTgca                             >  1:3921167/1‑58 (MQ=255)
                                                     aCATTGCTGCCAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCt                  >  1:2520160/1‑71 (MQ=255)
                                                     aCATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTg                               >  1:3673769/1‑58 (MQ=255)
                                                      cATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCt                                                     >  1:1182637/1‑35 (MQ=255)
                                                      cATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTg                 <  1:2655954/71‑1 (MQ=255)
                                                          gctgccAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCa              <  1:3094874/67‑1 (MQ=255)
                                                          gccgccAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGc                                       >  1:3002658/1‑45 (MQ=255)
                                                          gccgccAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGc                                       >  1:1904084/1‑45 (MQ=255)
                                                          gccgccAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCAt             >  1:2216343/1‑71 (MQ=255)
                                                          gccgccAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCAt             >  1:2221078/1‑71 (MQ=255)
                                                           ctgccAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCAt             <  1:2241161/68‑1 (MQ=255)
                                                           ctgccAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCAt             >  1:3379490/3‑70 (MQ=255)
                                                             gccAGCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCATCAg          >  1:1137738/1‑71 (MQ=255)
                                                             gccAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAg          >  1:2827844/1‑71 (MQ=255)
                                                              ccAGCCTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTg                 <  1:1938514/63‑1 (MQ=255)
                                                               cAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCa        <  1:1910719/71‑1 (MQ=255)
                                                               cAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCa        <  1:235307/71‑1 (MQ=255)
                                                               cAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCa        <  1:462589/71‑1 (MQ=255)
                                                               cAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCa        >  1:3621450/1‑71 (MQ=255)
                                                               cAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCa        <  1:1507521/71‑1 (MQ=255)
                                                               cAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCa        <  1:3841247/71‑1 (MQ=255)
                                                                aGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCAt       >  1:792378/1‑71 (MQ=255)
                                                                 gCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCATCAGCAtg      <  1:480872/71‑2 (MQ=255)
                                                                 gCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCATCAGCAtg      <  1:1881162/71‑2 (MQ=255)
                                                                  cGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATtg     >  1:1084140/1‑71 (MQ=255)
                                                                   gTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATtgt    <  1:415203/71‑1 (MQ=255)
                                                                   gTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATtg     <  1:371696/70‑1 (MQ=255)
                                                                     cACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGc  >  1:2726561/1‑71 (MQ=255)
                                                                     |                                                                      
GGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGC  >  minE/1259748‑1259887

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: