Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,569,817 T→C 8.9% P29P (CCT→CCC yibI → predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,569,8170TC 8.9% 251.8 / 12.5 90P29P (CCT→CCCyibIpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (48/34);  new base C (5/3);  total (53/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

CGTTGGGAGTGCTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCT  >  minE/2569750‑2569885
                                                                   |                                                                    
cgTTGGGAGTGCTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTAt                                                                   >  1:960427/1‑71 (MQ=255)
  ttGGGAGTGCTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCtt                                                                     <  1:101214/67‑1 (MQ=255)
   tGGGAGTGCTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCt                                                                 >  1:894829/1‑70 (MQ=255)
   tGGGAGTGCTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCt                                                                 >  1:587957/1‑70 (MQ=255)
   tGGGAGTGCTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCt                                                                 >  1:3114019/1‑70 (MQ=255)
   tGGGAGTGCTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCt                                                                 >  1:1007652/1‑70 (MQ=255)
      gAGTGCTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGAt                                                              >  1:467429/1‑70 (MQ=255)
       aGTGCTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTg                                                            <  1:3123798/71‑1 (MQ=255)
        gTGCTTATTTTTGTATTTCTGGTGATCTTCTATGGGATCATTATTTTACATGATATTCCCTACCTGATTg                                                            >  1:1295997/1‑70 (MQ=255)
        gTGCTTATTTTTGTATTTCTGGTGATCTTCTATGGGATCATTATTTTACATGATATTCCCTACCTGATTg                                                            >  1:1642966/1‑70 (MQ=255)
        gTGCTTATTTTTGTATTTCTGGTGATCTTCTATGGGATCATTATTTTACATGATATTCCCTACCTGATTGc                                                           >  1:3934834/1‑71 (MQ=255)
          gCTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGa                                                               <  1:1467725/65‑1 (MQ=255)
           cTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCa                                                         >  1:462838/1‑70 (MQ=255)
           cTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCa                                                         >  1:1497526/1‑70 (MQ=255)
            tGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGAtt                                                             >  1:3436943/1‑65 (MQ=255)
            tGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGAt                                                              <  1:3059743/64‑1 (MQ=255)
              aTCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCaaa                                                       <  1:2674449/69‑1 (MQ=255)
              aTCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCaaa                                                       <  1:3897787/69‑1 (MQ=255)
                cTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCaaaaa                                                     >  1:2114060/1‑69 (MQ=255)
                cTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAgc                                                   >  1:2045046/1‑71 (MQ=255)
                cTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAg                                                    >  1:3757642/1‑70 (MQ=255)
                  tCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGAt                                                              <  1:3266974/58‑1 (MQ=255)
                   cGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCt                                                                 >  1:1502056/1‑54 (MQ=255)
                   cGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCaaaaa                                                     <  1:995146/66‑1 (MQ=255)
                   cGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCaaaaa                                                     <  1:3618881/66‑1 (MQ=255)
                   cGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCaaaaa                                                     <  1:2397959/66‑1 (MQ=255)
                   cGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCaaaaa                                                     <  1:3864142/66‑1 (MQ=255)
                   cGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCaaaaa                                                     <  1:3868828/66‑1 (MQ=255)
                   cGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAgcgc                                                 <  1:3941678/70‑1 (MQ=255)
                   cGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAg                                                    >  1:3363620/1‑67 (MQ=255)
                    gTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCaa                                               <  1:478405/71‑1 (MQ=255)
                     tCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCt                                                                 <  1:2798126/52‑1 (MQ=255)
                     tCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCa                                                >  1:2025948/1‑69 (MQ=255)
                       tttCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCa                                            >  1:798974/1‑71 (MQ=255)
                        ttCTGGTGATCTTCTATGGGATCATTATTTTACATGATATTCCCTACCTGATTGCCAAAAAACGTaatc                                             >  1:2529465/1‑67 (MQ=255)
                        ttCTGGTGATCTTCTATGGGATCATTATTTTACATGATATTCCCTACCTGATTGCCAAAAAACGTaat                                              >  1:2307023/1‑67 (MQ=255)
                         tCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAAcc                                             <  1:1955310/68‑1 (MQ=255)
                         tCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCAt                                           <  1:995434/70‑1 (MQ=255)
                         tCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATc                                          >  1:2266703/1‑71 (MQ=255)
                         tCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATc                                          >  1:1074034/1‑71 (MQ=255)
                         tCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGACAAAAAgcg                                                  >  1:3170979/1‑63 (MQ=255)
                          cTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATcc                                         >  1:2047506/1‑71 (MQ=255)
                           tGGTGATCTTCTATGGGATCATTATTTTACATGATATTCCCTACCTGATTGCCaaaaa                                                     <  1:696601/58‑1 (MQ=255)
                            ggTGATCTTCTATGGGATCATTATTTTACATGATATTCCCTACCTGATTGCCAAAAAacgta                                                <  1:1893992/62‑3 (MQ=255)
                             gTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCaaaa                                                      <  1:1866585/55‑1 (MQ=255)
                             gTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCaa                                               <  1:748539/62‑1 (MQ=255)
                             gTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATccc                                        <  1:1767787/69‑1 (MQ=255)
                             gTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCa                                      <  1:2789212/71‑1 (MQ=255)
                             gTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCa                                      >  1:2533292/1‑71 (MQ=255)
                             gTGATCTTCTATGGGATCATTATTTTACATGATATTCCCTACCTGATTGCCAAAAAACGTAATCATcct                                        <  1:3379106/69‑2 (MQ=255)
                              tGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCaa                                                        <  1:3122978/52‑1 (MQ=255)
                                aTTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGAt                                                              >  1:1876429/1‑44 (MQ=255)
                                aTTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGAt                                                              >  1:3765969/1‑44 (MQ=255)
                                 ttttttATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGcc                                  >  1:2100605/1‑71 (MQ=255)
                                  tttttATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGcc                                  >  1:2744494/1‑70 (MQ=255)
                                  tttttATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGcc                                  >  1:3350040/1‑70 (MQ=255)
                                   ttttATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCg                                 <  1:249353/70‑1 (MQ=255)
                                    tttATGGAATCATCGCGATACATGACATCCCTTATCTg                                                                >  1:572273/1‑38 (MQ=255)
                                    tttATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGAc                               >  1:571351/1‑71 (MQ=255)
                                         ggAATCATCGCGATACATGACATCCCTTATCTGATTGCCaaaa                                                      <  1:3220887/43‑1 (MQ=255)
                                            atcatcGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTc                        <  1:1982544/70‑1 (MQ=255)
                                             tcatcGCGATACATGACATCCCTTATCTGATTGCCaaaaa                                                     >  1:2465508/1‑40 (MQ=255)
                                             tcatcGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCa                       >  1:75986/1‑70 (MQ=255)
                                             tcatcGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCa                       >  1:57758/1‑70 (MQ=255)
                                             tcatcGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCa                       >  1:1920899/1‑70 (MQ=255)
                                                  gcgATACATGACATCCCTTATCTGATTGCCAAAAAg                                                    >  1:1934114/1‑36 (MQ=255)
                                                  gcgATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATAcggc                 >  1:3639745/1‑71 (MQ=255)
                                                   cgATACATGACATCCCTTATCTGATTGCCAAAAAGCGCACCCATCCCCATGCCGACGCTATTCATAcggag                >  1:916638/1‑69 (MQ=255)
                                                   cgATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATAcggcg                <  1:1953051/71‑1 (MQ=255)
                                                     aTACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTAtt                         <  1:3783623/60‑1 (MQ=255)
                                                      tACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATc                                          <  1:1830793/42‑1 (MQ=255)
                                                        cATGACATCCCTTATCTGATTGCCAAAAAGCGCAAc                                              >  1:3943071/1‑36 (MQ=255)
                                                        cATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATa                     >  1:990938/1‑61 (MQ=255)
                                                        cATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATa                     >  1:1921031/1‑61 (MQ=255)
                                                        cATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCt            <  1:124014/70‑1 (MQ=255)
                                                          tGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACg                              >  1:1094361/1‑50 (MQ=255)
                                                           gACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGt        >  1:1269066/1‑71 (MQ=255)
                                                           gACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGt        >  1:2138577/1‑71 (MQ=255)
                                                             cttCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTa                           >  1:2725533/3‑50 (MQ=255)
                                                             cATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGcc                                  <  1:276661/43‑1 (MQ=255)
                                                             cATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGt        >  1:3620947/1‑69 (MQ=255)
                                                             cATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGa      <  1:3176992/71‑1 (MQ=255)
                                                             cATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGa      <  1:1503116/71‑1 (MQ=255)
                                                              aTCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCa                       <  1:663654/53‑1 (MQ=255)
                                                              aTCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGAg     >  1:2305542/1‑71 (MQ=255)
                                                               tCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGAGc    <  1:149024/71‑1 (MQ=255)
                                                                cccTTATCTGATTGCCAAAAAGCGCAACCATCCCCAt                                     >  1:1496976/1‑37 (MQ=255)
                                                                cccTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGAGcc   >  1:3491670/1‑71 (MQ=255)
                                                                cccTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGAGcc   >  1:842159/1‑71 (MQ=255)
                                                                cccTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGAGcc   >  1:843393/1‑71 (MQ=255)
                                                                 ccTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCt  <  1:1558395/71‑1 (MQ=255)
                                                                   ttATCTGATTGCCAAAAAGCGCAACCATCCCCATg                                    >  1:990655/1‑35 (MQ=255)
                                                                   |                                                                    
CGTTGGGAGTGCTGATCTTCGTCTTTCTGGTGATTTTTTATGGAATCATCGCGATACATGACATCCCTTATCTGATTGCCAAAAAGCGCAACCATCCCCATGCCGACGCTATTCATACGGCGGGCTGGGTGAGCCT  >  minE/2569750‑2569885

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: