Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,495,891 T→C 7.5% L107P (CTG→CCG)  mdeA → methionine gamma‑lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,495,8910TC 7.5% 145.8 / 2.4 53L107P (CTG→CCG) mdeAmethionine gamma‑lyase
Reads supporting (aligned to +/- strand):  ref base T (24/25);  new base C (2/2);  total (26/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.60e-01

CCGCATTTCCAACCCCACGCTGGCCTTGCTGGAGCAGCGCATGGCTTCGCTGGAGGGCGGTGAAGCGGGGTTGGCGTTGGCATCGGGGATGGGCGCAATTACCTCAACCCTCTGGACCCTGCT‑CAGGCCTGGGGATGAGCTGATTGTCGGGCGCACCCTGTATGGCTGTACCTTCGCTTTTCTGCACCACGGTATTGGTGAGTTCGGGGTGAAAGTTCGCCATGTCGATCTCAACGATGCCAAGGCCCTGAAAGCGGCGATCAGCAGCAAAACACGGATGATCTACTT  >  NC_002947/1495751‑1496038
                                                                                                                                             |                                                                                                                                                   
ccGCATTTCCAACCCCACGCTGGCCTTGCTGGAGCAGCGCATGGCTTCGCTGGAGGGCGGTGAAGCGGGGTTGGCGTTGGCATCGGGGATGGGCGCAATTACCTCAACCCTCTGGACCCTGCT‑CAGGCCTGGGGATGCGCTGATTGTCg                                                                                                                                             >  4:312289/1‑149 (MQ=255)
ccGCATTTCCAACCCCACGCTGGCCTTGCTGGAGCAGCGCATGGCTTCGCTGGAGGGCGGTGAAGCGGGGTTGGCGTTGGCATCGGGGATGGGCGCAATTACCTCAACCCTCTGGACCCTGCT‑CAGGCCTGGGGATGAGCTGATTGTCg                                                                                                                                             <  3:160920/149‑1 (MQ=255)
ccGCATTTCCAACCCCACGCTGGCCTTGCTGGAGCAGCGCATGGCTTCGCTGGAGGGCGGTGAAGCGGGGTTGGCGTTGGCATCGGGGATGGGCGCAATTACCTCAACCCTCTGGACCCTGCT‑CAGGCCTGGGGATGAGCTGATTGTCg                                                                                                                                             >  4:160920/1‑149 (MQ=255)
         cAACCCCACGCTGGCATTGCTGGAGCAGCGCATCGCTTCGCTGGAGGGCGGTGAAGCGGGGTTGGCTTTGGCATCGGGGATGGGCGCAATTACCTCAACCCTCTGGACCCTGCT‑CAGGCCTGGGGATGAGCTGATTGTCGGGCGCAccc                                                                                                                                    >  5:110924/1‑149 (MQ=255)
             cccACGCTGGCCTTGCTGGAGCAGCCCATGGCTTCGCTGGAGGGCGGTGAAGCGGGGTTGGCGTAGGCATCGGGGATGGGCGCAATTACCTCAACCCTCTGGACCCTGCT‑CAGGCCTGGGGATGAGCTGATTGTCGGGCGCACCCTGTa                                                                                                                                >  3:72546/1‑149 (MQ=255)
                  gCTGGCCTTGCTGGAGCAGCGCATGGCTTCGCTGGAGGGCGGTGAAGCGGGGTTGGCGTTGGCATCGGGGATGGGCGCAATTACCTCAACCCTCTGGACCCTGCT‑CAGGCCTGGGGATGAGCTGATTGTCGGGCGCACCCt                                                                                                                                   <  2:356168/141‑1 (MQ=255)
                  gCTGGCCTTGCTGGAGCAGCGCATGGCTTCGCTGGAGGGCGGTGAAGCGGGGTTGGCGTTGGCATCGGGGATGGGCGCAATTACCTCAACCCTCTCGACCCTGCT‑CAGGCCTGGGGATGAGCTGATTGTCGGGCGCACCCt                                                                                                                                   >  1:356168/1‑141 (MQ=255)
                      gCCTTGCTGGAGCAGCGCATGGCTTCGCTGGAGGGCGGTGAAGCGGGGTTGGCGTTGGCATCGGGGATGGGCGCAATTACCTCAACCCTCTGGACCCTGCT‑CAGGCCTGGGGATGCGCTGATTGTCGGGCGCACCCTGTATGGCTGTAc                                                                                                                       <  3:312289/149‑1 (MQ=255)
                       ccTTGCTGGAGCAGCGCATGGCTTCGCTGGAGGGCGGTGAAGCGGGGTTGGCGTTGGCATCGGGGATGGGCGCAATTACCTCAACCCTCTGGACCCTGCT‑CAGGCCTGGGGATGAGCTGATTCTCGGGCGCACCCTGTATGGCTGTAcc                                                                                                                      <  4:305190/149‑1 (MQ=255)
                        cTTGCTGGAGCAGCGCATGGCTTCGCTGGAGGGCGGTGAAGCGGGGTTGGCGTTGGCATCGGGGATGGGCGCAATTACCTCAACCCTCTGGACCCTGCT‑CAGGCCTGGGGATGAGCTGATTGTCGGGCGCACCCTGTATGGCTGTACCt                                                                                                                     >  3:60946/1‑149 (MQ=255)
                             tGGAGCAGCGCATGGCTTCGCTGGAGGGCGGTGAAGCGGGGTTGGCGTTGGCATCGGGGATGGGCGCAATTACCTCAACCCTCTGGACCCTGCT‑CAGGCCTGGGGATGAGCTGATTGTCGGGCGCACCCTGTATGGCTGTACCTTCGCt                                                                                                                >  2:266625/1‑149 (MQ=255)
                                     cgcATGGCTTCGCTGGAGGGCGGTGAAGCGGGGTTGGCGTTGGCATCGGGGATGGGCGCAATTACCTCAACCCTCTGGACCCTGCT‑CAGGCCTGGGGATGAGCTGATTGTCGGGCGCACCCTGTATGGCTGTACCTTCGCTTTTCTGca                                                                                                        >  4:193319/1‑149 (MQ=255)
                                      gcATGGCTTCGCTGGAGGGCGGTGAAGCGGGGTTGGCGTTGGCATCGGGGATGGGCGCAATTACCTCAACCCTCTGGACCCTGCT‑CAGGCCTGGGGATGAGCCGATTGTCGGGCGCACCCTGTATGGCTGTACCTTCGCTTTTCTGcac                                                                                                       >  1:179761/1‑149 (MQ=255)
                                      gcATGGCTTCGCTGGAGGGCGGTGAAGCGGGGTTGGCGTTGGCATCGGGGATGGGCGCAATTACCTCAACCCTCTGGACCCTGCT‑CAGGCCTGGGGATGAGCCGATTGTCGGGCGCACCCTGTATGGCTGTACCTTCGCTTTTCTGcac                                                                                                       >  1:174433/1‑149 (MQ=255)
                                           gCTTCGCTGGAGGGCGGTGAAGCGGGGTTGGCTTTGGCATGGGGGTCGGGCGCAATTACCTCAACCCTCTGGACCCTGCT‑CATGCCTGGGGATGAGCTGATTGTCGGGCGCACCCTGTATGTCTGTACCTTCGCTTTTCTGCACCACg                                                                                                   <  6:286013/148‑1 (MQ=255)
                                           gCTTCGCTGGAGGGCGGTGAAGCGGGGTTGGCGTTGGCATCGGGGATGGGCGCAATTACCTCAACCCTCTGGACCCTGCT‑CAGGCCTGGGGATGAGCTGATTGTCGGGCGCACCCTGTATgg                                                                                                                             >  4:129952/1‑122 (MQ=255)
                                           gCTTCGCTGGAGGGCGGTGAAGCGGGGTTGGCGTTGGCATCGGGGATGGGCGCAATTACCTCAACCCTCTGGACCCTGCT‑CAGGCCTGGGGATGAGCTGATTGTCGGGCGCACCCTGTATgg                                                                                                                             <  3:129952/122‑1 (MQ=255)
                                                gCTGGAGGGCGGTGAAGCGGGGTTGGCGTTGGCATCGGGGATGGGCGCAAATACCTCAAGCCTCTGGACCCTGCT‑GAGGCCTGGGGATGAGCTGATTGTCGGGCGCACCCTGTATGGCTGTACCTTCGCTTTTCTGCTCCACGGTATTg                                                                                             <  2:402865/149‑1 (MQ=255)
                                                   ggAGGGCGGTGAAGCGGGGTTGGCGTTGGCATCGGGGATGGGCGCAATTACCTCAACCCTCTGGACCCTGCT‑CAGGCCTGGGGATGAGCTGATTGTCGGGCGCACCCTGTATGGCTGTACCTTCGCTTTTCTGCACCACGGTATTGGt                                                                                           >  3:160827/1‑148 (MQ=255)
                                                    gAGGGCGGTGAAGCGGGGTTGGCGTTGGCATCGGGGATGGGCGCAATTACCTCAACCCTCTGGACCCTGCT‑CAGGCCTGGGGATGAGCTGATTGTCGGGCGCACCCTGTATGGCTGTACCTTCGCTTTTCTGCACCACGGTATTGGTGa                                                                                         <  8:4234/149‑1 (MQ=255)
                                                            tGAAGCGGGGTTGGCGTTGGCATCGGGGATGGGCGCAATTACCTCAACCCTCTGGACCCTGCT‑CAGGCCTGGGGATGAGCTGATTGTCGGGCGCACCCTGTATGGCTGTACCTTCGCTTTTCTGCACCACGGTATTGGTGAGTTCgggg                                                                                 >  8:117160/1‑149 (MQ=255)
                                                            tGAAGCGGGGTTGGCGTTGGCATCGGGGATGGGCGCAATTACCTCAACCCTCTGGACCCTGCT‑CAGGCCTGGGGATGAGCTGATTGTCGGGCGCACCCTGTATGGCTGTACCTTCGCTTTTCTGCACCACGGTATTGGTGAGTTCgggg                                                                                 <  3:186177/149‑1 (MQ=255)
                                                                gCGGGGTTGGCGTTGGCATCGGGGATGGGCGCAATTACCTCAACCCTCTGGACCCTGCT‑CAGGCCTGGGGATGAGCTGATTGTCgg                                                                                                                                            >  6:7112/1‑86 (MQ=255)
                                                                gCGGGGTTGGCGTTGGCATCGGGGATGGGCGCAATTACCTCAACCCTCTGGACCCTGCT‑CAGGCCTGGGGATGAGCTGATTGTCgg                                                                                                                                            <  5:7112/86‑1 (MQ=255)
                                                                gCGGGGTTGGCGTTGGCATCGGGGATGGGCGCAATTACCTCAACCCTCTGGACCCTGCT‑CAGGCCTGGGGATGAGCTGATTGTCGGGCGCACCCTGTATGGCTGTACCTTCGCTTTTCTGCACCACg                                                                                                   >  1:401325/1‑127 (MQ=255)
                                                                gCGGGGTTGGCGTTGGCATCGGGGATGGGCGCAATTACCTCAACCCTCTGGACCCTGCT‑CAGGCCTGGGGATGAGCTGATTGTCGGGCGCACCCTGTATGGCTGTACCTTCGCTTTTCTGCACCACg                                                                                                   <  2:401325/127‑1 (MQ=255)
                                                                 gcgggTTGGCCTAGGGTTCGGGGTTGGGGGCTATTACCTCAACCCTCTGGACCCTGCT‑CAGGCCTGGGGATGAGCCGATTGTCGGGCGCACCCTGTATGGCTGTACCTTCGCTTTTCTGCACCACGGTATTGGTGAGTTCGGGGTGaaa                                                                            <  2:174433/147‑1 (MQ=255)
                                                                 cGGGGTTGGCGTTGGCATCGGGGATGGGCGCAATTACCTCAACCCTCTGGACCCTGCT‑CAGGCCTGGGGATGAGCTGATTGTCGGGCGCACCCTGTATGGCTGt                                                                                                                         >  3:334625/1‑104 (MQ=255)
                                                                 cGGGGTTGGCGTTGGCATCGGGGATGGGCGCAATTACCTCAACCCTCTGGACCCTGCT‑CAGGCCTGGGGATGAGCTGATTGTCGGGCGCACCCTGTATGGCTGt                                                                                                                         <  4:334625/104‑1 (MQ=255)
                                                                 cGGGGTTGGCGTTGGCATCGGGGATGGGCGCAATTACCTCAACCCTCTGGACCCTGCT‑CAGGCCTGGGGATGAGCTGATTGTCGGGCGCACCCTGTATGGCTGTACCTTCGCTTTTCTGCACCACGGTATTGGTGAGTTCGGGGTGaaa                                                                            <  3:399637/149‑1 (MQ=255)
                                                                 cGGGGTTGGCGTTGGCATCGGGGATGGGCGCAATTACCTCAACCCTCTGGACCCTGCT‑CAGGCCTGGGGATGAGCCGATTGTCGGGCGCACCCTGTATGGCTGTACCTTCGCTTTTCTGCACCACGGTATTGGTGAGTTCGGGGTGaaa                                                                            <  2:179761/149‑1 (MQ=255)
                                                                       tggcgttggcATCGGGGATGGGCGCAATTACCTCAACCCTCTGGACCCTGCT‑CAGGCCTGGGGATGAGCTGATTGTCGGGCGCACCCTGTATGGCTGTACCTTCGCTTTTCTGCACCACGGTATTGGTGAGTTCGGGGTGAAAGTTCGc                                                                      >  3:41921/1‑149 (MQ=255)
                                                                         gcgttggcATCGGGGATGGGCGCAATTACCTCAACCCTCTGGACCCTGCTACAGGCCTGGGGATGAGCTGATTGTCGGGCGCACCCTGTATGGCTGTACCTTCGCTTTTGTGCACCACGGTATTGGTGAGTTCGGGGTGAAAGTTCGc                                                                      >  1:108845/1‑148 (MQ=255)
                                                                            ttggcATCGGGGATGGGCGCAATTACCTCAACCCTCTGGACCCTGCT‑CAGGCCTGGGGATGAGCTGATTGTCGGGCGCACCCTGTATGGCTGTACCTTCGCTTTTCTGCACCACGGTATTGGTGAGTTCGGGGTGAAAGTTCGCCATGt                                                                 >  5:97279/1‑149 (MQ=255)
                                                                                      ggATGGGCGCAATTACCTCAACCCTCTGGACCCTGCT‑CAGGCCTGGGGATGAGGTGATTGTCGGGCGCACCCTGTATGGCTGTACCTTCGCTTTTCTGCACCACGGTATTGGTGAGTTCGGGGTGAAAGTTCGCCATGTCGATCTCAAc                                                       <  6:110924/149‑1 (MQ=255)
                                                                                                     ccaccaCCCTCTGGCCCCTGCT‑CCGGCCGGGGGATGAGCTGATTGTCGGTCGCACCCTGTATTGCTGTTCCTTCGCTTTTCTGCACCTCTGTATTGGTGAGTTCTTGGTGAAAGTTCGCCATGTCGATCTCAACGATGCCAAGGCCCt                                         <  4:72546/143‑1 (MQ=255)
                                                                                                      cTCAACCCTCTGGACCCTGCT‑CAGGCCTGGGGATGAGCTGATTGTCGGGCGCACCCTGTATGGCTGTACCTTCGCTTTTCTGCACCACGGTATTGGTGAGTTCGGGGTGAAAGTTCGCCATGTCGATCTCAACGATGCCAAGGCCCTg                                        >  3:340253/1‑148 (MQ=255)
                                                                                                        cAACCCTCTGGACCCTGCT‑CAGGCCTGGGGATGAGCTGATTGTCGGGCGCACCCTGTATGGCTGTACCTTCGCTTTTCTGCACCACGGTATTGGTGAGTTCGGGGTGAAAGTTCGCCATGTCGATCTCAACGATGCCAAGGCCCTGa                                       >  3:186444/1‑147 (MQ=255)
                                                                                                             ctctGGACCCTGCT‑CAGGCCTGGGGATGAGCTGATTGTCGGGCGCACCCTGTATGGCTGTACCTTCGCTTTTCTGCACCACGGTATTGGTGAGTTCGGGGTGAAAGTTCGCCATGTCGATCTCAACGATGCCAAGGCCCTGAAAgcggc                                <  1:437309/149‑1 (MQ=255)
                                                                                                                  gACCCTGCT‑CAGGCCTGGGGATGAGCTGATTGTCGGGCGCACCCTGTATGGCTGTACCTTCGCTTTTCTGCACCACGGTATTGGTGAGTTCGGGGTGAAAGTTCGCCATGTCGATCTCAACGATGCCAAGGCCCTGAAAGCGGCGa                              <  3:193319/146‑1 (MQ=255)
                                                                                                                    cccTGCT‑CAGGCCTGGGGATGAGCTGATTGTCGGGCGCACCCTGTATGGCTGTACCTTCGCTTTTCTGCACCACGGTATTGGTGAGTTCGGGGTGAAAGTTCGCCATGTCGATCTCAACGATGCCAAGGc                                            <  4:396214/130‑1 (MQ=255)
                                                                                                                    cccTGCT‑CAGGCCTGGGGATGAGCTGATTGTCGGGCGCACCCTGTATGGCTGTACCTTCGCTTTTCTGCACCACGGTATTGGTGAGTTCGGGGTGAAAGTTCGCCATGTCGATCTCAACGATGCCAAGGc                                            >  3:396214/1‑130 (MQ=255)
                                                                                                                    cccTGCT‑CAGGCCTGGGGATGAGCTGATTGTCGGGCGCACCCTGTATGGCTGTACCTTCGCTTTTCTGCACCACGGTATTGGTGAGTTCGGGGTGAAAGTTCGCCATGTCGATCTCAACGATGCCAAGGCCCTGAAAGCGGCGATcagc                         <  4:41921/149‑1 (MQ=255)
                                                                                                                     ccTGCT‑CAGGCCTGGGGATGAGCTGATTGTCGGGCGCACCCTGTATGGCTGTACCTTCGCTTTTCTGCACCACGGTATTGGTGAGTTCGGGGTGAAAGTTCGCCATGTCGATCTCAACGATGCCAAGGCCCTGAAAGCGGCGATcagca                        >  4:366809/1‑149 (MQ=255)
                                                                                                                     ccTGCT‑CAGGCCTGGGGATGAGCTGATTGTCGGGCGCACCCTGTATGGCTGTACCTTCGCTTTTCTGCACCACGGTATTGGTGAGTTCAGGGTGAAAGTTCGCCATGTCGATCTCAACGATGCCAAGGCCCTGAAAGCGGCGATcagca                        <  1:266625/149‑1 (MQ=255)
                                                                                                                      cTGCT‑CAGGCCTGGGGATGAGCTGAGTGTCGGGCGCAGCCTGTATGGCTGTACCTTCGCTTTTCTGCACCACGGTATTGGTGAGTTCGGGGTGAAAGTTCGCGATGTCGATCTCAACGATGCCAAGGCCCTGAAAGCGGCGATcagcag                       <  8:181098/149‑1 (MQ=255)
                                                                                                                        gCT‑CAGGCCTGGGGATGAGCTGATTGTCGGGCGCACCCTGTATGGCTGTACCTTCGCTTTTCTGCACCACGGTATTGGTGAGTTCGGGGTGAAAGTTCGCCATGTCGATCTCAACGATGCCAAGGCCCTGAAAGCGGCGATCagcagc                      >  5:282082/1‑148 (MQ=255)
                                                                                                                                ccTGGGGATGAGCTGATTGTCGGGCGCACCCTGTATGGCTGTACCTTCGCTTTTCTGCACCACGGTATTGGTGAGTTCGGGGTGAAAGTTCGCCATGTCGATCTCAACGATGCCAAGGCCCTGAAAGCGGCGATCAGCAGCAAAACACg              <  4:60946/149‑1 (MQ=255)
                                                                                                                                ccTGGGGATGAGCTGATTGTCGGGCGCACCCTGTATGGCTGTACCTTCGCTTTTCTGCACCACGGTATTGGTGAGTTCGGGGTGAAAGTTCGCCATGTCGATCTCAACGATGCCAAGGCCCTGAAAGCGGCGATCAGCAGCAAAACACg              <  4:340253/149‑1 (MQ=255)
                                                                                                                                      gATGAGCTGATTGTCGGGCGCACCCTGTATGGCTGTACCTTCGCTTTTCTGCACCACGGTATTGGTGAGTTCGGGGTGAAAGTTCGCCATGTCGATCTCAACGATGCCAAGGCCCTGAAAGCGGCGATCAGCAGCAAAACACGgatgat        <  3:366809/149‑1 (MQ=255)
                                                                                                                                          aGCTGATTGTCGGGCGCACCCTGTATGGCTGTACCTTCGCTTTTCTGCACCACGGTATTGGTGAGTTCGGGGTGAAAGTTCGCCATGTCGATCTCAACGATGCCAAGGCCCTGAAAGCGGCGATCAGCAGCAAAACACGGATGATCTa     <  8:257545/148‑1 (MQ=255)
                                                                                                                                           gCTGATTGTCGGGCGCACCCTGTATGGCTGTATCTTCGCTTTTCTGCACAACGGTATTGGTGAGTTCGGGGTGAAAGTTCGCCATGTCGATCTCAACGATGCCAAGGCCCTGAAAGCGGCGATCAGCAGCAAAACACGGATGATCTACt   >  5:416492/1‑149 (MQ=255)
                                                                                                                                            cTGATTGTCGGGCGCACCCTGTATGGCTGTACCTTCGCTTTTCTGCACCACGGTATTGGTGAGTTCGGGGTGAAAGTTCGCCATGTCGATCTCAACGATGCCAAGGCCCTGAAAGCGGCGATCAGCAGCAAAACACGGATGATCTACtt  >  2:248418/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
CCGCATTTCCAACCCCACGCTGGCCTTGCTGGAGCAGCGCATGGCTTCGCTGGAGGGCGGTGAAGCGGGGTTGGCGTTGGCATCGGGGATGGGCGCAATTACCTCAACCCTCTGGACCCTGCT‑CAGGCCTGGGGATGAGCTGATTGTCGGGCGCACCCTGTATGGCTGTACCTTCGCTTTTCTGCACCACGGTATTGGTGAGTTCGGGGTGAAAGTTCGCCATGTCGATCTCAACGATGCCAAGGCCCTGAAAGCGGCGATCAGCAGCAAAACACGGATGATCTACTT  >  NC_002947/1495751‑1496038

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