Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,080,920 A→T 9.5% intergenic (‑144/+19) amaC ← / ← dpkA D‑lysine aminotransferase/delta 1‑piperideine‑2‑carboxylate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,080,9200AT 9.5% 143.8 / 4.0 50intergenic (‑144/+19)amaC/dpkAD‑lysine aminotransferase/delta 1‑piperideine‑2‑carboxylate reductase
Reads supporting (aligned to +/- strand):  ref base A (21/24);  new base T (3/2);  total (24/26)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

CACAGCGTGCAGCTCCTTGGATCAGGACGGTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATGCGCGTCAGCCCCACCGCCTGCATGTGC  >  NC_002947/4080774‑4081068
                                                                                                                                                  |                                                                                                                                                    
cacaGCGTGCAGCTCCTTGGATCAGGACGGTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAAc                                                                                                                                                    >  8:88404/1‑149 (MQ=255)
   aGCGTGCAGCTCCTTGGATCAGGACGGTACGGAAAGTCAGTGTTGAGCATAGCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCcg                                                                                                                                                 <  5:10051/149‑1 (MQ=255)
       tGCAGCTCCTTGGATCAGGACGGTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCcgcg                                                                                                                                               >  6:168965/1‑147 (MQ=255)
               cTTGGATCAGGACGGTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGcc                                                                                                                                       <  5:117312/147‑1 (MQ=255)
                        ggACGGTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGcc                                                                                                                              >  5:279955/1‑147 (MQ=255)
                           cGGTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAagca                                                                                                                         <  4:350086/149‑1 (MQ=255)
                             gTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTTTACCCGCGaaaaa                                                                                                                                          >  8:282953/1‑127 (MQ=255)
                             gTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTTTACCCGCGaaaaa                                                                                                                                          <  7:282953/130‑4 (MQ=255)
                             gTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAagcagc                                                                                                                       >  6:70392/1‑149 (MQ=255)
                              tACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGct                                                                                                                      <  7:88404/149‑1 (MQ=255)
                                  gAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCttt                                                                                                                  <  6:279955/149‑1 (MQ=255)
                                     aGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAg                                                                                                               <  5:286589/149‑1 (MQ=255)
                                               gAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAg                                                                                                               <  4:76643/139‑1 (MQ=255)
                                               gAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAg                                                                                                               >  3:76643/1‑139 (MQ=255)
                                               gAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAAc                                                                                                     >  1:321382/1‑149 (MQ=255)
                                                  cATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGCTGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTcc                                                                                                  <  4:75069/149‑1 (MQ=255)
                                                     aCCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGc                                                                                               >  4:197265/1‑149 (MQ=255)
                                                     aCCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGc                                                                                               >  4:110094/1‑149 (MQ=255)
                                                      cccGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGc                                                                                               <  3:110094/148‑1 (MQ=255)
                                                         ggAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCgg                                                                                           <  3:319713/149‑1 (MQ=255)
                                                                  gTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAAc                                                                                                     >  4:120497/1‑130 (MQ=255)
                                                                  gTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAAc                                                                                                     <  3:120497/130‑1 (MQ=255)
                                                                      aGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCAcccc                                                                              >  8:80167/1‑149 (MQ=255)
                                                                         aaCATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCcctc                                                                           <  2:304354/149‑1 (MQ=255)
                                                                               cACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCctcctcctc                                                                     >  1:166930/1‑149 (MQ=255)
                                                                                      gtgCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCAc                                                               <  5:168965/148‑1 (MQ=255)
                                                                                       tgCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCt                                                             >  2:155101/1‑149 (MQ=255)
                                                                                         cAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCcctc                                                                           >  1:83274/1‑133 (MQ=255)
                                                                                         cAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCcctc                                                                           <  2:83274/133‑1 (MQ=255)
                                                                                         cAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTcg                                                           >  7:218017/1‑149 (MQ=255)
                                                                                                gTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCa                                                     <  5:70392/148‑1 (MQ=255)
                                                                                                 tCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCAc                                                    <  2:166930/148‑1 (MQ=255)
                                                                                                   tATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAg                                                                                                               <  7:45037/87‑1 (MQ=255)
                                                                                                   tATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAg                                                                                                               >  8:45037/1‑87 (MQ=255)
                                                                                                    aTGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTcgc                                                          <  1:36072/139‑1 (MQ=255)
                                                                                                    aTGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTcgc                                                          >  2:36072/1‑139 (MQ=255)
                                                                                                     tGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACgg                                                  <  8:218017/146‑1 (MQ=255)
                                                                                                      gCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTcgc                                                          <  5:340080/137‑1 (MQ=255)
                                                                                                      gCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTcgc                                                          >  6:340080/1‑137 (MQ=255)
                                                                                                        aaCATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAgcgc                                            >  6:134779/1‑149 (MQ=255)
                                                                                                                  gACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAAGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCg                                  >  4:51557/1‑149 (MQ=255)
                                                                                                                        gATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATgc                            <  4:55433/149‑1 (MQ=255)
                                                                                                                                ctggctTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGc                                                                                                             >  3:332997/1‑60 (MQ=255)
                                                                                                                                ctggctTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGc                                                                                                             <  4:332997/60‑1 (MQ=255)
                                                                                                                                 tggctTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATGCGCGTCAGcc                   >  8:245130/1‑149 (MQ=255)
                                                                                                                                     cccttCGCGGGTTTACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCgg                                 <  6:29690/128‑1 (MQ=255)
                                                                                                                                     cccttCGCGGGTTTACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCgg                                 >  5:29690/4‑131 (MQ=255)
                                                                                                                                       tttCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATGCGCGTCAGCCCCAcc              <  5:132448/148‑1 (MQ=255)
                                                                                                                                        ttCGCGGGTTTACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATGCGCGTCAGCCCCACCGc            >  6:327881/1‑149 (MQ=255)
                                                                                                                                           gcgGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATGCGCGTCAGCCCCACCGCCTg         <  3:197265/149‑1 (MQ=255)
                                                                                                                                                 aaaCCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATGCGCGTCAGCCCCACCGCCTGCAtgtg   >  6:124601/1‑149 (MQ=255)
                                                                                                                                                  aaCCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATGCGCGTCAGCCCCACCGCCTGCATGTGc  >  5:368322/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CACAGCGTGCAGCTCCTTGGATCAGGACGGTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATGCGCGTCAGCCCCACCGCCTGCATGTGC  >  NC_002947/4080774‑4081068

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