Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 584,947 G→A 46.8% D316D (GAC→GAT cydC ← 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
*minE584,9470GA46.8% 13.0 / 79.5 94D316D (GAC→GATcydCfused cysteine transporter subunits and membrane component and ATP‑binding component of ABC superfamily
Reads supporting (aligned to +/- strand):  ref base G (24/26);  new base A (26/18);  total (50/44)
Fisher's exact test for biased strand distribution p-value = 3.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

AAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCC  >  minE/584883‑585013
                                                                |                                                                  
aaaCGCGATCGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGTTTTTGATCCGTTAAATCag                                                                <  1:2850113/69‑1 (MQ=255)
aaaCGCGATCGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGTTTTTGATCCGTTAAATCag                                                                <  1:636638/69‑1 (MQ=255)
aaaCGCGATCGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGTTTTTGATCCGTTAAATCag                                                                <  1:2108167/69‑1 (MQ=255)
aaaCGCGATCGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGTTTTTGATCCGTTAAATCag                                                                <  1:3191267/69‑1 (MQ=255)
aaaCGCGATCGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGTTTTTGATCCGTTAAATCag                                                                <  1:2220452/69‑1 (MQ=255)
aaaCGCGATCGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGTTTTTGATCCGTTAAATCag                                                                <  1:2320704/69‑1 (MQ=255)
aaaCGCGATCGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGTTTTTGATCCGTTAAATCag                                                                <  1:2521649/69‑1 (MQ=255)
aaaCGCGATCGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGTTTTTGATCCGTTAAATCag                                                                <  1:2986330/69‑1 (MQ=255)
aaaCGCGATCGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGTTTTTGATCCGTTAAATCag                                                                <  1:1391922/69‑1 (MQ=255)
aaaCGCGATCGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGTTTTTGATCCGTTAAATCag                                                                <  1:790720/69‑1 (MQ=255)
aaaCGCGATCGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGTTTTTGATCCGTTAAATCag                                                                <  1:2732936/69‑1 (MQ=255)
aaaCGCGATCGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGTTTTTGATCCGTTAAATCag                                                                <  1:1312758/69‑1 (MQ=255)
aaaCGCGATCGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGTTTTTGATCCGTTAAATCag                                                                <  1:953526/69‑1 (MQ=255)
aaaCGCGATCGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGTTTTTGATCCGTTAAATCag                                                                <  1:1095508/69‑1 (MQ=255)
 aaCGCGATCGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGTTTTTGATCCGTTAAATCag                                                                <  1:2728868/68‑1 (MQ=255)
 aaCGCGATCGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGTTTTTGATCCGTTAAATCag                                                                <  1:2762838/68‑1 (MQ=255)
 aaCGCGATCGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGTTTTTGATCCGTTAAATCag                                                                <  1:775721/68‑1 (MQ=255)
 aaCGCGATCGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGTTTTTGATCCGTTAAATCag                                                                <  1:71539/68‑1 (MQ=255)
 aaCGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCaga                                                               <  1:199720/69‑1 (MQ=255)
  aCGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCaga                                                               <  1:132196/68‑1 (MQ=255)
       aTCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGAtacg                                                         <  1:2837226/69‑1 (MQ=255)
       aTCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGAtacg                                                         <  1:2886885/69‑1 (MQ=255)
         cGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGTTTTTGATCCGTTAAATCAGTGAtccgta                                                       >  1:1957572/1‑69 (MQ=255)
         cGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGTTTTTGATCCGTTAAATCAGTGAtccgta                                                       >  1:501320/1‑69 (MQ=255)
         cGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGTTTTTGATCCGTTAAATCAGTGAtccgta                                                       >  1:1573157/1‑69 (MQ=255)
         cGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGTTTTTGATCCGTTAAATCAGTGAtccgta                                                       >  1:2417322/1‑69 (MQ=255)
         cGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGTTTTTGATCCGTTAAATCAGTGAtccgta                                                       >  1:34068/1‑69 (MQ=255)
         cGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGTTTTTGATCCGTTAAATCAGTGAtccgta                                                       >  1:1283573/1‑69 (MQ=255)
         cGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGTTTTTGATCCGTTAAATCAGTGAtccgta                                                       >  1:1712305/1‑69 (MQ=255)
         cGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGTTTTTGATCCGTTAAATCAGTGAtccgta                                                       >  1:3385605/1‑69 (MQ=255)
         cGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGTTTTTGATCCGTTAAATCAGTGAtccgta                                                       >  1:2605155/1‑69 (MQ=255)
         cGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGTTTTTGATCCGTTAAATCAGTGAtccgta                                                       >  1:1646098/1‑69 (MQ=255)
         cGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGTTTTTGATCCGTTAAATCAGTGAtccgta                                                       >  1:2668414/1‑69 (MQ=255)
         cGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGTTTTTGATCCGTTAAATCAGTGAtccgta                                                       >  1:2263063/1‑69 (MQ=255)
                   aGTTTGGGTATCAGGAAAGGTGACCTCCGTTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGc                                             <  1:3095659/69‑1 (MQ=255)
                   aGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGc                                             <  1:2559149/69‑1 (MQ=255)
                   aGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGc                                             <  1:181623/69‑1 (MQ=255)
                   aGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGc                                             <  1:2973010/69‑1 (MQ=255)
                   aGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGc                                             <  1:2652395/69‑1 (MQ=255)
                   aGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGc                                             <  1:922357/69‑1 (MQ=255)
                   aGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGc                                             <  1:2777570/69‑1 (MQ=255)
                   aGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGc                                             <  1:120003/69‑1 (MQ=255)
                   aGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGc                                             <  1:1015901/69‑1 (MQ=255)
                   aGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGc                                             <  1:2872876/69‑1 (MQ=255)
                   aGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGc                                             <  1:3088804/69‑1 (MQ=255)
                   aGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGc                                             <  1:2323975/69‑1 (MQ=255)
                   aGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGc                                             <  1:558691/69‑1 (MQ=255)
                   aGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGc                                             <  1:3118903/69‑1 (MQ=255)
                   aGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGc                                             <  1:3331875/69‑1 (MQ=255)
                   aGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGc                                             <  1:740767/69‑1 (MQ=255)
                   aGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGc                                             <  1:2024600/69‑1 (MQ=255)
                               aGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACtttc                                   >  1:2105606/1‑65 (MQ=255)
                               aGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGcc                                  >  1:2105191/1‑68 (MQ=255)
                               aGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGc                                   >  1:1447063/1‑67 (MQ=255)
                               aGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGc                                   >  1:1594560/1‑67 (MQ=255)
                               aGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGc                                   >  1:1355932/1‑67 (MQ=255)
                               aGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGc                                   >  1:2735116/1‑67 (MQ=255)
                                      gTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGc                                             <  1:2761266/50‑1 (MQ=255)
                                        gACTTCCGGTTTTTGATCCGTTAAATCAGTGAtccgt                                                        >  1:1090744/1‑37 (MQ=25)
                                        gACTTCCGGTTTTTGATCCGTTAAATCAGTGAtccgt                                                        >  1:2919907/1‑37 (MQ=25)
                                        gACTTCCGGTTTTTGATCCGTTAAATCAGTGAtccgt                                                        >  1:78441/1‑37 (MQ=25)
                                        gACTTCCGGTTTTTGATCCGTTAAATCAGTGAtccgt                                                        >  1:2989821/1‑37 (MQ=25)
                                        gACTTCCGGTTTTTGATCCGTTAAATCAGTGAtccgt                                                        >  1:2899658/1‑37 (MQ=25)
                                        gACTTCCGGTTTTTGATCCGTTAAATCAGTGAtccgt                                                        >  1:1987505/1‑37 (MQ=25)
                                        gACTTCCGGTTTTTGATCCGTTAAATCAGTGAtccgt                                                        >  1:2084554/1‑37 (MQ=25)
                                        gACTTCCGGTTTTTGATCCGTTAAATCAGTGAtccgt                                                        >  1:2150907/1‑37 (MQ=25)
                                        gACTTCCGGTTTTTGATCCGTTAAATCAGTGAtccgt                                                        >  1:1681681/1‑37 (MQ=25)
                                        gACTTCCGGTTTTTGATCCGTTAAATCAGTGAtccgt                                                        >  1:1590842/1‑37 (MQ=25)
                                        gACTTCCGGTTTTTGATCCGTTAAATCAGTGAtccgt                                                        >  1:2557881/1‑37 (MQ=25)
                                        gACTTCCGGTTTTTGATCCGTTAAATCAGTGAtccgt                                                        >  1:2568208/1‑37 (MQ=25)
                                        gACTTCCGGTTTTTGATCCGTTAAATCAGTGAtccgt                                                        >  1:2596846/1‑37 (MQ=25)
                                        gACTTCCGGTTTTTGATCCGTTAAATCAGTGAtccgt                                                        >  1:2618648/1‑37 (MQ=25)
                                                 tttttGATCCGTTAAGTCAGATATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCAcc               <  1:1334439/69‑1 (MQ=255)
                                                  ttttGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCAcc               <  1:3091760/68‑1 (MQ=255)
                                                  ttttGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCAcc               <  1:1387319/68‑1 (MQ=255)
                                                         ccGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACt        <  1:2984505/68‑1 (MQ=255)
                                                         ccGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTg       >  1:1714620/1‑69 (MQ=255)
                                                         ccGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTg       >  1:1885132/1‑69 (MQ=255)
                                                         ccGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTg       >  1:1958362/1‑69 (MQ=255)
                                                         ccGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTg       >  1:968100/1‑69 (MQ=255)
                                                         ccGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTg       >  1:526077/1‑69 (MQ=255)
                                                         ccGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTg       >  1:998555/1‑69 (MQ=255)
                                                         ccGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTg       >  1:3421050/1‑69 (MQ=255)
                                                         ccGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTg       >  1:3344868/1‑69 (MQ=255)
                                                         ccGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTg       >  1:2211364/1‑69 (MQ=255)
                                                         ccGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTg       >  1:2964187/1‑69 (MQ=255)
                                                         ccGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTg       >  1:2817261/1‑69 (MQ=255)
                                                              aaGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGcc  >  1:1975113/1‑69 (MQ=255)
                                                              aaGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGcc  >  1:2165285/1‑69 (MQ=255)
                                                              aaGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGcc  >  1:3309495/1‑69 (MQ=255)
                                                              aaGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGcc  >  1:3252085/1‑69 (MQ=255)
                                                              aaGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGcc  >  1:2310151/1‑69 (MQ=255)
                                                              aaGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGcc  >  1:1105474/1‑69 (MQ=255)
                                                              aaGTCAGAGAAACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGcc  >  1:2402443/1‑69 (MQ=255)
                                                                |                                                                  
AAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCC  >  minE/584883‑585013

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: