Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 894,155 T→G 100% 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,1550TG95.6% 227.3 / 0.5 68D124E (GAT→GAGpotFputrescine ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base G (31/34);  total (32/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.22e-01

TGACTGCGGGAGTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGATCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGAAAACGCGCCC  >  NC_000913/894061‑894251
                                                                                              |                                                                                                
tgccTGCGGGAGTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCAGAAGCTGGTCGCCAAACACGAGCCCGa                                                                                             >  2:487589/4‑100 (MQ=255)
   cTGCGGGATTGTTCCAGCCGCTGGACAAAAGCAAATTGTCGGAGTGGAAGAATCTCGACACGGAACTGCTGAAGCTGGCTGCCAAACACGAGCCCGa                                                                                             <  2:1261670/97‑1 (MQ=255)
   cTGCGGGAGTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAATCTGGTCGCTAAACACGAGCCCg                                                                                              <  2:2139841/96‑1 (MQ=255)
   cTGCGGGAGTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACaa                                                                                          >  2:119421/1‑100 (MQ=255)
     gCGGGAGTTTTTCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATa                                                                                        >  2:1913308/1‑100 (MQ=255)
     gCGGGAGTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATa                                                                                        >  1:515606/1‑100 (MQ=255)
        ggAGTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGTGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAAt                                                                                     <  2:347891/100‑1 (MQ=255)
        ggAGTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGACGCTGGTCGCCAAACACGAGCCCGACAATAAAt                                                                                     <  1:1646241/100‑1 (MQ=255)
         gAGTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAAtt                                                                                    >  2:620957/1‑100 (MQ=255)
           gTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAAttt                                                                                   >  1:1010023/1‑99 (MQ=255)
           gTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTg                                                                                  <  1:1833596/100‑1 (MQ=255)
            ttttCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGc                                                                                 <  1:1830242/100‑1 (MQ=255)
             tttCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCt                                                                                <  1:1837065/100‑1 (MQ=255)
              ttCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTa                                                                               <  1:1041848/100‑1 (MQ=255)
               tCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTAt                                                                              >  1:15716/1‑100 (MQ=255)
                   gCCGCTGTACAAAAGCAAATTGCCGGAGTTGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGc                                                                            <  1:18745/98‑1 (MQ=255)
                     cGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCta                                                                        >  2:788149/1‑100 (MQ=255)
                      gCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCtat                                                                       >  2:694125/1‑100 (MQ=255)
                       cTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCtata                                                                      <  1:358987/100‑1 (MQ=255)
                           aCAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATAtgtg                                                                  <  2:561027/100‑1 (MQ=255)
                           aCAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATAtgtg                                                                  <  2:382443/100‑1 (MQ=255)
                            cAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTgg                                                                 <  2:846465/100‑1 (MQ=255)
                             aaaaGCAAATTGCCTGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGATCCCGACAATAAATTTGCTATGCCCTATATGTggg                                                                <  2:767622/100‑1 (MQ=255)
                             aaaaGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTggg                                                                <  1:2245799/100‑1 (MQ=255)
                              aaaGCAAATTGTCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGc                                                               >  2:387950/1‑100 (MQ=255)
                              aaaGCAAATTGCCGGAGTGGAAGAATCTCGATCCGGAGCTGCTGAAGTTGGTCGCCAAACACGATCCCGACAATAAATTTGCTATGCCCTATATGTGGGc                                                               <  1:1833464/100‑1 (MQ=255)
                              aaaGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGc                                                               >  2:1716371/1‑100 (MQ=255)
                              aaaGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGc                                                               >  2:2090400/1‑100 (MQ=255)
                                aGCAAATTGCCGGAGTGGAAGAATCTCTACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGcga                                                             <  2:2180381/100‑1 (MQ=255)
                                aGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGcga                                                             >  1:509954/1‑100 (MQ=255)
                                        gCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGAAGACCggg                                                     <  1:1255485/100‑1 (MQ=255)
                                              gTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGc                                               <  2:789723/100‑1 (MQ=255)
                                              gTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGc                                               <  2:509951/100‑1 (MQ=255)
                                               tGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCt                                              <  1:1716391/100‑1 (MQ=255)
                                               tGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCt                                              >  1:421850/1‑100 (MQ=255)
                                                ggAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCta                                             >  1:230621/1‑100 (MQ=255)
                                                ggAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCta                                             >  1:2184415/1‑100 (MQ=255)
                                                ggAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGACCTATATGTGGGCGACGACCGGGATTGGCta                                             >  2:151023/1‑100 (MQ=255)
                                                 gaagaaTCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCtat                                            >  1:2381918/1‑100 (MQ=255)
                                                 gaagaaTCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCtat                                            >  2:2001249/1‑100 (MQ=255)
                                                     aaTCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATaa                                          >  1:1374310/1‑98 (MQ=255)
                                                     aaTCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACg                                        >  1:1371763/1‑100 (MQ=255)
                                                     aaTCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACg                                        <  1:831491/100‑1 (MQ=255)
                                                       tctcGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACg                                        >  1:441160/1‑98 (MQ=255)
                                                       tctcGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGtt                                      <  1:1671043/100‑1 (MQ=255)
                                                       tctcGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGtt                                      <  2:1737823/100‑1 (MQ=255)
                                                        ctcGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGAGAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTg                                     <  2:921087/100‑1 (MQ=255)
                                                         tcGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGa                                    >  2:1209171/1‑100 (MQ=255)
                                                         tcGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGa                                    <  2:2307881/100‑1 (MQ=255)
                                                         tcGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGa                                    <  1:486768/100‑1 (MQ=255)
                                                         tcGACATGGAACTGCTGAAGATGGTCGACAAAAACGAGCGCGACAATAAATTTGGTATGCCCTATATGTGGGCGACGACCGGGATTGGATATAAGGtt                                      >  2:446303/1‑98 (MQ=255)
                                                                     tgctgAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGCCTATAACGTTGATAAAGTTAAAGc                        >  2:341651/1‑100 (MQ=255)
                                                                      gctgAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCg                       >  1:1809849/1‑100 (MQ=255)
                                                                       ctgAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCgg                      >  2:1959534/1‑100 (MQ=255)
                                                                         gAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGGTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTg                    <  2:341256/100‑1 (MQ=255)
                                                                         gAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTg                    >  1:317365/1‑100 (MQ=255)
                                                                            gCTGGTCGCCAAGAAGGAGCCCGCCAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTg                    <  2:1595599/97‑1 (MQ=255)
                                                                             cTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTgg                <  2:1563550/100‑1 (MQ=255)
                                                                               ggTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGc              <  1:387953/100‑1 (MQ=255)
                                                                               ggTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGc              <  2:515603/100‑1 (MQ=255)
                                                                                gTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCg             >  1:765645/1‑100 (MQ=255)
                                                                                 tCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGa            >  1:1788165/1‑100 (MQ=255)
                                                                                   gCCAAACACGATCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTTATAAAGTTAAAGCGGTGCTGGGCGaaa          >  2:1838119/1‑100 (MQ=255)
                                                                                   gCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGaaa          <  2:2381890/100‑1 (MQ=255)
                                                                                         cacGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGAAAACgcgc    <  2:553783/100‑1 (MQ=255)
                                                                                         cacGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGAAAACgcgc    <  2:240126/100‑1 (MQ=255)
                                                                                          acGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGAAAACGCGcc   <  1:949717/100‑1 (MQ=255)
                                                                                           cGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGAAAACGCGccc  >  2:2051250/1‑100 (MQ=255)
                                                                                              |                                                                                                
TGACTGCGGGAGTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGATCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGAAAACGCGCCC  >  NC_000913/894061‑894251

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