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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,144,3450GA13.3% 118.5 / 9.1 45D206N (GAT→AAT) argGargininosuccinate synthetase
Reads supporting (aligned to +/- strand):  ref base G (17/21);  new base A (3/3);  total (21/24)
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

CTACAAAATGTCTGTCGAAAAAGCCTACTCCACAGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGAT  >  minE/2144275‑2144410
                                                                      |                                                                 
cTACAAAATGTCTGTCGAAAAAGCCTACTCCACAGACTCCAACATGCTTGGTGCAACGCATGAAGCTAAgg                                                                   <  1:1582816/71‑1 (MQ=255)
      aaTGTCTGTCGAAAAAGCCTACTCCACAGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGGATCTg                                                              >  1:1546749/1‑70 (MQ=255)
       aTGTCTGTCGAAAAAGCCTACTCCACAGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGGAt                                                                 >  1:944713/1‑66 (MQ=255)
              tcGAAAAAGCCTACTCCACAGCCTCCAACATGCTTGGTGCAACGCATGAAGCGAAGAATCTGGAATACCTc                                                     <  1:2184836/71‑1 (MQ=255)
              tcGAAAAAGCCTACTCCACAGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGGATCTgg                                                             <  1:289380/63‑1 (MQ=255)
                    aaGCCTACTCCACAGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGGATCTGGAATACCTCAACTcc                                               <  1:1648376/71‑1 (MQ=255)
                              cacaGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGGATCTGGaa                                                           >  1:1228614/1‑49 (MQ=255)
                              cacaGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGGATCTGGaa                                                           >  1:618925/1‑49 (MQ=255)
                              cacaGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGAATCTGGAATACCTCAACTCCAGCGTCaaaa                                     >  1:949909/1‑71 (MQ=255)
                              cacaGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGAATCTGGAATACCTCAACTCCAGCGTCaaa                                      <  1:1091290/70‑1 (MQ=255)
                                      ccAACATGCTTGGTGCAACGCATGAAGCGAAGAATCTGGAATACCTc                                                     <  1:2185299/47‑1 (MQ=255)
                                       cAACATGCTTGGTGCAACGCATGAAGCGAAGGATCTGGAATACCTCAACTCCa                                              >  1:1424315/1‑53 (MQ=255)
                                       cAACATGCTTGGTGCAACGCATGAAGCGAAGGATCTGGAATACCTCAACTCCa                                              >  1:1453132/1‑53 (MQ=255)
                                       cAACATGCTTGGTGCAACGCATGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTc                                >  1:280699/1‑67 (MQ=255)
                                        aaCATGCTTGGTGCAACGCATGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAAccc                           >  1:498122/1‑71 (MQ=255)
                                        aaCATGCTTGGTGCAACGCATGAAGCGAAGGATCAGGAATAc                                                        >  1:1287977/1‑42 (MQ=255)
                                         aCATGCTTGGTGCAACGCATGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAAccc                           >  1:742278/1‑70 (MQ=255)
                                          cATGCTTGGTGCAACGCATGAAGCGAAGAATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGa                         >  1:1097019/1‑71 (MQ=255)
                                           aTGCTTGGTGCAACGCATGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGAt                        <  1:2125537/71‑1 (MQ=255)
                                           aTGCTTGGTGCAACGCATGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGAt                        >  1:752234/1‑71 (MQ=255)
                                           aTGCTTGGTGCAACGCATGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGAt                        <  1:253271/71‑1 (MQ=255)
                                            tGCTTGGTGCAACGCATGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGAtt                       <  1:794976/71‑1 (MQ=255)
                                             gCTTGGTGCAACGCATGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTa                      >  1:616734/1‑71 (MQ=255)
                                              cTTGGTGCAACGCATGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTa                      <  1:1621110/70‑1 (MQ=255)
                                              cTTGGTGCAACGCATGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTa                      <  1:723416/70‑1 (MQ=255)
                                              cTTGGTGCAACGCATGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTa                      <  1:1099505/70‑1 (MQ=255)
                                              cTTGGTGCAACGCATGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTa                      <  1:52543/70‑1 (MQ=255)
                                              cTTGGTGCAACGCATGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTAt                     >  1:2138256/1‑71 (MQ=255)
                                                    gCAACGCATGAAGCGAAGTATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGt               >  1:754317/1‑71 (MQ=255)
                                                     cAACGCATGAAGCGAAGGATCTGGAATACCTCAACTCCAGGGTCAAAATCGTCAACCCGATTATGGGCGt               >  1:2080030/1‑70 (MQ=255)
                                                      aaCGCATGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTGa             <  1:2180675/71‑1 (MQ=255)
                                                      aaCGCATGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTGa             <  1:1383694/71‑1 (MQ=255)
                                                       aCGCATGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTGa             <  1:1014810/70‑1 (MQ=255)
                                                       aCGCATGAAGCGAAGAATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTGa             >  1:567252/1‑70 (MQ=255)
                                                        cGCATTAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCaaa                                      <  1:1066384/44‑1 (MQ=255)
                                                        cGCATGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCaaa                                      <  1:2211165/44‑1 (MQ=255)
                                                        cGCATGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTGaaa           >  1:1241471/1‑71 (MQ=255)
                                                          cATGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAAtt         <  1:767502/71‑1 (MQ=255)
                                                              aaGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGTGCGTGAAATTCTgg     >  1:1728584/1‑71 (MQ=255)
                                                              aaGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTgg     >  1:635623/1‑71 (MQ=255)
                                                              aaGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTgg     <  1:2133270/71‑1 (MQ=255)
                                                              aaGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTgg     <  1:1057634/71‑1 (MQ=255)
                                                              aaGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTg      <  1:482547/70‑1 (MQ=255)
                                                                 cGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGAt  <  1:1054398/71‑1 (MQ=255)
                                                                   aaGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTgg     <  1:726541/66‑1 (MQ=255)
                                                                   aaGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTgg     <  1:98755/66‑1 (MQ=255)
                                                                     ggATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCt       >  1:968904/1‑62 (MQ=255)
                                                                      |                                                                 
CTACAAAATGTCTGTCGAAAAAGCCTACTCCACAGACTCCAACATGCTTGGTGCAACGCATGAAGCGAAGGATCTGGAATACCTCAACTCCAGCGTCAAAATCGTCAACCCGATTATGGGCGTGAAATTCTGGGAT  >  minE/2144275‑2144410

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