Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,080,919 A→T 9.2% intergenic (‑143/+20) 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,9190AT 9.2% 149.7 / 4.3 54intergenic (‑143/+20)amaC/dpkAD‑lysine aminotransferase/delta 1‑piperideine‑2‑carboxylate reductase
Reads supporting (aligned to +/- strand):  ref base A (24/25);  new base T (3/2);  total (27/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 = 9.03e-01

TCCTTGGATCAGGACGGTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATGCGCGTCAGCCCCACCGCCTGCATGTG  >  NC_002947/4080787‑4081067
                                                                                                                                    |                                                                                                                                                    
tCCTTGGATCAGGACGGTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGc                                                                                                                                       >  2:114071/1‑148 (MQ=255)
     ggATCAGGACGGTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAAt                                                                                                                                  <  5:106729/148‑1 (MQ=255)
     ggATCAGGACGGTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAAt                                                                                                                                  >  5:119959/1‑148 (MQ=255)
          aGGACGGTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGcc                                                                                                                             <  2:268212/148‑1 (MQ=255)
            gACGGTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAag                                                                                                                          >  7:171106/1‑149 (MQ=255)
            ccgagctcGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCaa                                                                                                                           <  5:200511/141‑1 (MQ=255)
                gTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATACACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCAGCGAAGGGGACAAATCAGCCAagcagc                                                                                                                      >  3:359414/1‑149 (MQ=255)
                 tACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAag                                                                                                                          >  8:248644/1‑144 (MQ=255)
                    ggAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCaa                                                                                                                           >  6:200511/1‑140 (MQ=255)
                    ggAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCtt                                                                                                                  >  4:164813/1‑149 (MQ=255)
                     gAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGTCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCttt                                                                                                                 <  3:287610/149‑1 (MQ=255)
                         gTCAGTGTTGAGCATTGCCGGAGTGGGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCATGATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCGAAATCAGCCAAGCAGCTCTTTCAgg                                                                                                             <  2:177379/149‑1 (MQ=255)
                              tgtTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGTGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCttt                                                                                                        <  7:149773/149‑1 (MQ=255)
                              tgtTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCttt                                                                                                        >  8:275091/1‑149 (MQ=255)
                                ttGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCa                                                                                                                                     <  5:193723/118‑1 (MQ=255)
                                ttGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCa                                                                                                                                     >  6:193723/1‑118 (MQ=255)
                                ttGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTc                                                                                                       <  7:275091/148‑1 (MQ=255)
                                  gAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAAccc                                                                                                                                                 >  5:301940/1‑104 (MQ=255)
                                  gAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAAccc                                                                                                                                                 <  6:301940/104‑1 (MQ=255)
                                    gCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTc                                                                                                  >  1:366673/1‑149 (MQ=255)
                                           cGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCg                                                                                           <  7:248644/149‑1 (MQ=255)
                                           cGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCg                                                                                           <  4:263906/149‑1 (MQ=255)
                                                  cGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCAcc                                                                                     <  8:171106/148‑1 (MQ=255)
                                                   gggTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGCCCTTTCAACTCCTGCTCGGTCACCGc                                                                                   <  1:10199/149‑1 (MQ=255)
                                                         aGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCTAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTc                                                                                                  <  6:63331/128‑1 (MQ=255)
                                                         aGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCTAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTc                                                                                                  >  5:63331/1‑128 (MQ=255)
                                                         aGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAg                                                                                                                               <  2:205777/99‑1 (MQ=255)
                                                         aGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAg                                                                                                                               >  1:205777/1‑99 (MQ=255)
                                                           gAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCcct                                                                           >  7:302782/1‑149 (MQ=255)
                                                                tCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCcctcctc                                                                       >  8:175693/1‑148 (MQ=255)
                                                                 ccACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCcctcctcct                                                                     >  4:46892/1‑149 (MQ=255)
                                                                   aCTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCg                                                                   <  6:119959/149‑1 (MQ=255)
                                                                   aCTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCg                                                                   >  8:80268/1‑149 (MQ=255)
                                                                   aCTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCg                                                                   <  3:164813/149‑1 (MQ=255)
                                                                     tATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGc                                                                 >  1:103703/1‑149 (MQ=255)
                                                                     tATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGc                                                                 <  3:46892/149‑1 (MQ=255)
                                                                     tATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGc                                                                 >  7:52095/1‑149 (MQ=255)
                                                                       ttGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGc                                                                 <  2:103703/147‑1 (MQ=255)
                                                                         gtgCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCAc                                                              <  1:114071/148‑1 (MQ=255)
                                                                            cAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTcg                                                          >  5:210765/1‑149 (MQ=255)
                                                                                           aaCATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTACTCGGCCACCTCGCGCTCACGGTAgctc                                           >  8:149551/1‑147 (MQ=255)
                                                                                                 gCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTc                                      <  2:366673/148‑1 (MQ=255)
                                                                                                           gATGAGGGCTGGCTTTTTAGCGGGTAAACACGCGAAGGGGCCAAATAAGACAAGCAGCTCTTTCAGGCCTTTCAACTCATGCTAGGTCAACGCCACACCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCAGCTCGCCCGGCAttc                           >  6:59068/1‑147 (MQ=255)
                                                                                                                 ggctggctTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATGCGCGTCa                     <  7:175693/149‑1 (MQ=255)
                                                                                                                   ctggctTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATGCGCGTCAGc                   <  7:80268/149‑1 (MQ=255)
                                                                                                                   ctggctTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATGCGCGTCAGc                   <  8:52095/149‑1 (MQ=255)
                                                                                                                        cccttCGCGGGTTTACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATgcgc                         <  2:54981/135‑1 (MQ=255)
                                                                                                                        cccttCGCGGGTTTACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATgcgc                         >  1:54981/4‑138 (MQ=255)
                                                                                                                         ccttCGCGGGTTTACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATgcgc                         >  2:183238/3‑137 (MQ=255)
                                                                                                                         ccttCGCGGGTTTACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATgcgc                         <  1:183238/135‑1 (MQ=255)
                                                                                                                           ttCGCGGGTTTACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATGCGCGTCAGCCCCACCGc           >  1:340171/1‑149 (MQ=255)
                                                                                                                              gcgGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATGCGCGTCAGCCCCACCGCCTg        >  3:350099/1‑149 (MQ=255)
                                                                                                                                gggTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATGCGCGTCAGCCCCACCGCCTGCa      <  3:296501/149‑1 (MQ=255)
                                                                                                                                  gTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATGCGCGTCAGCCCCACCGCCTGCAtg    >  3:75665/1‑149 (MQ=255)
                                                                                                                                    aaaCCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATGCGCGTCAGCCCCACCGCCTGCAtgtg  >  6:339729/1‑149 (MQ=255)
                                                                                                                                    aaaCCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATGCGCGTCAGCCCCACCGCCTGCAt     >  3:330231/1‑146 (MQ=255)
                                                                                                                                    aaaCCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATGCGCGTCAGCCCCACCGCCTGCAt     <  4:330231/146‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                                    
TCCTTGGATCAGGACGGTACGGAAAGTCAGTGTTGAGCATACCCGGAGTGCGGGTTCAGGAACATCCACTATTGTGCAAGGTTGTCTATGCAACATTGCATGACCGGGATGAGGGCTGGCTTTTTCGCGGGTAAACCCGCGAAGGGGCCAAATCAGCCAAGCAGCTCTTTCAGGCCTTTCAACTCCTGCTCGGTCACCGCCACCCCCTCCTCCTCGGCCACCTCGCGCTCACGGTAGCGCCGCTCGCCCGGCATGCGCGTCAGCCCCACCGCCTGCATGTG  >  NC_002947/4080787‑4081067

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: ‑