Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,370,859 T→C 6.7% intergenic (+116/‑14) yohK → / → cdd predicted inner membrane protein/cytidine/deoxycytidine deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,370,8590TC 6.7% 293.0 / 9.4 104intergenic (+116/‑14)yohK/cddpredicted inner membrane protein/cytidine/deoxycytidine deaminase
Reads supporting (aligned to +/- strand):  ref base T (39/58);  new base C (4/3);  total (43/61)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.94e-01

Only 100 of 105 total aligned reads displayed.

TAATTAATGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCAATCT  >  minE/1370791‑1370926
                                                                    |                                                                   
taattaatGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTAt                                                                   <  1:2501191/71‑1 (MQ=255)
taattaatGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTAt                                                                   <  1:2746322/71‑1 (MQ=255)
taattaatGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTAt                                                                   >  1:143100/1‑71 (MQ=255)
 aattaatGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTAt                                                                   >  1:81105/1‑70 (MQ=255)
  attaatGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATg                                                                  <  1:334099/70‑1 (MQ=255)
      atGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGc                                                              <  1:3583566/70‑1 (MQ=255)
      atGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGc                                                              <  1:1317168/70‑1 (MQ=255)
       tGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGa                                                                 <  1:2995922/66‑1 (MQ=255)
         agaTTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACg                                                          >  1:3300635/1‑71 (MQ=255)
         agaTTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACg                                                          >  1:1426207/1‑71 (MQ=255)
          gaTTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGc                                                         <  1:317386/71‑1 (MQ=255)
            ttCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGcc                                                        <  1:1181131/70‑1 (MQ=255)
            ttCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCa                                                       >  1:3784143/1‑71 (MQ=255)
              cAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATg                                                     >  1:10785/1‑71 (MQ=255)
               aGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATg                                                     <  1:152029/70‑1 (MQ=255)
               aGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGc                                                    <  1:1189717/71‑1 (MQ=255)
                  tCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATc                                                 <  1:1816825/71‑1 (MQ=255)
                  tCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATc                                                 >  1:704160/1‑71 (MQ=255)
                  tCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATc                                                 >  1:1122918/1‑71 (MQ=255)
                   cacaTATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATc                                                 >  1:215249/1‑70 (MQ=255)
                      atatAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAgg                                                               <  1:26645/53‑1 (MQ=255)
                      atatAAAGCCAAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGt                                            <  1:2095069/71‑1 (MQ=255)
                       tataAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAgg                                                               >  1:593599/1‑52 (MQ=255)
                       tataAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCGCGt                                            >  1:1953288/1‑71 (MQ=255)
                        ataAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACg                                                          >  1:2520542/1‑56 (MQ=255)
                        ataAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGcc                                                        <  1:87491/58‑1 (MQ=255)
                        ataAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGtt                                           >  1:1943721/1‑71 (MQ=255)
                          aaaGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAgg                                                               <  1:1594945/49‑1 (MQ=255)
                          aaaGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAgg                                                               <  1:2330119/49‑1 (MQ=255)
                            aGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTc                                        >  1:1618794/1‑70 (MQ=255)
                            aGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCa                                       <  1:3279408/71‑1 (MQ=255)
                            aGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCa                                       <  1:3296937/71‑1 (MQ=255)
                            aGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCa                                       <  1:3605270/71‑1 (MQ=255)
                                acaACGGGTTCGTAAACTGTTATCCCATTACATGATTATGa                                                                 >  1:303457/1‑41 (MQ=255)
                                acaACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCAt                                                      <  1:1613061/52‑1 (MQ=255)
                                acaACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGTATCCACGTTTTCAAAcc                                   <  1:578236/71‑1 (MQ=255)
                                acaACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAAcc                                   <  1:1988627/71‑1 (MQ=255)
                                 caACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGc                                                         <  1:2803365/48‑1 (MQ=255)
                                 caACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGc                                                         <  1:689547/48‑1 (MQ=255)
                                  aaCGGGTTCGTAAACTGTTATCCCATTACATGATTa                                                                    >  1:3512963/1‑36 (MQ=255)
                                  aaCGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCaa                                                            <  1:3350369/44‑1 (MQ=255)
                                  aaCGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGc                                 <  1:1303317/71‑1 (MQ=255)
                                  aaCGGGTTCGTAAACTGTTATCCCATTACATGATCATGAGGCAACGCCATGCATCCACGTTTTCa                                       >  1:2229880/1‑65 (MQ=255)
                                  aaCGGGTTCGTAAACTGTTATCCCATTACATGATCATGAGGCAACGCCATGCATCCACGTTTTCa                                       >  1:3649683/1‑65 (MQ=255)
                                   aCGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCAt                                                      >  1:3748363/1‑49 (MQ=255)
                                   aCGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGc                                 <  1:1459726/70‑1 (MQ=255)
                                    cGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCtt                               >  1:3918730/1‑71 (MQ=255)
                                    cGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCtt                               <  1:2132094/71‑1 (MQ=255)
                                     gggTTCGTAAACTGTTATCCCATTACATGATTATTAGGCAACGCCATGCATCCACGTTTTCAAACCGCttt                              <  1:2006276/71‑1 (MQ=255)
                                     gggTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGc                                                         <  1:68363/44‑1 (MQ=255)
                                     gggTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCtt                               <  1:803736/70‑1 (MQ=255)
                                      ggTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCtttt                             <  1:264524/71‑1 (MQ=255)
                                      ggTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCtttt                             <  1:3735249/71‑1 (MQ=255)
                                       gTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTg                            <  1:972134/71‑1 (MQ=255)
                                         tCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCaaa                                     >  1:1606603/1‑60 (MQ=255)
                                          cGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGtttt                                         >  1:3108275/1‑55 (MQ=255)
                                          cGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGcc                          >  1:3430834/1‑70 (MQ=255)
                                          cGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGcc                          >  1:3105297/1‑70 (MQ=255)
                                            tAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGc                                                    >  1:3106378/1‑42 (MQ=255)
                                            tAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGc                                                    <  1:81305/42‑1 (MQ=255)
                                            tAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCaa                       >  1:770377/1‑71 (MQ=255)
                                            tAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCaa                       <  1:1231122/71‑1 (MQ=255)
                                              aaCTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGttt                                          <  1:483338/50‑1 (MQ=255)
                                              aaCTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACt                     >  1:2788784/1‑71 (MQ=255)
                                               aCTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACtt                    <  1:811517/71‑1 (MQ=255)
                                               aCTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACtt                    <  1:2660740/71‑1 (MQ=255)
                                               aCTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACtt                    <  1:2080647/71‑1 (MQ=255)
                                               aCTGTTATCCCATTACATGATCATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACtt                    <  1:2413661/71‑1 (MQ=255)
                                                cTGTTATCCCATTACATGATCATGAGGCAACGCCATGCATCCa                                               >  1:3127054/1‑43 (MQ=255)
                                                cTGTTATCCCATAACATGATCATGAGGCAACGCCATGCATCCACGATTTCAAACCGCTATTGCCCAACTTg                   >  1:3838197/1‑71 (MQ=255)
                                                 tGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGtt                                           <  1:1053505/46‑1 (MQ=255)
                                                 tGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGtt                                           <  1:495659/46‑1 (MQ=255)
                                                  gTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGATTTTGc                           <  1:512942/61‑1 (MQ=255)
                                                   ttATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGc                  <  1:3319171/69‑1 (MQ=255)
                                                   ttATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCg                 <  1:1623341/70‑1 (MQ=255)
                                                   ttATCCCATTACATGATTATGAGGCAACCCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCgg                <  1:3904087/71‑1 (MQ=255)
                                                      tCCCATTACATGATTATGAGGCAACGCCATGCATCCACGttt                                          <  1:2022629/42‑1 (MQ=255)
                                                      tCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGAt              <  1:1459177/70‑1 (MQ=255)
                                                        ccATTACATGATTATGAGGCAACGCCATGCATCCAc                                              >  1:3952541/1‑36 (MQ=255)
                                                        ccATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACtt                    <  1:1124580/62‑1 (MQ=255)
                                                        ccATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAAc           <  1:3773826/71‑1 (MQ=255)
                                                        ccATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAAc           >  1:2760523/1‑71 (MQ=255)
                                                         cATTACATGATTATGAGGCAACGTCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACt          >  1:31761/1‑71 (MQ=255)
                                                         cATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACt          >  1:3509086/1‑71 (MQ=255)
                                                         cATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACt          >  1:2956425/1‑71 (MQ=255)
                                                         cATTACATGATCATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACt          <  1:711767/71‑1 (MQ=255)
                                                           ttACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTg        >  1:33375/1‑71 (MQ=255)
                                                           ttACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTg        <  1:1492510/71‑1 (MQ=255)
                                                           ttACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTg        >  1:2080999/1‑71 (MQ=255)
                                                             aCATGATCATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGc       <  1:1822022/70‑1 (MQ=255)
                                                              cATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCaa     <  1:692356/71‑1 (MQ=255)
                                                               aTGATTATGAGGCAACGCCATGCATCCACGTTTTCa                                       <  1:668470/36‑1 (MQ=255)
                                                               aTGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCAAt    >  1:2290867/1‑71 (MQ=255)
                                                               aTGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCAAt    <  1:1642469/71‑1 (MQ=255)
                                                               aTGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCAAt    <  1:1915914/71‑1 (MQ=255)
                                                                tGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCAAt    <  1:2442651/70‑1 (MQ=255)
                                                                tGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCAATc   >  1:1136871/1‑71 (MQ=255)
                                                                tGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCAATc   <  1:3467809/71‑1 (MQ=255)
                                                                 gATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCAATCt  >  1:2965457/1‑71 (MQ=255)
                                                                   ttATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTg        >  1:1784000/1‑63 (MQ=255)
                                                                    |                                                                   
TAATTAATGAGATTCAGATCACATATAAAGCCACAACGGGTTCGTAAACTGTTATCCCATTACATGATTATGAGGCAACGCCATGCATCCACGTTTTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCAATCT  >  minE/1370791‑1370926

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