Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,630,986	A→G	intergenic (‑73/+16)	ydfI ← / ← ydfJ	putative NAD‑dependent D‑mannonate oxidoreductase/pseudogene, MFS transporter family, interrupted by Qin prophage,Phage or Prophage Related, putative transport protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,630,986	0	A	G	100.0%	207.7 / NA	60	intergenic (‑73/+16)	ydfI/ydfJ	putative NAD‑dependent D‑mannonate oxidoreductase/pseudogene, MFS transporter family, interrupted by Qin prophage,Phage or Prophage Related, putative transport protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (31/29);  total (31/29)

AGCCAGGTTATTAAGATCATAAACAGGGAGTGTCGCTTTTGCTGATAACAAATTATTTCCCATAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCA  >  NC_000913/1630851‑1631110                                                                                                                                        |                                                                                                                             aGCCAGGTTATTAAGATCATAAACAGGGAGTGTCGCTTTTGCTGATAACAAATTATTTCCCATAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGaa                                                                                                                           <  1:701981/139‑1 (MQ=255)                aTCATAAACAGGGAGTGTCGCTTTTGCTGATAACAAATTATTTCCCATAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGAAg                                                                                                                          >  1:691529/1‑125 (MQ=255)                aTCATAAACAGGGAGTGTCGCTTTTGCTGATAACAAATTATTTCCCATAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGAAg                                                                                                                          <  2:691529/125‑1 (MQ=255)                aTCATAAACAGGGAGTGTCGCTTTTGCTGATAACAAATTATTTCCCATAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGAAAGGGAAGGGATGAAGGGTTTa                                                                                                            >  2:655630/1‑139 (MQ=255)                 tCATAAACAGGGAGTGTCGCTTTTGCTGATAACAAATTATTTCCCATAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTAAgg                                                                                                                 <  1:55392/133‑1 (MQ=255)                 tCATAAACAGGGAGTGTCGCTTTTGCTGATAACAAATTATTTCCCATAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTAc                                                                                                           >  2:497352/1‑139 (MQ=255)                 tCATAAACAGGGAGTGTCGCTTTTGCTGATAACAAATTATTTCCCATAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTAc                                                                                                           >  1:493664/1‑139 (MQ=255)                 tCATAAACAGGGAGTGTCGCTTTTGCTGATAACAAATTATTTCCCATAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCAGGTAAGGGAAGGGATTAAgg                                                                                                                 >  2:55392/1‑133 (MQ=255)                  cATAAACAGGGAGTGTCGCTTTTGCTGATAACAAATTATTTCCCATAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTACa                                                                                                          <  2:493664/139‑1 (MQ=255)                          gggAGTGTCGCTTTTGCTGATAACAAATTATTTCCCATAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTa                                                                                                            >  2:817324/1‑129 (MQ=255)                          gggAGTGTCGCTTTTGCTGATAACAAATTATTTCCCATAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTa                                                                                                            <  1:817324/129‑1 (MQ=255)                           ggAGTGTCGCTTTTGCTGATAACAAATTATTTCCCATAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGaa                                                                                                                           >  1:307448/1‑113 (MQ=255)                           ggAGTGTCGCTTTTGCTGATAACAAATTATTTCCCATAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGaa                                                                                                                           <  2:307448/113‑1 (MQ=255)                           ggACTGTCGCATATGCTGATCACAAATTATTTTCCGTAACAACTCCTGAAAGATAAATATGGCAGGCTATATGTTTTGCTATATGAATAAAAATCCCATCTCCGGTTAGGGAAGGGATTAAGGCTGTCCAGACTGCTgg                                                                                                 >  2:805430/1‑139 (MQ=255)                                gtCGCTTTTGCTGATAACAAATTATTTCCCATAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAgtt                                                                                            >  2:58569/1‑137 (MQ=255)                                gtCGCTTTTGCTGATAACAAATTATTTCCCATAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGt                                                                                            >  2:443838/1‑139 (MQ=255)                                 tCGCTTTTGCTGATAACAAATTATTTCCCATAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGaa                                                                                                                           <  2:587640/107‑1 (MQ=255)                                 tCGCTTTTGCTGATAACAAATTATTTCCCATAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGaa                                                                                                                           >  1:587640/1‑107 (MQ=255)                                 tCGCTTTTGCTGATAACAAATTATTTCCCATAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGtt                                                                                           >  1:145437/1‑139 (MQ=255)                                     ttttGCTGATAACAAATTATTTCCCATAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTgcgc                                                                                       >  2:811495/1‑139 (MQ=255)                                      tttGCTGATAACAAATTATTTCCCATAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCa                                                                                      <  1:811495/139‑1 (MQ=255)                                                    attattTCCCATAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCa                                                                                      >  1:1083951/1‑125 (MQ=255)                                                    attattTCCCATAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCa                                                                                      <  2:1083951/125‑1 (MQ=255)                                                              aTAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCt                                                              >  2:323075/1‑139 (MQ=255)                                                                aaCAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTc                                                            >  1:760652/1‑139 (MQ=255)                                                                  cAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCg                                                          <  1:323075/139‑1 (MQ=255)                                                                    aTTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCa                                                        <  2:290048/139‑1 (MQ=255)                                                                      tCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAg                                                                                     <  2:998260/108‑1 (MQ=255)                                                                      tCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAg                                                                                     >  1:998260/1‑108 (MQ=255)                                                                        cTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGaa                                                                                                                           >  2:148159/1‑68 (MQ=255)                                                                        cTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGaa                                                                                                                           <  1:148159/68‑1 (MQ=255)                                                                              tataaatatGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCa                                              >  1:1061694/1‑139 (MQ=255)                                                                                taaatatGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAgc                                            <  1:655630/139‑1 (MQ=255)                                                                                 aaatatGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAgcg                                           <  2:760652/139‑1 (MQ=255)                                                                                       ggCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTg                                                                                                  >  1:309964/1‑78 (MQ=255)                                                                                       ggCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTg                                                                                                  <  2:309964/78‑1 (MQ=255)                                                                                       ggCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTa                                     >  2:420047/1‑139 (MQ=255)                                                                                        gCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTaa                                    >  1:240783/1‑139 (MQ=255)                                                                                            gCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCa                                <  1:53709/139‑1 (MQ=255)                                                                                             cTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCAc                               <  1:58569/139‑1 (MQ=255)                                                                                                    ttttGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCGCGTACAACACGCGGTTGATCGTCCGCAGCGTCGTCCAGCGCACTTAAATCACGGTCttt                        <  2:363598/139‑1 (MQ=255)                                                                                                       tGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTACa                                                                                                          <  2:809651/54‑1 (MQ=255)                                                                                                       tGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTACa                                                                                                          >  1:809651/1‑54 (MQ=255)                                                                                                       tGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCAc                     >  2:243800/1‑139 (MQ=255)                                                                                                       tGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCAc                     >  2:885065/1‑139 (MQ=255)                                                                                                          tataTGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTgcg                                                                                        <  1:847921/69‑1 (MQ=255)                                                                                                          tataTGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTgcg                                                                                        >  2:847921/1‑69 (MQ=255)                                                                                                            tatGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGt                                                                    <  1:1052676/87‑1 (MQ=255)                                                                                                            tatGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGt                                                                    >  2:1052676/1‑87 (MQ=255)                                                                                                                aaTAAAAATCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCAc                               <  1:994044/120‑1 (MQ=255)                                                                                                                aaTAAAAATCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCAc                               >  2:994044/1‑120 (MQ=255)                                                                                                                   aaaaaTCCCCTCTCCGGTAAGGGTAGGGGTTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGACTTTCACCTCTGGCAtttt         >  1:264614/1‑139 (MQ=255)                                                                                                                   aaaaaTCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCAtttt         <  2:164201/139‑1 (MQ=255)                                                                                                                     aaaTCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCa       <  1:885065/139‑1 (MQ=255)                                                                                                                         ccccTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAgcgc   <  2:240783/139‑1 (MQ=255)                                                                                                                          cccTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTgcgc                                                                                       >  1:172503/1‑54 (MQ=255)                                                                                                                          cccTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTgcgc                                                                                       <  2:172503/54‑1 (MQ=255)                                                                                                                          cccTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTgcg                                                                                        <  2:1009461/53‑1 (MQ=255)                                                                                                                          cccTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTgcg                                                                                        >  1:1009461/1‑53 (MQ=255)                                                                                                                          cccTCTCCGGTAAGGGAAGGGATTAAGGGTCTACAGCATTCTGGAAGGTTGCGAAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTAACAGCGAACTTAAATCACGGTCTTTCACATCTGGCATTTTCAGCGCa  <  2:264614/139‑1 (MQ=255)                                                                                                                                        |                                                                                                                             AGCCAGGTTATTAAGATCATAAACAGGGAGTGTCGCTTTTGCTGATAACAAATTATTTCCCATAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCA  >  NC_000913/1630851‑1631110

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