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% 147.2 / ‑3.8 43D124E (GAT→GAGpotFputrescine ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (22/20);  total (23/20)
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.53e-01

TGACTGCGGGAGTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGATCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGAAAACGCG  >  NC_000913/894061‑894248
                                                                                              |                                                                                             
tGACTGCGGGAGTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGa                                                                                          >  1:1511996/1‑100 (MQ=255)
 gACTGCGGGAGTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGAc                                                                                         <  1:6864/100‑1 (MQ=255)
      cGGGAGTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATaa                                                                                    >  2:1179244/1‑100 (MQ=255)
        ggAGTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACAAGCCCGACAATAAAt                                                                                  >  1:809045/1‑100 (MQ=255)
          aGTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAAttt                                                                                >  2:1252274/1‑100 (MQ=255)
            ttttCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGc                                                                              <  1:1252282/100‑1 (MQ=255)
            ttttCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCCCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGCGCCCGACAATAAATTTGc                                                                              >  2:39224/1‑100 (MQ=255)
             tttCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCt                                                                             >  1:212166/1‑100 (MQ=255)
             tttCCAGCCGCTGGAAAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCt                                                                             <  1:1545961/100‑1 (MQ=255)
                  aGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGcc                                                                        <  1:729401/100‑1 (MQ=255)
                    ccGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCt                                                                      >  1:291062/1‑100 (MQ=255)
                        tGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTatat                                                                  >  2:547770/1‑100 (MQ=255)
                         ggACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTAAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATAtg                                                                 <  2:185793/100‑1 (MQ=255)
                           aCAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATAtgtg                                                               <  2:1539420/100‑1 (MQ=255)
                           aCAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATAtgtg                                                               >  1:1308867/1‑100 (MQ=255)
                            cAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTgg                                                              <  2:400439/100‑1 (MQ=255)
                            cAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTgg                                                              >  1:679674/1‑100 (MQ=255)
                             aaaaGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTggg                                                             <  2:98443/100‑1 (MQ=255)
                               aaGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGcg                                                           <  1:611378/100‑1 (MQ=255)
                                       tGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCgg                                                   <  1:487296/100‑1 (MQ=255)
                                       tGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCgg                                                   <  1:1403055/100‑1 (MQ=255)
                                        gCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCggg                                                  <  2:212165/100‑1 (MQ=255)
                                         ccGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGa                                                 >  1:801942/1‑100 (MQ=255)
                                               tGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCt                                           >  1:1033812/1‑100 (MQ=255)
                                                ggAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCta                                          >  1:355283/1‑100 (MQ=255)
                                                 gaagaaTCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCtat                                         >  1:1494365/1‑100 (MQ=255)
                                                      aTCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGt                                    <  1:1410612/100‑1 (MQ=255)
                                                        ctcGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTg                                  <  2:526384/100‑1 (MQ=255)
                                                             cccGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATaa                              <  2:603228/99‑1 (MQ=255)
                                                                    ctgctgAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAg                      <  1:1408191/100‑1 (MQ=255)
                                                                      gctgAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCg                    >  1:81341/1‑100 (MQ=255)
                                                                      gctgAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCg                    <  2:1511986/100‑1 (MQ=255)
                                                                       ctgAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCgg                   >  1:1271928/1‑100 (MQ=255)
                                                                         gAAGCTGGTCGCCAAACACGATCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTg                 >  2:1149789/1‑100 (MQ=255)
                                                                            gCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTg              <  2:1494355/100‑1 (MQ=255)
                                                                              tGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTggg            <  1:870165/100‑1 (MQ=255)
                                                                              tGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTggg            <  1:334339/100‑1 (MQ=255)
                                                                                   gCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGaaa       >  2:1447428/1‑100 (MQ=255)
                                                                                    ccAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGaaaa      >  1:114019/1‑100 (MQ=255)
                                                                                    ccAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATGGGCTATAACGTTGCTCAAGTTAAAGCGGTGCTGGGCg          >  2:215626/1‑96 (MQ=255)
                                                                                     cAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGAAAAc     >  1:1053712/1‑100 (MQ=255)
                                                                                     cAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGAAAAc     >  1:162060/1‑100 (MQ=255)
                                                                                        acacGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGAAAAcgcg  >  2:925878/1‑100 (MQ=255)
                                                                                              |                                                                                             
TGACTGCGGGAGTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGATCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGAAAACGCG  >  NC_000913/894061‑894248

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