Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 584,917 A→G 55.9% F326F (TTT→TTC 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,9170AG55.9% 27.7 / 79.6 102F326F (TTT→TTCcydCfused cysteine transporter subunits and membrane component and ATP‑binding component of ABC superfamily
Reads supporting (aligned to +/- strand):  ref base A (22/23);  new base G (28/29);  total (50/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.23e-01

Only 100 of 103 total aligned reads displayed.

ATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCC  >  minE/584851‑584981
                                                                  |                                                                
aTAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGaaa                                                                <  1:1195908/69‑1 (MQ=255)
     tgaactgaacATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTg                                                            <  1:498568/68‑1 (MQ=255)
     tgaactgaacATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTg                                                            <  1:1435148/68‑1 (MQ=255)
      gaactgaacATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGAc                                                          >  1:1950797/1‑69 (MQ=255)
      gaactgaacATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGAc                                                          >  1:1695749/1‑69 (MQ=255)
      gaactgaacATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGAc                                                          >  1:1575694/1‑69 (MQ=255)
      gaactgaacATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGAc                                                          >  1:645326/1‑69 (MQ=255)
              cATCCCGTAATGTCTGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCgg                                                    <  1:3039935/67‑1 (MQ=255)
              cATCCCGTAATGTCTGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCgg                                                    <  1:2058359/67‑1 (MQ=255)
              cATCCCGTAATGTCTGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCgg                                                    <  1:2027247/67‑1 (MQ=255)
              cATCCCGTAATGTCTGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCgg                                                    <  1:2854218/67‑1 (MQ=255)
              cATCCCGTAATGTCTGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCgg                                                    <  1:419817/67‑1 (MQ=255)
              cATCCCGTAATGTCTGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCgg                                                    <  1:2894927/67‑1 (MQ=255)
              cATCCCGTAATGTCTGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCgg                                                    <  1:2759708/67‑1 (MQ=255)
              cATCCCGTAATGTCTGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCgg                                                    <  1:2868622/67‑1 (MQ=255)
              cATCCCGTAATGTCTGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCgg                                                    <  1:944109/67‑1 (MQ=255)
              cATCCCGTAATGTCTGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCgg                                                    <  1:1220645/67‑1 (MQ=255)
                    gTAATGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGttt                                                 >  1:176744/1‑64 (MQ=255)
                    gTAATGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGttt                                                 >  1:1987394/1‑64 (MQ=255)
                    gTAATGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGttt                                                 >  1:1905352/1‑64 (MQ=255)
                    gTAATGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGttt                                                 >  1:427320/1‑64 (MQ=255)
                    gTAATGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGttt                                                 >  1:77616/1‑64 (MQ=255)
                    gTAATGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGttt                                                 >  1:49244/1‑64 (MQ=255)
                    gTAATGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGttt                                                 >  1:1285841/1‑64 (MQ=255)
                    gTAATGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGttt                                                 >  1:2878521/1‑64 (MQ=255)
                    gTAATGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGttt                                                 >  1:152905/1‑64 (MQ=255)
                    gTAATGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGttt                                                 >  1:1525206/1‑64 (MQ=255)
                    gTAATGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGtt                                                  >  1:1532091/1‑63 (MQ=255)
                    gTAATGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGtt                                                  >  1:1203064/1‑63 (MQ=255)
                    gTAATGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGTTTTTGag                                            >  1:859107/1‑68 (MQ=255)
                    gTAATGTCAGCGAAACGCGATAGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGttt                                                 >  1:2467483/1‑64 (MQ=255)
                      aaCGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATcc                                          >  1:1825163/1‑69 (MQ=255)
                      aaCGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATcc                                          >  1:931577/1‑69 (MQ=255)
                      aaCGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATcc                                          >  1:1500884/1‑69 (MQ=255)
                      aaCGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATcc                                          >  1:334420/1‑69 (MQ=255)
                       aCGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATcc                                          >  1:642648/1‑68 (MQ=255)
                           ctgCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCgg                                                    <  1:848727/52‑1 (MQ=255)
                            aGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTaa                                     >  1:2309371/1‑68 (MQ=255)
                            aGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTaa                                     >  1:2386915/1‑68 (MQ=255)
                            aGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTaa                                     >  1:1689968/1‑68 (MQ=255)
                            aGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTaa                                     >  1:37985/1‑68 (MQ=255)
                            aGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTaa                                     >  1:425729/1‑68 (MQ=255)
                            aGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTaa                                     >  1:1149133/1‑68 (MQ=255)
                                aaaCGCGATCGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGTTTTTGATCCGTTAAATCag                                <  1:2732936/69‑1 (MQ=255)
                                aaaCGCGATCGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGTTTTTGATCCGTTAAATCag                                <  1:2850113/69‑1 (MQ=255)
                                aaaCGCGATCGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGTTTTTGATCCGTTAAATCag                                <  1:1095508/69‑1 (MQ=255)
                                aaaCGCGATCGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGTTTTTGATCCGTTAAATCag                                <  1:2986330/69‑1 (MQ=255)
                                aaaCGCGATCGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGTTTTTGATCCGTTAAATCag                                <  1:3191267/69‑1 (MQ=255)
                                aaaCGCGATCGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGTTTTTGATCCGTTAAATCag                                <  1:636638/69‑1 (MQ=255)
                                aaaCGCGATCGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGTTTTTGATCCGTTAAATCag                                <  1:790720/69‑1 (MQ=255)
                                aaaCGCGATCGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGTTTTTGATCCGTTAAATCag                                <  1:953526/69‑1 (MQ=255)
                                aaaCGCGATCGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGTTTTTGATCCGTTAAATCag                                <  1:1391922/69‑1 (MQ=255)
                                aaaCGCGATCGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGTTTTTGATCCGTTAAATCag                                <  1:2521649/69‑1 (MQ=255)
                                aaaCGCGATCGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGTTTTTGATCCGTTAAATCag                                <  1:2108167/69‑1 (MQ=255)
                                aaaCGCGATCGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGTTTTTGATCCGTTAAATCag                                <  1:2220452/69‑1 (MQ=255)
                                aaaCGCGATCGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGTTTTTGATCCGTTAAATCag                                <  1:1312758/69‑1 (MQ=255)
                                aaaCGCGATCGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGTTTTTGATCCGTTAAATCag                                <  1:2320704/69‑1 (MQ=255)
                                 aaCGCGATCGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGTTTTTGATCCGTTAAATCag                                <  1:2728868/68‑1 (MQ=255)
                                 aaCGCGATCGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGTTTTTGATCCGTTAAATCag                                <  1:775721/68‑1 (MQ=255)
                                 aaCGCGATCGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGTTTTTGATCCGTTAAATCag                                <  1:2762838/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)
                                         cGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGt                                                   >  1:3074663/1‑41 (MQ=38)
                                         cGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGTTTTTGATCCGTTAAATCAGTGAtccgta                       >  1:3385605/1‑69 (MQ=255)
                                         cGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGTTTTTGATCCGTTAAATCAGTGAtccgta                       >  1:501320/1‑69 (MQ=255)
                                         cGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGTTTTTGATCCGTTAAATCAGTGAtccgta                       >  1:34068/1‑69 (MQ=255)
                                         cGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGTTTTTGATCCGTTAAATCAGTGAtccgta                       >  1:1573157/1‑69 (MQ=255)
                                         cGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGTTTTTGATCCGTTAAATCAGTGAtccgta                       >  1:1957572/1‑69 (MQ=255)
                                         cGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGTTTTTGATCCGTTAAATCAGTGAtccgta                       >  1:1712305/1‑69 (MQ=255)
                                         cGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGTTTTTGATCCGTTAAATCAGTGAtccgta                       >  1:1646098/1‑69 (MQ=255)
                                         cGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGTTTTTGATCCGTTAAATCAGTGAtccgta                       >  1:2263063/1‑69 (MQ=255)
                                         cGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGTTTTTGATCCGTTAAATCAGTGAtccgta                       >  1:2417322/1‑69 (MQ=255)
                                         cGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGTTTTTGATCCGTTAAATCAGTGAtccgta                       >  1:1283573/1‑69 (MQ=255)
                                         cGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGTTTTTGATCCGTTAAATCAGTGAtccgta                       >  1:2605155/1‑69 (MQ=255)
                                         cGGCAACACGAGTTTGGGTATCAGGGAAGGTGACTTCCGGTTTTTGATCCGTTAAATCAGTGAtccgta                       >  1:2668414/1‑69 (MQ=255)
                                                   aGTTTGGGTATCAGGAAAGGTGACCTCCGTTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGc             <  1:3095659/69‑1 (MQ=255)
                                                   aGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGc             <  1:740767/69‑1 (MQ=255)
                                                   aGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGc             <  1:922357/69‑1 (MQ=255)
                                                   aGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGc             <  1:120003/69‑1 (MQ=255)
                                                   aGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGc             <  1:2652395/69‑1 (MQ=255)
                                                   aGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGc             <  1:558691/69‑1 (MQ=255)
                                                   aGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGc             <  1:1015901/69‑1 (MQ=255)
                                                   aGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGc             <  1:181623/69‑1 (MQ=255)
                                                   aGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGc             <  1:2024600/69‑1 (MQ=255)
                                                   aGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGc             <  1:3331875/69‑1 (MQ=255)
                                                   aGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGc             <  1:3118903/69‑1 (MQ=255)
                                                   aGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGc             <  1:3088804/69‑1 (MQ=255)
                                                   aGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGc             <  1:2872876/69‑1 (MQ=255)
                                                   aGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGc             <  1:2323975/69‑1 (MQ=255)
                                                   aGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGc             <  1:2777570/69‑1 (MQ=255)
                                                   aGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGc             <  1:2559149/69‑1 (MQ=255)
                                                               aGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACtttc   >  1:2105606/1‑65 (MQ=255)
                                                               aGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGcc  >  1:2105191/1‑68 (MQ=255)
                                                               aGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGc   >  1:1594560/1‑67 (MQ=255)
                                                               aGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGc   >  1:1447063/1‑67 (MQ=255)
                                                               aGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGc   >  1:1355932/1‑67 (MQ=255)
                                                               aGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGc   >  1:2735116/1‑67 (MQ=255)
                                                                  |                                                                
ATAAGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCC  >  minE/584851‑584981

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: