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_000913	894,155	0	T	G	100.0%	207.8 / NA	58	D124E (GAT→GAG)	potF	putrescine ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (27/31);  total (27/31)

CTGCGGGAGTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGATCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGAAAACGCGCC  >  NC_000913/894064‑894250                                                                                            |                                                                                                cTGCGGGAGTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACaa                                                                                         >  2:19015/1‑100 (MQ=255)    cGGGAGTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGACAAACACGAGCCCGACAATaa                                                                                      >  2:343059/1‑100 (MQ=255)          ttttCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGc                                                                                <  2:697085/100‑1 (MQ=255)            ttCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACACTAACTTTGCTa                                                                              >  2:1546124/1‑100 (MQ=255)            ttCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTa                                                                              <  2:597221/100‑1 (MQ=255)            ttCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTa                                                                              >  1:537983/1‑100 (MQ=255)             tCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTAt                                                                             >  2:1533742/1‑100 (MQ=255)              ccAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATg                                                                            >  2:694642/1‑100 (MQ=255)               cAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGc                                                                           >  2:281597/1‑100 (MQ=255)                agcgtcTGGATTAAAGCAAATTGCCGGAGCGGAAGAATCTCGCCCCGAAACTGCTGAAGCTTGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGcc                                                                          <  2:1274817/95‑1 (MQ=255)                aGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGcc                                                                          >  2:1400221/1‑100 (MQ=255)                aGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGcc                                                                          <  2:1515708/100‑1 (MQ=255)                 gCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGccc                                                                         <  2:760774/100‑1 (MQ=255)                 gCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGccc                                                                         <  1:427090/100‑1 (MQ=255)                  ccGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCt                                                                        >  1:1293266/1‑100 (MQ=255)                     cTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCtata                                                                     <  2:1550303/100‑1 (MQ=255)                     cTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCtata                                                                     <  2:532141/100‑1 (MQ=255)                       ggACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTatat                                                                    <  2:1496250/99‑1 (MQ=255)                        taCAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATAtgt                                                                  >  2:230631/2‑100 (MQ=255)                         aCAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATAtgtg                                                                 <  2:1049070/100‑1 (MQ=255)                            aaaGCAAATTGCCGGAGTGGAAGAATCTCGACCGGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGc                                                              >  2:1028638/1‑100 (MQ=255)                            aaaGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGc                                                              >  1:1519533/1‑100 (MQ=255)                            aaaGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGc                                                              >  1:1458897/1‑100 (MQ=255)                             aaGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGcg                                                             >  1:215953/1‑100 (MQ=255)                                   aTTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGAcc                                                       <  2:567389/100‑1 (MQ=255)                                     tGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGGGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCgg                                                     <  1:1328893/100‑1 (MQ=255)                                      gCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCggg                                                    <  2:1003770/100‑1 (MQ=255)                                        cGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGAt                                                  >  1:42322/1‑100 (MQ=255)                                          gAGTGGAAGAATCTCGACCCGGAACTGCTCAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTTCTATGCCCTATATGTGGGCGACGACCGGGATTg                                                <  2:951800/100‑1 (MQ=255)                                             tGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCt                                             <  2:1409572/100‑1 (MQ=255)                                             tGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCt                                             <  2:1419546/100‑1 (MQ=255)                                              ggAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCta                                            >  1:424714/1‑100 (MQ=255)                                               gaagaaTCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCtat                                           >  2:1580011/1‑100 (MQ=255)                                                aagaaTCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCtata                                          <  1:1182920/100‑1 (MQ=255)                                                 agaaTCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATaa                                         >  2:403981/1‑100 (MQ=255)                                                      ctcGACCCGGAACTGCCGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTg                                    <  2:924450/100‑1 (MQ=255)                                                        cGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGAt                                  <  2:537982/100‑1 (MQ=255)                                                        cGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGAt                                  <  1:1354113/100‑1 (MQ=255)                                                         gACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATa                                 >  2:899852/1‑100 (MQ=255)                                                         gACCCGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATa                                 >  1:1299007/1‑100 (MQ=255)                                                            ccGGAACTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAg                              <  2:424713/100‑1 (MQ=255)                                                                aaCTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTaa                          <  2:1240107/100‑1 (MQ=255)                                                                aaCTGCTGAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTaa                          >  2:880738/1‑100 (MQ=255)                                                                     ctgAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCgg                     >  2:1114569/1‑100 (MQ=255)                                                                      tgAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGt                    >  1:1519112/1‑100 (MQ=255)                                                                       gAAGCTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTg                   >  2:887901/1‑100 (MQ=255)                                                                           cTGGTCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGGGCTgg               >  2:1141736/1‑100 (MQ=255)                                                                               tCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGa           <  1:1446692/100‑1 (MQ=255)                                                                               tCGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGa           <  1:214931/100‑1 (MQ=255)                                                                                cGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGaa          >  1:654580/1‑100 (MQ=255)                                                                                cGCCAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGaa          <  1:777457/100‑1 (MQ=255)                                                                                  ccAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGaaaa        <  1:281597/100‑1 (MQ=255)                                                                                   cAAACACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGAAAAc       <  2:1458889/100‑1 (MQ=255)                                                                                     aaCACGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGAAAAcgc     <  1:1255171/100‑1 (MQ=255)                                                                                      acacGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGAAAAcgcg    >  1:377419/1‑100 (MQ=255)                                                                                        acGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGAAAACGCGcc  <  1:1580021/100‑1 (MQ=255)                                                                                        acGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGAAAACGCGcc  <  1:1337810/100‑1 (MQ=255)                                                                                        acGAGCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGAAAACGCGcc  <  1:694644/100‑1 (MQ=255)                                                                                            |                                                                                                CTGCGGGAGTTTTCCAGCCGCTGGACAAAAGCAAATTGCCGGAGTGGAAGAATCTCGACCCGGAACTGCTGAAGCTGGTCGCCAAACACGATCCCGACAATAAATTTGCTATGCCCTATATGTGGGCGACGACCGGGATTGGCTATAACGTTGATAAAGTTAAAGCGGTGCTGGGCGAAAACGCGCC  >  NC_000913/894064‑894250

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: ‑