Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,001,235 A→G 8.1% F356L (TTC→CTC)  amaA ← L‑pipecolate oxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,001,2350AG 8.1% 132.7 / 2.8 49F356L (TTC→CTC) amaAL‑pipecolate oxidase
Reads supporting (aligned to +/- strand):  ref base A (21/24);  new base G (2/2);  total (23/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.45e-01

ACACCGCTCATCGGCGAAAACGGAATCGGCAGGCTGTCCGCCTCGCCCTTGAGCAGGAACTCGGCAAACGCCCGGCCAATCACGGTGCCGGTCGTGTTACCCCGCCCGTTGTACCCAGTCACCGCCACCAACCCAGGCGCCGGCTCGAACAGGCGCATCAGGTGGTCAGGGGTGAAGTCGATGCAGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGATGCGATCGGCCCAGCTGCGCACGAACCAGGCTGGCTTGTTGTCGACCCGGCC  >  NC_002947/6001087‑6001382
                                                                                                                                                    |                                                                                                                                                   
acacCGCTCATCGGCGAAAACGGAATCGGCAGGCTGTCCGCCTCGCCCTTGAGCAGGAACTCGGCAAACGCCCGGCCAATCACGGTGCCGGTCGTGTTACCCCGCCCGTTGTACCCAGTCACCGCCACCAACCCAGGCGCCGGCTCGaa                                                                                                                                                     <  2:127231/149‑1 (MQ=255)
   ccGCTCATCGGCGAAAACGGAATCGGCAGGCTGTCCGCCTCGCCCTTGAGCAGGAACTCGGCAAACGCCCGGCCAATCACGGTGCCGGTCGTGTTACCCCGCCCGTTGTACCCAGTCACCGCCACCAACCCAGGCGCCGGCTCGAACAg                                                                                                                                                  <  2:42687/149‑1 (MQ=255)
          tCGGCGAAAACTGAATCGGCAGGCGGTCCGCCTGGCCCTTGAGCAGGAACTCGGCAAACGCCCGGCCAATAACGGTGCCGGTCGTGTTACCCCGACCGTTGTACCCAGTCACCGCCACCAACCCAGGCGCCTGTTCGAGCCGgcgcctg                                                                                                                                           <  7:164202/149‑4 (MQ=255)
          tCGGCGAAAACGGAATCGGCAGGCTGTCCGCCTCGCCCTTGAGCAGGAACTCGGCAAACGCCCGGCCAATCACGGTGCCGGTCGTGTTACCCCGCCCGTTGTACCCAGTCACCGCCACCAACCCAGGCGCCGGCTCGAACAGGCGCAt                                                                                                                                            >  2:285895/1‑148 (MQ=255)
            ggCGAAAACGGAATCGGCAGGCTGTCCGCCTCGCCCTTGAGCAGGAACTCGGCAAACGCCCGGCCAATCACGGTGCCGGTCGTGTTACCCCGCCCGTTGTACCCAGTCACCGCCACCAACCCAGGCGCCGGCTCGAACAGGCGCATCAg                                                                                                                                         >  6:1316/1‑149 (MQ=255)
            ggCGAAAACGGAATCGGCAGGCTGTCCGCCTCGCCCTTGAGCAGGAACTCGGCAAACGCCCGGCCAATCACGGTGCCGGTCGTGTTACCCCGCCCGTTGTACCCAGTCACCGCCACCAACCCAGGCGCCGGCTCGAACAGGCGCATCAg                                                                                                                                         <  3:221581/149‑1 (MQ=255)
             gCGAAAACGGAATCGGCAGGCTGTCCGCCTCGCCCTTGAGCAGGAACTCGGCAAACGCCCGGCCAATCACGGTGCCGGTCGTGTTACCCCGCCCGTTGTACCCAGTCACCGCCACCAACCCAGGCGCCGGCTCGAACAGGCGCATCAgg                                                                                                                                        <  4:263114/149‑1 (MQ=255)
               gAAAACGGAATCGGCAGGCTGTCCGCCTCGCCCTTGAGCAGGAACTCGGCAAACGCCCGGCCACTCACGGTGCCGGTCGTGTTACCCCGCCCGTTGTACCCAGTCACCGCCACCAACCCAGGCGCCGGCTCGAACAGGCGCATCAgggg                                                                                                                                      >  7:347703/1‑147 (MQ=255)
               gAAAACGGAATCGGCAGGCTGTCCGCCTCGCCCTTGAGCAGGAACTCGGCAAACGCCCGGCCAATCACCGTGCCGGTCGTCTTACCCCGCCCGTTGTACCCAGTCACCGCCACCACCCCAGGCGCCGGATCGAACAGGCGCATCAggtg                                                                                                                                      >  3:30596/1‑149 (MQ=255)
                 aaaCGGAATCGGCAGGCTGTCCGCCTCGCCCTTGAGCAGGAACTCGGCAAACGCCCGGCCAATCACGGTGCCGGTCGTGTTACCCCGCCCGTTGTACCCAGTCACCGCCACCAACCCAGGCGCCGGCTCGAACAGGCGCATCAggtggt                                                                                                                                    <  6:300695/149‑1 (MQ=255)
                  aaCGGAATCGGCAGGCTGTCCGCCTCGCCCTTGAGCAGGAACTCGGCAAACGCCCGGCCAATCACGGTGCCGGTCGTGTTACCCCGCCCGTTGTACCCAGTCACCGCCACCAACCCAGGCGCCGGCTCGAACAGGCGCATCAggtggt                                                                                                                                    >  1:9405/1‑148 (MQ=255)
                    cGGAATCGGCAGGCTGTCCGCCTCGCCCTTGAGCAGGAACTCGGCAAACGCCCGGCCAATCACGGTGCCGGTCGTGTTACCCCGCCCGTTGTACCCAGTCACCGCCACCAACCCAGGCGCCGGCTCGAACAGGCGCATCAGGTGGTCAg                                                                                                                                 <  6:352441/149‑1 (MQ=255)
                     ggAATCGGCAGGCTGTCCGCCTCGCCCTTGAGCAGGAACTCGGCAAACGCCCGGCCAATCACGGTGCCGGTCGTGTTACCCCGCCCGTTGTACCCAGTCACCGCCACCAACCCAGGCGCCTGTTCGAGCCGgc                                                                                                                                                >  1:190431/1‑133 (MQ=255)
                     ggAATCGGCAGGCTGTCCGCCTCGCCCTTGAGCAGGAACTCGGCAAACGCCCGGCCAATCACGGTGCCGGTCGTGTTACCCCGCCCGTTGTACCCAGTCACCGCCACCAACCCAGGCGCCTGTTCGAGCCGgc                                                                                                                                                <  2:190431/133‑1 (MQ=255)
                      gAATCGGCAGGCTGTCCGCCTCGCCCTTGAGCAGGAACTCGGCAAACGCCCGGCCAATCACGGTGCCGGTCGTGTTACCCCGCCCGTTGTACCCAGTCACCGCCACCAACCCAGGCGCCGGCTCGAACAGGCGCATCAGGTGGTCAggg                                                                                                                               <  2:332944/149‑1 (MQ=255)
                      gAATCGGCAGGCTGTCCGCCTCGCCCTTGAGCAGGAACTCGGCAAACGCCCGGCCAATCACGGTGCCGGTCGTGTTACCCCGCCCGTTGTACCCAGTCACCGCCACCAACCCAGGCGCCGGCTCGAACAGGCGCATCAGGTGGTCAggg                                                                                                                               <  7:62085/149‑1 (MQ=255)
                       aaTCGGCAGGCTGTCCGCCTCGCCCTTGAGCAGGAACTCGGCAAACGCCCGGCCAATCACGGTGCCGGTCGTGTTACCCCGCCCGTTGTACCCAGTCACCGCCACCAACCCAGGCGCCGGCTCGAACAGGCGCATCAGGTGGTCAgggg                                                                                                                              >  3:131568/1‑149 (MQ=255)
                       aaTCGGCAGGCTGTCCGCCTCGCCCTTGAGCAGGAACTCGGCAAACGCCAGGCCAATCACGGTGCCGGTCGTGTTACCCCGCCCGTTGTACCCAGTCACCGCCACCAACCCAGGCGCCGGCTCGAACAGGCGCATCAGGTGGTCAgggg                                                                                                                              >  5:57501/1‑149 (MQ=255)
                           ggCAGGCTGTCCGCCTCGCCCTTGAGCAGGAACTCGGCAAACGCCCGGCCAATCACGGTGCCGGTCGTGTTACCCCGCCCGTTGTACCCAGTCACCGCCACCAACCCAGGCGCCGGCTCGAACAGGCGCATCAGGTGGTCAGGGGTGaa                                                                                                                          <  8:170432/149‑1 (MQ=255)
                             cAGGCTGTCCGCCTCGCCCTTGAGCAGGAACTCGGCAAACGCCCGGCCAATCACGGTGCCGGTCGTGTTACCCCGCCCGTTGTACCCAGTCACCGCCACCAACCCAGGCGCCGGCTCGAACAGGCGCATCAGGTGGTCAGGGGTGAAGt                                                                                                                        <  8:116908/149‑1 (MQ=255)
                             cAGGCTGTCCGCCTCGCCCTTGAGCAGGAACTCGGCAAACGCCCGGCCAATCACGGTGCCGGTCGTGTTACCCCGCCCGTTGTACCCAGTCACCGCCACCAACCCAGGCGCCGGCTCGAACAGGCGCATCAGGTGGTCAGGGGTGAAGt                                                                                                                        >  7:84864/1‑149 (MQ=255)
                                              ccTTGAGCAGGAACTCGGCAAACGCCCGGCCAATCACGGTGCCGGTCGTGTTACCCCGCCCGTTGTACCCAGTCACCGCCACCAACCCAGGCGCCGGCTCGAACAGGCGCATCAGGTGGTCAGGGGTGAAGTCGATGCAGCCGGTCCAg                                                                                                       >  4:16673/1‑149 (MQ=255)
                                                    gCAGGAACTCGGCAAACGCCCGGCCAATCACGGTGCCGGTCGTGTTACCCCGCCCGTTGTACCCAGTCACCGCCACCAACCCAGGCGCCGGCTCGAACAGGCGCATCAGGTGGTCAGGGGTGAAGTCGATGCAGCCGGTCCAGTGCAtt                                                                                                 >  4:48050/1‑149 (MQ=255)
                                                    gCAGGAACTCGGCAAACGCCCGGCCAATCACGGTGCCGGTCGTGTTACCCCGCCCGTTGTACCCAGTCACCGCCACCAACCCAGGCGCCGGCTCGAACAGGCGCATCAGGTGGTCAGGGGTGAAGTCGATGCAGCCGGTCCAGTGCAtt                                                                                                 <  4:131568/149‑1 (MQ=255)
                                                      aGGAACTCGGCAAACGCCCGGCCAATCACGGTGCCGGTCGTGTTACCCCGCCCGTTGTACCCAGTCACCGCCACCAACCCAGGCGCCGGCTCGAACAGGCGCATCAGGTGGTCAGGGGTGAAGTCGATGCAGCCGGTCCAGTGCATTTc                                                                                               <  5:1316/149‑1 (MQ=255)
                                                              ggCAAACGCCCGGCCAATCACGGTGCCGGTCGTGTTACCCCGCCCGTTGTACCCAGTCACCGCCACCAACCCAGGCGCCGGCTCGAACAGGCGCATCAGGTGGTCAGGGGTGAAGTCGATGCAGCCGGTCCAGTGCAt                                                                                                  <  4:160690/138‑1 (MQ=255)
                                                              ggCAAACGCCCGGCCAATCACGGTGCCGGTCGTGTTACCCCGCCCGTTGTACCCAGTCACCGCCACCAACCCAGGCGCCGGCTCGAACAGGCGCATCAGGTGGTCAGGGGTGAAGTCGATGCAGCCGGTCCAGTGCAt                                                                                                  >  3:160690/1‑138 (MQ=255)
                                                              ggCAAACGCCCGGCCAATCACGGTGCCGGTCGTGTTACCCCGCCCGTTGTACCCAGTCACCGCCACCAACCCAGGCGCCGGCTCGAACAGGCGCATCAGGTGGTCAGGGGTGAAGTCGATGCAGCCGGTCCAGTGCATTTCCCACTc                                                                                         <  6:260088/147‑1 (MQ=255)
                                                                 aaaCGCCCGGCCAATCACGGTGCCGGTCGTGTTACCCCGCCCGTTGTACCCAGTCACCGCCACCAACCCAGGCGCCGGCTCGAACAGGCGCATCAGGGGGTCAGGGGTGAAGTCGATGCAGCCGGTCCAGTGCCTTTCCCCCGCGCCCt                                                                                    >  6:327775/1‑149 (MQ=255)
                                                                            cAATCACGGTGCCGGTCGTGTTACCCCGCCCGTTGTACCCAGTCACCGCCACCAACCCAGGCGCCGGCTCGAACAGGCGCATCAGGTGGTCAGGGGTGAAGTCGATGCAGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTcc                                                                         <  4:359162/149‑1 (MQ=255)
                                                                               tCACGGTGCCGGTCGTGTTACCCCGCCCGTTGTACCCAGTCACCGCCACCAACCCAGGCGCCGGCTCGAACAGGCGCATCAGGTGGTCAGGGGTGAAGTCGATGCAGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCggg                                                                      <  4:185816/149‑1 (MQ=255)
                                                                               tCACGGTGCCGGTCGTGTTACCCCGCCCGTTGTACCCAGTCACCGCCACCAACCCAGGCGCCGGCTCGAACAGGCGCATCAGGTGGTCAGGGGTGAAGTCGATGCAGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCggg                                                                      >  6:9210/1‑149 (MQ=255)
                                                                                  cGGTGCCGGTCGTGTTACCCCGCCCGTTGTACCCAGTCACCGCCACCAACCCAGGCGCCGGCTCGAACAGGCGCATCAGGTGGTCAGGGGTGAAGTCGATGCAGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGTaa                                                                   <  2:272969/149‑1 (MQ=255)
                                                                                  cGGTGCCGGTCGTGTTACCCCGCCCGTTGTACCCAGTCACCGCCACCAACCCAGGCGCCGGCTCGAACAGGCGCATCAGGTGGTCAGGGGTGAAGTCGATGCAGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGTaa                                                                   >  1:272969/1‑149 (MQ=255)
                                                                                       ccGGTCGTGTTACCCCGCCCGTTGTACCCAGTCACCGCCACCAACCCAGGCGCCGGCTCGAACAGGCGCATCAGGTGGTCAGGGGTGAAGTCGATGCAGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGTAATAGt                                                               <  8:84864/148‑1 (MQ=255)
                                                                                        cGGTCGTGTTACCCCGCCCGTTGTACCCAGTCACCGCCACCAACCCAGGCGCCGGCTCGAACAGGCGCATCAGGTGGTCAGGGGTGAAGTCGATGCAGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGTAATAGTgg                                                             <  2:9405/149‑1 (MQ=255)
                                                                                          gTCGTGTTACCCCGCCCGTTGTACCCAGTCACCGCCACCAACCCAGGCGCCGGCTCGAACAGGCGCATCAGGTGGTCAGGGGTGAAGTCGATGCAGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAg                                                                             <  1:198956/131‑1 (MQ=255)
                                                                                          gTCGTGTTACCCCGCCCGTTGTACCCAGTCACCGCCACCAACCCAGGCGCCGGCTCGAACAGGCGCATCAGGTGGTCAGGGGTGAAGTCGATGCAGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAg                                                                             >  2:198956/1‑131 (MQ=255)
                                                                                             gtgtTACCCCGCCCGTTGTACCCAGTCACCGCCACCAACCCAGGCGCCGGCTCGAACAGGCGCATCAGGTGGTCAGGGGTGAAGTCGATGCAGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGTAATAGTGGctct                                                         <  3:16673/148‑1 (MQ=255)
                                                                                                   ccccGCCCGTTGTACCCAGTCACCGCCACCAACCCAGGCGCCGGCTCGAACAGGCGCATCAGGTGGTCAGGGGTGAAGTCGATGCAGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGATGc                                                   >  6:131759/1‑148 (MQ=255)
                                                                                                    cccgcccgTTGTACCCAGTCACCGCCACCAACCCAGGCGCCGGCTCGAACAGGCGCATCAGGTGGTCAGGGGTGAAGTCGATGCAGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGATGCGa                                                 >  4:149653/1‑149 (MQ=255)
                                                                                                     ccgcccgTTGTACCCAGTCACCGCCACCAACCCAGGCGCCGGCTCGAACAGGCGCATCAGGTGGTCAGGGGTGAAGTCGATGCAGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGATGCGAt                                                >  7:302671/1‑149 (MQ=255)
                                                                                                                  ccAGTCACCGCCACCAACCCAGGCGCCGGCTCGAACAGGCGCATCAGGTGGTCAGGGGTGAAGTCGATGCAGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCggg                                                                      >  4:361910/1‑114 (MQ=255)
                                                                                                                  ccAGTCACCGCCACCAACCCAGGCGCCGGCTCGAACAGGCGCATCAGGTGGTCAGGGGTGAAGTCGATGCAGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCggg                                                                      <  3:361910/114‑1 (MQ=255)
                                                                                                                   cAGTCACCGCCACCAACCCAGGCGCCGGCTCGAACAGGCGCATCAGGTGGTCAGGGGTGAAGTCGATGCAGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCgg                                                                       >  6:135528/1‑112 (MQ=255)
                                                                                                                   cAGTCACCGCCACCAACCCAGGCGCCGGCTCGAACAGGCGCATCAGGTGGTCAGGGGTGAAGTCGATGCAGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCgg                                                                       <  5:135528/112‑1 (MQ=255)
                                                                                                                             caccaACCCAGGCGCCGGCTCGAACAGGCGCATCAGGTGGTCAGGGGTGAAGTCGATGCAGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGATGCGATCGGCCCAGCTGCGCACGAACCAgg                        >  6:85671/1‑149 (MQ=255)
                                                                                                                                   ctgatgcgcCTGTTCGAGCCGGCGCATCAGGTGGTCAGGGGTGAAGTCGATGCAGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGATGCGATCGGCCCAGCTGCGCACGAACCAggctggct                  >  4:132262/6‑149 (MQ=255)
                                                                                                                                       ggCGCCGGCTCGCACAGGCGCATCAGGTGGGCAGGGGGGAAGTCGCTGCAGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGATGCGAGCGGCCCAGCTGCGCACGAACCAGGCTGGCttgtt              <  8:347703/149‑1 (MQ=255)
                                                                                                                                          gcCGGCTCGAACAGGCGCATCAGGTGGTCAGGGGTGCAGTCGATGCAGCCGTTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGATGCGATCGGCCCAGCTGCGCACGAACCAGGCTGGCTTGTTGTc           <  4:30596/149‑1 (MQ=255)
                                                                                                                                                   aaCAGGCGCATCAGGTGGTCAGGGGTGAAGTCGATGCAGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGATGCGATCGGCCCAGCTGCGCACGAACCAGGCTGGCTTGTTGTCGACCCGGcc  <  5:131759/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ACACCGCTCATCGGCGAAAACGGAATCGGCAGGCTGTCCGCCTCGCCCTTGAGCAGGAACTCGGCAAACGCCCGGCCAATCACGGTGCCGGTCGTGTTACCCCGCCCGTTGTACCCAGTCACCGCCACCAACCCAGGCGCCGGCTCGAACAGGCGCATCAGGTGGTCAGGGGTGAAGTCGATGCAGCCGGTCCAGTGCATTTCCCACTCGACCTTGCCCAGCTCCGGGTAATAGTGGCTCTGGATGCGATCGGCCCAGCTGCGCACGAACCAGGCTGGCTTGTTGTCGACCCGGCC  >  NC_002947/6001087‑6001382

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