Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,123,008	G→A	intergenic (‑29/+76)	cpsG ← / ← cpsB	phosphomannomutase/mannose‑1‑phosphate guanylyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,123,008	0	G	A	100.0%	157.8 / NA	46	intergenic (‑29/+76)	cpsG/cpsB	phosphomannomutase/mannose‑1‑phosphate guanylyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (22/24);  total (22/24)

TGAGAAATTCGCCATAGGCGCGACCAATGCGCCAGGCGATATCTTCATTCAGTTCTTCGCCTAATTTCCCGCGAATATCATAGGCTTTAAAGCAGGTTAATTTTTTCATATCGTTACCCTTTTTCAGGCAATAGTTGGCCCTGACCGAAGCGGCCAATTTTATTTTTGTCGTTATTTACACCGCGGTTTCGCATTCATTGCCTGATGCGACGTTTACACCCGTCCGTAGCGATCCGCGAAACGCACCACATCATCCTCTTCGAGATAAGAG  >  NC_000913/2122871‑2123141                                                                                                                                          |                                                                                                                                      tGAGAAATTCGCCATAGGCGCGACCAATGCGCCAGGCGATATCTTCATTCAGTTCTTCGCCTAATTTCCCGCGAATATCATAGGCTTTAAAGCAGGTTAATTTTTTCATATCGTTACCCTTTTTCAGGCAATAGTTGAcc                                                                                                                                     <  2:229488/140‑1 (MQ=255)                       aCCAATGCGCCAGGCGATATCTTCATTCAGTTCTTCGCCTAATTTCCCGCGAATATCATAGGCTTTAAAGCAGGTTAATTTTTTCATATCGTTACCCTTTTTCAGGCAATAGTTGACCCTGACCGAAGCGGCCAatttta                                                                                                               >  1:371353/1‑140 (MQ=255)                       aCCAATGCGCCAGGCGATATCTTCATTCAGTTCTTCGCCTAATTTCCCGCGAATATCATAGGCTTTAAAGCAGGTTAATTTTTTCATATCGTTACCCTTTTTCAGGCAATAGTTGACCCTGACCGAAGCGGCCAatttta                                                                                                               >  1:298586/1‑140 (MQ=255)                          aaTGCGCCAGGCGATATCTTCATTCAGTTCTTCGCCTAATTTCCCGCGAATATCATAGGCTTTAAAGCAGGTTAATTTTTTCATATCGTTACCCTTTTTCAGGCAATAGTTGACCCTGACCGAAGCGGCCAattttattt                                                                                                            <  1:707383/140‑1 (MQ=255)                            tGCGCCAGGCGATATCTTCATTCAGTTCTTCGCCTAATTTCCCGCGAATATCATAGGCTTTAAAGCAGGTTAATTTTTTCATATCGTTACCCTTTTTCAGGCAATAGTTGACCCTGACCGAAGCGGCCAATTTTAttttt                                                                                                          >  2:48670/1‑140 (MQ=255)                            tGCGCCAGGCGATATCTTCATTCAGTTCTTCGCCTAATTTCCCGCGAATATCATAGGCTTTAAAGCAGGTTAATTTTTTCATATCGTTACCCTTTTTCAGGCAATAGTTGACCCTGACCGAAGCGGCCAATTTTAttttt                                                                                                          >  2:274123/1‑140 (MQ=255)                             gcgcCAGGCGATATCTTCATTCAGTTCTTCGCCTAATTTCCCGCGAATATCATAGGCTTTAAAGCAGGTTAATTTTTTCATATCGTTACCCTTTTTCAGGCAATAGTTGACCCTGACCGAAGCGGCCAATTTTATTTTTg                                                                                                         <  2:840360/140‑1 (MQ=255)                               gcCAGGCGATATCTTCATTCAGTTCTTCGCCTAATTTCCCGCGAATATCATAGGCTTTAAAGCAGGTTAATTTTTTCATATCGTTACCCTTTTTCAGGCAATAGTTGACCCTGACCGAAGCGGCCAATTTTATTTTTGTc                                                                                                       <  2:251371/140‑1 (MQ=255)                               gcCAGGCGATATCTTCATTCAGTTCTTCGCCTAATTTCCCGCGAATATCATAGGCTTTAAAGCAGGTTAATTTTTTCATATCGTTACCCTTTTTCAGGCAATAGTTGACCCTGACCGAAGCGGCCAATTTTATTTTTGTc                                                                                                       <  2:163287/140‑1 (MQ=255)                                  aGGCGATATCTTCATTCAGTTCTTCGCCTAATTTCCCGCGAATATCATAGGCTTTAAAGCAGGTTAATTTTTTCATATCGTTACCCTTTTTCAGGCAATAGTTGACCCTGACCGAAGCGGCCAATTTTATTTTTGTCGtt                                                                                                    >  1:831283/1‑140 (MQ=255)                                   ggCGATATCTTCATTCAGTTCTTCGCCTAATTTCCCGCGAATATCATAGGCTTTAAAGCAGGTTAATTTTTTCATATCGTTACCCTTTTTCAGGCAATAGTTGACCCTGACCGAAGCGGCCAATTTTATTTTTGTCGTTa                                                                                                   >  1:80643/1‑140 (MQ=255)                                    gCGATATCTTCATTCAGTTCTTCGCCTAATTTCCCGCGAATATCATAGGCTTTAAAGCAGGTTAATTTTTTCATATCGTTACCCTTTTTCAGGCAATAGTTGACCCTGACCGAAGCGGCCAATTTTATTTTTGTCGTTAt                                                                                                  >  1:401570/1‑140 (MQ=255)                                              cattcaGTTCTTCGCCTAATTTCCCGCGAATATCATAGGCTTTAAAGCAGGTTAATTTTTTCATATCGTTACCCTTTTTCAGGCAATAGTTGACCCTGACCGAAGCGGCCAATTTTATTTTTGTCGt                                                                                                     <  2:820018/127‑1 (MQ=255)                                              cattcaGTTCTTCGCCTAATTTCCCGCGAATATCATAGGCTTTAAAGCAGGTTAATTTTTTCATATCGTTACCCTTTTTCAGGCAATAGTTGACCCTGACCGAAGCGGCCAATTTTATTTTTGTAGTTATTTACACcgcg                                                                                        <  1:1034531/140‑1 (MQ=255)                                              cattcaGTTCTTCGCCTAATTCCCCGCGAATATCATAGGCTTTAAAGCAGGTTAATTTTTTCATATCGTTACCCTTTTTCAGGCAATAGTTGACCCTGACCGAAGCGGCCAATTTTATTTTTGTCGt                                                                                                     >  1:820018/1‑127 (MQ=255)                                                ttcaGTTCTTCGCCTAATTTCCCGCGAATATCATAGGCTTTAAAGCAGGTTAATTTTTTCATATCGTTACCCTTTTTCAGGCAATAGTTGACCCTGACCGAAGCGGCCAATTTTATTTTTGTCGTTATTTACACCGCGGt                                                                                      <  1:204710/140‑1 (MQ=255)                                                   aGTTCTTCGCCTAATTTCCCGCGAATATCATAGGCTTTAAAGCAGGTTAATTTTTTCATATCGTTACCCTTTTTCAGGCAATAGTTGACCCTGACCGAAGCGGCCAATTTTATTTTTGTCGTTATTTACACCGCGGTTTc                                                                                   <  2:1000742/140‑1 (MQ=255)                                                    gTTCTTCGCCTAATTTCCCGCGAATATCATAGGCTTTAAAGCAGGTTAATTTTTTCATATCGTTACCCTTTTTCAGGCAATAGTTGACCCTGACCGAAGCGGCCAATTTTATTTTTGTCGTTATTTACACCGCGGTTTCg                                                                                  <  2:371353/140‑1 (MQ=255)                                                     ttcttcGCCTAATTTCCCGCGAATATCATAGGCTTTAAAGCAGGTTAATTTTTTCATATCGTTACCCTTTTTCAGGCAATAGTTGACCCTGACCGAAGCGGCCAATTTTATTTTTGTCGTTATTTACACCGCGGTTTCGc                                                                                 >  1:1085059/1‑140 (MQ=255)                                                     ttcttcGCCTAATTTCCCGCGAATATCATAGGCTTTAAAGCAGGTTAATTTTTTCATATCGTTACCCTTTTTCAGGCAATAGTTGACCCTGACCGAAGCGGCCAATTTTATTTTTGTCGTTATTTACACCGCGGTTTCGc                                                                                 >  2:690053/1‑140 (MQ=255)                                                       cttcGCCTAATTTCCCGCGAATATCATAGGCTTTAAAGCAGGTTAATTTTTTCATATCGTTACCCTTTTTCAGGCAATAGCTGACCCTGACCGAAGCGGCCAATTTTATTTTTGTCGTTATTTACACCGCGGTTTCGcat                                                                               >  1:796281/1‑140 (MQ=255)                                                             cTAATTTCCCGCGAATATCATAGGCTTTAAAGCAGGTTAATTTTTTCATATCGTTACCCTTTTTCAGGCAATAGTTGACCCTGACCGAAGCGGCCAATTTTATTTTTGTCGTTATTTACACCGCGGTTTCGCATTCATTg                                                                         >  1:431734/1‑140 (MQ=255)                                                              tAATTTCCCGCGAATATCATAGGCTTTAAAGCAGGTTAATTTTTTCATATCGTTACCCTTTTTCAGGCAATAGTTGACCCTGACCGAAGCGGCCAATTTTATTTTTGTCGTTATTTACACCGCGGTTTCGCATTCATTGc                                                                        <  1:360946/140‑1 (MQ=255)                                                                      cgcgAATATCATAGGCTTTAAAGCAGGTTAATTTTTTCATATCGTTACCCTTTTTCAGGCAATAGTTGACCCTGACCGAAGCGGCCAATTTTATTTTTGTCGTTATTTACACCGCGGTTTCGCATTCATTGCCTGATGCg                                                                <  1:804432/140‑1 (MQ=255)                                                                      cgcgAATATCATAGGCTTTAAAGCAGGTTAATTTTTTCATATCGTTACCCTTTTTCAGGCAATAGTTGACCCTGACCGAAGCGGCCAATTTTATTTTTGTCGTTATTTACACCGCGGTTTCGCATTCATTGCCTGATGCg                                                                <  2:80643/140‑1 (MQ=255)                                                                               cATAGGCTTTAAAGCAGGTTAATTTTTTCATATCGTTACCCTTTTTCAGGCAATAGTTGACCCTGACCGAAGCGGCCAATTTTATTTTTGTCGTTATTTACACCGCGGTTTCGCATTCATTGCCTGATGCGACGTTTaca                                                       <  2:1085059/140‑1 (MQ=255)                                                                               cATAGGCTTTAAAGCAGGTTAATTTTTTCATATCGTTACCCTTTTTCAGGCAATAGTTGACCCTGACCGAAGCGGCCAATTTTATTTTTGTCGTTATTTACACCGCGGTTTCGCATTCATTGCCTGATGCGACGTTTaca                                                       <  1:690053/140‑1 (MQ=255)                                                                               cATAGGCTTTAAAGCAGGTTAATTTTTTCATATCGTTACCCTTTTTCAGGCAATAGCTGACCCTGACCGAAGCGGCCAATTTTATTTTTGTCGTTATTTACACCGCGGTTTCGCATTCATTGCCTGATGCGACGTTTaca                                                       <  2:796281/140‑1 (MQ=255)                                                                                aTAGGCTTTAAAGCAGGTTAATTTTTTCATATCGTTACCCTTTTTCAGGCAATAGTTGACCCTGACCGAAGCGGCCAATTTTATTTTTGTCGTTATTTACACCGCGGTTTCGCATTCATTGCCTGATGCGACGTTTacac                                                      <  2:831283/140‑1 (MQ=255)                                                                                  aGGCTTTAAAGCAGGTTAATTTTTTCATATCGTTACCCTTTTTCAGGCAATAGTTGACCCTGACCGAAGCGGCCAATTTTATTTTTGTCGTTATTTACACCGCGGTTTCGCATTCATTGCCTGATGCGACGTTTACAccc                                                    <  2:159111/140‑1 (MQ=255)                                                                                       ttAAAGCAGGTTAATTTTTTCATATCGTTACCCTTTTTCAGGCAATAGTTGACCCTGACCGAAGCGGCCAATTTTATTTTTGTCGTTATTTACACCGCGGTTTCGCATTCATTGCCTGATGCGACGTTTACACccgtccg                                               >  1:676876/1‑140 (MQ=255)                                                                                        tAAAGCAGGTTAATTTTTTCATATCGTTACCCTTTTTCAGGCAATAGTTGACCCTGACCGAAGCGGCCAATTTTATTTTTGTCGTTATTTACACCGCGGTTTCGCATTCATTGCCTGATGCGACGTTTACACccgtccgt                                              >  2:727865/1‑140 (MQ=255)                                                                                        tAAAGCAGGTTAATTTTTTCATATCGTTACCCTTTTTCAGGCAATAGTTGACCCTGACCGAAGCGGCCAATTTTATTTTTGTCGTTATTTACACCGCGGTTTCGCATTCATTGCCTGATGCGACGTTTACACccgtccgt                                              >  2:768234/1‑140 (MQ=255)                                                                                        tAAAGCAGGTTAATTTTTTCATATCGTTACCCTTTTTCAGGCAATAGTTGACCCTGACCGAAGCGGCCAATTTTATTTTTGTCGTTATTTACACCGCGGTTTCGCATTCATTGCCTGATGCGACGTTTACACccgtccgt                                              >  1:211049/1‑140 (MQ=255)                                                                                         aaaGCAGGTTAATTTTTTCATATCGTTACCCTTTTTCAGGCAATAGTTGACCCTGACCGAAGCGGCCAATTTTATTTTTGTCGTTATTTACACCGCGGTTTCGCATTCATTGCCTGATGCGACGTTTACACCCGTCCGTa                                             >  2:796442/1‑140 (MQ=255)                                                                                            gCAGGTTAATTTTTTCATATCGTTACCCTTTTTCAGGCAATAGTTGACCCTGACCGAAGCGGCCAATTTTATTTTTGTCGTTATTTACACCGCGGTTTCGCATTCATTGCCTGATGCGACGTTTACACCCGTCCGTAGCg                                          <  2:996512/140‑1 (MQ=255)                                                                                                          tCATATCGTTACCCTTTTTCAGGCAATAGTTGACCCTGACCGAAGCGGCCAATTTTATTTTTGTCGTTATTTACACCGCGGTTTCGCATTCATTGCCTGATGCGACGTTTACACCCGTCCGTAGCGATCCGCGAAACGca                            <  2:211049/140‑1 (MQ=255)                                                                                                          tCATATCGTTACCCTTTTTCAGGCAATAGTTGACCCTGACCGAAGCGGCCAATTTTATTTTTGTCGTTATTTACACCGCGGTTTCGCATTCATTGCCTGATGCGACGTTTACACCCGTCCGTAGCGATCCGCGAAACGca                            <  1:48670/140‑1 (MQ=255)                                                                                                                          tttCAGGCAATAGTTGACCCTGACCGAAGCGGCCAATTTTATTTTTGTCGTTATTTACACCGCGGTTTCGCATTCATTGCCTGATGCGACGTTTACACCCGTCCGTAGCGATCCGCGAAACGCACCAcatca                    >  2:105625/1‑132 (MQ=255)                                                                                                                          tttCAGGCAATAGTTGACCCTGACCGAAGCGGCCAATTTTATTTTTGTCGTTATTTACACCGCGGTTTCGCATTCATTGCCTGATGCGACGTTTACACCCGTCCGTAGCGATCCGCGAAACGCACCAcatca                    <  1:105625/132‑1 (MQ=255)                                                                                                                          tttCAGGCAATAGTTGACCCTGACCGAAGCGGCCAATTTTATTTTTGTCGTTATTTACACCGCGGTTTCGCATTCATTGCCTGATGCGACGTTTACACCCGTCCGTAGCGATCCGCGAAACGCACCACATCATCCTCTTc            <  1:274123/140‑1 (MQ=255)                                                                                                                          tttCAGGCAATAGTTGACCCTGACCGAAGCGGCCAATTTTATTTTTGTCGTTATTTACACCGCGGTTTCGCATTCATTGCCTGATGCGACGTTTACACCCGTCCGTAGCGATCCGCGAAACGCACCACATCATCCTCTTc            >  1:1051595/1‑140 (MQ=255)                                                                                                                          tttCAGGCAATAGTTGACCCTGACCGAAGCGGCAAATTTTATTTTTGTCGTTATTTACACCGCGGTTTCGCATTCATTGCCTGATGCGACGTTTACACCCGTCCGTAGCGATCCGCGAAACGCACCACATCATCCTCTTc            >  2:933180/1‑140 (MQ=255)                                                                                                                           ttCAGGCAATAGTTGACCCTGACCGAAGCGGCCAATTTTATTTTTGTCGTTATTTACACCGCGGTTTCGCATTCATTGCCTGATGCGACGTTTACACCCGTCCGTAGCGATCCGCGAAACGCACCACATCATCCTCTTCg           >  2:394141/1‑140 (MQ=255)                                                                                                                                gCAATAGTTGACCCTGACCGAAGCGGCCAATTTTATTTTTGTCGTTATTTACACCGCGGTTTCGCATTCATTGCCTGATGCGACGTTTACACCCGTCCGTAGCGATCCGCGAAACGCACCACATCATCCTCTTCGAGATa      >  1:873778/1‑140 (MQ=255)                                                                                                                                    tAGTTGACCCTGACCGAAGCGGCCAATTTTATTTTTGTCGTTATTTACACCGCGGTTTCGCATTCATTGCCTGATGCGACGTTTACACCCGTCCGTAGCGATCCGCGAAACGCACCACATCATCCTCTTCGAGATAagag  <  2:873778/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                      TGAGAAATTCGCCATAGGCGCGACCAATGCGCCAGGCGATATCTTCATTCAGTTCTTCGCCTAATTTCCCGCGAATATCATAGGCTTTAAAGCAGGTTAATTTTTTCATATCGTTACCCTTTTTCAGGCAATAGTTGGCCCTGACCGAAGCGGCCAATTTTATTTTTGTCGTTATTTACACCGCGGTTTCGCATTCATTGCCTGATGCGACGTTTACACCCGTCCGTAGCGATCCGCGAAACGCACCACATCATCCTCTTCGAGATAAGAG  >  NC_000913/2122871‑2123141

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