Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 894,155 T→G D124E (GAT→GAG potF → putrescine ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913894,1550TG97.7% 145.2 / ‑3.5 43D124E (GAT→GAGpotFputrescine ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (19/23);  total (19/24)
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.14e-01

GCGGGAGTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGATCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGAAAACGCGC  >  NC_000913/894066‑894249
                                                                                         |                                                                                              
gCGGGAGTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATa                                                                                      >  1:484166/1‑100 (MQ=255)
   ggAGTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAAt                                                                                   >  1:1011126/1‑100 (MQ=255)
   ggAGTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAAt                                                                                   >  1:773724/1‑100 (MQ=255)
      gTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTg                                                                                >  1:229876/1‑100 (MQ=255)
       ttttCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGc                                                                               >  1:310144/1‑100 (MQ=255)
        tttCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCt                                                                              <  2:1011117/100‑1 (MQ=255)
         ttCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTa                                                                             >  2:830116/1‑100 (MQ=255)
            cAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGc                                                                          >  2:1178714/1‑100 (MQ=255)
              gCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGccc                                                                        <  1:211923/100‑1 (MQ=255)
               ccGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCt                                                                       <  2:1297429/100‑1 (MQ=255)
                    ggACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATAtg                                                                  <  1:451040/100‑1 (MQ=255)
                       cAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTgg                                                               <  2:470046/100‑1 (MQ=255)
                       cAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTgg                                                               >  1:953375/1‑100 (MQ=255)
                         aaaGCAGATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGc                                                             <  2:313560/100‑1 (MQ=255)
                           aGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACCCGAGCCCGACAATAAATTTGCTATGCCCTATCTGTGGGcga                                                           >  1:1035989/1‑100 (MQ=255)
                               aaTTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCgacgac                                                       <  2:493963/100‑1 (MQ=255)
                                 ttGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGATCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCg                                                     <  1:907697/100‑1 (MQ=255)
                                       gAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTg                                               <  2:527179/100‑1 (MQ=255)
                                           ggAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCta                                           <  1:901933/100‑1 (MQ=255)
                                            gaagacTCTCGTACCGGCACTGCTGACGCTGGTTGCCAAACACGCGCCCGAGACTCAATTTGCTCTGCTCTATATTTGGGCGACGACCGGGATTGGc                                             <  2:107518/97‑1 (MQ=255)
                                            gaagaaTCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCtat                                          <  1:618769/100‑1 (MQ=255)
                                              agaaTCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATaa                                        >  2:306711/1‑100 (MQ=255)
                                              agaaTCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCAACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTCGCtata                                         >  2:614316/1‑99 (MQ=255)
                                                   ctcGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTg                                   <  1:1178725/100‑1 (MQ=255)
                                                     cGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGAt                                 <  1:781000/100‑1 (MQ=255)
                                                         ccGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAg                             <  2:474544/100‑1 (MQ=255)
                                                          cGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGt                            <  2:551804/100‑1 (MQ=255)
                                                            gAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTa                          <  1:68563/100‑1 (MQ=255)
                                                              aCTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTaaa                        <  1:314603/100‑1 (MQ=255)
                                                              aCTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTaaa                        <  1:173640/100‑1 (MQ=255)
                                                                  ctgAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCgg                    >  2:645166/1‑100 (MQ=255)
                                                                   tgAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGt                   >  2:829135/1‑100 (MQ=255)
                                                                    gAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTg                  >  1:186418/1‑100 (MQ=255)
                                                                        cTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCAGTGCTgg              >  2:484024/1‑100 (MQ=255)
                                                                          ggTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGc            >  1:640250/1‑100 (MQ=255)
                                                                          ggTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGc            >  1:840776/1‑100 (MQ=255)
                                                                           gTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCg           <  1:320059/100‑1 (MQ=255)
                                                                            tCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGa          >  2:451320/1‑100 (MQ=255)
                                                                            tCGCCAAACACGAGCCCGACAATAAAGTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGa          <  2:846424/100‑1 (MQ=255)
                                                                             cGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGa          <  2:415734/99‑1 (MQ=255)
                                                                              gCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGaaa        >  2:512661/1‑100 (MQ=255)
                                                                                 aaaCACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCgcgtttt     <  2:240823/100‑7 (MQ=255)
                                                                                    cacGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGAAAACgcgc  <  1:924016/100‑1 (MQ=255)
                                                                                         |                                                                                              
GCGGGAGTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGATCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGAAAACGCGC  >  NC_000913/894066‑894249

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