Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 982,587 T→C 11.0% S11S (TCT→TCC ydfA → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE982,5870TC11.0% 335.4 / 24.7 127S11S (TCT→TCCydfAhypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (56/57);  new base C (8/6);  total (64/63)
Fisher's exact test for biased strand distribution p-value = 7.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 142 total aligned reads displayed.

TCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAA  >  minE/982519‑982654
                                                                    |                                                                   
tCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATctct                                                                   >  1:2057105/1‑71 (MQ=255)
tCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATctct                                                                   >  1:4062144/1‑71 (MQ=255)
 ccGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTg                                                                  >  1:1394761/1‑71 (MQ=255)
 ccGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTg                                                                  >  1:417821/1‑71 (MQ=255)
    cTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGgtg                                                               <  1:1415266/71‑1 (MQ=255)
         ccAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACgg                                                          >  1:4572432/1‑71 (MQ=255)
          cAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGc                                                         >  1:4226460/1‑71 (MQ=255)
           aGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCg                                                        >  1:917044/1‑71 (MQ=255)
            gCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCgt                                                       >  1:2822184/1‑71 (MQ=255)
            gCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCgt                                                       <  1:878853/71‑1 (MQ=255)
             cGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAAtctc                                                                    <  1:5628697/57‑1 (MQ=255)
                 aaGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGgt                                                                >  1:3648786/1‑57 (MQ=255)
                 aaGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGgt                                                                >  1:2688634/1‑57 (MQ=255)
                 aaGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCgt                                                       <  1:143826/66‑1 (MQ=255)
                 aaGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCgt                                                       <  1:5353364/66‑1 (MQ=255)
                 aaGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCgt                                                       <  1:3120689/66‑1 (MQ=255)
                    gggAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCaa                                               <  1:4421474/71‑1 (MQ=255)
                         aTGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGgtg                                                               <  1:1262847/50‑1 (MQ=255)
                         aTGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCgt                                                       <  1:993457/58‑1 (MQ=255)
                         aTGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCgt                                                       <  1:1151288/58‑1 (MQ=255)
                         aTGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCgt                                                       <  1:2769632/58‑1 (MQ=255)
                         aTGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAg                                          <  1:272388/71‑1 (MQ=255)
                          tGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAgg                                         <  1:4847719/71‑1 (MQ=255)
                          tGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAg                                          >  1:5218018/1‑70 (MQ=255)
                           gAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGa                                        >  1:3614632/1‑71 (MQ=255)
                           gAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGa                                        >  1:3253683/1‑71 (MQ=255)
                           gAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGAGTGTTTCCAAACCAGGa                                        >  1:3639645/1‑71 (MQ=255)
                               tCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCgtgt                                                     <  1:2040631/54‑1 (MQ=255)
                                 aGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCAc                                  >  1:4293902/1‑71 (MQ=255)
                                 aGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCAc                                  >  1:4859633/1‑71 (MQ=255)
                                  gCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACg                                 <  1:5078245/71‑1 (MQ=255)
                                  gCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACg                                 <  1:1306891/71‑1 (MQ=255)
                                  gCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTAAGGCGTGTTTCCAAACCAGGACGGCACg                                 <  1:3760797/71‑1 (MQ=255)
                                   cATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGt                                >  1:1875911/1‑71 (MQ=255)
                                    aTGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGtt                               <  1:1904906/71‑1 (MQ=255)
                                     tGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTc                              <  1:347761/71‑1 (MQ=255)
                                       gATACTATCGATCTTGGCAACAACGAATCTCTGGTGTa                                                             >  1:2077663/1‑38 (MQ=255)
                                       gATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCAc                            <  1:2359212/71‑1 (MQ=255)
                                       gATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCAc                            >  1:5263158/1‑71 (MQ=255)
                                        aTACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGAc                                       >  1:217861/1‑59 (MQ=255)
                                        aTACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGAc                                       >  1:3013321/1‑59 (MQ=255)
                                        aTACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCAcc                           >  1:93754/1‑71 (MQ=255)
                                         tACTATCGATCTTGGCAACAACGAATCTCTGGTGTACTGCGTGTTTCCAAACCAGGACGGCACGTTCACCg                          <  1:1293860/71‑1 (MQ=255)
                                         tACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCg                          <  1:1743708/71‑1 (MQ=255)
                                           cTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTc                              >  1:51302/1‑65 (MQ=255)
                                           cTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGc                         <  1:5578374/70‑1 (MQ=255)
                                           cTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGc                         <  1:165066/70‑1 (MQ=255)
                                           cTATCGATCTTGGCAACAACGAATCTCTGGTGTACAGCGTGTTTCCAAACCAGGACGGCACGTTCACCGc                         <  1:4647135/70‑1 (MQ=255)
                                            tATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACgg                                     <  1:115781/57‑1 (MQ=255)
                                            tATCGATCTTGGCAACAACGAATCTCTGGTGTACGGAGTGTTTc                                                  >  1:161236/1‑44 (MQ=255)
                                            aaTCGACCTCGGCAACAACGAATCCCTGGTGTGCg                                                           >  1:3371419/2‑35 (MQ=16)
                                            aaTCGACCTCGGCAACAACGAATCCCTGGTGTGCGGCGTGTTCCCCAACCAGGATGTAACGTTCACTGcc                        >  1:2488450/2‑69 (MQ=255)
                                             aTCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAAt                      >  1:4801204/1‑71 (MQ=255)
                                              tCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGtt                               >  1:103855/1‑61 (MQ=255)
                                              tCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAAt                      >  1:2985504/1‑70 (MQ=255)
                                              tCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAAt                      >  1:2598837/1‑70 (MQ=255)
                                               cGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATg                     >  1:2933886/1‑70 (MQ=255)
                                               cGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGa                    >  1:1089754/1‑71 (MQ=255)
                                               cGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGa                    >  1:2444799/1‑71 (MQ=255)
                                                gATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGAc                   >  1:1627117/1‑71 (MQ=255)
                                                   ctcggCAACAACGAATCCCTGGTGTGCGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCATGACg                  <  1:2015484/66‑1 (MQ=255)
                                                   cTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGta                >  1:3541868/1‑71 (MQ=255)
                                                   cTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGta                >  1:1714743/1‑71 (MQ=255)
                                                   cTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGta                >  1:2658389/1‑71 (MQ=255)
                                                     cggCAACAACGAATCCCTGGTGTGCGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCATGa                    >  1:1535482/2‑65 (MQ=255)
                                                      tgCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGa                                        >  1:2325243/2‑44 (MQ=255)
                                                      ggCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCAGGTTCACCGCAATGACGTATAc             >  1:5165757/1‑71 (MQ=255)
                                                      ggCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACg                  <  1:1106693/66‑1 (MQ=255)
                                                      ggCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACg                  <  1:5432214/66‑1 (MQ=255)
                                                       gCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGta                <  1:1698663/67‑1 (MQ=255)
                                                       gCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGta                <  1:3138491/67‑1 (MQ=255)
                                                       gCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATAcc            <  1:5574020/71‑1 (MQ=255)
                                                       gCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATAcc            <  1:4958940/71‑1 (MQ=255)
                                                        caacaaCGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGc                         <  1:342520/57‑1 (MQ=255)
                                                        caacaaCGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATAcc            >  1:4368369/1‑70 (MQ=255)
                                                        caacaaCGAATCCCTGGTGTGCGGCGTGTTCCCCAACCAGGATGGAACGTTc                              >  1:1368259/1‑52 (MQ=255)
                                                         aacaacGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCaa          <  1:4636744/71‑1 (MQ=255)
                                                         aacaacGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCaa          <  1:2344521/71‑1 (MQ=255)
                                                         aacaacGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCaa          <  1:5425847/71‑1 (MQ=255)
                                                         aacaacGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCaa          <  1:587474/71‑1 (MQ=255)
                                                         aacaacGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCaa          <  1:327097/71‑1 (MQ=255)
                                                         aacaacGAATCCCTGGTGTGCGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCATGACGTATACCaa          <  1:5612030/71‑1 (MQ=255)
                                                          acaacGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCaaa         <  1:3543836/71‑1 (MQ=255)
                                                           caacGAATCTGTGGTGTACGGCGTGTTTCCAAAca                                            <  1:3063435/35‑2 (MQ=255)
                                                           caacGAATCTCTGGTGTACGGCGTGTTTCCAAAcc                                            >  1:2182536/1‑35 (MQ=255)
                                                            aacGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCa                        <  1:2332322/54‑1 (MQ=255)
                                                               gAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCaa    <  1:348569/71‑1 (MQ=255)
                                                               gAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCaa    <  1:411584/71‑1 (MQ=255)
                                                               gAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCaa    <  1:3691746/71‑1 (MQ=255)
                                                                aaTCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCaa    >  1:5408693/1‑70 (MQ=255)
                                                                 aTCCCTGGTGTGCGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCaaaa  <  1:4875973/71‑1 (MQ=255)
                                                                  tctcTGGTGTACGGCGTGTTTCCAAACCAGGACGGc                                    <  1:672586/36‑1 (MQ=255)
                                                                  tctcTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGc      <  1:4505428/66‑1 (MQ=255)
                                                                  tctcTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGc      <  1:4885688/66‑1 (MQ=255)
                                                                   ctctGGTGTACGTCGTGTTTCCAAACCAGGACGGCACGTTCACCGCaa                       >  1:693582/1‑48 (MQ=255)
                                                                   ctctGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCaaa   >  1:1740954/1‑68 (MQ=255)
                                                                   ctctGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCaaa   >  1:2412398/1‑68 (MQ=255)
                                                                   ctctGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCaaa   >  1:2656630/1‑68 (MQ=255)
                                                                    tctGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGtt                               >  1:4400558/1‑39 (MQ=255)
                                                                    tctGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGtata              >  1:1999018/1‑56 (MQ=255)
                                                                    |                                                                   
TCCGCTACCCCAGCGTTAAGGGGAAATGAGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAA  >  minE/982519‑982654

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: