Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,313,825	A→C	intergenic (+26/‑43)	rpiB → / → yjdP	ribose 5‑phosphate isomerase B/allose 6‑phosphate isomerase/putative periplasmic protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,313,825	0	A	C	100.0%	169.1 / NA	47	intergenic (+26/‑43)	rpiB/yjdP	ribose 5‑phosphate isomerase B/allose 6‑phosphate isomerase/putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (25/22);  total (25/22)

CGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCG  >  NC_000913/4313677‑4313972                                                                                                                                                     |                                                                                                                                                    cGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGCGt                                                                                                                                                   <  2:236261/151‑1 (MQ=255)      ggTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGCGTAATg                                                                                                                                               <  2:172491/150‑1 (MQ=255)        ttGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGCGTAATGAt                                                                                                                                             <  2:538331/150‑1 (MQ=255)          ggCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGCGTAATGATTaa                                                                                                                                          <  1:856953/151‑1 (MQ=255)            ccTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGCGTAATGATTAAtg                                                                                                                                        >  1:860868/1‑151 (MQ=255)             cTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGCGTAATGATTAAtgt                                                                                                                                       >  1:16655/1‑151 (MQ=255)             cTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGCGTAATGATTAAtgt                                                                                                                                       >  1:314517/1‑151 (MQ=255)                gAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGCGTAATGATTAATGTgg                                                                                                                                     >  1:475733/1‑150 (MQ=255)                 aaCTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGCGTAATGATTAATGTGGAt                                                                                                                                   <  1:554846/151‑1 (MQ=255)                 aaCTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGCGTAATGATTAATGTGGAt                                                                                                                                   >  2:598284/1‑151 (MQ=255)                      gCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGCGTAATGATTAATGTGGATAGAGt                                                                                                                              >  2:669631/1‑151 (MQ=255)                        aaaaaTGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGCGTAATGATTAATGTGGATAGAGtt                                                                                                                             <  2:568110/150‑1 (MQ=255)                         aaaaTGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGCGTAATGATTAATGTGGATAGAGTTTc                                                                                                                           <  2:6669/151‑1 (MQ=255)                               aTTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGCGTAATGATTAATGTGGATAGAGTTTCTTTTTg                                                                                                                     <  2:905926/151‑1 (MQ=255)                                    ggATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGCGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGt                                                                                                                 <  2:96243/150‑1 (MQ=255)                                       tGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGCGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTgg                                                                                                              >  1:640502/1‑150 (MQ=255)                                                gggCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGCGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAAc                                                                                                    <  2:915358/151‑1 (MQ=255)                                                      acaGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGCGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTc                                                                                              >  1:339751/1‑151 (MQ=255)                                                       caGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGCGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTcc                                                                                             <  1:254588/151‑1 (MQ=255)                                                           aCGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGCGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCt                                                                                          >  1:375207/1‑150 (MQ=255)                                                             gAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGCGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTAtt                                                                                       <  2:67660/151‑1 (MQ=255)                                                             gAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGCGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTAtt                                                                                       >  2:780510/1‑151 (MQ=255)                                                             gAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGCGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTAtt                                                                                       >  2:117593/1‑151 (MQ=255)                                                              aaGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGCGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTAttt                                                                                      >  1:680946/1‑151 (MQ=255)                                                                   ggTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGCGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCtttt                                                                                 >  1:789708/1‑151 (MQ=255)                                                                    gtcgtcATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGCGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCtttt                                                                                 >  1:292403/1‑150 (MQ=255)                                                                       gtcATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGCGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTc                                                                             <  1:131182/151‑1 (MQ=255)                                                                       gtcATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGCGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTc                                                                             <  1:569161/151‑1 (MQ=255)                                                                           tcaACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGCGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCAccc                                                                         >  2:393234/1‑151 (MQ=255)                                                                            caacaaCGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGCGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCt                                                                        <  1:669612/151‑1 (MQ=255)                                                                              acaacGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGCGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGc                                                                      <  2:38241/151‑1 (MQ=255)                                                                               caacGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGCGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCt                                                                     <  2:48374/151‑1 (MQ=255)                                                                                     gTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGCGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCt                                                               >  2:144680/1‑151 (MQ=255)                                                                                                tACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGCGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTc                                                    >  2:813585/1‑151 (MQ=255)                                                                                                 aCGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGCGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTcc                                                   <  1:213452/151‑1 (MQ=255)                                                                                                    gCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGCGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGc                                                >  1:100476/1‑151 (MQ=255)                                                                                                    gCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGCGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGc                                                >  2:525500/1‑151 (MQ=255)                                                                                                     cAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGCGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGcc                                               <  2:307087/151‑1 (MQ=255)                                                                                                            gcagcGGAGAAATTGAGATTCATCCACTACTTGCATGGATGCGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCg                                        >  1:111828/1‑151 (MQ=255)                                                                                                                 ggAGAAATTGAGATTCATCCACTACTTGCATGGATGCGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGAtatt                                   <  1:633301/151‑1 (MQ=255)                                                                                                                        ttGAGATTCATCCACTACTTGCATGGATGCGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATcgacg                            >  1:189982/1‑151 (MQ=255)                                                                                                                          gagaTTCATCCACTACTTGCATGGATGCGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCgacgac                          <  1:777741/151‑1 (MQ=255)                                                                                                                           agaTTCATCCACTACTTGCATGGATGCGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCgacgac                          >  1:885670/1‑150 (MQ=255)                                                                                                                                cATCCACTACTTGCATGGATGCGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCAtt                    >  1:491403/1‑151 (MQ=255)                                                                                                                                 aTCCACTACTTGCATGGATGCGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTg                   <  1:505870/151‑1 (MQ=255)                                                                                                                                       tactTGCATGGATGCGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAat              >  1:262108/1‑150 (MQ=255)                                                                                                                                                  aTGCGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCg  >  1:46027/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                    CGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCG  >  NC_000913/4313677‑4313972

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