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%	200.6 / NA	57	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 (34/23);  total (34/23)

TCATAAACAGGGAGTGTCGCTTTTGCTGATAACAAATTATTTCCCATAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCAGA  >  NC_000913/1630867‑1631112                                                                                                                        |                                                                                                                               tCATAAACAGGGAGTGTCGCTTTTGCTGATAACAAATTATTTCCCATAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTAc                                                                                                             >  1:673545/1‑139 (MQ=255) tCATAAACAGGGAGTGTCGCTTTTGCTGATAACAAATTATTTCCCATAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTAc                                                                                                             >  1:1007349/1‑139 (MQ=255)  cATAAACAGGGAGTGTCGCTTTTGCTGATAACAAATTATTTCCCATAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTata                                                                                                            <  1:574940/139‑3 (MQ=255)  cATAAACAGGGAGTGTCGCTTTTGCTGATAACAAATTATTTCCCATAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTata                                                                                                            >  2:574940/1‑137 (MQ=255)   aTAAACAGGGAGTGTCGCTTTTGCTGATAACAAATTATTTCCCATAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAggga                                                                                                                              >  2:56334/1‑120 (MQ=255)   aTAAACAGGGAGTGTCGCTTTTGCTGATAACAAATTATTTCCCATAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAggga                                                                                                                              <  1:56334/120‑1 (MQ=255)          gggAGTGTCGCTTTTGCTGATAACAAATTATTTCCCATAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGaa                                                                                                                             <  2:815910/114‑1 (MQ=255)          gggAGTGTCGCTTTTGCTGATAACAAATTATTTCCCATAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGaa                                                                                                                             >  1:815910/1‑114 (MQ=255)           ggAGTGTCGCTTTTGCTGATAACAAATTATTTCCCATATCAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTa                                                                                                                      <  2:775080/120‑1 (MQ=255)           ggAGTGTCGCTTTTGCTGATAACAAATTATTTCCCATAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTa                                                                                                                      >  1:775080/1‑120 (MQ=255)           ggAGTGTCGCTTTTGCTGATAACAAATTATTTCCCATAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTgg                                                                                                   >  1:91724/1‑139 (MQ=255)               tgtCGCTTTTGCTGATAACAAATTATTTCCCATAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTATACAGACTTCTGGAAgg                                                                                               >  2:108055/1‑139 (MQ=255)                gtCGCTTTTGCTGATAACAAATTATTTCCCATAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGt                                                                                              >  1:863707/1‑139 (MQ=255)                    cTTTTGCTGATAACAAATTATTTCCCATAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTgcg                                                                                          <  2:673545/139‑1 (MQ=255)                     ttttGCTGATAACAAATTATTTCCCATAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTgcgc                                                                                         <  2:1007349/139‑1 (MQ=255)                     ttttGCTGATAACAAATTATTTCCCATAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTgcgc                                                                                         >  2:1156786/1‑139 (MQ=255)                                     ttattTCCCATAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACAcgc                                                                         <  1:1156786/139‑1 (MQ=255)                                         ttCCCATAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTa                                                                                                                      <  2:510517/90‑1 (MQ=255)                                         ttCCCATAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTa                                                                                                                      >  1:510517/1‑90 (MQ=255)                                               tAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCtt                                                               >  2:1264793/1‑139 (MQ=255)                                                      tCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGc                                                        >  2:960752/1‑139 (MQ=255)                                                          ttaatataaatatGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACAcg                                                                          <  2:481969/115‑1 (MQ=255)                                                          taaatataaataTGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACAcg                                                                          >  1:481969/1‑117 (MQ=255)                                                             atataaatatGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGaa                                                                                                                             >  2:234961/1‑63 (MQ=255)                                                             atataaatatGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGaa                                                                                                                             <  1:234961/63‑1 (MQ=255)                                                                taaatatGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAgc                                              <  1:1264793/139‑1 (MQ=255)                                                                       ggCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTacaaca                                                                            >  2:1148776/1‑102 (MQ=255)                                                                       ggCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTacaaca                                                                            <  1:1148776/102‑1 (MQ=255)                                                                       ggCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTa                                       >  2:770223/1‑139 (MQ=255)                                                                       ggCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTa                                       >  2:357383/1‑139 (MQ=255)                                                                       ggCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCCTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTa                                       >  2:1157741/1‑139 (MQ=255)                                                                       ggCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTAGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTa                                       >  2:831226/1‑139 (MQ=255)                                                                       ggCAAGCTATATGTTCTGTTATATGAATAAAAATCACCTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTa                                       >  2:877483/1‑139 (MQ=255)                                                                                  tGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCAc                                 >  2:462349/1‑134 (MQ=255)                                                                                  tGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCAc                                 <  1:462349/134‑1 (MQ=255)                                                                                    ttttGTTATATGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCttt                          <  1:770223/139‑1 (MQ=255)                                                                                          tataTGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGCTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACctc                    >  1:694555/1‑139 (MQ=255)                                                                                            tatGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTg                  >  1:995037/1‑139 (MQ=255)                                                                                              tGAATAAAAATCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGc                >  2:737919/1‑139 (MQ=255)                                                                                                   ttcccttaccTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCa                                                          <  1:1085514/84‑1 (MQ=255)                                                                                                   ttcccttaccTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCa                                                          >  2:1085514/9‑92 (MQ=255)                                                                                                   aaaaaTCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACg                                <  1:502143/118‑1 (MQ=255)                                                                                                   aaaaaTCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACg                                >  2:502143/1‑118 (MQ=255)                                                                                                      aaTCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAg        >  1:1145876/1‑139 (MQ=255)                                                                                                       aTCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCg                                                            >  1:848274/1‑86 (MQ=255)                                                                                                       aTCCCCTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCg                                                            <  2:848274/86‑1 (MQ=255)                                                                                                         ccccTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAgcgc     <  1:960752/139‑1 (MQ=255)                                                                                                          cccTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTgc                                                                                           >  2:226049/1‑52 (MQ=255)                                                                                                          cccTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTgc                                                                                           <  1:226049/52‑1 (MQ=255)                                                                                                          cccTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCa    <  1:831226/139‑1 (MQ=255)                                                                                                          cccTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCa    >  1:581214/1‑139 (MQ=255)                                                                                                          cccTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCa    <  2:694555/139‑1 (MQ=255)                                                                                                          cccTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCa    >  1:365801/1‑139 (MQ=255)                                                                                                          cccTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCa    <  1:357383/139‑1 (MQ=255)                                                                                                          cccTCTCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCa    <  2:863707/139‑1 (MQ=255)                                                                                                            ctctCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGt                                                                      <  1:409297/71‑1 (MQ=255)                                                                                                            ctctCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGt                                                                      >  2:409297/1‑71 (MQ=255)                                                                                                            ctctCCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCaga  >  1:250689/1‑139 (MQ=255)                                                                                                            ccctcCGGTAAGGGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCagc  >  2:24382/3‑138 (MQ=255)                                                                                                                        |                                                                                                                               TCATAAACAGGGAGTGTCGCTTTTGCTGATAACAAATTATTTCCCATAACAATTCCTTAAATATAAATATGGCAAGCTATATGTTTTGTTATATGAATAAAAATCCCCTCTCCGGTAAGAGAAGGGATTAAGGGTTTACAGACTTCTGGAAGGTTGCGCAGCTCTTACAACACGCGGTTGATCTTCCGCAGCGTCTTCCAGCGCACTTAAATCACGGTCTTTCACCTCTGGCATTTTCAGCGCAGA  >  NC_000913/1630867‑1631112

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