Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,259,808 C→T 16.4% A372A (GCG→GCA 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,8080CT16.4% 149.9 / 18.0 61A372A (GCG→GCAtorZtrimethylamine N‑oxide reductase system III, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base C (23/28);  new base T (6/4);  total (29/32)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

GCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCAT  >  minE/1259739‑1259876
                                                                     |                                                                    
gCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATcgccg                                                                     >  1:1662212/1‑71 (MQ=255)
gCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATcgccg                                                                     >  1:3842005/1‑71 (MQ=255)
    ccGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAg                                                                 <  1:744832/71‑1 (MQ=255)
         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)
                   tAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCa                                                                  >  1:1484681/1‑55 (MQ=255)
                      gagaAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGc                                                       <  1:1773659/63‑1 (MQ=255)
                           aCCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAgt                                                >  1:1623018/1‑65 (MQ=255)
                            ccAAAACCGCCGCCCGGTGTACCAATCTGCCCTATCATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTg                                         <  1:1866248/71‑1 (MQ=255)
                             cAAAACCGCCGCCGGGTGTGCCGATTTGCCCCAACATTGCTGCCAg                                                                 <  1:1463004/46‑1 (MQ=16)
                             cAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGt                                        >  1:237533/1‑71 (MQ=255)
                             cAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGt                                        >  1:1587592/1‑71 (MQ=255)
                              aaaaCCGCCGCCCGGTGTACCAATCTGCCCTAACATcgccgc                                                                    >  1:855298/1‑42 (MQ=255)
                               aaaCCGCCGCCGGGTGTGCCGATTTGCCCCAACATTGCTGCCAGCGTGACGATCATCCAGTGTTTTTgct                                       >  1:1019307/1‑68 (MQ=255)
                               aaaCCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCa                                                                  >  1:3627077/1‑43 (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:3115781/35‑1 (MQ=255)
                                                             aaCATCGCCGCCAGCGTCACCAGCATCCAGTGttt                                            <  1:1163994/35‑1 (MQ=255)
                                                             aaCATCGCCGCCAGCGTCACCAGCATCCAGTGttt                                            <  1:1159233/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)
                                                                   gccgccAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGc                            >  1:3002658/1‑45 (MQ=255)
                                                                   gccgccAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGc                            >  1:1904084/1‑45 (MQ=255)
                                                                   gccgccAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCAt  >  1:2221078/1‑71 (MQ=255)
                                                                   gccgccAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCAt  >  1:2216343/1‑71 (MQ=255)
                                                                     |                                                                    
GCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCAT  >  minE/1259739‑1259876

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: