Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,144,345 G→A 45.0% D206N (GAT→AAT)  argG → argininosuccinate synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,144,3450GA45.0% 35.4 / 123.6 91D206N (GAT→AAT) argGargininosuccinate synthetase
Reads supporting (aligned to +/- strand):  ref base G (20/30);  new base A (17/24);  total (37/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATGTCTGTCGAAAAAGCCTACTCCACAGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGAT  >  minE/2144281‑2144410
                                                                |                                                                 
aaTGTCTGTCGAAAAAGCCTACTCCACAGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGAATCTgg                                                             <  1:2671956/71‑1 (MQ=255)
    tctgtcGAAAAAGCCTACTCCACAGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGGATCTGGAATa                                                         <  1:3617213/71‑1 (MQ=255)
        tcGAAAAAGCCTACTCCACAGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGGATCTgg                                                             <  1:2603356/63‑1 (MQ=255)
        tcGAAAAAGCCTACTCCACAGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGAATCTgg                                                             <  1:2975273/63‑1 (MQ=255)
        tcGAAAAAGCCTACTCCACAGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGAATCTgg                                                             <  1:459278/63‑1 (MQ=255)
        tcGAAAAAGCCTACTCCACAGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGAATCTgg                                                             <  1:3809518/63‑1 (MQ=255)
          gAAAAAGCCTACTCCACAGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGGATCTGGaa                                                           >  1:2020302/1‑63 (MQ=255)
           aaaaaGCCTACTCCACAGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGGATc                                                                <  1:211703/57‑1 (MQ=255)
           aaaaaGCCTACTCCACAGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGGATc                                                                <  1:3759460/57‑1 (MQ=255)
                gCCTACTCCACAGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGAATCTgg                                                             <  1:922399/55‑1 (MQ=255)
                   tACTCCACAGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGGATCTGGAATACCTCAACTCCa                                              <  1:3520407/67‑1 (MQ=255)
                   tACTCCACAGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGGATCTGGAATACCTCAACTCCa                                              <  1:539660/67‑1 (MQ=255)
                   tACTCCACAGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGGATCTGGAATACCTCAACTCCa                                              <  1:158298/67‑1 (MQ=255)
                      tCCACAGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGAATCTGGAATACCTCAACTCCAGCGTCa                                        >  1:4065092/1‑70 (MQ=255)
                      tCCACAGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGAATCTGGAATACCTCAACTCCAGCGTCa                                        >  1:2349126/1‑70 (MQ=255)
                      tCCACAGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGAATCTGGAATACCTCAACTCCAGCGTCa                                        >  1:3036346/1‑70 (MQ=255)
                        cacaGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGGa                                                                  <  1:2869574/42‑1 (MQ=255)
                        cacaGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCaaa                                      <  1:424499/70‑1 (MQ=255)
                        cacaGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCaaa                                      <  1:1210859/70‑1 (MQ=255)
                        cacaGACTACAACATGCTTGGTGCAACGCATGAAGCGAAGGATCTGGaa                                                           >  1:2342585/1‑49 (MQ=255)
                         acaGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGGATCTGGAATa                                                         >  1:1464837/1‑50 (MQ=255)
                         acaGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGGATCTGGAATAc                                                        >  1:1815530/1‑51 (MQ=255)
                         acaGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGAATCTGGAATACCTCAACTCCAGCGTCAAAAt                                    >  1:3469066/1‑71 (MQ=255)
                          caGACTCCAACATGCTTGGTGCAACGCATGAAGCTAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATc                                   <  1:610878/71‑1 (MQ=255)
                          caGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATc                                   <  1:3779486/71‑1 (MQ=255)
                             aCTCCAACATGCTTGGTGCAACGCATGAAGCGAAGAATCTGGAATACCTCAACTCCAGCGTCaaa                                      <  1:116041/65‑1 (MQ=255)
                               tCCAACATGCTTGGTGCAACGCATGAAGCGAAGGATCTgg                                                             >  1:2699081/1‑40 (MQ=255)
                               tCCAACATGCTTGGTGCAACGCATGAAGCGAAGGATCTgg                                                             >  1:2643810/1‑40 (MQ=255)
                               tCCAACATGCTTGGTGCAACGCATGAAGCGAAGGATCTgg                                                             >  1:980873/1‑40 (MQ=255)
                               tCCAACATGCTTGGTGCAACGCATGAAGCGAAGGATCTgg                                                             >  1:1177503/1‑40 (MQ=255)
                               tCCAACATGCTTGGTGCAACGCATGAAGCGAAGGATCTgg                                                             >  1:1279475/1‑40 (MQ=255)
                               tCCAACATGCTTGGTGCAACGCATGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATc                                   <  1:843391/66‑1 (MQ=255)
                               tCCAACATGCTTGGTGCAACGCATGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATc                                   <  1:2567271/66‑1 (MQ=255)
                               tCCAACATGCTTGGTGCAACGCATGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATc                                   <  1:160864/66‑1 (MQ=255)
                               tCCAACATGCTTGGTGCAACGCATGAAGCGAAGAATCTGGAATACCTCAACTCCAGCGTCaaaa                                     >  1:2524624/1‑64 (MQ=255)
                               tCCAACATGCTTGGTGCAACGCATGAAGCGAAGAATCTGGAATACCTCAACTCCAGCGTCaaaa                                     >  1:2540427/1‑64 (MQ=255)
                               tCCAACATGCTTGGTGCAACGCATGAAGCGAAGAATCTGGAATACCTCAACTCCAGCGTCaaaa                                     >  1:2935372/1‑64 (MQ=255)
                               tCCAACATGCTTGGTGCAACGCATGAAGCGAAGAATCTGGAATACCTCAACTCCAGCGTCAAAAt                                    >  1:1109877/1‑65 (MQ=255)
                                ccAACATGCTTGGTGCAACGCATGAAGCGAAGAATCTgg                                                             <  1:140316/39‑1 (MQ=255)
                                ccAACATGCTTGGTGCAACGCATGAAGCGAAGAATCTgg                                                             <  1:3730138/39‑1 (MQ=255)
                                ccAACATGCTTGGTGCAACGCATGAAGCGAAGAATCTgg                                                             <  1:358343/39‑1 (MQ=255)
                                ccAACATGCTTGGTGCAACGCATGAAGCGAAGAATCTgg                                                             <  1:2759927/39‑1 (MQ=255)
                                 cAACATGCTTGGTGCAACGCATGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAAcc                            >  1:123452/1‑71 (MQ=255)
                                  aaCATGCTTGGTGCAACGCATGAAGCGAAGAATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAAccc                           <  1:3319367/71‑1 (MQ=255)
                                   aCATGCTTGGTGCAACGCATGAAGCGAAGAATCTGGAATACCTCaa                                                   >  1:2670236/1‑46 (MQ=255)
                                     aTGCTTGGTGCAACGCATGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGAt                        <  1:505541/71‑1 (MQ=255)
                                     aTGCTTGGTGCAACGCATGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGAt                        <  1:1616168/71‑1 (MQ=255)
                                      tGCTTGGTGCAACGCATGAAGCGAAGAATCTGGAATa                                                         >  1:413737/1‑37 (MQ=255)
                                      tGCTTGGTGCAACGCATGAAGCGAAGAATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGAtt                       <  1:3639102/71‑1 (MQ=255)
                                        cTTGGTGCAACGCATGAAGCGAAGAATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTAt                     <  1:3135097/71‑1 (MQ=255)
                                         ttGGTGCAACGCATGAAGCGAAGAATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATg                    >  1:1360774/1‑71 (MQ=255)
                                         ttGGTGAAACGCATGAAGCGAAGAATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATg                    >  1:2854659/1‑71 (MQ=255)
                                          tGGTGCAACGCATGAAGCGAAGAATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATgg                   >  1:3176158/1‑71 (MQ=255)
                                          tGGTGCAACGCATGAAGCGAAGAATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATgg                   >  1:3248112/1‑71 (MQ=255)
                                           ggTGCAACGCATGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCg                                  <  1:3281305/55‑1 (MQ=255)
                                           ggTGCAACGCATGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCg                                  <  1:3670049/55‑1 (MQ=255)
                                           ggTGCAACGCATGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCg                                  <  1:2208329/55‑1 (MQ=255)
                                           ggTGCAACGCATGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCg                                  <  1:2651325/55‑1 (MQ=255)
                                           ggTGCAACGCATGAAGCGAAGGATCTAGAATACCTCAACTCCAGCGTCAAAATCg                                  <  1:811956/55‑1 (MQ=255)
                                           ggTGCAACGCATGAAGCGAAGAATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGAtt                       >  1:3633810/1‑66 (MQ=255)
                                           ggTGCAACGCATGAAGCGAAGAATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATgg                   <  1:253528/70‑1 (MQ=255)
                                           ggTGCAACGCATGAAGCGAAGAATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATgg                   <  1:2348378/70‑1 (MQ=255)
                                           gatGCAACGCATGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCg                                  <  1:752698/53‑1 (MQ=255)
                                            gTGCAACGCATGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGc                 <  1:1472408/71‑1 (MQ=255)
                                            gTGCAACGCATGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGc                 <  1:1544375/71‑1 (MQ=255)
                                             tGCAACGCATGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTAt                     <  1:4146528/66‑1 (MQ=255)
                                             tGCAACGCATGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTAt                     <  1:3655767/66‑1 (MQ=255)
                                             tGCAACGCATGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTAt                     <  1:1368956/66‑1 (MQ=255)
                                             tGCAACGCATGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTAt                     <  1:3791085/66‑1 (MQ=255)
                                             tGCAACGCATGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTAt                     <  1:387807/66‑1 (MQ=255)
                                              gCAACGCATGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGt               >  1:862464/1‑71 (MQ=255)
                                              gCAACGCATGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGt               >  1:4188067/1‑71 (MQ=255)
                                                 aCGCATGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTg              >  1:1368665/1‑69 (MQ=255)
                                                 aCGCATGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTGa             >  1:293130/1‑70 (MQ=255)
                                                    cATGAAGCGAAGAATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAAtt         <  1:2684475/71‑1 (MQ=255)
                                                    cATGAAGCGAAGAATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAAtt         <  1:208133/71‑1 (MQ=255)
                                                    cATGAAGCGAAGAATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAAtt         <  1:2247840/71‑1 (MQ=255)
                                                    cATGAAGCGAAGAATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAAtt         <  1:622440/71‑1 (MQ=255)
                                                    cATGAAGCGAAGAATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAAtt         <  1:29189/71‑1 (MQ=255)
                                                      tGAAGCGAAGAATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCt       <  1:1926860/71‑1 (MQ=255)
                                                      tGAAGCGAAGAATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCt       <  1:4166984/71‑1 (MQ=255)
                                                        aaGCGAAGAATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAAccc                           <  1:2299597/49‑1 (MQ=255)
                                                        aaGCGAAGAATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAAccc                           <  1:1448471/49‑1 (MQ=255)
                                                          gCGAAGAATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGa                         >  1:133168/1‑49 (MQ=255)
                                                           cGAAGAATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGAt  >  1:3157783/1‑71 (MQ=255)
                                                            gaagaaTCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGAt                        <  1:2124468/48‑1 (MQ=255)
                                                            gaagaaTCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGAt                        <  1:495088/48‑1 (MQ=255)
                                                            gaagaaTCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGAt                        <  1:2962342/48‑1 (MQ=255)
                                                            gaagaaTCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGAt                        <  1:2698116/48‑1 (MQ=255)
                                                            gaagaaTCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGAt                        <  1:1159370/48‑1 (MQ=255)
                                                            gaagaaTCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGAt                        <  1:1836574/48‑1 (MQ=255)
                                                             aaGGATCTGGAATACCTCAACTCCAGCGTCAAAAt                                    >  1:1265812/1‑35 (MQ=255)
                                                             aaGGATCTGGAATACCTCAACTCCAGCGTCAAAAt                                    >  1:434778/1‑35 (MQ=255)
                                                             aaGGATCTGGAATACCTCAACTCCAGCGTCAAAAt                                    >  1:3647628/1‑35 (MQ=255)
                                                             aaGGATCTGGAATACCTCAACTCCAGCGTCAAAAt                                    >  1:3642615/1‑35 (MQ=255)
                                                             aaGGATCTGGAATACCTCAACTCCAGCGTCAAAAt                                    >  1:1167476/1‑35 (MQ=255)
                                                             aaGGATCTGGAATACCTCAACTCCAGCGTCAAAAt                                    >  1:2916938/1‑35 (MQ=255)
                                                                |                                                                 
AATGTCTGTCGAAAAAGCCTACTCCACAGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGAT  >  minE/2144281‑2144410

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: