Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,940,923 C→A 11.8% G180C (GGC→TGC)  fecC ← iron‑dicitrate transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,940,9230CA11.8% 308.9 / 16.6 101G180C (GGC→TGC) fecCiron‑dicitrate transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (44/45);  new base A (4/8);  total (48/53)
Fisher's exact test for biased strand distribution p-value = 3.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01

Only 100 of 102 total aligned reads displayed.

AAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGAT  >  minE/2940859‑2940992
                                                                |                                                                     
aaGAGCTGCCAGACATCCTGCCAGCGGGTGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGt                                                                     <  1:1680369/67‑1 (MQ=255)
   agCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCa                                                                >  1:3928100/1‑69 (MQ=255)
   agCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCAt                                                               >  1:2804368/1‑70 (MQ=255)
   agCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATg                                                              <  1:3311750/71‑1 (MQ=255)
    gCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGa                                                             >  1:5677543/1‑71 (MQ=255)
      tGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTa                                                                    <  1:866669/62‑1 (MQ=255)
        ccAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCtt                                                         >  1:4696969/1‑71 (MQ=255)
        ccAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCtt                                                         <  1:5536495/71‑1 (MQ=255)
        ccAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCt                                                          <  1:851021/70‑1 (MQ=255)
          aGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCg                                                       >  1:2440963/1‑71 (MQ=255)
             cATCCTGCCAGCGGGCGTTGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGcc                                                    <  1:728398/71‑1 (MQ=255)
              aTCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATc                                                           >  1:4623525/1‑63 (MQ=255)
                   gCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCcagcag                                               >  1:4538324/1‑70 (MQ=255)
                    ccAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCcagcagca                                             <  1:899833/71‑1 (MQ=255)
                     cAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCcagcagcag                                            <  1:4445747/71‑1 (MQ=255)
                     cAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCcagcagcag                                            >  1:2600287/1‑71 (MQ=255)
                     cAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCcagcagcag                                            >  1:15255/1‑71 (MQ=255)
                       gCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCgg                                                      <  1:1869783/59‑1 (MQ=255)
                       gCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAggg                                          <  1:4979176/71‑1 (MQ=255)
                        cGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGt                                         <  1:2335046/71‑1 (MQ=255)
                         gggCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGt                                         >  1:1040648/1‑70 (MQ=255)
                          taccggttggtttagcgagaagAGCCAGTAAAAGACGCAGt                                                                     <  1:4281735‑M2/19‑1 (MQ=255)
                          taccggttgggttagcgagaagAGCCAGTAAAAGACGCAGt                                                                     <  1:315000‑M2/19‑1 (MQ=255)
                           accggttgggttagcgagaagAGCCAGTAAAAGACGCAGt                                                                     >  1:221378‑M2/22‑40 (MQ=255)
                           accggttgggttagcgagaagAGCCAGTAAAAGACGCAGt                                                                     >  1:4638053‑M2/22‑40 (MQ=255)
                           accggttgggttagcgagaagAGCCAGTAAAAGACGCAGt                                                                     >  1:311424‑M2/22‑40 (MQ=255)
                           accggttgggttagcgagaagAGCCAGTAAAAGACGCAGt                                                                     >  1:5203469‑M2/22‑40 (MQ=255)
                            cGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATg                                     <  1:1323210/71‑1 (MQ=255)
                             gTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGc                                    >  1:1657209/1‑71 (MQ=255)
                             gTGGCACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGc                                    <  1:4619452/71‑1 (MQ=255)
                             cggttgggttagcgagaagAGCCAGTAAAAGACGCAGt                                                                     <  1:3476691‑M2/19‑1 (MQ=255)
                             cggttgggttagcgagaagAGCCAGTAAAAGACGCAGt                                                                     <  1:296584‑M2/19‑1 (MQ=255)
                               gttgggttagcgagaagAGCCAGTAAAAGACGCAGt                                                                     <  1:5605238‑M2/19‑1 (MQ=255)
                               gttgggttagcgagaagAGCCAGTAAAAGACGCAGt                                                                     <  1:4304800‑M2/19‑1 (MQ=255)
                               gggACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCtt                                                         >  1:3187980/1‑48 (MQ=255)
                               gggACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCgg                                  >  1:1920005/1‑71 (MQ=255)
                                ttgggttagcgagaagAGCCAGTAAAAGACGCAGt                                                                     <  1:5464082‑M2/19‑1 (MQ=255)
                                ttgggttagcgagaagAGCCAGTAAAAGACGCAGt                                                                     <  1:3815422‑M2/19‑1 (MQ=255)
                                ggACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGa                                                             >  1:396076/1‑43 (MQ=255)
                                 gACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGt                                <  1:5315348/71‑1 (MQ=255)
                                 gACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGt                                <  1:5000159/71‑1 (MQ=255)
                                  acacTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCcagcagc                                              >  1:2260745/1‑56 (MQ=255)
                                  acacTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCcagcagc                                              >  1:2478118/1‑56 (MQ=255)
                                  acacTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGt                                >  1:2037066/1‑70 (MQ=255)
                                  acacTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTa                               <  1:1299109/71‑1 (MQ=255)
                                  acacTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTa                               <  1:3417407/71‑1 (MQ=255)
                                    acTCCGCCTGCGAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGt                                <  1:2841071/68‑1 (MQ=255)
                                    acTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAg                             >  1:1212340/1‑71 (MQ=255)
                                     cTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGt                                <  1:844676/67‑1 (MQ=255)
                                       ccGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTc                                                        >  1:2409728/1‑41 (MQ=255)
                                       ccGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGcc                          >  1:4873813/1‑71 (MQ=255)
                                       ccGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGcc                          >  1:806127/1‑71 (MQ=255)
                                        cGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGAt                                                            >  1:2516497/1‑36 (MQ=255)
                                        cGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGAt                                                            >  1:155178/1‑36 (MQ=255)
                                        cGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGcc                          >  1:4559467/1‑70 (MQ=255)
                                        cGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGcc                          >  1:4806669/1‑70 (MQ=255)
                                        cGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGcc                          >  1:3849812/1‑70 (MQ=255)
                                        cGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGcc                          <  1:3160664/70‑1 (MQ=255)
                                            tGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATa                      >  1:1464831/1‑70 (MQ=255)
                                            tGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATa                      >  1:3405125/1‑70 (MQ=255)
                                             gccagccaGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATAc                     >  1:979894/1‑70 (MQ=255)
                                              ccagccagTAAAAGATGCCGTAAGCATGATCTTCGGCcagcag                                               <  1:965132/43‑1 (MQ=255)
                                              ccagccagTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAggg                                          <  1:3761631/48‑1 (MQ=255)
                                              ccagccagTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAggg                                          <  1:4773106/48‑1 (MQ=255)
                                              ccagccagTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAggg                                          <  1:3583650/48‑1 (MQ=255)
                                              ccagccagTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACaa                   <  1:2575713/71‑1 (MQ=255)
                                              ccagccagTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACaa                   <  1:63706/71‑1 (MQ=255)
                                              ccagccagTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACa                    <  1:3190775/70‑1 (MQ=255)
                                                agccagTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGa                                       >  1:4863534/1‑49 (MQ=255)
                                                  ccagTAAAAGATGCCGTAAGCATGATCTTCGGCcagcagc                                              >  1:3908302/1‑40 (MQ=255)
                                                    agTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATa                      >  1:3454589/1‑62 (MQ=255)
                                                     gTAAAAGATGCCGTAAGCATGATCTTCGGCcagcagca                                             <  1:5151891/38‑1 (MQ=255)
                                                     gTAAAAGATGCCGTAAGCATGATCTTCGGCcagcagca                                             <  1:5023045/38‑1 (MQ=255)
                                                     gTAAAAGATGCCGTAAGCATGATCTTCGGCcagcagca                                             <  1:2304844/38‑1 (MQ=255)
                                                     gTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGt                                         <  1:3997697/42‑1 (MQ=255)
                                                     gTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGcc             >  1:1506570/1‑70 (MQ=255)
                                                      tAAAAGATGCCGTAAGCATGATCTTCGGCcagcagcag                                            <  1:1887268/38‑1 (MQ=255)
                                                      tAAAAGATGCCGTAAGCATGATCTTCGGCcagcag                                               <  1:2654600/35‑1 (MQ=255)
                                                      tAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGt                                         <  1:2759335/41‑1 (MQ=255)
                                                      tAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGt                                         <  1:2891433/41‑1 (MQ=255)
                                                      tAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGt                                         <  1:3398288/41‑1 (MQ=255)
                                                      tAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACa                    <  1:5138802/62‑1 (MQ=255)
                                                       aaaaGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGt                                         <  1:815855/40‑1 (MQ=255)
                                                       aaaaGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGa           >  1:184489/1‑70 (MQ=255)
                                                       aaaaGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGa           <  1:3715774/70‑1 (MQ=255)
                                                        aaaGATGCCGTAAGCATGATCTTCGGCcagcagca                                             >  1:4631440/1‑35 (MQ=255)
                                                        aaaGATGCCGTAAGCATGATCTTCGGCcagcagca                                             >  1:51199/1‑35 (MQ=255)
                                                        aaaGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGa                                       >  1:110389/1‑41 (MQ=255)
                                                        aaaGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGa           >  1:908964/1‑69 (MQ=255)
                                                        aaaGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGa           >  1:1571628/1‑69 (MQ=255)
                                                        aaaGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGa           >  1:4838999/1‑69 (MQ=255)
                                                        aaaGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGa           >  1:4021089/1‑69 (MQ=255)
                                                          aGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAgc       <  1:2018120/71‑1 (MQ=255)
                                                            aTGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGCGTAAGGcc                          <  1:1874828/50‑1 (MQ=255)
                                                             tGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGc                                    <  1:167498/39‑1 (MQ=255)
                                                             tGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGc                                    <  1:4846957/39‑1 (MQ=255)
                                                             tGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGcgcg    <  1:5632986/71‑1 (MQ=255)
                                                              gCCGTAAGCATGATCTTCGGCCGGCAGCAGGGTGATGCGGGTaa                              >  1:5216156/1‑44 (MQ=255)
                                                              gCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTaa                              >  1:751354/1‑44 (MQ=255)
                                                               ccGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGAt  <  1:3962875/71‑1 (MQ=255)
                                                                |                                                                     
AAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGAT  >  minE/2940859‑2940992

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: