Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 585,997 A→T 37.0% V555V (GTT→GTA cydD ← fused cysteine transporter subunits and membrane component and ATP‑binding component of ABC superfamily

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE585,9970AT37.0% 69.4 / 75.2 89V555V (GTT→GTAcydDfused cysteine transporter subunits and membrane component and ATP‑binding component of ABC superfamily
Reads supporting (aligned to +/- strand):  ref base A (27/29);  new base T (18/15);  total (45/44)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATGGGCCACCAGCCACACTTAATTCCGCGTAACGTCCTTGCTCAATAATCCGGCCATCCTGCATAACCCAAATGACATCCCAGTCAGCAAGATCTTCTAACTGGTGGGTGACCATTAACGTTGTCTG  >  minE/585932‑586059
                                                                 |                                                              
aaTGGGGCGCCAGCCACACTTAATTCCGCGTAACGTCCTTGCTCAATAATCTGACCATCCTGCATTAc                                                              >  1:48847/1‑68 (MQ=255)
aaTGGGCCGCCAGCCACACTTAATTCCGCGTAACGTCCTTGCTCAATAATCTGACCATCCTGCATTAc                                                              >  1:2328049/1‑68 (MQ=255)
aaTGGGCCGCCAGCCACACTTAATTCCGCGTAACGTCCTTGCTCAATAATCTGACCATCCTGCATTAc                                                              >  1:2091397/1‑68 (MQ=255)
aaTGGGCCGCCAGCCACACTTAATTCCGCGTAACGTCCTTGCTCAATAATCTGACCATCCTGCATTAc                                                              >  1:1942871/1‑68 (MQ=255)
aaTGGGCCGCCAGCCACACTTAATTCCGCGTAACGTCCTTGCTCAATAATCTGACCATCCTGCATTAc                                                              >  1:193236/1‑68 (MQ=255)
aaTGGGCCGCCAGCCACACTTAATTCCGCGTAACGTCCTTGCTCAATAATCTGACCATCCTGCATTAc                                                              >  1:1894418/1‑68 (MQ=255)
aaTGGGCCGCCAGCCACACTTAATTCCGCGTAACGTCCTTGCTCAATAATCTGACCATCCTGCATTAc                                                              >  1:2129693/1‑68 (MQ=255)
aaTGGGCCGCCAGCCACACTTAATTCCGCGTAACGTCCTTGCTCAATAATCTGACCATCCTGCATTAc                                                              >  1:2669593/1‑68 (MQ=255)
aaTGGGCCGCCAGCCACACTTAATTCCGCGTAACGTCCTTGCTCAATAATCTGACCATCCTGCATTAc                                                              >  1:2697749/1‑68 (MQ=255)
aaTGGGCCGCCAGCCACACTTAATTCCGCGTAACGTCCTTGCTCAATAATCTGACCATCCTGCATTAc                                                              >  1:493621/1‑68 (MQ=255)
aaTGGGCCGCCAGCCACACTTAATTCCGCGTAACGTCCTTGCTCAATAATCTGACCATCCTGCATTAc                                                              >  1:1588206/1‑68 (MQ=255)
aaTGGGCCGCCAGCCACACTTAATTCCGCGTAACGTCCTTGCTCAATAATCTGACCATCCTGCATTAc                                                              >  1:2889982/1‑68 (MQ=255)
aaTGGGCCGCCAGCCACACTTAATTCCGCGTAACGTCCTTGCTCAATAATCTGACCATCCTGCATTAc                                                              >  1:3047904/1‑68 (MQ=255)
aaTGGGCCGCCAGCCACACTTAATTCCGCGTAACGTCCTTGCTCAATAATCTGACCATCCTGCATTAc                                                              >  1:1522732/1‑68 (MQ=255)
aaTGGGCCGCCAGCCACACTTAATTCCGCGTAACGTCCTTGCTCAATAATCTGACCATCCTGCATTAc                                                              >  1:603908/1‑68 (MQ=255)
aaTGGGCCGCCAGCCACACTTAATTCCGCGTAACGTCCTTGCTCAATAATCTGACCATCCTGCATTAc                                                              >  1:3097572/1‑68 (MQ=255)
aaTGGGCCGCCAGCCACACTTAATTCCGCGTAACGTCCTTGCTCAATAATCTGACCATCCTGCATTAc                                                              >  1:3145816/1‑68 (MQ=255)
aaTGGGCCGCCAGCCACACTTAATTCCGCGTAACGTCCTTGCTCAATAATCTGACCATCCTGCATTAc                                                              >  1:2465614/1‑68 (MQ=255)
  tGGGCCACCAGCCACACTTAATTCCGCGTAACGTCCTTGCTCAATAATCCGGCCATCCTGCATAACCCa                                                           >  1:458935/1‑69 (MQ=255)
  tGGGCCACCAGCCACACTTAATTCCGCGTAACGTCCTTGCTCAATAATCCGGCCATCCTGCATAACCCa                                                           >  1:2483220/1‑69 (MQ=255)
  tGGGCCACCAGCCACACTTAATTCCGCGTAACGTCCTTGCTCAATAATCCGGCCATCCTGCATAACCCa                                                           >  1:3211149/1‑69 (MQ=255)
  tGGGCCACCAGCCACACTTAATTCCGCGTAACGTCCTTGCTCAATAATCCGGCCATCCTGCATAACCCa                                                           >  1:3183625/1‑69 (MQ=255)
           aGCCACACTTAATTCCGCGTAACGTCCTTGCTCAATAATCCGGCCATCCTGCATAACCCAAATGACAt                                                   >  1:2703576/1‑68 (MQ=255)
              cacaCTTAATTCCGCGTAACGTCCTTGCTCAATAATCCGGCCATCCTGCATAACCCAAATGACATCCCa                                               <  1:426582/69‑1 (MQ=255)
              cacaCTTAATTCCGCGTAACGTCCTTGCTCAATAATCCGGCCATCCTGCATAACCCAAATGACATCCCa                                               <  1:2149091/69‑1 (MQ=255)
              cacaCTTAATTCCGCGTAACGTCCTTGCTCAATAATCCGGCCATCCTGCATAACCCAAATGACATCCCa                                               <  1:1807029/69‑1 (MQ=255)
              cacaCTTAATTCCGCGTAACGTCCTTGCTCAATAATCCGGCCATCCTGCATAACCCAAATGACATCCCa                                               <  1:850117/69‑1 (MQ=255)
              cacaCTTAATTCCGCGTAACGTCCTTGCTCAATAATCCGGCCATCCTGCATAACCCAAATGACATCCCa                                               <  1:1107010/69‑1 (MQ=255)
              cacaCTTAATTCCGCGTAACGTCCTTGCTCAATAATCCGGCCATCCTGCATAACCCAAATGACATCCCa                                               <  1:2160879/69‑1 (MQ=255)
               acacTTAATTCCGCGTAACGTCCTTGCTCAATAATCCGGCCATCCTGCATAACCCAAATGACATCCCa                                               <  1:3343583/68‑1 (MQ=255)
                cacTTAATTCCGCGTAACGTCCTTGCTCAATAATCCGGCCATCCTGCATAACCCAAATGACATccaa                                               <  1:2260333/67‑3 (MQ=255)
                        tCCGCGTAACGTCCTTGCTCAATAATCCGGCCATCCTGCATAACCCAAATGACATCCCAGTCAGCAAg                                      <  1:1329655/68‑1 (MQ=255)
                        tCCGCGTAACGTCCTTGCTCAATAATCCGGCCATCCTGCATAACCCAAATGACATCCCAGTCAGCAAGa                                     <  1:291656/69‑1 (MQ=255)
                              tAACGTCCTTGCTCAATAATCCGGCCATCCTGCATAACCCAAATGACATCCCAGTCAGCAAGAtcttct                               >  1:525651/1‑69 (MQ=255)
                              tAACGTCCTTGCTCAATAATCCGGCCATCCTGCATAACCCAAATGACATCCCAGTCAGCAAGAtcttct                               >  1:262611/1‑69 (MQ=255)
                              tAACGTCCTTGCTCAATAATCCGGCCATCCTGCATAACCCAAATGACATCCCAGTCAGCAAGAtcttct                               <  1:4103/69‑1 (MQ=255)
                              tAACGTCCTTGCTCAATAATCCGGCCATCCTGCATAACCCAAATGACATCCCAGTCAGCAAGAtcttct                               <  1:2241409/69‑1 (MQ=255)
                              tAACGTCCTTGCTCAATAATCCGGCCATCCTGCATAACCCAAATGACATCCCAGTCAGCAAGAtcttct                               >  1:617126/1‑69 (MQ=255)
                              tAACGTCCTTGCTCAATAATCCGGCCATCCTGCATAACCCAAATGACATCCCAGTCAGCAAGAtcttct                               >  1:1769293/1‑69 (MQ=255)
                              tAACGTCCTTGCTCAATAATCCGGCCATCCTGCATAACCCAAATGACATCCCAGTCAGCAAGAtcttc                                >  1:2740864/1‑68 (MQ=255)
                                     cTTGCTCAATAATCCGGCCATCCTGCATAACCCAAATGACATCCCAGTCAGCAAGATCTTCTAACtggt                        <  1:2326978/69‑1 (MQ=255)
                                     cTTGCTCAATAATCCGGCCATCCTGCATAACCCAAATGACATCCCAGTCAGCAAGATCTTCTAACtggt                        <  1:862042/69‑1 (MQ=255)
                                     cTTGCTCAATAATCCGGCCATCCTGCATAACCCAAATGACATCCCAGTCAGCAAGATCTTCTAACtggt                        <  1:331706/69‑1 (MQ=255)
                                     cTTGCTCAATAATCCGGCCATCCTGCATAACCCAAATGACATCCCAGTCAGCAAGATCTTCTAACtggt                        <  1:805937/69‑1 (MQ=255)
                                     cTTGCTCAATAATCCGGCCATCCTGCATAACCCAAATGACATCCCAGTCAGCAAGATCTTCTAACtggt                        <  1:803052/69‑1 (MQ=255)
                                     cTTGCTCAATAATCCGGCCATCCTGCATAACCCAAATGACATCCCAGTCAGCAAGATCTTCTAACtggt                        >  1:700610/1‑69 (MQ=255)
                                     cTTGCTCAATAATCCGGCCATCCTGCATAACCCAAATGACATCCCAGTCAGCAAGATCTTCTAACtggt                        >  1:1508284/1‑69 (MQ=255)
                                     cTTGCTCAATAATCCGGCCATCCTGCATAACCCAAATGACATCCCAGTCAGCAAGATCTTCTAACtggt                        <  1:1528722/69‑1 (MQ=255)
                                     cTTGCTCAATAATCCGGCCATCCTGCATAACCCAAATGACATCCCAGTCAGCAAGATCTTCTAACtggt                        >  1:1574930/1‑69 (MQ=255)
                                     cTTGCTCAATAATCCGGCCATCCTGCATAACCCAAATGACATCCCAGTCAGCAAGATCTTCTAACtggt                        <  1:1587616/69‑1 (MQ=255)
                                     cTTGCTCAATAATCCGGCCATCCTGCATAACCCAAATGACATCCCAGTCAGCAAGATCTTCTAACtggt                        <  1:1600252/69‑1 (MQ=255)
                                     cTTGCTCAATAATCCGGCCATCCTGCATAACCCAAATGACATCCCAGTCAGCAAGATCTTCTAACtggt                        <  1:594172/69‑1 (MQ=255)
                                     cTTGCTCAATAATCCGGCCATCCTGCATAACCCAAATGACATCCCAGTCAGCAAGATCTTCTAACtggt                        <  1:2267019/69‑1 (MQ=255)
                                     cTTGCTCAATAATCCGGCCATCCTGCATAACCCAAATGACATCCCAGTCAGCAAGATCTTCTAACtggt                        <  1:2653324/69‑1 (MQ=255)
                                     cTTGCTCAATAATCCGGCCATCCTGCATAACCCAAATGACATCCCAGTCAGCAAGATCTTCTAACtggt                        <  1:1781287/69‑1 (MQ=255)
                                     cTTGCTCAATAATCCGGCCATCCTGCATAACCCAAATGACATCCCAGTCAGCAAGATCTTCTAACtggt                        <  1:563233/69‑1 (MQ=255)
                                     cTTGCTCAATAATCCGGCCATCCTGCATAACCCAAATGACATCCCAGTCAGCAAGATCTTCTAACtggt                        <  1:1787670/69‑1 (MQ=255)
                                     cTTGCTCAATAATCCGGCCATCCTGCATAACCCAAATGACATCCCAGTCAGCAAGATCTTCTAACtggt                        >  1:2546377/1‑69 (MQ=255)
                                     cTTGCTCAATAATCCGGCCATCCTGCATAACCCAAATGACATCCCAGTCAGCAAGATCTTCTAACtggt                        <  1:2541508/69‑1 (MQ=255)
                                     cTTGCTCAATAATCCGGCCATCCTGCATAACCCAAATGACATCCCAGTCAGCAAGATCTTCTAACtggt                        >  1:251216/1‑69 (MQ=255)
                                      ttGCTCAATTTTCCGGCGATCCTGCATAACCCAAATGACATCCCAGTCAGCAAGATCTTCTAACtggt                        <  1:1740590/68‑1 (MQ=255)
                                              taatCCGGCCATCCTGCATAACCCAAATGACATCCCAGTCAGCAAGATCTTCTAACTGGTGGGTGAcc                >  1:469805/1‑68 (MQ=255)
                                               aatCCGGCCATCCTGCATAACCCAAATGACATCCCAGTCAGCAAGa                                     >  1:3331269/1‑46 (MQ=255)
                                               aatCCGGCCATCCTGCATAACCCAAATGACATCCCAGTCAGCAAGATCTTCTAACTGGTGGGTGACCa               >  1:2534428/1‑68 (MQ=255)
                                               aatCCGGCCATCCTGCATAACCCAAATGACATCCCAGTCAGCAAGATCTTCTAACTGGTGGGTGACCa               >  1:2000678/1‑68 (MQ=255)
                                               aatCCGGCCATCCTGCATAACCCAAATGACATCCCAGTCAGCAAGATCTTCTAACTGGTGGGTGACCa               >  1:1341974/1‑68 (MQ=255)
                                               aatCCGGCCATCCTGCATAACCCAAATGACATCCCAGTCAGCAAGATCTTCTAACTGGTGGGTGACCa               >  1:1662179/1‑68 (MQ=255)
                                               aatCCGGCCACCCTGCATAACCCAAATGACATCCCAGTCAGCAAGATCTTCTAACTGGTGGGTGACCa               >  1:899866/1‑68 (MQ=255)
                                                     gCCATCCTGCATAACCCAAATGACATCCCAGTCAGCAAGATCTTCTAACTGGTGGGTGACCATTAACGt        <  1:2175436/69‑1 (MQ=255)
                                                           cTGCATTACCCAAATGACATCCCAGTCAGCAAGATCTTCTATCTGGTGGGTGACCATTAACGTTGTCTg  <  1:995805/69‑1 (MQ=255)
                                                           cTGCATTACCCAAATGACATCCCAGTCAGCAAGATCTTCTAACTGGTGGGTGACCATTAACGTTGTCtt  <  1:2699330/69‑2 (MQ=255)
                                                           cTGCATTACCCAAATGACATCCCAGTCAGCAAGATCTTCTAACTGGTGGGTGACCATTAACGTTGTCTg  <  1:1440130/69‑1 (MQ=255)
                                                           cTGCATTACCCAAATGACATCCCAGTCAGCAAGATCTTCTAACTGGTGGGTGACCATTAACGTTGTCTg  <  1:1034359/69‑1 (MQ=255)
                                                           cTGCATTACCCAAATGACATCCCAGTCAGCAAGATCTTCTAACTGGTGGGTGACCATTAACGTTGTCTg  <  1:1076358/69‑1 (MQ=255)
                                                           cTGCATTACCCAAATGACATCCCAGTCAGCAAGATCTTCTAACTGGTGGGTGACCATTAACGTTGTCTg  <  1:1216411/69‑1 (MQ=255)
                                                           cTGCATTACCCAAATGACATCCCAGTCAGCAAGATCTTCTAACTGGTGGGTGACCATTAACGTTGTCTg  <  1:1220606/69‑1 (MQ=255)
                                                           cTGCATTACCCAAATGACATCCCAGTCAGCAAGATCTTCTAACTGGTGGGTGACCATTAACGTTGTCTg  <  1:708869/69‑1 (MQ=255)
                                                           cTGCATTACCCAAATGACATCCCAGTCAGCAAGATCTTCTAACTGGTGGGTGACCATTAACGTTGTCTg  <  1:1445357/69‑1 (MQ=255)
                                                           cTGCATTACCCAAATGACATCCCAGTCAGCAAGATCTTCTAACTGGTGGGTGACCATTAACGTTGTCTg  <  1:1657892/69‑1 (MQ=255)
                                                           cTGCATTACCCAAATGACATCCCAGTCAGCAAGATCTTCTAACTGGTGGGTGACCATTAACGTTGTCTg  <  1:2577803/69‑1 (MQ=255)
                                                           cTGCATTACCCAAATGACATCCCAGTCAGCAAGATCTTCTAACTGGTGGGTGACCATTAACGTTGTCTg  <  1:2826833/69‑1 (MQ=255)
                                                           cTGCATTACCCAAATGACATCCCAGTCAGCAAGATCTTCTAACTGGTGGGTGACCATTAACGTTGTCTg  <  1:2861560/69‑1 (MQ=255)
                                                           cTGCATTACCCAAATGACATCCCAGTCAGCAAGATCTTCTAACTGGTGGGTGACCATTAACGTTGTCTg  <  1:3053068/69‑1 (MQ=255)
                                                           cTGCATAACCCAAATGACATCCCAGTCAGCAAGATCTTCTAACt                           >  1:3161930/1‑44 (MQ=255)
                                                           cTGCATAACCCAAATGACATCCCAGTCAGCAAGATCTTCTAACTggtgggt                    >  1:2373168/1‑51 (MQ=255)
                                                           cTGCATAACCCAAATGACATCCCAGTCAGCAAGATCTTCTAACTGGTGGGTGACCATTAACGTTGTCTg  >  1:1372855/1‑69 (MQ=255)
                                                           cTGCATAACCCAAATGACATCCCAGTCAGCAAGATCTTCTAACTGGTGGGTGACCATTAACGTTGTCTg  >  1:127710/1‑69 (MQ=255)
                                                           cTGCATAACCCAAATGACATCCCAGTCAGCAAGATCTTCTAACTGGTGGGTGACCATTAACGTTGTCTg  >  1:814837/1‑69 (MQ=255)
                                                            tGCATTACCCAAATGACATCCCAGTCAGCAAGATCTTCTAACTGGTGGGTGACCATTAACGTTGTCTg  <  1:1144558/68‑1 (MQ=255)
                                                                 |                                                              
AATGGGCCACCAGCCACACTTAATTCCGCGTAACGTCCTTGCTCAATAATCCGGCCATCCTGCATAACCCAAATGACATCCCAGTCAGCAAGATCTTCTAACTGGTGGGTGACCATTAACGTTGTCTG  >  minE/585932‑586059

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