Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,396,394 A→G 11.7% V265A (GTA→GCA) ‡ corA ← magnesium/nickel/cobalt transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,396,3940AG11.7% 181.2 / 12.9 68V265A (GTA→GCA) ‡corAmagnesium/nickel/cobalt transporter
Reads supporting (aligned to +/- strand):  ref base A (42/18);  new base G (3/5);  total (45/23)
Fisher's exact test for biased strand distribution p-value = 1.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCACTTCAGTTCTGGCATAAACTCAAAGTTCATGCCATAGCTGGAAGCAACGAGCGTCGGCGGCAGGAATACCACGGATACCACCGAGAAGATTTTGATGATGCGGTTCTGCTCGATGTTGATAAAACCCATTGCCGCCT  >  minE/2396324‑2396463
                                                                      |                                                                     
ccACTTCAGTTCTGGCATAAACTCAAAGTTCATGCCATAGCTGGAAGCAACGAGCGTCGGCGGCAGGAATa                                                                       >  1:2384133/1‑71 (MQ=255)
ccACTTCAGTTCTGGCATAAACTCAAAGTTCATGCCATAGCTGGAAGCAACGAGCGTCGGCGGCAGGAATa                                                                       >  1:4074728/1‑71 (MQ=255)
ccACTTCAGTTCTGGCATAAACTCAAAGTTCATGCCATAGCTGGAAGCAACGAGCGTCGGCGGCAGGAATa                                                                       >  1:3599123/1‑71 (MQ=255)
  aCTTCAGTTCTGGCATAAACTCAAAGTTCATGCCATAGCTGGAAGCAACGAGCGTCGGCGGCAGGAATAc                                                                      >  1:3149922/1‑70 (MQ=255)
  aCTTCAGTTCTGGCATAAACTCAAAGTTCATGCCATAGCTGGAAGCAACGAGCGTCGGCGGCAGGAATAc                                                                      >  1:5086693/1‑70 (MQ=255)
  aCTTCAGTTCTGGCATAAACTCAAAGTTCATGCCATAGCTGGAAGCAACGAGCGTCGGCGGCAGGAATAc                                                                      >  1:2302459/1‑70 (MQ=255)
  aCTTCAGTTCTGGCATAAACTCAAAGTTCATGCCATAGCTGGAAGCAACGAGCGTCGGCGGCAGGAATAc                                                                      >  1:1227991/1‑70 (MQ=255)
  aCTTCAGTTCTGGCATAAACTCAAAGTTCATGCCATAGCTGGAAGCAACGAGCGTCGGCGGCAGGAATAc                                                                      >  1:5036079/1‑70 (MQ=255)
    ttCAGTTCTGGCATAAACTCAAAGTTCATGCCATAGCTGGAAGCAACGAGCGTCGGCGGCAGGAATACCAc                                                                   >  1:4681309/1‑71 (MQ=255)
    ttCAGTTCTGGCATAAACTCAAAGTTCATGCCATAGCTGGAAGCAACGAGCGTCGGCGGCAGGAATACCAc                                                                   >  1:1489538/1‑71 (MQ=255)
    ttCAGTTCTGGCATAAACTCAAAGTTCATGCCATAGCTGGAAGCAACGAGCGTCGGCGGCAGGAATACCAc                                                                   >  1:2626403/1‑71 (MQ=255)
    ttCAGTTCTGGCATAAACTCAAAGTTCATGCCATAGCTGGAAGCAACGAGCGTCGGCGGCAGGAATACCAc                                                                   >  1:5075900/1‑71 (MQ=255)
              gCATAAACTCAAAGTTCATGCCATAGCTGGAAGCAACGAGCGTCGGCGGCAGGAATACCACGGATacca                                                           <  1:4314796/69‑1 (MQ=255)
                   aaCTCAAAGTTCATGCCATAGCTGGAAGCAACGAGCGTCGGCGGCAGGAATACCACGGATACCACCgag                                                      >  1:3349645/1‑69 (MQ=255)
                   aaCTCAAAGTTCATGCCATAGCTGGAAGCAACGAGCGTCGGCGGCAGGAATACCACGGATACCACCgag                                                      >  1:4509401/1‑69 (MQ=255)
                   aaCTCAAAGTTCATGCCATAGCTGGAAGCAACGAGCGTCGGCGGCAGGAATACCACGGATACCACCgag                                                      >  1:4151875/1‑69 (MQ=255)
                      tCAAAGTTCATGCCATAGCTGGAAGCAACGAGCGTCGGCGGCAGGAATACCACGGAt                                                               <  1:4545789/57‑1 (MQ=255)
                          aGTTCATGCCATAGCTGGAAGCAACGAGCGTCGGCGGCAGGAATACCACGGATag                                                             >  1:1389904/1‑54 (MQ=255)
                          aGTTCATGCCATAGCTGGAAGCAACGAGCGTCGGCGGCAGGAATACCACGGATac                                                             >  1:1920749/1‑55 (MQ=255)
                          aGTTCATGCCATAGCTGGAAGCAACGAGCGTCGGCGGCAGGAATACCACGGATac                                                             >  1:2692619/1‑55 (MQ=255)
                          aGTTCATGCCATAGCTGGAAGCAACGAGCGTCGGCGGCAGGAATACCACGGATac                                                             >  1:4255323/1‑55 (MQ=255)
                          aGTTCATGCCATAGCTGGAAGCAACGAGCGTCGGCGGCAGGAATACCACGGATac                                                             >  1:4434747/1‑55 (MQ=255)
                          aGTTCATGCCATAGCTGGAAGCAACGAGCGTCGGCGGCAGGAATACCACGGATac                                                             >  1:175046/1‑55 (MQ=255)
                          aGTTCATGCCATAGCTGGAAGCAACGAGCGTCGGCGGCAGGAATACCACGGATac                                                             >  1:163453/1‑55 (MQ=255)
                          aGTTCATGCCATAGCTGGAAGCAACGAGCGTCGGCGGCAGGAATACCACGGATac                                                             >  1:477543/1‑55 (MQ=255)
                          aGTTCATGCCATAGCTGGAAGCAACGAGCGTCGGCGGCAGGAATACCACGGATACCACCGAGAAGAtttt                                              >  1:2085929/1‑70 (MQ=255)
                          aGTTCATGCCATAGCTGGAAGCAACGAGCGTCGGCGGCAGGAATACCACGGATACCACCGAGAAGAtttt                                              >  1:3837953/1‑70 (MQ=255)
                          aGTTCATGCCATAGCTGGAAGCAACGAGCGTCGGCGGCAGGAATACCACGGATACCACCGAGAAGAtttt                                              >  1:3351686/1‑70 (MQ=255)
                          aGTTCATGCCATAGCTGGAAGCAACGAGCGTCGGCGGCAGGAATACCACGGATACCACCGAGAAGAtttt                                              >  1:3221235/1‑70 (MQ=255)
                          aGTTCATGCCATAGCTGGAAGCAACGAGCGTCGGCGGCAGGAATACCACGGATACCACCGAGAAGAtttt                                              >  1:2932107/1‑70 (MQ=255)
                          aGTTCATGCCATAGCTGGAAGCAACGAGCGTCGGCGGCAGGAATACCACGGATACCACCGAGAAGAtttt                                              >  1:5139216/1‑70 (MQ=255)
                          aGTTCATGCCATAGCTGGAAGCAACGAGCGTCGGCGGCAGGAATACCACCGATac                                                             >  1:4440803/1‑55 (MQ=255)
                            ttCATGCCATAGCTGGAAGCAACGAGCGTCGGCGGCAGGAATACCACGGATACCACCGAGAAGATTTtgat                                           <  1:401522/71‑1 (MQ=255)
                            ttCATGCCATAGCTGGAAGCAACGAGCGTCGGCGGCAGGAATACCACGGATACCACCGAGAAGATTTtgat                                           <  1:1142886/71‑1 (MQ=255)
                                tGCCATAGCTGGAAGCAACGAGCGTCGGCGGCAGGAATACCACGGATACCACTGAGAAGATTTTgatgat                                        >  1:3294223/1‑70 (MQ=255)
                                tGCCATAGCTGGAAGCAACGAGCGTCGGCGGCAGGAATACCACGGATACCACCGAGAAGATTTTgatgat                                        >  1:98226/1‑70 (MQ=255)
                                tGCCATAGCTGGAAGCAACGAGCGTCGGCGGCAGGAATACCACGGATACCACCGAGAAGATTTTGatgatg                                       >  1:2458664/1‑71 (MQ=255)
                                       gCTGGAAGCAACGAGCGTCGGCGGCAGGAATACCACGGATACCACCGAGAAGATTTTGATGATGCGGTTCt                                >  1:3331190/1‑71 (MQ=255)
                                       gCTGGAAGCAACGAGCGTCGGCGGCAGGAATACCACGGATACCACCGAGAAGATTTTGATGATGCGGTTCt                                >  1:2531082/1‑71 (MQ=255)
                                           gAAGCAACGAGCGTCGGCGGCAGGAATACCACGGATACCACCGAGAAGATTTTGATGATGCGGTTCTGCTc                            <  1:2526692/71‑1 (MQ=255)
                                           gAAGCAACGAGCGTCGGCGGCAGGAATACCACGGATACCACCGAGAAGATTTTGAAGATGCGGTTCTGCTc                            <  1:4773092/71‑1 (MQ=255)
                                            aaGCAACGAGCGTCGGCGGCAGGAATACCACGGATACCACCGAGAAGATTTtgat                                           <  1:1574528/55‑1 (MQ=255)
                                             aGCAACGAGCGTCGGCGGCAGGAATACCACGGATACCACCGAGAAGATTTTGATGATGCGGTTCTGCTCGa                          <  1:2167132/71‑1 (MQ=255)
                                             aGCAACGAGCGTCGGCGGCAGGAATACCACGGATACCACCGAGAAGATTTTGATGATGCGGCTCTGCTCGa                          <  1:3130823/71‑1 (MQ=255)
                                              gCAACGAGCGTCGGCGGCAGGAATACCACGGATACCACCGAGAAGATTTTGATGATGCGGTTCTGCTCGAt                         <  1:3354515/71‑1 (MQ=255)
                                                       gTCGGCGGCAGGAAGGCGACCGAGGCCACCGAGAAGATCTTGATGATCTTGTTCTGGTTGATGTTGATgaa                <  1:1813713/71‑1 (MQ=255)
                                                       gTCGGCGGCAGGAAGGCGACCGAGGCCACCGAGAAGATCTTGATGATCTTGTTCTGGTTGATGTTGATgaa                <  1:2447158/71‑1 (MQ=255)
                                                       gTCGGCGGCAGGAAGGCGACCGAGGCCACCGAGAAGATCTTGATGATCTTGTTCTGGTTGATGTTGATgaa                <  1:5057904/71‑1 (MQ=255)
                                                       gTCGGCGGCAGGAAGGCGACCGAGGCCACCGAGAAGATCTTGATGATCTTGTTCTGGTTGATGTTGATgaa                <  1:5072763/71‑1 (MQ=255)
                                                       gTCGGCGGCAGGAAGGCGACCGAGGCCACCGAGAAGATCTTGATGATCTTGTTCTGGTTGATGTTGATgaa                <  1:549962/71‑1 (MQ=255)
                                                        tCGGCGGCAGGAATACCACGGATACCACCGAGAAGATTTTGATGATGCGGTTCTGCTCGATGTTGATaaaa               >  1:258502/1‑71 (MQ=255)
                                                          ggcggcAGGAAGGCGACCGAGGCCACCGAGAAGATCTTGATGATCTTGTTCTGGTTGATGTTGATgaagc              >  1:4217374/1‑68 (MQ=255)
                                                          ggcggcAGGAAGGCGACCGAGGCCACCGAGAAGATCTTGATGATCTTGTTCTGGTTGATGTTGATGAAGcc             >  1:1692191/1‑71 (MQ=255)
                                                          ggcggcAGGAAGGCGACCGAGGCCACCGAGAAGATCTTGATGATCTTGTTCTGGTTGATGTTGATGAAGcc             >  1:4184570/1‑71 (MQ=255)
                                                              gcAGGAATACCACGGATACCACCGAGAAGATTTTGATGATGCGGTTCTGCTCGATGTTGAt                   <  1:1228375/61‑1 (MQ=255)
                                                              gcAGGAATACCACGGATACCACCGAGAAGATTTTGATGATGCGGTTCTGCTCGATGTTGAt                   <  1:1777772/61‑1 (MQ=255)
                                                              gcAGGAATACCACGGATACCACCGAGAAGATTTTGATGATGCGGTTCTGCTCGATGTTGAt                   <  1:1047516/61‑1 (MQ=255)
                                                              gcAGGAATACCACGGATACCACCGAGAAGATTTTGATGATGCGGTTCTGCTCGATGTTGAt                   <  1:3384193/61‑1 (MQ=255)
                                                              gcAGGAATACCACGGATACCACCGAGAAGATTTTGATGATGCGGTTCTGCTCGATGTTGAt                   <  1:4141456/61‑1 (MQ=255)
                                                              gcAGGAATACCACGGATACCACCGAGAAGATTTTGATGATGCGGTTCTGCTCGATGTTGAt                   <  1:3557998/61‑1 (MQ=255)
                                                              gcAGGAATACCACGGATACCACCGAGAAGATTTTGATGATGCGGTTCTGCTCGATGTTGATAAAACCCAtt         <  1:1087850/71‑1 (MQ=255)
                                                              gcAGGAATACCACGGATACCACCGAGAAGATTTTGATGATGCGGTTCTGCTCGATGTTGATAAAACCCAtt         <  1:2477986/71‑1 (MQ=255)
                                                               cAGGAATACCACGGATACCACCGAGAAGATTTTGATGATGCGGTTCTGCTCGATGTTGATAAAACCCAt          >  1:812618/1‑69 (MQ=255)
                                                               cAGGAATACCACGGATACCACCGAGAAGATTTTGATGATGCGGTTCTGCTCGATGTTGATAAAACCCAt          >  1:3800745/1‑69 (MQ=255)
                                                                aGGAATACCACGGATACCACCGAGAAGATTTTGATGATGCGGTTCTGCTCGATGTTGATaa                 >  1:251984/1‑61 (MQ=255)
                                                                  gAATACCACGGATACCACCGAGAAGATTTTGATGATGCGGTTCTGCTCGATGTTGATAAAACCCAt          >  1:2162886/1‑66 (MQ=255)
                                                                  gAATACCACGGATACCACCGAGAAGATTTTGATGATGCGGTTCTGCTCGATGTTGATAAAACCCAt          >  1:1675692/1‑66 (MQ=255)
                                                                  gAATACCACGGATACCACCGAGAAGATTTTGATGATGCGGTTCTGCTCGATGTTGATAAAACCCAt          >  1:1695809/1‑66 (MQ=255)
                                                                    ataccacggataccacggAGAAGATTTTGATGCTGCGGTTCTGCTCGATGTTGATAAAAccc            >  1:2891766/1‑62 (MQ=255)
                                                                    aTACCACGGATACCACCGAGAAGATTTTGATGATGCGGTTCTGCTCGATGTTGATa                  >  1:4469958/1‑56 (MQ=255)
                                                                    aTACCACGGATACCACCGAGAAGATTTTGATGATGCGGTTCTGCTCGATGTTGATAAAAccc            >  1:4330856/1‑62 (MQ=255)
                                                                     tACCACGGATACCACCGAGAAGATTTTGATGATGCGGTTCTGCTCGATGTTGATAAAACCCATTGCCGCCt  >  1:1174252/1‑71 (MQ=255)
                                                                      aCCACGGATACCACCGAGAAGATTTTGATGATGCGGTTCTGc                              <  1:1383812/42‑1 (MQ=255)
                                                                      aCCACGGATACCACCGAGAAGATTTTGATGATGCGGTTCTGc                              <  1:3525462/42‑1 (MQ=255)
                                                                      aCCACGGATACCACCGAGAAGATTTTGATGATGCGGTTCTGc                              <  1:1085906/42‑1 (MQ=255)
                                                                      aCCACGGATACCACCGAGAAGATTTTGATGATGCGGTTCTGc                              <  1:699964/42‑1 (MQ=255)
                                                                      |                                                                     
CCACTTCAGTTCTGGCATAAACTCAAAGTTCATGCCATAGCTGGAAGCAACGAGCGTCGGCGGCAGGAATACCACGGATACCACCGAGAAGATTTTGATGATGCGGTTCTGCTCGATGTTGATAAAACCCATTGCCGCCT  >  minE/2396324‑2396463

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: