Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 897,409 C→T 8.3% P109L (CCG→CTG)  potI → putrescine ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913897,4090CT 8.3% 139.7 / 3.1 49P109L (CCG→CTG) potIputrescine ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (20/24);  new base T (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 = 6.32e-01

GATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTATGATCACCGCGCCGCTGGTAATGCCAGATGTCATCACGGGCTTGTCGCTGTTGTTATTATTCGTCGCGCTTGCTCATGCCATTGGCTGGCCTGCGGACCGCGGTATGCTCACCATCTGGCTGGCGCATGTCACTTTTTGTACGGCTTATGTGGCGGT  >  NC_000913/897266‑897553
                                                                                                                                               |                                                                                                                                                
gatgaGTGCGGTTGGTTTAAGGCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTATGATCACCGCGCCGCTg                                                                                                                                              <  4:25242/148‑1 (MQ=255)
 atgaGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTATGATCACCGCGCCGCTgg                                                                                                                                             <  8:117770/148‑1 (MQ=255)
     gTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTATGATCACCGCGCCGCTGGTAATg                                                                                                                                        >  7:132807/1‑149 (MQ=255)
      tGCGGTTGGTTGAAGTCTGACAATTGCGGCCTGTGCTGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCGTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTATGATCACCGCGCCGCTGGTAATGc                                                                                                                                       <  4:282580/149‑1 (MQ=255)
       gCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTATGATCACCGCGCCGCTGGTAATGc                                                                                                                                       >  6:193544/1‑148 (MQ=255)
          gttggttTAAGTCTTACAATTGCGTCCTGTGCGGCGACGTCGTCGGCGATCCTCTGAACTAGTGCGTCGGTGGTGCTGGTGAGCTTTGGCAGGTTTCGAGGATCAAATGGCTTTTCCTTTATGATCACCGCGCCGCTGGTAATGCCAGa                                                                                                                                   <  1:194898/149‑1 (MQ=255)
           ttggttTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTATGATCACCGCGCCGCTGGTAATGCCAGAt                                                                                                                                  <  4:138076/149‑1 (MQ=255)
              gttTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTATGATCACCGCGCCGCTGGTAATGCCAGATGt                                                                                                                                >  1:178274/1‑148 (MQ=255)
                     tCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGCGGATCAAACGGCTTTGCCTTTATGATCACCGCGCCGCTGGTAATGCCAGATGTCATCCCgg                                                                                                                        >  7:358986/1‑149 (MQ=255)
                         aCAATTGCGTCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTATGATCACCGCGCCGCTGGTAATGCCAGAt                                                                                                                                  <  2:324252/135‑1 (MQ=255)
                         aCAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTATGATCACCGCGCCGCTGGTAATGCCAGAt                                                                                                                                  >  1:324252/1‑135 (MQ=255)
                          cAATTGCGGCCTGTGCGGCGACGGCGGCTGCGATCCTCGGAACTATAGCGGCGGAGGTGCTGGTGCGCTTTGGCAGTTTTCGCGGATCAAATGGCTTTGCCATAATGATCACCGCGCAGCTCGTAATGCCAGATGTCAACACGGGCTTg                                                                                                                   >  4:247514/1‑149 (MQ=255)
                             ttGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTATGATCACCGCGCCGCTGGTAATGCCAGATGTCATCACGGGCTTGTCg                                                                                                                >  7:1863/1‑149 (MQ=255)
                                       tgCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTATGATCACCGCGCCGCTGGTAATGCCAGATGTCATCACGGGCTTGTCGCTGTTGttat                                                                                                      <  8:1863/149‑1 (MQ=255)
                                         cggcgACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGCGGATCAAATGGGTTTGCCTTTATGATCACCGCGCCGCTGGTAATGCCAGATGTCATCACGGGCTTGTCGCTGTTGttatta                                                                                                    <  2:178274/149‑1 (MQ=255)
                                         cggcgACGGCGGAGGCGATCCTCGGAACAATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTATCGCGGATCAAATGGCTTTGCCTTTATGATCACCGCGCCGCTGGTAATGCCAGATGTCATCACTGGCTTGTCGCTGTTGttatta                                                                                                    >  8:52472/1‑149 (MQ=255)
                                               cggcggcggCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTATGATCACCGCGCCGCTGGTAATGCCAGATGTCATCACGGGCTTGTCGCTGTTGTTATTATtcgtc                                                                                              <  8:132807/149‑1 (MQ=255)
                                                      ggcgATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTATGATCACCGCGCCGCTGGTAATGCCAGATGTCATCACGGGCTTGTCGCTGTTGTTATTATTCGTCGCGCTTg                                                                                       >  1:168010/1‑149 (MQ=255)
                                                             cTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTATGATCACCGCGCCGCTGGTAATGCCAGATGTCATCACGGGCTTGTCGCTGTTGTTATTATTCGTCGCGCTTGCTCATGc                                                                                >  3:50101/1‑149 (MQ=255)
                                                                       ttGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTATGATCACCGCGCCGCTGGTAATGCCAGATGTCATCACGGGCTTGTCGCTGTTGTTATTATTCGTCGCGCTTGCTCATg                                                                                 <  1:385742/138‑1 (MQ=255)
                                                                       ttGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTATGATCACCGCGCCGCTGGTAATGCCAGATGTCATCACGGGCTTGTCGCTGTTGTTATTATTCGTCGCGCTTGCTCATg                                                                                 >  2:385742/1‑138 (MQ=255)
                                                                       ttGCGGCGGGGGTGCGGGTGCGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTATGATCACCGCGCCGCTGGTAATGCCAGATGTCATCACGGGCTTGTCGCTGTTGTTATTATTCGTCGCGCTTGCTCATGCCATtggctg                                                                       <  8:353769/148‑1 (MQ=255)
                                                                         gcggcgGTGGTGCTGGTGCGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTATGATCACCGCGCCGCTGGTAATGCCAGATGTCATCACGGGCTTGTCGCTGTTGTTATTATTCGTCGCGCTTGCTCATGCCATtggctggc                                                                     <  7:37622/148‑1 (MQ=255)
                                                                         gcggcgGTGGTGCTGGTGCGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTATGATCACCGCGCCGCTGGTAATGCCAGATGTCATCACGGGCTTGTCGCTGTTGTTATTATTCGTCGCGCTTGCTCATGCCATTGGCTGGcc                                                                    >  7:292293/1‑149 (MQ=255)
                                                                             cggTGGTGCTGGTGCGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTATGATCACCGCGCCGCTGGTAATGCCAGATGTCATCACGGGCTTGTCGCTGTTGTTATTATTCGTCGCGCTTGCTCATGCCATTGGCTGGCCTGCg                                                                <  8:107871/149‑1 (MQ=255)
                                                                                     ctggtgcGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTATGATCACCGCGCCGCTGGTAATGCCAGATGTCATCACGGGCTTGTCGCTGTTGTTATTATTCGTCGCGCTTGCTCATGCCATTGGCTGGCCTGCGGACcgcg                                                         <  8:274723/148‑1 (MQ=255)
                                                                                     ctggtgcGCTTTGGCAGGTTTCGCGGATCAAATGGCTGTGCCTTTATGATCACCGCGCCGCGGGTAATGCCAGATGTCATCACGGGCTTGTCGCTGTTGTTATTATTCGTCGCGCTTGCTCATGCCATTGGCTGGCCTGCGGACcgcg                                                         <  3:80795/148‑1 (MQ=255)
                                                                                     ctggtgaGCTTTGGCAGGTTTCGCGGATCAAATGGCTGTGCCTTTATGATCACCGCGCCGCTGGTAATGCCAGATGTCATCACGGGCTTGTCGCTGTTGTTATTATTCGTCGCGCTTGCTCATGCCATTGGCTGGCCTGCGGACcgcg                                                         <  3:247514/148‑1 (MQ=255)
                                                                                       ggtgcGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTATGATCACCGCGCCGCTGGTAATGCCAGATGTCATCACTGGc                                                                                                                      <  7:321261/85‑1 (MQ=255)
                                                                                       ggtgcGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTATGATCACCGCGCCGCTGGTAATGCCAGATGTCATCACGGGc                                                                                                                      >  8:321261/1‑85 (MQ=255)
                                                                                         tgcGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTATGATCACCGCGCCGCTGGTAATGCCAGATGTCATCACGGGCTTGTCGCTGTTGTTATTATTCGTCGCGCTTGCTCATGCCATTGGCTGGCCTGCGGACCGCGGTATg                                                    <  4:194860/149‑1 (MQ=255)
                                                                                            gcTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTATGATCACCGCGCCGCTGGTAATGCCAGATGTCATCACGGGCTTGTCGCTGTTGTTATTATTCGTCGCGCTTGCTCATGCCATTGGCTGGCCTGCGGACCGCGGTATGCTc                                                 <  7:52472/149‑1 (MQ=255)
                                                                                                  gCAGGTTTCGCGGATCAAATGGCTTTGCCTTTATGATCACCGCGCCGCTGGTAATGCCAGATGTCATCACggg                                                                                                                       <  8:292767/73‑1 (MQ=255)
                                                                                                  gCAGGTTTCGCGGATCAAATGGCTTTGCCTTTATGATCACCGCGCCGCTGGTAATGCCAGATGTCATCACggg                                                                                                                       >  7:292767/1‑73 (MQ=255)
                                                                                                  gCAGGTTTCGCGGATCAAATGGCTTTGCCTTTATGATCACCGCGCCGCTGGTAATGCCAGATGTCATCACGGGCTTGTCGCTGTTGTTATTATTCGTCGCGCTTGCTCATGCCATTGGCTGGCCTGCGGACCGCGGTATGCTCACCATc                                           >  7:306460/1‑149 (MQ=255)
                                                                                                       tttCGCGGATCAAATGGCTTTGCCTTTATCATCACCGAGACGCTGGTAATGACAGATGTCATCACTGGCTTGTCGCTGTTGTTATTATTCGTCGCTCTTGCTCATGCCATTGGCTGGCCTGCGGACCTCGGTAAGCTCACCATctggct                                      >  3:163373/1‑149 (MQ=255)
                                                                                                           gctgATCAAATGGCTTTGCCTTTATGATCACCGCGCCGCTGGTAATGCCAGATGTCATCACGGGCTTGTCGCTGTTGTTATTATTCGACGCGCTTGCTCATGCCATTGGCTGGCCTGCGGACCGCGGTATGCTCACCATCtggctggc                                   <  8:292293/145‑1 (MQ=255)
                                                                                                              gATCAAATGGCTTTGCCTTTATGATCACCGCGCCGCTGGTAATGCCAGATGTCATCACGGGCTTGTCGCTGTTGTTATTATTCGTCGCGCTTGCTCATGCCATTGGCTGGCCTGCGGACCGCGGTATGCTCACCATCTGGCTGGCGCAt                               >  8:216139/1‑149 (MQ=255)
                                                                                                                 cAAATGGCTTTGCCTTTATGATCACCGCGCCGCTGGTAATGCCAGATGTCATCACGGGCTTGTCGCTGTTGTTATTATTCGTCGCGCTTGCTCATGCCATTGGCTGGCCTGCGGACCGCGGTATGCTCACCATCTGGCTGGCGCATGt                             >  7:360677/1‑148 (MQ=255)
                                                                                                                  aaaTGGCTATGCCTTTATGATCACCGCGCCGCTGGTAATGCCAGATGTCATCACGGGCTTGTCGCTGTTGTTATTATTCGTCGCGCTTGCTCATGCCATTGGCTGGCCTGCGGACCGCGGTATGCTCACCATCTGGCTGGCGCATGTCa                           >  2:251645/1‑149 (MQ=255)
                                                                                                                     tGGCTTTGCCTTTATGATCACCGCGCCGCTGGTAATGCCAGATGTCATCACGGGCTTGTCGCTGTTGTTATTATTCGTCGCGCTTGCTCATGCCATTGGCTGGCCTGCGGACCGCGGTATGCTCACCATCTGGCTGGCGCATGTCACtt                        >  4:305708/1‑149 (MQ=255)
                                                                                                                       gCTTTGCCTTTATGATCACCGCGCCGCTGGTAATGCCAGATGTCATCACGGGCTTGTCGCTGTTGTTATTATTCGTCGCGCTTGCTCATGCCATTGGCTGGCCTGCGGACCGCGGTATGCTCACCATCTGGCTGGCGCATGTCACttt                       <  5:370534/148‑1 (MQ=255)
                                                                                                                                    gATCACCGCGCCGCTGGTAATGCCAGATGTCATCACGGGCTTGTCGCTGTTGTTATTATTCGTCGCGCTTGCTCATGCCATTGGCTGGCCTGCGGACCGCGGTATGCTCACCATCTGGCTGGCGCATGTCACTTTTTGTACGGCTTAtg         >  3:112049/1‑149 (MQ=255)
                                                                                                                                        aCCGCGCTGCTGGTAATGCCAGATGTCATCACGGGCTTGTCGCTGTTGTTATTATTCGTCGCGCTTGCTCATGCCATTGGCTGGCCTGCGGACCGCGGTATGCTCACCATCTGGCTGgcgc                                 <  1:105265/121‑1 (MQ=255)
                                                                                                                                        aCCGCGCTGCTGGTAATGCCAGATGTCATCACGGGCTTGTCGCTGTTGTTATTATTCGTCGCGCTTGCTCATGCCATTGGCTGGCCTGCGGACCGCGGTATGCTCACCATCTGGCTGgcgc                                 >  2:110427/1‑121 (MQ=255)
                                                                                                                                        aCCGCGCTGCTGGTAATGCCAGATGTCATCACGGGCTTGTCGCTGTTGTTATTATTCGTCGCGCTTGCTCATGCCATTGGCTGGCCTGCGGACCGCGGTATGCTCACCATCTGGCTGgcgc                                 >  2:105265/1‑121 (MQ=255)
                                                                                                                                        aCCGCGCTGCTGGTAATGCCAGATGTCATCACGGGCTTGTCGCTGTTGTTATTATTCGTCGCGCTTGCTCATGCCATTGGCTGGCCTGCGGACCGCGGTATGCTCACCATCTGGCTGgcgc                                 <  1:110427/121‑1 (MQ=255)
                                                                                                                                          cgcgCCGCTGGTAATGCCAGATGTCATCACGGGCTTGTCGCTGTTGTTATTATTCGTCGCGCTTGCTCATGCCATTGGCTGGCCTGCGGACCGCGGTATGCTCACCATCTGGCTGGCGCATGTCACTTTTTGTACGGCTTATGTGGCgg   <  4:50101/149‑1 (MQ=255)
                                                                                                                                           gcgcCGCTGGTAATGCCAGATGTCATCACGGGCTTGTCGCTGTTGTTATTATTCGTCGCGCTTGCTCATGCCATTGGCTGGCCTGCGGACCGCGGTATGCTCACCATCTGGCTGGCGCATGTCACTTTTTGTACGGCTTATGTGGCGGt  <  2:168010/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                
GATGAGTGCGGTTGGTTTAAGTCTGACAATTGCGGCCTGTGCGGCGACGGCGGCGGCGATCCTCGGAACTATTGCGGCGGTGGTGCTGGTGCGCTTTGGCAGGTTTCGCGGATCAAATGGCTTTGCCTTTATGATCACCGCGCCGCTGGTAATGCCAGATGTCATCACGGGCTTGTCGCTGTTGTTATTATTCGTCGCGCTTGCTCATGCCATTGGCTGGCCTGCGGACCGCGGTATGCTCACCATCTGGCTGGCGCATGTCACTTTTTGTACGGCTTATGTGGCGGT  >  NC_000913/897266‑897553

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