Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,364,922 A→G I231V (ATT→GTT)  puuB → gamma‑glutamylputrescine oxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,364,9220AG97.6% 133.7 / ‑4.2 41I231V (ATT→GTT) puuBgamma‑glutamylputrescine oxidase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (15/25);  total (15/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 = 9.06e-01

ATTGGCGAAGCGGACGCCATCCGCCTCAACGGCGGGCGCGTGTATGAACTTTCTGCCGTGACGCAGATCCAGCACACCACGCCAGCCGTTGTGCGAACTGCCAAAGGTCAGGTGACGGCGAAGTATGTGATTGTCGCCGGGAATGCGTATCTGGGCGATAAAGTAGAGCCGGAACTGGCGAAACGCAGCATGCCGTGCGGCACGCAGGTGATCACCACCGAACGGCTGTCGGAAGATTTAGCCCGTTCGCTGATCCCGAAAAACTA  >  NC_000913/1364793‑1365058
                                                                                                                                 |                                                                                                                                         
aTTGGCGAAGCGGACGCCATCCGCCTCAACGGCGGGCGCGTGTATGAACTTTCTGCCGTGACGCAGATCCAGCACACCACGCCAGCCGTTGTGCGAACTGCCAAAGGTCAGGTGACGGCGAAGTATGTGGTTGTCGCCgg                                                                                                                                <  1:410132/140‑1 (MQ=255)
 ttGGCGAAGCGGACGCCATCCGCCTCAACGGCGGGCGCGTGTATGAACTTTCTGCCGTGACGCAGATCCAGCACACCACGCCAGCCGTTGTGCGAACTGCCAAAGGTCAGGTGACGGCGAAGTATGTGGTTGTCGCCggg                                                                                                                               <  1:500301/140‑1 (MQ=255)
    gCGAAGCGGACGCCATCCGCCTCAACGGCGGGCGCGTGTATGAACTTTCTGCCGTGGCGCAGATCCAGCACACCACGCCCGCCGTTGTGCGAACTGGCCAAGGTCAGGTGACGGCGGAGTATTTGGTTGTTGCCCGGGaa                                                                                                                             >  2:473163/1‑140 (MQ=255)
    gCGAAGCGGACGCCATCCGCCTCAACGGCGGGCGCGTGTATGAACTTTCTGCCGTGACGCAGATCCAGCACACCACGCCAGCCGTTGTGCGAACTGCCAAAGGTCAGGTGACGGCGAAGTATGTGGTTGTCGCCGGGAAt                                                                                                                            <  1:767396/140‑1 (MQ=255)
             aCGCCATCCGCCTCAACGGCGGGCGCGTGTATGAACTTTCTGCCGTGACGCAGATCCAGCACACCACGCCAGCCGTTGTGCGAACTGCCAAAGGTCAGGTGACGGCGAAGTATGTGGTTGTCGCCGGGAATGCGTATCTg                                                                                                                   <  1:155722/140‑1 (MQ=255)
             aCGCCATCCGCCTCAACGGCGGGCGCGTGTATGAACTTTCTGCCGTGACGCAGATCCAGCACACCACGCCAGCCGTTGTGCGAACTGCCAAAGGTCAGGTGACGGCGAAGTATGTGGTTGTCGCCGGGAATGCGTATCTg                                                                                                                   <  2:393775/140‑1 (MQ=255)
             aCGCCATCCGCCTCAACGGCGGGCGCGTGAATGAACTTTCTGCCGTGACGCAGATCCAGCACACCACGCCAGCCGTTGTGCGAACTGCCAAAGGTCAGGTGACGGCGAAGTATGTGGTTGTCGCCGGGAATGCGTATCTg                                                                                                                   <  1:152062/140‑1 (MQ=255)
               gcaaTCCGCCTCAACGGCGGGCGCGTGTATGAACTTTCTGCCGTGACGCAGATCCAGCACACCCCGCCAGCCGTTGTGCGAACTGCCAAAGGTCAGGTGACGGCGAAGTATGTGGTTGTCGCCGGGAATGCGTATCTggg                                                                                                                 <  2:864528/137‑1 (MQ=255)
                     cGCCTCAACGGCGGGCGCGTGAATGAACTTTCTGCCGTGACGCAGATCCAGCACACCACGCCAGCCGTTGTGCGAACTGCCAAAGGGCAGGTGACGGCGAAGGATGTGGTTGTCGCCGGGGATGCGGATCTGGGCGGTaa                                                                                                           >  2:152062/1‑140 (MQ=255)
                           aaCGGCGGGCGCGTGTATGAACTTTCTGCCGTGACGCAGATCCAGCACACCACGCCAGCCGTTGTGCGAACTGCCAAAGGTCAGGTGACGGCGAAGTATGTGGTTGTCGCCGGGAATGCGTATCTggg                                                                                                                 >  1:349714/1‑128 (MQ=255)
                           aaCGGCGGGCGCGTGTATGAACTTTCTGCCGTGACGCAGATCCAGCACACCACGCCAGCCGTTGTGCGAACTGCCAAAGGTCAGGTGACGGCGAAGTATGTGGTTGTCGCCGGGAATGCGTATCTggg                                                                                                                 <  2:349714/128‑1 (MQ=255)
                                          tATGAACTTTCTGCCGTGACGCAGATCCAGCACACCACGCCAGCCGTTGTGCGAACTGCCAAAGGTCAGGTGACGGCGAAGTATGTGGTTGTCGCCGGGAATGCGTATCTGGGCGATAAAGTAGAGCCGGAACTGGCGaa                                                                                      <  2:88515/140‑1 (MQ=255)
                                                   tCTGCCGTGACGCAGATCCAGCACACCACGCCAGCCGTTGTGCGAACTGCCAAAGGTCAGGTGACGGCGAAGTATGTGGTTGTCGCCGGGAATGCGTATCTGGGCGATAAAGTAGAGCCGGAACTGGCGAAACgcagc                                                                               >  1:541166/1‑138 (MQ=255)
                                                   tCTGCCGTGACGCAGATCCAGCACACCACGCCAGCCGTTGTGCGAACTGCCAAAGGTCAGGTGACGGCGAAGTATGTGGTTGTCGCCGGGAATGCGTATCTGGGCGATAAAGTAGAGCCGGAACTGGCGAAACgcagc                                                                               <  2:541166/138‑1 (MQ=255)
                                                   tCTGCCGTGACGCAGATCCAGCACACCACGCCAGCCGTTGTGCGAACTGCCAAAGGTCAGGTGACGGCGAAGTATGTGGTTGTCGCCGGGAATGCGTATCTGGGCGATAAAGTAGAGCCGGAACTGGCGAAACGCAGCAt                                                                             >  1:98243/1‑140 (MQ=255)
                                                   tCTGCCGTGACGCAGATCCAGCACACCACGCCAGCCGTTGTGCGAACTGCCAAAGGTCAGGTGACGGCGAAGTATGTGGTTGTCGCCGGGAATGCGTATCTGGGCGATAAAGTAGAGCCGGAACTGGCGAAACGCAGCAt                                                                             >  1:692551/1‑140 (MQ=255)
                                                         gTGACGCAGATCCAGCACACCACGCCAGCCGTTGTGCGAACTGCCAAAGGTCAGGTGACGGCGAAGTATGTGATTGTCGCCGGGAATGCGTATCTGGGCGATAAAGTAGAGCCGGAACTGGCGAAACGCAGCATGCCGTg                                                                       <  1:473163/140‑1 (MQ=255)
                                                           gACGCAGATCCAGCACACCACGCCAGCCGTTGTGCGAACTGCCAAAGGTCAGGTGACGGCGAAGTATGTGGTTGTCGCCGGGAATGCGTATCTGGGCGATAAAGTAGAGCCGGAACTGGCGAAACGCAGCATGCCGTGCg                                                                     <  2:98243/140‑1 (MQ=255)
                                                              gCAGATCCAGCACACCACGCCAGCCGTTGTGCGAACTGCCAAAGGTCAGGTGACGGCGAAGTATGTGGTTGTCGCCGGGAATGCGTATCTGGGCGATAAAGTAGAGCCGGAACTGGCGAAACGCAGCATGCCGTGCGGCa                                                                  <  2:597524/140‑1 (MQ=255)
                                                                 gATCCAGCACACCACGCCAGCCGTTGTGCGAACTGCCAAAGGTCAGGTGACGGCGAAGTATGTGGTTGTCGCCGGGAATGCGTATCTGGGCGATAAAGTAGAGCCGGAACTGGCGAAACGCAGCATGCCGTGCGGCACGc                                                               <  1:17108/140‑1 (MQ=255)
                                                                 gATCCAGCACACCACGCCAGCCGTTGTGCGAACTGCCAAAGGTCAGGTGACGGCGAAGTATGTGGTTGTCGCCGGGAATGCGTATCTGGGCGATAAAGTAGAGCCGGAACTGGCGAAACGCAGCATGCCGTGCGGCACGc                                                               <  1:650012/140‑1 (MQ=255)
                                                                    ccAGCACACCACGCCAGCCGTTGTGCGAACTGCCAAAGGTCAGGTGACGGCGAAGTATGTGGTTGTCGCCGGGAATGCGTATCTGGGCGATAAAGTAGAGCCGGAACTGGCGAAACGCAGCATGCCGTGCGGCACGCAgg                                                            >  1:109627/1‑140 (MQ=255)
                                                                          caccacGCCAGCCGTTGTGCGAACTGCCAAAGGTCAGGTGACGGCGAAGTATGTGGTTGTCGCCGGGAATGCGTATCTGGGCGATAAAGTAGAGCCGGAACTGGCGAAACGCAGCATGCCGTGCGGCACGCAGGTGATca                                                      >  1:61327/1‑140 (MQ=255)
                                                                          caccacGCCAGCCGTTGTGCGAACTGCCAAAGGTCAGGTGACGGCGAAGTATGTGGTTGTCGCCGGGAATGCGTATCTGGGCGATAAAGTAGAGCCGGAACTGGCGAAACGCAGCATGCCGTGCGGCACGCAGGTGATca                                                      <  1:306520/140‑1 (MQ=255)
                                                                            ccacGCCAGCCGTTGTGCGAACTGCCAAAGGTCAGGTGACGGCGAAGTATGTGGTTGTCGCCGGGAATGCGTATCTGGGCGATAAAGTAGAGCCGGAACTGGCGAAACGCAGCATGCCGTGCGGCACGCAGGTGATcacc                                                    >  2:567638/1‑140 (MQ=255)
                                                                             cacGCCAGCCGTTGTGCGAACTGCCAAAGGTCAGGTGACGGCGAAGTATGTGGTTGTCGCCGGGAATGCGTATCTGGGCGATAAAGTAGAGCCGGAACTGGCGAAACGCAGCATGCCGTGCGGCa                                                                  <  1:825720/125‑1 (MQ=255)
                                                                             cacGCCAGCCGTTGTGCGAACTGCCAAAGGTCAGGTGACGGCGAAGTATGTGGTTGTCGCCGGGAATGCGTATCTGGGCGATAAAGTAGAGCCGGAACTGGCGAAACGCAGCATGCCGTGCGGCa                                                                  >  2:825720/1‑125 (MQ=255)
                                                                             cacGCCAGCCGTTGTGCGAACTGCCAAAGGTCAGGTGACGGCGAAGTATGTGGTTGTCGCCGGGAATGCGTATCTGGGCGATAAAGTAGAGCCGGAACTGGCGAAACGCAGCATGCCGTGCGGCACGc                                                               <  2:841985/128‑1 (MQ=255)
                                                                             cacGCCAGCCGTTGTGCGAACTGCCAAAGGTCAGGTGACGGCGAAGTATGTGGTTGTCGCCGGGAATGCGTATCTGGGCGATAAAGTAGAGCCGGAACTGGCGAAACGCAGCATGCCGTGCGGCACGc                                                               >  1:841985/1‑128 (MQ=255)
                                                                             cacGCCAGCCGTTGTGCGAACTGCCAAAGGTCAGGTGACGGCGAAGTATGTGGTTGTCGCCGGGAATGCGTATCTGGGCGATAAAGTAGAGCCGGAACTGGCGAAACGCAGCATGCCGTGCGGCACGCAGGTGATcacca                                                   >  2:809789/1‑140 (MQ=255)
                                                                              acGCCAGCCGTTGTGCGAACTGCCAAAGGTCAGGTGACGGCGAAGTATGTGGTTGTCGCCGGGAATGCGTATCTGGGCGATAAAGTAGAGCCGGAACTGGCGAAACGCAGCATGCCGTGCGGCACGCAGGTGATcaccac                                                  >  2:642418/1‑140 (MQ=255)
                                                                              acGCCAGCCGTTGTGCGAACTGCCAAAGGTCAGGTGACGGCGAAGTATGTGGTTGTCGCCGGGAATGCGTATCTGGGCGATAAAGTAGAGCCGGAACTGGCGAAACGCAGCATGCCGTGCGGCACGCAGGTGATcac                                                     <  1:679131/137‑1 (MQ=255)
                                                                              acGCCAGCCGTTGTGCGAACTGCCAAAGGTCAGGTGACGGCGAAGTATGTGGTTGTCGCCGGGAATGCGTATCTGGGCGATAAAGTAGAGCCGGAACTGGCGAAACGCAGCATGCCGTGCGGCACGCAGGTGATcac                                                     >  2:679131/1‑137 (MQ=255)
                                                                                     ccGTTGTGCGAACTGCCAAAGGTCAGGTGACGGCGAAGTATGTGGTTGTCGCCGGGAATGCGTATCTGGGCGATAAAGTAGAGCCGGAACTGGCGAAACGCAGCATGCCGTGCGGCACGCAGGTGATCACCACCGAACgg                                           <  1:690361/140‑1 (MQ=255)
                                                                                        ttGTGCGAACTGCCAAAGGTCAGGTGACGGCGAAGTATGTGGTTGTCGCCGGGAATGCGTATCTGGGCGATAAAGTAGAGCCGGAACTGGCGAAACGCAGCATGCCGTGCGGCACGCAGGTGATCACCACCGAACGGCTg                                        <  1:127970/140‑1 (MQ=255)
                                                                                            gCGAACTGCCAAAGGTCAGGTGACGGCGAAGTATGTGGTTGTCGCCGGGAATGCGTATCTGGGCGATAAAGTAGAGCCGGAACTGGCGAAACGCAGCATGCCGTGCGGCACGCAGGTGATCACCACCGAACGGCTGTCgg                                    <  1:50807/140‑1 (MQ=255)
                                                                                               aaCTGCCAAAGGTCAGGTGACGGCGAAGTATGTGGTTGTCGCCGGGAATGCGTATCTGGGCGATAAAGTAGAGCCGGAACTGGCGAAACGCAGCATGCCGTGCGGCACGCAGGTGATCaccacc                                                 <  1:519318/124‑1 (MQ=255)
                                                                                               aaCTGCCAAAGGTCAGGTGACGGCGAAGTATGTGGTTGTCGCCGGGAATGCGTATCTGGGCGATAAAGTAGAGCCGGAACTGGCGAAACGCAGCATGCCGTGCGGCACGCAGGTGATCaccacc                                                 >  2:519318/1‑124 (MQ=255)
                                                                                               aaCTGCCAAAGGTCAGGTGACGGCGAAGTATGTGGTTGTCGCCGGGAATGCGTATCTGGGCGATAAAGTAGAGCCGGAACTGGCGAAACGCAGCATGCCGTGCGGCACGCAGGTGATCACCACCGAACGGCTGTCGGAAg                                 >  2:72888/1‑140 (MQ=255)
                                                                                                  tGCCAAAGGTCAGGTGACGGCGAAGTATGTGGTTGTCGCCGGGAATGCGTATCTGGGCGATAAAGTAGAGCCGGAACTGGCGAAACGCAGCATGCCGTGCGGCACGCAGGTGATCACCACCGAACGGCTGTCgg                                    >  2:829286/1‑134 (MQ=255)
                                                                                                  tGCCAAAGGTCAGGTGACGGCGAAGTATGTGGTTGTCGCCGGGAATGCGTATCTGGGCGATAAAGTAGAGCCGGAACTGGCGAAACGCAGCATGCCGTGCGGCACGCAGGTGATCACCACCGAACGGCTGTCgg                                    <  1:829286/134‑1 (MQ=255)
                                                                                                      aaaGGTCAGGTGACGGCGAAGTATGTGGTTGTCGCCGGGAATGCGTATCTGGGCGATAAAGTAGAGCCGGAACTGGCGAAACGCAGCATGCCGTGCGGCACGCAGGTGATCACCACCGAACGGCTGTCGGAAGATTTAGc                          <  2:82861/140‑1 (MQ=255)
                                                                                                                            aTGTGGTTGTCGCCGGGAATGCGTATCTGGGCGATAAAGTAGAGCCGGAACTGGCGAAACGCAGCATGCCGTGCGGCACGCAGGTGATCACCACCGAACGGCTGTCGGAAGATTTAGCCCGTTCGCTGATCCCGAAAAac    <  2:611866/140‑1 (MQ=255)
                                                                                                                              gtgGTTGTCGCCGGGAATGCGTATCTGGGCGATAAAGTAGAGCCGGAACTGGCGAAACGCAGCATGCCGTGCGGCACGCAGGTGATCACCACCGAACGGCTGTCGGAAGATTTAGCCCGTTCGCTGATCCCGAAAAacta  <  2:692551/140‑1 (MQ=255)
                                                                                                                                 |                                                                                                                                         
ATTGGCGAAGCGGACGCCATCCGCCTCAACGGCGGGCGCGTGTATGAACTTTCTGCCGTGACGCAGATCCAGCACACCACGCCAGCCGTTGTGCGAACTGCCAAAGGTCAGGTGACGGCGAAGTATGTGATTGTCGCCGGGAATGCGTATCTGGGCGATAAAGTAGAGCCGGAACTGGCGAAACGCAGCATGCCGTGCGGCACGCAGGTGATCACCACCGAACGGCTGTCGGAAGATTTAGCCCGTTCGCTGATCCCGAAAAACTA  >  NC_000913/1364793‑1365058

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: