Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,259,820 C→G 20.0% L368L (CTG→CTC) ‡ 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,8200CG20.0% 158.5 / 32.9 75L368L (CTG→CTC) ‡torZtrimethylamine N‑oxide reductase system III, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base C (28/32);  new base G (8/7);  total (36/39)
Fisher's exact test for biased strand distribution p-value = 7.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.83e-01

GATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGC  >  minE/1259756‑1259887
                                                                |                                                                   
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)
                  cgccgcCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCa                                             >  1:956166/1‑71 (MQ=255)
                  cgccgcCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCa                                             <  1:2071491/71‑1 (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:3416274/1‑65 (MQ=255)
                                      tGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTg                               >  1:3121549/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: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:2221078/1‑71 (MQ=255)
                                                  gccgccAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCAt             >  1:2216343/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:462589/71‑1 (MQ=255)
                                                       cAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCa        <  1:1507521/71‑1 (MQ=255)
                                                       cAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCa        <  1:1910719/71‑1 (MQ=255)
                                                       cAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCa        <  1:3841247/71‑1 (MQ=255)
                                                       cAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCa        >  1:3621450/1‑71 (MQ=255)
                                                       cAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCa        <  1:235307/71‑1 (MQ=255)
                                                        aGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCAt       >  1:792378/1‑71 (MQ=255)
                                                         gCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCATCAGCAtg      <  1:1881162/71‑2 (MQ=255)
                                                         gCGTGACGATCATCCAGTGTTTTTGCTCACCAAACTGTTGGCGCTGCATTCCCCAGCCTGCCATCAGCAtg      <  1:480872/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)
                                                                |                                                                   
GATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTGC  >  minE/1259756‑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: