Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,155,561 G→A 52.2% A130A (GCC→GCT astB ← succinylarginine dihydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,155,5610GA52.2% ‑2.0 / 61.1 69A130A (GCC→GCTastBsuccinylarginine dihydrolase
Reads supporting (aligned to +/- strand):  ref base G (14/19);  new base A (22/14);  total (36/33)
Fisher's exact test for biased strand distribution p-value = 1.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

AATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTA  >  minE/1155495‑1155628
                                                                  |                                                                   
aaTCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGc                                                                   <  1:1028302/69‑1 (MQ=255)
aaTCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGc                                                                   <  1:2864478/69‑1 (MQ=255)
aaTCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGc                                                                   <  1:2629808/69‑1 (MQ=255)
aaTCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGc                                                                   <  1:2489060/69‑1 (MQ=255)
aaTCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGc                                                                   <  1:2431201/69‑1 (MQ=255)
          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:536039/1‑69 (MQ=255)
                           gggCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAg                                        >  1:280575/1‑69 (MQ=255)
                           gggCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAg                                        >  1:1699092/1‑69 (MQ=255)
                           gggCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAg                                        >  1:237653/1‑69 (MQ=255)
                           gggCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAg                                        >  1:869265/1‑69 (MQ=255)
                            ggCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAg                                        >  1:847210/1‑68 (MQ=255)
                              cgcTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGt                                     >  1:1083001/1‑69 (MQ=255)
                              cgcTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGt                                     >  1:2692988/1‑69 (MQ=255)
                                          aCGGTGAAATTTATTGTTCAGGTTAGCAACAGTGAGATg                                                       >  1:1319573/1‑39 (MQ=37)
                                          aCGGTGAAATTTATTGTTCAGGTTAGCAACAGTGAGATGCACTTTGCCATCCAGc                                       >  1:927697/1‑55 (MQ=255)
                                          aCGGTGAAATTTATTGTTCAGGTTAGCAACAGTGAGATGCACTTTGCCATCCAGc                                       >  1:1740743/1‑55 (MQ=255)
                                          aCGGTGAAATTTATTGTTCAGGTTAGCAACAGTGAGATGCACTTTGCCATCCAGc                                       >  1:3123713/1‑55 (MQ=255)
                                          aCGGTGAAATTTATTGTTCAGGTTAGCAACAGTGAGATGCACTTTGCCATCCAGc                                       >  1:3306962/1‑55 (MQ=255)
                                          aCGGTGAAATTTATTGTTCAGGTTAGCAACAGTGAGATGCACTTTGCCATCCAGc                                       >  1:1295548/1‑55 (MQ=255)
                                          aCGGTGAAATTTATTGTTCAGGTTAGCAACAGTGAGATGCACTTTGCCATCCAGc                                       >  1:1741090/1‑55 (MQ=255)
                                          aCGGTGAAATTTATTGTTCAGGTTAGCAACAGTGAGATGCACTTTGCCATCCAGc                                       >  1:2789781/1‑55 (MQ=255)
                                          aCGGTGAAATTTATTGTTCAGGTTAGCAACAGTGAGATGCACTTTGCCATCCAGc                                       >  1:2668898/1‑55 (MQ=255)
                                          aCGGTGAAATTTATTGTTCAGGTTAGCAACAGTGAGATGCACTTTGCCATCCAGc                                       >  1:53117/1‑55 (MQ=255)
                                          aCGGTGAAATTTATTGTTCAGGTTAGCAACAGTGAGATGCACTTTGCCATCCAGc                                       >  1:263734/1‑55 (MQ=255)
                                          aCGGTGAAATTTATTGTTCAGGTTAGCAACAGTGAGATGCACTTTGCCATCCAGc                                       >  1:620687/1‑55 (MQ=255)
                                          aCGGTGAAATTTATTGTTCAGGTTAGCAACAGTGAGATGCACTTTGCCATCCAGc                                       >  1:2139938/1‑55 (MQ=255)
                                          aCGGTGAAATTTATTGTTCAGGTTAGCAACAGTGAGATGCACTTTGCCATCCAGc                                       >  1:85290/1‑55 (MQ=255)
                                          aCGGTGAAATTTATTGTTCAGGTTAGCAACAGTGAGATGCACTTTGCCATCCAGc                                       >  1:1885805/1‑55 (MQ=255)
                                          aCGGTGAAATTTATTGTTCAGGTTAGCAACAGTGAGATGCACTTTGCCATCCAGc                                       >  1:753893/1‑55 (MQ=255)
                                          aCGGTGAAATTTATTGTTCAGGTTAGCAACAGTGAGATGCACTTTGCCATCCAGc                                       >  1:2459009/1‑55 (MQ=255)
                                          aCGGTGAAATTTATTGTTCAGGTTAGCAACAGTGAGATGCACTTTGCCATCCAGc                                       >  1:2461741/1‑55 (MQ=255)
                                          aCGGTGAAATTTATTGTTCAGGTTAGCAACAGTGAGATGCACTTTGCCATCCAGc                                       >  1:1562547/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:510346/1‑69 (MQ=255)
                                              tGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGcgcg                     >  1:2681862/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:2761874/69‑1 (MQ=255)
                                                    ttATTGTTCAGGTTAGCAACAGTGAGATGCACTTTGCCATCCAGCGTATCGGCGGATGGCGCGATCGtt               <  1:66746/69‑1 (MQ=255)
                                                    ttATTGTTCAGGTTAGCAACAGTGAGATGCACTTTGCCATCCAGCGTATCGGCGGATGGCGCGATCGtt               <  1:3244358/69‑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:2256453/69‑1 (MQ=255)
                                                    ttATTGTTCAGGTTAGCAACAGTGAGATGCACTTTGCCATCCAGCGTATCGGCGGATGGCGCGATCGtt               <  1:2313554/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:266340/69‑1 (MQ=255)
                                                    ttATTGTTCAGGTTAGCAACAGTGAGATGCACTTTGCCATCCAGCGTATCGGCGGATGGCGCGATCGtt               <  1:288434/69‑1 (MQ=255)
                                                         gTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGc          <  1:256481/69‑1 (MQ=255)
                                                                  ggCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGTCAGATGGCGCGATCGTTGCCGCATTGGCt   >  1:2675510/1‑67 (MQ=255)
                                                                  ggCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTa  >  1:1620976/1‑68 (MQ=255)
                                                                  ggCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTa  >  1:2874055/1‑68 (MQ=255)
                                                                  ggCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTa  >  1:3124965/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)
                                                                  |                                                                   
AATCGCTTTTAACAGCGATTCAGTGACGGGCGCTTCCAGCGAACGGTGAAATTTATTGTTCAGGTTGGCAACGGTGAGATGCACTTTGCCATCCAGCGTATCGGCAGATGGCGCGATCGTTGCCGCATTGGCTA  >  minE/1155495‑1155628

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: