Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,805,493	G→T	E227* (GAA→TAA)	proV →	glycine betaine/proline ABC transporter periplasmic binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,805,493	0	G	T	100.0%	179.8 / NA	45	E227* (GAA→TAA)	proV	glycine betaine/proline ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (26/19);  total (26/19)

CGTTAGCGATTAATCCGGATATATTATTAATGGACGAAGCCTTCTCGGCGCTCGATCCATTAATTCGCACCGAGATGCAGGATGAGCTGGTAAAATTACAGGCGAAACATCAGCGCACCATTGTCTTTATTTCCCACGATCTTGATGAAGCCATGCGTATTGGCGACCGAATTGCCATTATGCAAAATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAGT  >  NC_000913/2805347‑2805639                                                                                                                                                   |                                                                                                                                                   cgTTAGCGATTAATCCGGATATATTATTAATGGACGAAGCCTTCTCGGCGCTCGATCCATTAATTCGCACCGAGATGCAGGATGAGCTGGTAAAATTACAGGCGAAACATCAGCGCACCATTGTCTTTATTTCCCACGATCTTGATTAAGc                                                                                                                                                >  2:511043/1‑151 (MQ=255)    tAGCGATTAATCCGGATATATTATTAATGGACGAAGCCTTCTCGGCGCTCGATCCATTAATTCGCACCGAGATGCAGGATGAGCTGGTAAAATTACAGGCGAAACATCAGCGCACCATTGTCTTTATTTCCCACGATCTTGATTAAGCCAt                                                                                                                                             >  1:944887/1‑151 (MQ=255)     aGCGATTAATCCGGATATATTATTAATGGACGAAGCCTTCTCGGCGCTCGATCCATTAATTCGCACCGAGATGCAGGATGAGCTGGTAAAATTACAGGCGAAACATCAGCGCACCATTGTCTTTATTTCCCACGATCTTGATTAAGCCATg                                                                                                                                            >  1:949915/1‑151 (MQ=255)      gCGATTAATCCGGATATATTATTAATGGACGAAGCCTTCTCGGCGCTCGATCCATTAATTCGCACCGAGATGCAGGATGAGCTGGTAAAATTACAGGCGAAACATCAGCGCACCATTGTCTTTATTTCCCACGATCTTGATTAAGCCATGc                                                                                                                                           >  2:172039/1‑151 (MQ=255)       cGATTAATCCGGATATATTATTAATGGACGAAGCCTTCTCGGCGCTCGATCCATTAATTCGCACCGAGATGCAGGATGAGCTGGTAAAATTACAGGCGAAACATCAGCGCACCATTGTCTTTATTTCCCACGATCTTGATTAAGCCATGCg                                                                                                                                          >  2:111672/1‑151 (MQ=255)       cGATTAATCCGGATATATTATTAATGGACGAAGCCTTCTCGGCGCTCGATCCATTAATTCGCACCGAGATGCAGGATGAGCTGGTAAAATTACAGGCGAAACATCAGCGCACCATTGTCTTTATTTCCCACGATCTTGATTAAGCCATGCg                                                                                                                                          <  1:688913/151‑1 (MQ=255)               ccGGATATATTATTAATGGACGAAGCCTTCTCGGCGCTCGATCCATTAATTCGCACCGAGATGCAGGATGAGCTGGTAAAATTACAGGCGAAACATCAGCGCACCATTGTCTTTATTTCCCACGATCTTGATTAAGCCATGCGTATTGGc                                                                                                                                   <  1:861522/150‑1 (MQ=255)                  gATATATTATTAATGGACGAAGCCTTCTCGGCGCTCGATCCATTAATTCGCACCGAGATGCAGGATGAGCTGGTAAAATTACAGGCGAAACATCAGCGCACCATTGTCTTTATTTCCCACGATCTTGATTAAGCCATGCGTATTGGCGAcc                                                                                                                               <  1:434563/151‑1 (MQ=255)                  gATATATTATTAATGGACGAAGCCTTCTCGGCGCTCGATCCATTAATTCGCACCGAGATGCAGGATGAGCTGGTAAAATTACAGGCGAAACATCAGCGCACCATTGTCTTTATTTCCCACGATCTTGATTAAGCCATGCGTATTGGCGAcc                                                                                                                               <  1:674813/151‑1 (MQ=255)                   atatatTATTAATGGACGAAGCCTTCTCGGCGCTCGATCCATTAATTCGCACCGAGATGCAGGATGAGCTGGTAAAATTACAGGCGAAACATCAGCGCACCATTGTCTTTATTTCCCACGATCTTGATTAAGCCATGCGTATTGGCGACCg                                                                                                                              >  1:356297/1‑151 (MQ=255)                    tatatTATTAATGGACGAAGCCTTCTCGGCGCTCGATCCATTAATTCGCACCGAGATGCAGGATGAGCTGGTAAAATTACAGGCGAAACATCAGCGCACCATTGTCTTTATTTCCCACGATCTTGATTAAGCCATGCGTATTGGCGACCg                                                                                                                              >  1:25649/1‑150 (MQ=255)                     atatTATTAATGGACGAAGCCTTCTCGGCGCTCGATCCATTAATTCGCACCGAGATGCAGGATGAGCTGGTAAAATTACAGGCGAAACATCAGCGCACCATTGTCTTTATTTCCCACGATCTTGATTAAGCCATGCGTATTGGCGACCGa                                                                                                                             >  1:170960/1‑150 (MQ=255)                      tattatTAATGGACGAAGCCTTCTCGGCGCTCGATCCATTAATTCGCACCGAGATGCAGGATGAGCTGGTAAAATTACAGGCGAAACATCAGCGCACCATTGTCTTTATTTCCCACGATCTTGATTAAGCCATGCGTATTGGCGACCGAAt                                                                                                                           >  1:644527/1‑151 (MQ=255)                        ttattaATGGACGAAGCCTTCTCGGCGCTCGATCCATTAATTCGCACCGAGATGCAGGATGAGCTGGTAAAATTACAGGCGAAACATCAGCGCACCATTGTCTTTATTTCCCACGATCTTGATTAAGCCATGCGTATTGGCGACCGAAtt                                                                                                                          <  1:771193/150‑1 (MQ=255)                         tattaATGGACGAAGCCTTCTCGGCGCTCGATCCATTAATTCGCACCGAGATGCAGGATGAGCTGGTAAAATTACAGGCGAAACATCAGCGCACCATTGTCTTTATTTCCCACGATCTTGATTAAGCCATGCGTATTGGCGACCGAATTg                                                                                                                         <  1:645732/150‑1 (MQ=255)                            taATGGACGAAGCCTTCTCGGCGCTCGATCCATTAATTCGCACCGAGATGCAGGATGAGCTGGTAAAATTACAGGCGAAACATCAGCGCACCATTGTCTTTATTTCCCACGATCTTGATTAAGCCATGCGTATTGGCGACCGAATTGCCAt                                                                                                                     >  1:551762/1‑151 (MQ=255)                              aTGGACGAAGCCTTCTCGGCGCTCGATCCATTAATTCGCACCGAGATGCAGGATGAGCTGGTAAAATTACAGGCGAAACATCAGCGCACCATTGTCTTTATTTCCCACGATCTTGATTAAGCCATGCGTATTGGCGACCGAATTGCCAtt                                                                                                                    >  2:722054/1‑150 (MQ=255)                                  aCGAAGCCTTCTCGGCGCTCGATCCATTAATTCGCACCGAGATGCAGGATGAGCTGGTAAAATTACAGGCGAAACATCAGCGCACCATTGTCTTTATTTCCCACGATCTTGATTAAGCCATGCGTATTGGCGACCGAATTGCCATTATGCa                                                                                                               <  1:393676/151‑1 (MQ=255)                                       gCCTTCTCGGCGCTCGATCCATTAATTCGCACCGAGATGCATGATGAGCTGGTAAAATTACAGGCGAAACATCAGCGCACCATTGTCTTTATTTCCCACGATCTTGATTAAGCCATGCGTATTGGCGACCGAATTGCCATTATGCAAAAt                                                                                                           <  1:90011/150‑1 (MQ=255)                                                     cGATCCATTAATTCGCACCGAGATGCAGGATGAGCTGGTAAAATTACAGGCGAAACATCAGCGCACCATTGTCTTTATTTCCCACGATCTTGATTAAGCCATGCGTATTGGCGACCGAATTGCCATTATGCAAAATGGTGAAGTGGTACAg                                                                                            <  1:456475/151‑1 (MQ=255)                                                      gATCCATTAATTCGCACCGAGATGCAGGATGAGCTGGTAAAATTACAGGCGAAACATCAGCGCACCATTGTCTTTATTTCCCACGATCTTGATTAAGCCATGCGTATTGGCGACCGAATTGCCATTATGCAAAATGGTGAAGTGGTACAgg                                                                                           >  2:674955/1‑151 (MQ=255)                                                       aTCCATTAATTCGCACCGAGATGCAGGATGAGCTGGTAAAATTACAGGCGAAACATCAGCGCACCATTGTCTTTATTTCCCACGATCTTGATTAAGCCATGCGTATTGGCGACCGAATTGCCATTATGCAAAATGGTGAAGTGGTACAGGt                                                                                          >  2:158169/1‑151 (MQ=255)                                                         ccATTAATTCGCACCGAGATGCAGGATGAGCTGGTAAAATTACAGGCGAAACATCAGCACACCATTGTCTTTATTTCCCACGATCTTGATTAAGCCATGCGTATTGGCGACCGAATTGCCATTATGCAAAATGGTGAAGTGGTACAGGTCg                                                                                        >  1:610524/1‑151 (MQ=255)                                                               aTTCGCACCGAGATGCAGGATGAGCTGGTAAAATTACAGGCGAAACATCAGCGCACCATTGTCTTTATTTCCCACGATCTTGATTAAGCCATGCGTATTGGCGACCGAATTGCCATTATGCAAAATGGTGAAGTGGTACAGGTCGGCacac                                                                                  <  2:841139/151‑1 (MQ=255)                                                                    cACCGAGATGCAGGATGAGCTGGTAAAATTACAGGCGAAACATCAGCGCACCATTGTCTTTATTTCCCACGATCTTGATTAAGCCATGCGTATTGGCGACCGAATTGCCATTATGCAAAATGGTGAAGTGGTACAGGTCGGCACACCGGAt                                                                             <  1:336483/151‑1 (MQ=255)                                                                      ccGAGATGCAGGATGAGCTGGTAAAATTACAGGCGAAACATCAGCGCACCATTGTCTTTATTTCCCACGATCTTGATTAAGCCATGCGTATTGGCGACCGAATTGCCATTATGCAAAATGGTGAAGTGGTACAGGTCGGCACACCGGATg                                                                            >  1:189737/1‑150 (MQ=255)                                                                          gaTGCAGGATGAGCTGGTAAAATTACAGGCGAAACATCAGCGCACCATTGTCTTTATTTCCCACGATCTTGATTAAGCCATGCGTATTGGCGACCGAATTGCCATTATGCAAAATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAAtt                                                                       <  1:876010/151‑1 (MQ=255)                                                                           aTGCAGGATGAGCTGGTAAAATTACAGGCGAAACATCAGCGCACCATTGTCTTTATTTCCCACGATCTTGATTAAGCCATGCGTATTGGCGATCGAATTGCCATTATGCAAAATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATtc                                                                      >  1:151854/1‑151 (MQ=255)                                                                                       cTGGTAAAATTACAGGCGAAACATCAGCGCACCATTGTCTTTATTTCCCACGATCTTGATTAAGCCATGCGTATTGGCGACCGAATTGCCATTATGCAAAATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCgg                                                          >  2:584033/1‑151 (MQ=255)                                                                                         ggTAAAATTACAGGCGAAACATCAGCGCACCATTGTCTTTATTTCCCACGATCTTGATTAAGCCATGCGTATTGGCGACCGAATTGCCATTATGCAAAATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCg                                                        >  1:572689/1‑151 (MQ=255)                                                                                           tAAAATTACAGGCGAAACATCAGCGCACCATTGTCTTTATTTCCCACGATCTTGATTAAGCCATGCGTATTGGCGACCGAATTGCCATTATGCAAAATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGa                                                       >  2:635301/1‑150 (MQ=255)                                                                                            aaaaTTACAGGCGAAACATCAGCGCACCATTGTCTTTATTTCCCACGATCTTGATTAAGCCATGCGTATTGGCGACCGAATTGCCATTATGCAAAATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGaa                                                      <  2:87066/150‑1 (MQ=255)                                                                                                 tACAGGCGAAACATCAGCGCACCATTGTCTTTATTTCCCACGATCTTGATTAAGCCATGCGTATTGGCGACCGAATTGCCATTATGCAAAATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGAtt                                                 >  2:807063/1‑150 (MQ=255)                                                                                                  aCAGGCGAAACATCAGCGCACCATTGTCTTTATTTCCCACGATCTTGATTAAGCCATGCGTATTGGCGACCGAATTGCCATTATGCAAAATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTAt                                               >  2:234397/1‑151 (MQ=255)                                                                                                  aCAGGCGAAACATCAGCGCACCATTGTCTTTATTTCCCACGATCTTGATTAAGCCATGCGTATTGGCGACCGAATTGCCATTATGCAAAATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTAt                                               <  2:817518/151‑1 (MQ=255)                                                                                                           aCATCAGCGCACCATTGTCTTTATTTCCCACGATCTTGATTAAGCCATGCGTATTGGCGACCGAATTGCCATTATGCAAAATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTAcc                                      >  1:573940/1‑151 (MQ=255)                                                                                                               cAGCGCACCATTGTCTTTATTTCCCACGATCTTGATTAAGCCATGCGTATTGGCGACCGAATTGCCATTATGCAAAATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACcttct                                  >  2:401462/1‑151 (MQ=255)                                                                                                               cAGCGCACCATTGTCTTTATTTCCCACGATCTTGATTAAGCCATGCGTATTGGCGACCGAATTGCCATTATGCAAAATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACcttct                                  <  1:399709/151‑1 (MQ=255)                                                                                                                   gcACCATTGTCTTTATTTCCCACGATCTTGATTAAGCCATGCGTATTGGCGACCGAATTGCCATTATGCAAAATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCAGGCGAATGATTATGTCCGTACCTTCTTCCg                              >  2:453701/1‑151 (MQ=255)                                                                                                                        aTTGTCTTTATTTCCCACGATCTTGATTAAGCCATGCGTATTGGCGACCGAATTGCCATTATGCAAAATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCg                         <  2:584302/151‑1 (MQ=255)                                                                                                                        aTTGTCTTTATTTCCCACGATCTTGATTAAGCCATGCGTATTGGCGACCGAATTGCCATTATGCAAAATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCg                         >  2:813822/1‑151 (MQ=255)                                                                                                                           gTCTTTATTTCCCACGATCTTGATTAAGCCATGCGTATTGGCGACCGAATTGCCATTATGCAAAATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTg                      <  1:436492/151‑1 (MQ=255)                                                                                                                           gTCTTTATTTCCCACGATCTTGATTAAGCCATGCGTATTGGCGACCGAATTGCCATTATGCAAAATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTg                      <  1:767260/151‑1 (MQ=255)                                                                                                                           gTCTTTATTGCCCACGATCTTGATTAAGCCATGCGTATTGGCGACCGAATTGCCATTATGCAAAATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGtt                       <  2:758128/150‑1 (MQ=255)                                                                                                                                       cACGATCTTGATTAAGCCATGCGTATTGGCGACCGAATTGCCATTATGCAAAATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAgg          >  1:924692/1‑151 (MQ=255)                                                                                                                                               tgattaAGCCATGCGTATTGGCGACCGAATTGCCATTATGCAAAATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAGt  >  1:954700/1‑151 (MQ=255)                                                                                                                                                   |                                                                                                                                                   CGTTAGCGATTAATCCGGATATATTATTAATGGACGAAGCCTTCTCGGCGCTCGATCCATTAATTCGCACCGAGATGCAGGATGAGCTGGTAAAATTACAGGCGAAACATCAGCGCACCATTGTCTTTATTTCCCACGATCTTGATGAAGCCATGCGTATTGGCGACCGAATTGCCATTATGCAAAATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAGT  >  NC_000913/2805347‑2805639

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑