Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,785,421	C→A	V54V (GTC→GTA)	gpmM →	2,3‑bisphosphoglycerate‑independent phosphoglycerate mutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,785,421	0	C	A	100.0%	110.3 / NA	33	V54V (GTC→GTA)	gpmM	2,3‑bisphosphoglycerate‑independent phosphoglycerate mutase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (18/15);  total (18/15)

CTGGTGATTCTGGATGGCTATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTAATCGACGCTTCCGGTCTGGAAGTCGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGC  >  NC_000913/3785287‑3785543                                                                                                                                       |                                                                                                                           cTGGTGATTCTGGATGGCTATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTAATCGACGCTTCCGGTCTGGAAGTAGGTCt                                                                                                                       >  1:255364/1‑140 (MQ=255)   ggTGATTCTGGATGGCTATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTAATCGACGCTTCCGGTCTGGAAGTAGGTCTGc                                                                                                                     >  1:435840/1‑140 (MQ=255)    gTGATTCTGGATGGCTATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTAATCGACGCTTCCGGTCTGGAAGTAGGTCTGcc                                                                                                                    >  1:127121/1‑140 (MQ=255)         tCTGGATGGCTATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTAATCGACGCTTCCGGTCTGGAAGTAGGTCTGCCTGAcc                                                                                                               >  1:226488/1‑140 (MQ=255)           tGGATGGCTATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTAATCGACGCTTCCGGTCTGGAAGTAGGTCTGCCTGACCGt                                                                                                             <  1:31739/140‑1 (MQ=255)                    atggctatCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTAATCGACGCTTCCGGTCTGGAAGTAGGTCTGCCTGACCGTCAGATGGGt                                                                                                    <  1:344129/140‑1 (MQ=255)                       gctatCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGCCCGCATACCCTAATCGACGCTTCCGGTCTGGAAGTAGGTCTGCCTGACCGTCAGATGGGTAAc                                                                                                 >  1:153760/1‑140 (MQ=255)                           tCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTAATCGACGCTTCCGGTCTGGAAGTAGGTCTGCCTGACCGTCAGATGGGTAACTCCg                                                                                             >  1:449836/1‑140 (MQ=255)                            cgcgAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTAATCGACGCTTCCGGTCTGGAAGTAGGTCTGCCTGACCGTCAGATGGGTAACTCCGa                                                                                            >  1:432187/1‑140 (MQ=255)                            cgcgAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTAATCGACGCTTCCGGTCTGGAAGTAGGTCTGCCTGACCGTCAGATGGGTAACTCCGa                                                                                            >  2:2790/1‑140 (MQ=255)                                   aaCAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTAATCGACGCTTCCGGTCTGGAAGTAGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGc                                                                                     <  2:86401/140‑1 (MQ=255)                                                cGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTAATCGACGCTTCCGGTCTGGAAGTAGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAgg                                                                                      >  1:162106/1‑126 (MQ=255)                                                cGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTAATCGACGCTTCCGGTCTGGAAGTAGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAgg                                                                                      <  2:162106/126‑1 (MQ=255)                                                          aGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTAATCGACGCTTCCGGTCTGGAAGTAGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTgccggccg                                                              <  2:127121/140‑1 (MQ=255)                                                           gTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTAATCGACGCTTCCGGTCTGGAAGTAGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGc                                                             <  1:367870/140‑1 (MQ=255)                                                               tAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTAATCGACGCTTCCGGTCTGGAAGTAGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCg                                                         >  2:443251/1‑140 (MQ=255)                                                                     cccGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTAATCGACGCTTCCGGTCTGGAAGTAGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATc                                                   <  1:377477/140‑1 (MQ=255)                                                                         gTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTAATCGACGCTTCCGGTCTGGAAGTAGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGa                                               >  1:276780/1‑140 (MQ=255)                                                                               gATGCACTGTGGGCCAATCGTCCGCATACCCTAATCGACGCTTCCGGTCTGGAAGTAGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGAc                                         <  2:286737/140‑1 (MQ=255)                                                                                    aCTGTGGGCCAATCGTCCGCATACCCTAATCGACGCTTCCGGTCTGGAAGTAGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTc                                    >  1:372999/1‑140 (MQ=255)                                                                                       gtgGGCCAATCGTCCGCATACCCTAATCGACGCTTCCGGTCTGGAAGTAGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTgg                                 <  1:443251/140‑1 (MQ=255)                                                                                        tgGGCCAATCGTCCGCATACCCTAATCGACGCTTCCGGTCTGGAAGTAGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGa                                >  2:416891/1‑140 (MQ=255)                                                                                        tgGGCCAATCGTCCGCATACCCTAATCGACGCTTCCGGTCTGGAAGTAGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGa                                >  1:454603/1‑140 (MQ=255)                                                                                          ggCCAATCGTCCGCATACCCTAATCGACGCTTCCGGTCTGGAAGTAGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACg                              >  2:292/1‑140 (MQ=255)                                                                                             cAATCGTCCGCATACCCTAATCGACGCTTCCGGTCTGGAAGTAGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACGTTg                           <  2:449836/140‑1 (MQ=255)                                                                                                  gTCCGCATACCCTAATCGACGCTTCCGGTCTGGAAGTAGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATc                      >  1:287755/1‑140 (MQ=255)                                                                                                            ccTAATCGACGCTTCCGGTCTGGAAGTAGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATcgcg            <  2:153105/140‑1 (MQ=255)                                                                                                             cTAATCGACGCTTCCGGTCTGGAAGTAGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCgcgc           >  1:465555/1‑140 (MQ=255)                                                                                                                aTCGACGCTTCCGGTCTGGAAGTAGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGcttt        <  2:235851/140‑1 (MQ=255)                                                                                                                 tCGACGCTTCCGGTCTGGAAGTAGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGctttc       <  2:372999/140‑1 (MQ=255)                                                                                                                   gACGCTTCCGGTCTGGAAGTAGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGctttctt     <  2:465555/140‑1 (MQ=255)                                                                                                                      gCTTCCGGTCTGGAAGTAGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGc  <  1:44588/140‑1 (MQ=255)                                                                                                                      gCTTCCGGTCTGGAAGTAGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGc  >  2:8815/1‑140 (MQ=255)                                                                                                                                       |                                                                                                                           CTGGTGATTCTGGATGGCTATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTAATCGACGCTTCCGGTCTGGAAGTCGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGGCCACGTTAACCTGGGTGCCGGCCGCATCGTGTATCAGGACCTGACTCGTCTGGACGTTGAAATCAAAGATCGCGCTTTCTTTGC  >  NC_000913/3785287‑3785543

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