Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,785,332	A→G	N25D (AAC→GAC)	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,332	0	A	G	100.0%	127.7 / NA	37	N25D (AAC→GAC)	gpmM	2,3‑bisphosphoglycerate‑independent phosphoglycerate mutase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (19/18);  total (19/18)

ATCGATTCAGCTGTAGTAAAATTACGCAAATTTTGACTCTTGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCTATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTAATCGACGCTTCCGGTCTGGAAGTCGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGG  >  NC_000913/3785201‑3785459                                                                                                                                    |                                                                                                                                aTCGATTCAGCTGTAGTAAAATTACGCAAATTTTGACTCTTGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCTATGGCTATCGCGAAGAACAGCAGGATGAc                                                                                                                               >  1:217998/1‑134 (MQ=255)   cGATTCAGCTGTAGTAAAATTACGCAAATTTTGACTCTTGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCTATGGCTATCGCGAAGAACAGCAGGATGACGCCAtttt                                                                                                                       >  1:243805/1‑140 (MQ=255)      ttCAGCTGTAGTAAAATTACGCAAATTTTGACTCTTGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCTATGGCTATCGCGAAGAACAGCAGGATGACGCCATTTTTAg                                                                                                                    <  2:243805/140‑1 (MQ=255)               agtaAAATTACGCAAATTTTGACTCTTGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCTATGGCTATCGCGAAGAACAGCAGGATGACGCCATTTTTAGTGCTAAAAc                                                                                                           <  2:925386/140‑1 (MQ=255)                   aaaTTACGCAAATTTTGACTCTTGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCTATGGCTATCGCGAAGAACAGCAGGATGACGCCATTTTTAGTGCTAAAACCCCg                                                                                                       <  1:615618/140‑1 (MQ=255)                     aTTACGCAAATTTTGACTCTTGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCTATGGCTATCGCGAAGAACAGCAGGATGACGCCATTTTTAGTGCTAAAACCCCGGt                                                                                                     >  1:14494/1‑140 (MQ=255)                        aCGCAAATTTTGACTCTTGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCTATGGCTATCGCGAAGAACAGCAGGATGACGCCATTTTTAGTGCTAAAACCCCGGTAAt                                                                                                  >  1:542095/1‑140 (MQ=255)                        aCGCAAATTTTGACTCTTGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCTATGGCTATCGCGAAGAACAGCAGGATGACGCCATTTTTAGTGCTAAAACCCCGGTAAt                                                                                                  >  2:438430/1‑140 (MQ=255)                         cGCAAATTTTGACTCTTGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCTATGGCTATCGCGAAGAACAGCAGGATGACGCCATTTTTAGTGCTAAAACCCCGGTAATg                                                                                                 >  1:317661/1‑140 (MQ=255)                             aaTTTTGACTCTTGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCTATGGCTATCGCGAAGAACAGCAGGATGACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATg                                                                                             <  1:339048/140‑1 (MQ=255)                             aaTTTTGACTCTTGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCTATGGCTATCGCGAAGAACAGCAGGATGACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATg                                                                                             <  2:317661/140‑1 (MQ=255)                                        ttGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCTATGGCTATCGCGAAGAACAGCAGGATGACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCa                                                                                           <  1:407289/131‑1 (MQ=255)                                        ttGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCTATGGCTATCGCGAAGAACAGCAGGATGACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCa                                                                                           >  2:407289/1‑131 (MQ=255)                                              aTGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCTATGGCTATCGCGAAGAACAGCAGGATGACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGt                                                                            <  2:14494/140‑1 (MQ=255)                                                   gTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCTATGGCTATCGCGAAGAACACCAGGATGACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCa                                                                       <  1:861391/140‑1 (MQ=255)                                                     tGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCTATGGCTATCGCGAAGAACAGCAGGATGACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATa                                                                     >  2:885517/1‑140 (MQ=255)                                                                  gTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCTATGGCTATCGCGAAGAACAGCAGGATGACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTAATCGACGc                                                        >  1:145573/1‑140 (MQ=255)                                                                  gTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCTATGGCTATCGCGAAGAACAGCAGGATGACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTAATCGACGc                                                        >  2:97268/1‑140 (MQ=255)                                                                     tCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCTATGGCTATCGCGAAGAACAGCAGGATGACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTAATCGACGCTTc                                                     <  1:438430/140‑1 (MQ=255)                                                                     tCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCTATGGCTATCGCGAAGAACAGCAGGATGACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTAATCGACGCTTc                                                     >  2:553477/1‑140 (MQ=255)                                                                       tAAAAAACCTATGGTACTGGTGATTCTGGATGGCTATGGCTATCGCGAAGAACAGCAGGATGACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTAATCGACGCTTCCg                                                   <  2:886903/140‑1 (MQ=255)                                                                           aaaCCTATGGTACTGGTGATTCTGGATGGCTATGGCTATCGCGAAGAACAGCAGGATGACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTAATCGACGCTTCCGGTCt                                               <  1:74798/140‑1 (MQ=255)                                                                                tATGGTACTGGTGATTCTGGATGGCTATGGCTATCGCGAAGAACAGCAGGATGACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTAATCGa                                                           >  1:918611/1‑123 (MQ=255)                                                                                tATGGTACTGGTGATTCTGGATGGCTATGGCTATCGCGAAGAACAGCAGGATGACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTAATCGa                                                           <  2:918611/123‑1 (MQ=255)                                                                                        tGGTGATTCTGGATGGCTATGGCTATCGCGAAGAACAGCAGGATGACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTAATCGACGCTTCCGGTCTGGAAGTCGGTCTg                                  >  1:280650/1‑140 (MQ=255)                                                                                          gTGATTCTGGATGGCTATGGCTATCGCGAAGAACAGCAGGATGACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTAATCGACGCTTCCGGTCTGGAAGTCGGTCTGcc                                <  1:894976/140‑1 (MQ=255)                                                                                                        ctatggctatCGCGAAGAACAGCAGGATGACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTAATCGACGCTTCCGGTCTGGAAGTCGGTCTGCCTGAc                            <  1:824012/130‑1 (MQ=255)                                                                                                        ctatggctatCGCGAAGAACAGCAGGATGACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTAATCGACGCTTCCGGTCTGGAAGTCGGTCTGCCTGAc                            >  2:824012/1‑130 (MQ=255)                                                                                                         tatggctatCGCGAAGAACAGCAGGATGACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTAATCGACGCTTCCGGTCTGGAAGTCGGTCTGCCTGACCGTCAGATggg                 >  2:933902/1‑140 (MQ=255)                                                                                                          atggctatCGCGAAGAACAGCAGGATGACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTAATCGACGCTTCCGGTCTGGAAGTCGGTCTGCCTGACCGTCAGATGGGt                <  2:188151/140‑1 (MQ=255)                                                                                                          aaggctaaCGCGGAGAACAGCAGGATGGCGCCAATTTTTGTGGTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTAATCGACGCTTCCGGTCTGGAAGTCGGTCTGCCTGACCGTCAGATGGGt                <  1:636777/138‑1 (MQ=255)                                                                                                                 tCGCGCAGAACAGCAGGATGCCGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTAATCGACGCTTCCGGTCTGGAAGTCGGTCTGCCTGACCGTCAGATGGGTAACTCCg         >  2:85257/1‑140 (MQ=255)                                                                                                                     gaagaaCAGCAGGATGACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTAATCGACGCTTCCGGTCTGGAAGTCGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGt     >  2:512848/1‑140 (MQ=255)                                                                                                                      aagaaCAGCAGGATGACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTAATCGACGCTTCCGGTCTGGAAGTCGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTa    >  2:387358/1‑140 (MQ=255)                                                                                                                        gaaCAGCAGGATGACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTAATCGACGCTTCCGGTCTGGAAGTCGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAgg  <  1:553477/140‑1 (MQ=255)                                                                                                                        gaaCAGCAGGATGACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTAATCGACGCTTCCGGTCTGGAAGTCGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAgg  <  1:933902/140‑1 (MQ=255)                                                                                                                           cagcagGATGACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTAATCGACGCTTCCGGTCTGGAAGTCGGTCTGCCTGACCGTCAGATGGGTAACTc           <  2:840438/128‑1 (MQ=255)                                                                                                                           cagcagGATGACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTAATCGACGCTTCCGGTCTGGAAGTCGGTCTGCCTGACCGTCAGATGGGTAACTc           >  1:840438/1‑128 (MQ=255)                                                                                                                                    |                                                                                                                                ATCGATTCAGCTGTAGTAAAATTACGCAAATTTTGACTCTTGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCTATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACCCTAATCGACGCTTCCGGTCTGGAAGTCGGTCTGCCTGACCGTCAGATGGGTAACTCCGAAGTAGG  >  NC_000913/3785201‑3785459

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