Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,577,522	A→C	M33R (ATG→AGG)	yddB ←	putative TonB‑dependent receptor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,577,522	0	A	C	100.0%	163.7 / NA	47	M33R (ATG→AGG)	yddB	putative TonB‑dependent receptor
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (24/23);  total (24/23)

ACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCATCTCTTCAAAAAAA‑TGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTA  >  NC_000913/1577399‑1577650                                                                                                                            |                                                                                                                                  aCCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCCTCTCTTCaa                                                                                                                          >  1:538007/1‑133 (MQ=255)      aGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCCTCTCTTCaaaaaaa                                                                                                                     <  1:309754/133‑1 (MQ=255)             gttCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCCTCTCTTCAAAAAAA‑TGCATGt                                                                                                             >  1:303275/1‑133 (MQ=255)                  cagcagATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCCTCTCTTCAAAAAAA‑TGCATGTACATg                                                                                                        >  2:94442/1‑133 (MQ=255)                        aTCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCCTCTCTTCAAAAAAA‑TGCATGTACATGTTTTTg                                                                                                  >  1:367076/1‑133 (MQ=255)                              gATATTGCCATTCCCGGTCGGTAAACGCTC‑ATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCCTCTCTTCAAAAAAA‑TGCATGTACATGTTTTTGTCATCGa                                                                                           <  1:38818/133‑1 (MQ=255)                                    ggcATTCCCCGTCGGGAAACGCTCGATGTTTTCACTGGTGTGATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCCTCTCTTCAAAAAAA‑TGCATGTACATGTTTTTGTCATCGACGGcc                                                                                      <  2:303275/131‑1 (MQ=255)                                    gCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCCTCTCTTCAAAAAAA‑TGCATGTACATGTTTTTGTCATCGACGGcc                                                                                      <  2:544529/133‑1 (MQ=255)                                    gCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCCTCTCTTCAAAAAAA‑TGCATGTACATGTTTTTGTCATCGACGGcc                                                                                      <  1:441762/133‑1 (MQ=255)                                      cATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCCTCTCTTCAAAAAAA‑TGCATGTACATGTTTTTGTCATCGACGGCCTg                                                                                    <  2:117346/133‑1 (MQ=255)                                         tCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCCTCTCTTCAAAAAAA‑TGCATGTACATGTTTTTGTCATCGACGGCCTgcgc                                                                                 >  1:441909/1‑133 (MQ=255)                                         tCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCCTCTCTTCAAAAAAA‑TGCATGTACATGTTTTTGTCATCGACGGCCTgcgc                                                                                 >  2:245220/1‑133 (MQ=255)                                              gtcggtAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCCTCTCTTCAAAAAAA‑TGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACAt                                                                            >  2:39255/1‑133 (MQ=255)                                                 ggtAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCCTCTCTTCAAAAAAA‑TGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAg                                                                         <  1:577328/133‑1 (MQ=255)                                                              tGCTTTCACTGGTGTAATGCGTGTTGCCGTTTGCGGGTACAGGGACAGGAGCATAGACCGTCCTCTCTTCAAAAAAA‑TGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAAt                                                            <  1:571791/133‑1 (MQ=255)                                                              tGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCCTCTCTTCAAAAAAA‑TGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAAt                                                            <  1:245220/133‑1 (MQ=255)                                                                   tCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCCTCTCTTCAAAAAAA‑TGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGc                                                       <  1:94442/133‑1 (MQ=255)                                                                    cACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCCTCTCTTCAAAAAAA‑TGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGcc                                                      <  2:538007/133‑1 (MQ=255)                                                                    cACTGGTGTAATGCGTGTTGCCGTTTACCGGTACAGGGACAGGAGCATAGACCGTCCTCTCTTCAAAAAAA‑TGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGcc                                                      >  1:391434/1‑133 (MQ=255)                                                                        ggTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCCTCTCTTCAAAAAAA‑TGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGa                                                  >  1:65973/1‑133 (MQ=255)                                                                          tgtAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCCTCTCTTCAAAAAAA‑TGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAAt                                                >  2:52144/1‑133 (MQ=255)                                                                           gtAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCCTCTCTTCAAAAAAA‑TGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATa                                               <  2:441909/133‑1 (MQ=255)                                                                            tAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCCTCTCTTCAAAAAAA‑TGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATaa                                              >  1:61735/1‑133 (MQ=255)                                                                               tGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCCTCTCTTCAAAAAAA‑TGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGaa                                           <  2:356314/133‑1 (MQ=255)                                                                               tGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCCTCTCTTCAAAAAAA‑TGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGaa                                           <  1:39255/133‑1 (MQ=255)                                                                                 cGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCCTCTCTTCAAAAAAA‑TGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACt                                         >  2:26237/1‑133 (MQ=255)                                                                                 cGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCCTCTCTTCAAAAAAA‑TGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACt                                         >  2:532648/1‑133 (MQ=255)                                                                                   tgtTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCCTCTCTTCAAAAAAA‑TGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCg                                       >  2:78705/1‑133 (MQ=255)                                                                                   tgtTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCCTCTCTTCAAAAAAA‑TGCATGTACATGTTTTTGTCATCGACGGCCTGCGCAACATCAGCGCCACATAGAATGACGCGAGGAATAAGAACTCg                                       >  2:260997/1‑133 (MQ=255)                                                                                     ttGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCCTCTCTTCAAAAAAA‑TGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCt                                     <  1:26237/133‑1 (MQ=255)                                                                                     ttGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCCTCTCTTCAAAAAAA‑TGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCt                                     >  1:325177/1‑133 (MQ=255)                                                                                      tGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCCTCTCTTCAAAAAAA‑TGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCtt                                    >  1:588038/1‑133 (MQ=255)                                                                                        ccGTTTACGGGTACAGGGACAGGAGCATAGACCGTCCTCTCTTCAAAAAAA‑TGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCa                                  <  1:52144/133‑1 (MQ=255)                                                                                         cGTTTACGGGTACAGGGACAGGAGCATAGACCGTCCTCTCTTCAAAAAAA‑TGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCat                                 <  2:65973/133‑1 (MQ=255)                                                                                         cGTTTACGGGTACAGGGACAGGAGCATAGACCGTCCTCTCTTCAAAAAAA‑TGCATGTACATGTTTTTGTCATCGACGGCATGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCat                                 <  2:61735/133‑1 (MQ=255)                                                                                                gggTACAGGGCCAGGTGCATAGACCGTCCTCTCTTCAAAAAGA‑TGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCAtttattta                          <  2:290317/133‑1 (MQ=255)                                                                                                gggTACAGGGACAGGAGCATAGACCGTCCTCTCTTCAAAAAAA‑TGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCAtttattta                          <  2:325177/133‑1 (MQ=255)                                                                                                gggTACAGGGACAGGAGCATAGACCGTCCTCTCTTCAAAAAAA‑TGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCAtttattta                          <  2:588038/133‑1 (MQ=255)                                                                                                     cAGGGACAGGAGCATAGACCGTCCTCTCTTCAAAAAAA‑TGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAGATGACGCCGGGAATAAGAACTCGCTCCATTGATTTATAtct                     >  2:193676/1‑133 (MQ=255)                                                                                                     cAGGGACAGGAGCATAGACCGTCCTCTCTTCAAAAAAA‑TGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATAtct                     >  2:41238/1‑133 (MQ=255)                                                                                                      aGGGACAGGAGCATAGACCGTCCTCTCTTCAAAAAAATTGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATAtct                     <  1:463729/133‑1 (MQ=255)                                                                                                        ggACAGGAGCATAGACCGTCCTCTCTTCAAAAAAA‑TGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCg                  >  1:473724/1‑133 (MQ=255)                                                                                                        ggACAGGAGCATAGACCGTCCTCTCTTCAAAAAAA‑TGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCg                  >  1:115716/1‑133 (MQ=255)                                                                                                        ggACAGGAGCATAGACCGTCCTCTCTTCAAAAAAA‑TGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCg                  >  1:601131/1‑133 (MQ=255)                                                                                                                      aCCGTCCTCTCTTCAAAAAAA‑TGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCttt    >  2:6375/1‑133 (MQ=255)                                                                                                                        cGTCCTCTCTTCAAAAAAA‑TGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCGTTTGTTTTTATCTCCGGACAGCAGGGCtttca  >  2:282656/1‑131 (MQ=255)                                                                                                                        cGTCCTCTCTTCAAAAAAA‑TGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTa  <  2:473724/133‑1 (MQ=255)                                                                                                                            |                                                                                                                                  ACCGCAGGGTTGGTTCTCAGCAGATCGCTGATATTGCCATTCCCGGTCGGTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCATCTCTTCAAAAAAA‑TGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTA  >  NC_000913/1577399‑1577650

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