Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	894,989	C→T	intergenic (+93/‑2)	potF → / → potG	putrescine ABC transporter periplasmic binding protein/putrescine ABC transporter ATP binding subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	894,989	0	C	T	100.0%	174.6 / NA	50	intergenic (+93/‑2)	potF/potG	putrescine ABC transporter periplasmic binding protein/putrescine ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (25/25);  total (25/25)

TGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCGGTGGATGATGTCAGCCTGACCATCTACAA  >  NC_000913/894854‑895121                                                                                                                                        |                                                                                                                                     tgtgCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGTCGTg                                                                                                                                  >  1:748415/1‑140 (MQ=255)   ggcgcCCCCGCGCGTGGACAAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCGGGGGGGGCGCCCCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGTCGTGaa                                                                                                                                <  2:118754/139‑1 (MQ=255)     cccccccGCGGGGGGACCAAAGTGAAGAGCGAAAAATAACCCGCGGTCGTAGATGCGGGGGGGGCGCCCCACCCCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGTCGTGAATg                                                                                                                              <  2:748415/136‑1 (MQ=255)         cccGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCCCCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGTCGTGAATGACGc                                                                                                                          <  2:37720/140‑1 (MQ=255)         cccGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGTCGTGAATGACGc                                                                                                                          <  1:221116/140‑1 (MQ=255)          ccGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGTCGTGAATGACGCt                                                                                                                         >  1:549295/1‑140 (MQ=255)           cgcgcgTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGTCGTGAATGACGCTa                                                                                                                        >  1:587314/1‑140 (MQ=255)             cgcgTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGTCGTGAATGACGCTATc                                                                                                                      >  2:436365/1‑140 (MQ=255)                      cAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGTCGTGAATGACGCTATCCCTCGCCCg                                                                                                             <  1:436365/140‑1 (MQ=255)                       aaaGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCCCCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGTCGTGAATGACGCTATCCCTCGCCCGc                                                                                                            <  1:53553/140‑1 (MQ=255)                            gAAGAGCGAAAAATAACCCGCAGTCGTAGAGGCGGGGGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGTCGTGAATGACGCTATACCTCGCCCGCAGGCg                                                                                                       <  2:549295/140‑1 (MQ=255)                              agagCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGTCGTGAATGACGCTATCCCTCGCCCGCAGGCGaa                                                                                                     >  2:218724/1‑140 (MQ=255)                                                 aGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGTCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGAcg                                                                                  >  2:677447/1‑140 (MQ=255)                                                   tCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGTCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGAcgcc                                                                                <  2:1158943/140‑1 (MQ=255)                                                    ggTAGATGCGGGGGGGGCGCCCCACCCCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGTCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGAcgccg                                                                               <  2:423259/139‑1 (MQ=255)                                                    cGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGTCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGAcgccg                                                                               <  1:1100181/140‑1 (MQ=255)                                                            ccGGAGGGGCGCACCACCCCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGTCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGa                                                                       <  1:548498/140‑1 (MQ=255)                                                             cGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGTCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGaa                                                                      >  1:555065/1‑140 (MQ=255)                                                             cGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGTCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGaa                                                                      >  2:1003033/1‑140 (MQ=255)                                                             cGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGTCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGaa                                                                      >  1:423259/1‑140 (MQ=255)                                                             cGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGTCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGaa                                                                      >  1:1171149/1‑140 (MQ=255)                                                             cGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATCCCGGAGAGCAGTCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGATGCCGCTATTAGaa                                                                      >  2:1204039/1‑140 (MQ=255)                                                                       caccacACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGTCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAAcc                                                            >  2:185999/1‑140 (MQ=255)                                                                             ccccGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGTCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCa                                                      <  1:677447/139‑1 (MQ=255)                                                                                  ccgGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGTCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATcc                                                 >  2:894289/1‑140 (MQ=255)                                                                                   cgGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGTCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCt                                                >  1:953815/1‑140 (MQ=255)                                                                                    ggCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGTCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTa                                               >  2:113818/1‑140 (MQ=255)                                                                                      cAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGTCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACg                                             <  1:239341/140‑1 (MQ=255)                                                                                       aaTTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGTCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGa                                            <  2:441733/140‑1 (MQ=255)                                                                                             cACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGTCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCa                                      <  2:135287/140‑1 (MQ=255)                                                                                                  ttATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGTCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATg                                 <  1:218724/140‑1 (MQ=255)                                                                                                   tATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGTCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGc                                <  2:206440/140‑1 (MQ=255)                                                                                                    aTGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGTCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACAcgcg                               <  2:1171149/140‑1 (MQ=255)                                                                                                    aTGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGTCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCg                               <  2:587314/140‑1 (MQ=255)                                                                                                         gcgcTTGCACACATTCAATGCCGGAGAGCAGTCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCGGTGga                          <  1:1003033/140‑1 (MQ=255)                                                                                                          cgcTTGCACACATTCAATGCCGGAGAGCAGTCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCGGTGgat                         >  2:800178/1‑140 (MQ=255)                                                                                                            cTTGCACACATTCAATGCCGGAGAGCAGTCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAAc                                    <  2:350961/127‑1 (MQ=255)                                                                                                            cTTGCACACATTCAATGCCGGAGAGCAGTCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAAc                                    >  1:350961/1‑127 (MQ=255)                                                                                                            cTTGCACACATTCAATGCCGGAGAGCAGTCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCgg                              <  1:504759/133‑1 (MQ=255)                                                                                                            cTTGCACACATTCAATGCCGGAGAGCAGTCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCgg                              >  2:504759/1‑133 (MQ=255)                                                                                                            cTTGCACACATTCAATGCCGGAGAGCAGTCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCGGTGgatga                       >  2:67349/1‑140 (MQ=255)                                                                                                            cTTGCACACATTCAATGCCGGAGAGCAGTCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCGGTGgatga                       >  1:165290/1‑140 (MQ=255)                                                                                                                       tCAATGCCGGAGAGCAGTCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCGGTGga                          <  1:400719/126‑1 (MQ=255)                                                                                                                       tCAATGCCGGAGAGCAGTCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCGGTGga                          >  2:400719/1‑126 (MQ=255)                                                                                                                               ggAGAGCAGTCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCGGTGGATGATGTCAGCCTGACCATCTAc    >  2:509467/1‑140 (MQ=255)                                                                                                                               ggAGAGCAGTCGTGAATGACGCTATCCCTCACCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCGGTGGATGATGTCAGCCTGACCATCTAc    >  1:392606/1‑140 (MQ=255)                                                                                                                                gagaGCAGTCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCGGTGGATGATGTCAGCCTGACCATCTACa   <  1:894289/140‑1 (MQ=255)                                                                                                                                gagaGCAGTCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCGGTGGATGATGTCAGCCTGACCATCTACa   <  1:509467/140‑1 (MQ=255)                                                                                                                                 agagCAGTCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGTGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCGGTGGATGATGTCAGCCTGACCATCTACaa  >  1:133862/1‑140 (MQ=255)                                                                                                                                 agagCAGTCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCGGTGGATGATGTCAGCCTGACCATCTACaa  >  1:1031119/1‑140 (MQ=255)                                                                                                                                 agagCAGTCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCGGTGGATGATGTCAGCCTGACCATCTACaa  >  1:751073/1‑140 (MQ=255)                                                                                                                                     cagtCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCGGTGGATGATGTCAGCCTGAcc          <  1:1006866/128‑1 (MQ=255)                                                                                                                                        |                                                                                                                                     TGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAGGCGAAAACCCGTAAGGCGCTGACGCCGCTATTAGAAATCCGCAACCTGACCAAATCCTACGATGGTCAACATGCGGTGGATGATGTCAGCCTGACCATCTACAA  >  NC_000913/894854‑895121

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