Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,155,567 G→A 48.6% T128T (ACC→ACT astB ← succinylarginine dihydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,155,5670GA48.6% 1.5 / 64.5 74T128T (ACC→ACTastBsuccinylarginine dihydrolase
Reads supporting (aligned to +/- strand):  ref base G (21/17);  new base A (22/14);  total (43/31)
Fisher's exact test for biased strand distribution p-value = 6.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACAT  >  minE/1155505‑1155635
                                                              |                                                                    
aaCAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTga                                                                  <  1:976745/67‑1 (MQ=255)
aaCAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTga                                                                  <  1:1090037/67‑1 (MQ=255)
aaCAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTga                                                                  <  1:1094052/67‑1 (MQ=255)
aaCAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTga                                                                  <  1:287788/67‑1 (MQ=255)
aaCAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTga                                                                  <  1:287675/67‑1 (MQ=255)
aaCAGCGATTCAGTGACGGGCGCTTCCAGCGAAAGGTGAAATTTATTGTTCAGGTTGGCAACGGTga                                                                  <  1:3139770/67‑1 (MQ=255)
   aGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATg                                                              <  1:1002439/68‑1 (MQ=255)
   aGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATg                                                              <  1:170202/68‑1 (MQ=255)
             tGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCa                                                    <  1:100318/68‑1 (MQ=255)
             tGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCa                                                    <  1:3377223/68‑1 (MQ=255)
             tGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCa                                                    <  1:1085094/68‑1 (MQ=255)
                 gggCGCTTTCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAg                                               >  1:2151030/1‑69 (MQ=255)
                 gggCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAg                                               >  1:237653/1‑69 (MQ=255)
                 gggCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAg                                               >  1:1699092/1‑69 (MQ=255)
                 gggCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAg                                               >  1:280575/1‑69 (MQ=255)
                 gggCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAg                                               >  1:536039/1‑69 (MQ=255)
                 gggCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAg                                               >  1:869265/1‑69 (MQ=255)
                  ggCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAg                                               >  1:847210/1‑68 (MQ=255)
                    cgcTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGt                                            >  1:2692988/1‑69 (MQ=255)
                    cgcTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGt                                            >  1:1083001/1‑69 (MQ=255)
                                aCGGTGAAATTTATTGTTCAGGTTAGCAACAGTGAGATg                                                              >  1:1319573/1‑39 (MQ=37)
                                aCGGTGAAATTTATTGTTCAGGTTAGCAACAGTGAGATGCACTTTGCCATCCAGc                                              >  1:2461741/1‑55 (MQ=255)
                                aCGGTGAAATTTATTGTTCAGGTTAGCAACAGTGAGATGCACTTTGCCATCCAGc                                              >  1:2459009/1‑55 (MQ=255)
                                aCGGTGAAATTTATTGTTCAGGTTAGCAACAGTGAGATGCACTTTGCCATCCAGc                                              >  1:85290/1‑55 (MQ=255)
                                aCGGTGAAATTTATTGTTCAGGTTAGCAACAGTGAGATGCACTTTGCCATCCAGc                                              >  1:1295548/1‑55 (MQ=255)
                                aCGGTGAAATTTATTGTTCAGGTTAGCAACAGTGAGATGCACTTTGCCATCCAGc                                              >  1:753893/1‑55 (MQ=255)
                                aCGGTGAAATTTATTGTTCAGGTTAGCAACAGTGAGATGCACTTTGCCATCCAGc                                              >  1:263734/1‑55 (MQ=255)
                                aCGGTGAAATTTATTGTTCAGGTTAGCAACAGTGAGATGCACTTTGCCATCCAGc                                              >  1:2789781/1‑55 (MQ=255)
                                aCGGTGAAATTTATTGTTCAGGTTAGCAACAGTGAGATGCACTTTGCCATCCAGc                                              >  1:2668898/1‑55 (MQ=255)
                                aCGGTGAAATTTATTGTTCAGGTTAGCAACAGTGAGATGCACTTTGCCATCCAGc                                              >  1:3123713/1‑55 (MQ=255)
                                aCGGTGAAATTTATTGTTCAGGTTAGCAACAGTGAGATGCACTTTGCCATCCAGc                                              >  1:3306962/1‑55 (MQ=255)
                                aCGGTGAAATTTATTGTTCAGGTTAGCAACAGTGAGATGCACTTTGCCATCCAGc                                              >  1:2139938/1‑55 (MQ=255)
                                aCGGTGAAATTTATTGTTCAGGTTAGCAACAGTGAGATGCACTTTGCCATCCAGc                                              >  1:53117/1‑55 (MQ=255)
                                aCGGTGAAATTTATTGTTCAGGTTAGCAACAGTGAGATGCACTTTGCCATCCAGc                                              >  1:1885805/1‑55 (MQ=255)
                                aCGGTGAAATTTATTGTTCAGGTTAGCAACAGTGAGATGCACTTTGCCATCCAGc                                              >  1:1741090/1‑55 (MQ=255)
                                aCGGTGAAATTTATTGTTCAGGTTAGCAACAGTGAGATGCACTTTGCCATCCAGc                                              >  1:1740743/1‑55 (MQ=255)
                                aCGGTGAAATTTATTGTTCAGGTTAGCAACAGTGAGATGCACTTTGCCATCCAGc                                              >  1:927697/1‑55 (MQ=255)
                                aCGGTGAAATTTATTGTTCAGGTTAGCAACAGTGAGATGCACTTTGCCATCCAGc                                              >  1:1562547/1‑55 (MQ=255)
                                aCGGTGAAATTTATTGTTCAGGTTAGCAACAGTGAGATGCACTTTGCCATCCAGc                                              >  1:620687/1‑55 (MQ=255)
                                aCGGTGAAATTTATTGTTCAGGTTAGCAACAGTGAGATGCACTTTGCCATCAAGc                                              >  1:865424/1‑55 (MQ=255)
                                aCGGTGAAATTTATTGTTCAGGTTAGCAACAGTGAGATGCACTTTGCCATACAGc                                              >  1:1996615/1‑55 (MQ=255)
                                aCGGAGAAATTTATTGTTCAGGTAAGCAACAGTgag                                                                 >  1:1438659/1‑36 (MQ=19)
                                    tGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTGTCGGCAGATGGcgcg                            <  1:3044006/69‑1 (MQ=255)
                                    tGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATgg                                >  1:1799290/1‑65 (MQ=255)
                                    tGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGcgcg                            >  1:2681862/1‑69 (MQ=255)
                                    tGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGcgcg                            >  1:510346/1‑69 (MQ=255)
                                    tGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGcgcg                            >  1:1411669/1‑69 (MQ=255)
                                    tGAAATATATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGcgcg                            >  1:2638643/1‑69 (MQ=255)
                                     gAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGcgcg                            <  1:539100/68‑1 (MQ=255)
                                          ttATTGTTCAGGTTAGCAACAGTGAGATGCACTTTGCCATCCAGCGTATCGGCGGATGGCGCGATCGtt                      <  1:3082433/69‑1 (MQ=255)
                                          ttATTGTTCAGGTTAGCAACAGTGAGATGCACTTTGCCATCCAGCGTATCGGCGGATGGCGCGATCGtt                      <  1:2014640/69‑1 (MQ=255)
                                          ttATTGTTCAGGTTAGCAACAGTGAGATGCACTTTGCCATCCAGCGTATCGGCGGATGGCGCGATCGtt                      <  1:2168274/69‑1 (MQ=255)
                                          ttATTGTTCAGGTTAGCAACAGTGAGATGCACTTTGCCATCCAGCGTATCGGCGGATGGCGCGATCGtt                      <  1:2241927/69‑1 (MQ=255)
                                          ttATTGTTCAGGTTAGCAACAGTGAGATGCACTTTGCCATCCAGCGTATCGGCGGATGGCGCGATCGtt                      <  1:3244358/69‑1 (MQ=255)
                                          ttATTGTTCAGGTTAGCAACAGTGAGATGCACTTTGCCATCCAGCGTATCGGCGGATGGCGCGATCGtt                      <  1:2256453/69‑1 (MQ=255)
                                          ttATTGTTCAGGTTAGCAACAGTGAGATGCACTTTGCCATCCAGCGTATCGGCGGATGGCGCGATCGtt                      <  1:66746/69‑1 (MQ=255)
                                          ttATTGTTCAGGTTAGCAACAGTGAGATGCACTTTGCCATCCAGCGTATCGGCGGATGGCGCGATCGtt                      <  1:2313554/69‑1 (MQ=255)
                                          ttATTGTTCAGGTTAGCAACAGTGAGATGCACTTTGCCATCCAGCGTATCGGCGGATGGCGCGATCGtt                      <  1:288434/69‑1 (MQ=255)
                                          ttATTGTTCAGGTTAGCAACAGTGAGATGCACTTTGCCATCCAGCGTATCGGCGGATGGCGCGATCGtt                      <  1:2395137/69‑1 (MQ=255)
                                          ttATTGTTCAGGTTAGCAACAGTGAGATGCACTTTGCCATCCAGCGTATCGGCGGATGGCGCGATCGtt                      <  1:2423392/69‑1 (MQ=255)
                                          ttATTGTTCAGGTTAGCAACAGTGAGATGCACTTTGCCATCCAGCGTATCGGCGGATGGCGCGATCGtt                      <  1:2535482/69‑1 (MQ=255)
                                          ttATTGTTCAGGTTAGCAACAGTGAGATGCACTTTGCCATCCAGCGTATCGGCGGATGGCGCGATCGtt                      <  1:2761874/69‑1 (MQ=255)
                                          ttATTGTTCAGGTTAGCAACAGTGAGATGCACTTTGCCATCCAGCGTATCGGCGGATGGCGCGATCGtt                      <  1:266340/69‑1 (MQ=255)
                                               gTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGc                 <  1:256481/69‑1 (MQ=255)
                                                        ggCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGTCAGATGGCGCGATCGTTGCCGCATTGGCt          >  1:2675510/1‑67 (MQ=255)
                                                        ggCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTa         >  1:3124965/1‑68 (MQ=255)
                                                        ggCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTa         >  1:2874055/1‑68 (MQ=255)
                                                        ggCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTa         >  1:1620976/1‑68 (MQ=255)
                                                        ggCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTa         >  1:945170/1‑68 (MQ=255)
                                                        ggCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTa         >  1:1014584/1‑68 (MQ=255)
                                                        ggCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGAGATCGTTGCCGCATTGGCTa         >  1:930720/1‑68 (MQ=255)
                                                           aaCGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCca     <  1:1169052/69‑1 (MQ=255)
                                                           aaCGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCca     <  1:2517697/69‑1 (MQ=255)
                                                           aaCGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCca     <  1:2135023/69‑1 (MQ=255)
                                                              ggTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACAt  >  1:2106191/1‑69 (MQ=255)
                                                              |                                                                    
AACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTACCCACAT  >  minE/1155505‑1155635

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: