Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,805,532	+T	coding (718/1203 nt)	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,532	1	.	T	100.0%	162.2 / NA	42	N240Y (AAT→TAT)	proV	glycine betaine/proline ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base T (20/22);  total (20/22)

AAGCCTTCTCGGCGCTCGATCCATTAATTCGCACCGAGATGCAGGATGAGCTGGTAAAATTACAGGCGAAACATCAGCGCACCATTGTCTTTATTTCCCACGATCTTGATGAAGCCATGCGTATTGGCGACCGAATTGCCATTATGCAAA‑ATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTTA  >  NC_000913/2805383‑2805675                                                                                                                                                       |                                                                                                                                                aaGCCTTCTCGACCCTCGATCCATTAATTAGCACCGAGATTCAGGATTAGCTGGTAAAATTACAGGCGAAACTGCATCGAACCATTGTCTTTATTTCCCACGATCTTGATGAAGCCATGCGGATTGGCGTCCGAATTGCCATTATGCaaat                                                                                                                                                 <  2:374513/151‑2 (MQ=255)              gcTCGATCCATTAATTCGCACCGAGATGCAGGATGAGCTGGTAAAATTACAGGCGAAACATCAGCGCACCATTGTCTTTATTTCCCACGATCTTGATGAAGCCATGCGTATTGGCGACCGAATTGCCATTATGCAAATATGGTGAAGTGGt                                                                                                                                    >  2:543603/1‑151 (MQ=255)                 cGATCCATTAATTCGCACCGAGATGCAGGATGAGCTGGTAAAATTACAGGCGAAACATCAGCGCACCATTGTCTTTATTTCCCACGATCTTGATGAAGCCATGCGTATTGGCGACCGAATTGCCATTATGCAAATATGGTGAAGTGGTACa                                                                                                                                 <  2:7370/151‑1 (MQ=255)                  gATCCATTAATTCGCACCGAGATGCAGGATGAGCTGGTAAAATTACAGGCGAAACATCAGCGCACCATTGTCTTTATTTCCCACGATCTTGATGAAGCCATGCGTATTGGCGACCGAATTGCCATTATGCAAATATGGTGAAGTGGTAc                                                                                                                                  >  2:433235/1‑149 (MQ=255)                    tCCATTAATTCGCACCGAGATGCAGGATGAGCTGGTAAAATTACAGGCGAAACATCAGCGCACCATTGTCTTTATTTCCCACGATCTTGATGAAGCCATGCGTATTGGCGACCGAATTGCCATTATGCAAATATGGTGAAGTGGTACAGGt                                                                                                                              >  1:520834/1‑151 (MQ=255)                        ttAATTCGCACCGAGATGCAGGATGAGCTGGTAAAATTACAGGCGAAACATCAGCGCACCATTGTCTTTATTTCCCACGATCTTGATGAAGCCATGCGTATTGGCGACCGAATTGCCATTATGCAAATATGGTGAAGTGGTACAGGTCGGc                                                                                                                          >  2:135479/1‑151 (MQ=255)                         tAATTCGCACCGAGATGCAGGATGAGCTGGTAAAATTACAGGCGAAACATCAGCGCACCATTGTCTTTATTTCCCACGATCTTGATGAAGCCATGCGTATTGGCGACCGAATTGCCATTATGCAAATATGGTGAAGTGGTACAGGTCGGc                                                                                                                          >  2:307237/1‑150 (MQ=255)                          aaTTCGCACCGAGATGCAGGATGAGCTGGTAAAATTACAGGCGAAACATCAGCGCACCATTGTCTTTATTTCCCACGATCTTGATGAAGCCATGCGTATTGGCGACCGAATTGCCATTATGCAAATATGGTGAAGTGGTACAGGTCGGcac                                                                                                                        <  1:670374/151‑1 (MQ=255)                          aaTTCGCACCGAGATGCAGGATGAGCTGGTAAAATTACAGGCGAAACATCAGCGCACCATTGTCTTTATTTCCCACGATCTTGATGAAGCCATGCGTATTGGCGACCGAATTGCCATTATGCAAATATGGTGAAGTGGTACAGGTCGGca                                                                                                                         <  2:90866/150‑1 (MQ=255)                                 aCCGAGATGCAGGATGAGCTGGTAAAATTACAGGCGAAACATCAGCGCACCATTGTCTTTATTTCCCACGATCTTGATGAAGCCATGCGTATTGGCGACCGAATTGCCATTATGCAAATATGGTGAAGTGGTACAGGTCGGCACACCGGAt                                                                                                                 >  1:152189/1‑151 (MQ=255)                                 aCCGAGATGCAGGATGAGCTGGTAAAATTACAGGCGAAACATCAGCGCACCATTGTCTTTATTTCCCACGATCTTGATGAAGCCATGCGTATTGGCGACCGAATTGCCATTATGCAAATATGGTGAAGTGGTACAGGTCGGCACACCGGAt                                                                                                                 <  1:466144/151‑1 (MQ=255)                                           aGGATGAGCTGGTAAAATTACAGGCGAAACATCAGCGCACCATTGTCTTTATTTCCCACGATCTTGATGAAGCCATGCGTATTGGCGACCGAATTGCCATTATGCAAATATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATtctc                                                                                                        >  2:372642/1‑150 (MQ=255)                                              aTGAGCTGGTAAAATTACAGGCGAAACATCAGCGCACCATTGTCTTTATTTCCCACGATCTTGATGAAGCCATGCGTATTGGCGACCGAATTGCCATTATGCAAATATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCaata                                                                                                    <  1:575761/151‑1 (MQ=255)                                               tGAGCTGGTAAAATTACAGGCGAAACATCAGCGCACCATTGTCTTTATTTCCCACGATCTTGATGAAGCCATGCGTATTGGCGACCGAATTGCCATTATGCAAATATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCaata                                                                                                    >  1:480319/1‑150 (MQ=255)                                                     ggTAAAATTACAGGCGAAACATCAGCGCACCATTGTCTTTATTTCCCACGATCTTGATGAAGCCATGCGTATTGGCGACCGAATTGCCATTATGCAAATATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGc                                                                                             <  2:321865/151‑1 (MQ=255)                                                         aaaTTACAGGCGAAACATCAGCGCACCATTGTCTTTATTTCCCACGATCTTGATGAAGCCATGCGTATTGGCGACCGAATTGCCATTATGCAAATATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAAt                                                                                         <  2:687422/151‑1 (MQ=255)                                                           aTTACAGGCGAAACATCAGCGCACCATTGTCTTTATTTCCCACGATCTTGATGAAGCCATGCGTATTGGCGACCGAATTGCCATTATGCAAATATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGa                                                                                       <  2:279277/151‑1 (MQ=255)                                                              aCAGGCGAAACATCAGCGCACCATTGTCTTTATTTCCCACGATCTTGATGAAGCCATGCGTATTGGCGACCGAATTGCCATTATGCAAATATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGAtt                                                                                     >  2:159088/1‑150 (MQ=255)                                                               cAGGCGAAACATCAGCGCACCATTGTCTTTATTTCCCACGATCTTGATGAAGCCATGCGTATTGGCGACCGAATTGCCATTATGCAAATATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTAt                                                                                   <  1:119388/151‑1 (MQ=255)                                                                    gAAACATCGGCGCACCATTGTCTTTATTTCCCACGATCTTGATGAAGCCATGCGTATTGGCGACCGAATTGCCATTATGCAAATATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCg                                                                              >  2:329373/1‑151 (MQ=255)                                                                      aaCATCAGCGCACCATTGTCTTTATTTCCCACGATCTTGATGAAGCCATGCGTATTGGCGACCGAATTGCCATTATGCAAATATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGt                                                                             >  2:226553/1‑150 (MQ=255)                                                                       aCATCAGCGCACCATTGTCTTTATTTCCCACGATCTTGATGAAGCCATGCGTATTGGCGACCGAATTGCCATTATGCAAATATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTAc                                                                           <  1:48330/151‑1 (MQ=255)                                                                             gcgcACCATTGTCTTTATTTCCCACGATCTTGATGAAGCCATGCGTATTGGCGACCGAATTGCCATTATGCAAATATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACcttctt                                                                     <  2:140007/151‑1 (MQ=255)                                                                                       gTCTTTATTTCCCACGATCTTGATGAAGCCATGCGTATTGGCGACCGAATTGCCATTATGCAAATATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGtt                                                           <  2:137476/151‑1 (MQ=255)                                                                                       gTCTTTATTTCCCACGATCTTGATGAAGCCATGCGTATTGGCGACCGAATTGCCATTATGCAAATATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGtt                                                           <  1:317300/151‑1 (MQ=255)                                                                                           ttATTTCCCACGATCTTGATGAAGCCATGCGTATTGGCGACCGAATTGCCATTATGCAAATATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGa                                                         <  2:119938/149‑1 (MQ=255)                                                                                               ttCCCACGATCTTGATGAAGCCATGCGTATTGGCGACCGAATTGCCATTATGCAAATATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATAtt                                                     >  2:623318/1‑149 (MQ=255)                                                                                                   cACGATCTTGATGAAGCCATGCGTATTGGCGACCGAATTGCCATTATGCAAATATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTc                                                 >  1:19147/1‑149 (MQ=255)                                                                                                    aCGATCTTGATGAAGCCATGCGTATTGGCGACCGAATTGCCATTATGCAAATATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAgg                                              <  1:520142/151‑1 (MQ=255)                                                                                                        tCTTGATGAAGCCATGCGTATTGGCGACCGAATTGCCATTATGCAAATATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTAtt                                          >  2:271476/1‑151 (MQ=255)                                                                                                          ttGATGAAGCCATGCGTATTGGCGACCGAATTGCCATTATGCAAATATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTc                                         >  1:536609/1‑150 (MQ=255)                                                                                                           tgatgaAGCCATGCGTATTGGCGACCGAATTGCCATTATGCAAATATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTAtt                                          <  1:307267/148‑1 (MQ=255)                                                                                                             atgaAGCCATGCGTATTGGCGACCGAATTGCCATTATGCAAATATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAGTg                                     <  2:92343/151‑1 (MQ=255)                                                                                                                     aTGCGTATTGGCGACCGAATTGCCATTATGCAAATATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGat                             <  1:5795/151‑1 (MQ=255)                                                                                                                            ttGGCGACCGAATTGCCATTATGCAAATATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGATATTGccc                      >  2:27355/1‑151 (MQ=255)                                                                                                                               gCGACCGAATTGCCATTATGCAAATATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGATATTGCccgcc                   <  1:396133/151‑1 (MQ=255)                                                                                                                                cGACCGAATTGCCATTATGCAAATATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGATATTGCccgccg                  >  2:500692/1‑151 (MQ=255)                                                                                                                                       aTTGCCATTATGCAAATATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTAAGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCg           <  1:372811/151‑1 (MQ=255)                                                                                                                                         tGCCATTATGCAAATATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGa          >  2:403336/1‑150 (MQ=255)                                                                                                                                          gCCATTATGCAAATATGGTGAAGTGGTACAGGTCGGCACGCCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGaa         <  2:665757/150‑1 (MQ=255)                                                                                                                                          gCCATTATGCAAATATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAAt        <  2:13661/151‑1 (MQ=255)                                                                                                                                           ccATTATGCAAATATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATg       >  1:472801/1‑151 (MQ=255)                                                                                                                                                 tGCAAATATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTTa  >  2:55685/1‑150 (MQ=255)                                                                                                                                                       |                                                                                                                                                AAGCCTTCTCGGCGCTCGATCCATTAATTCGCACCGAGATGCAGGATGAGCTGGTAAAATTACAGGCGAAACATCAGCGCACCATTGTCTTTATTTCCCACGATCTTGATGAAGCCATGCGTATTGGCGACCGAATTGCCATTATGCAAA‑ATGGTGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTTA  >  NC_000913/2805383‑2805675

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