Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,259,799 T→C 12.5% L375L (TTA→TTG 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,7990TC12.5% 145.8 / 11.2 56L375L (TTA→TTGtorZtrimethylamine N‑oxide reductase system III, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base T (24/25);  new base C (2/5);  total (26/30)
Fisher's exact test for biased strand distribution p-value = 4.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

TCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCA  >  minE/1259732‑1259867
                                                                   |                                                                    
tCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCt                                                                      <  1:2559601/68‑1 (MQ=255)
tCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAAc                                                                   <  1:2414199/71‑1 (MQ=255)
tCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAAc                                                                   <  1:1932896/71‑1 (MQ=255)
    gTTACCACCATTGGCAAAATGGTAAGAAAGACCAAAACCGCCGCCGGGTGTGCCGATTTGCCCCAACAttg                                                               <  1:372648/71‑3 (MQ=255)
       gCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATcgccg                                                            >  1:3842005/1‑71 (MQ=255)
       gCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATcgccg                                                            >  1:1662212/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)
                   aaaaTGGTAAGAAAGACCAAAACCGCCGCCGGGTGTGCCGATTTGCCCCAACAt                                                                 <  1:1836766/52‑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:1587592/1‑71 (MQ=255)
                                    cAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGt                               >  1:237533/1‑71 (MQ=255)
                                     aaaaCCGCCGCCCGGTGTACCAATCTGCCCTAACAt                                                                 >  1:634121/1‑36 (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:1156444/1‑71 (MQ=255)
                                                  ggTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTg                 <  1:1504740/71‑1 (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:3576895/44‑2 (MQ=12)
                                                               gCCCCAACATTGCTGCCAGCGTGACGATCATCCAGTGTTTTTgc                               <  1:3303665/44‑2 (MQ=12)
                                                                cccTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATAtt                     <  1:3882319/53‑1 (MQ=255)
                                                                 ccTAACATCGCCGCCAGCGTCACCAGCATCCAgtgt                                     <  1:1739990/36‑1 (MQ=255)
                                                                 ccTAACATCGCCGCCAGCGTCACCAGCATCCAgtgt                                     <  1:2709702/36‑1 (MQ=255)
                                                                 ccTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGt                               <  1:3317346/42‑1 (MQ=255)
                                                                 ccTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGacgc             >  1:2492214/1‑60 (MQ=255)
                                                                 ccTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCa  <  1:3514365/71‑1 (MQ=255)
                                                                 ccTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCa  <  1:843594/71‑1 (MQ=255)
                                                                   tAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTg                    <  1:1012511/51‑1 (MQ=255)
                                                                   |                                                                    
TCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCA  >  minE/1259732‑1259867

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: