Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,649,041	G→A	*73* (TAG→TAA)	ydfC →	uncharacterized protein, Qin prophage

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,649,041	0	G	A	100.0%	199.6 / NA	61	*73* (TAG→TAA)	ydfC	uncharacterized protein, Qin prophage
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (31/30);  total (31/30)

CGCAGGATCTGAAATGGGCTAATGAAGCCAGGGCTAAGCAGGCTGAAAATCTGTCCGCTCTGGAGCAGGATATCACCAATCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCA  >  NC_000913/1648908‑1649174                                                                                                                                      |                                                                                                                                      cgcAGGATCTGAAATGGGCTAATGAAGCCAGGGCTAAGCGGGCTGAAAATCTGTCCGCTCTGGAGCAGGATATCACCAATCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCt                                                                                                                                  >  1:680357/1‑139 (MQ=255)              aTGGGCTAATGAAGCCAGGGCTAAGCAGGCTGAAAATCTGTCCGCTCTGGAGCAGGATATCACCAATCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCt                                                                                                                     <  2:1953065/139‑1 (MQ=255)                 ggCTAATGAAGCCAGGGCTAAGCAGGCTGAAAATCTGTCCGCTCTGGAGCAGGATATCACCAATCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATcaa                                                                                                                            >  2:240004/1‑129 (MQ=255)                 ggCTAATGAAGCCAGGGCTAAGCAGGCTGAAAATCTGTCCGCTCTGGAGCAGGATATCACCAATCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATcaa                                                                                                                            <  1:240004/129‑1 (MQ=255)                   cTAATGAAGCCAGGGCTAAGCAGGCTGAAAATCTGTCCGCTCTGGAGCAGGATATCACCAATCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCTAAGGt                                                                                                                >  1:955524/1‑139 (MQ=255)                      aTGAAGCCAGGGCTAAGCAGGCTGAAAATCTGTCCGCTCTGGAGCAGGATATCACCAATCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAAt                                                                                                                               <  1:543850/121‑1 (MQ=255)                      aTGAAGCCAGGGCTAAGCAGGCTGAAAATCTGTCCGCTCTGGAGCAGGATATCACCAATCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAAt                                                                                                                               >  2:543850/1‑121 (MQ=255)                      aTGAAGCCAGGGCTAAGCAGGCTGAAAATCTGTCCGCTCTGGAGCAGGATATCACCAATCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCTAAGGTTAg                                                                                                             >  2:1479590/1‑139 (MQ=255)                       tGAAGCCAGGGCTAAGCAGGCTGAAAATCTGTCCGCTCTGGAGCAGGATATCACCAATCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCTAAGGTTAGt                                                                                                            <  2:390545/139‑1 (MQ=255)                           gCCAGGGCTAAGCAGGCTGAAAATCTGTCCGCTCTGGAGCAGGATATCACCAATCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCTAAGGTTAGTAAtt                                                                                                        >  1:817677/1‑139 (MQ=255)                                ggCTAAGCAGGCTGAAAATCTGTCCGCTCTGGAGCAGGATATCACCAATCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGa                                                                                                   >  2:493204/1‑139 (MQ=255)                                     aGCAGGCTGAAAATCTGTCCGCTCTGGAGCAGGATATCACCAATCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCTAAGGtt                                                                                                               >  1:131020/1‑122 (MQ=255)                                     aGCAGGCTGAAAATCTGTCCGCTCTGGAGCAGGATATCACCAATCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCTAAGGtt                                                                                                               <  2:131020/122‑1 (MQ=255)                                     aGCAGGCTGAAAATCTGTCCGCTCTGGAGCAGGATATCACCAATCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTtct                                                                                              >  2:875817/1‑139 (MQ=255)                                        aGGCTGAAAATCTGTCCGCTCTGGAGCAGGATATCACCAATCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCa                                                                                           <  1:1479590/139‑1 (MQ=255)                                           cTGAAAATCTGTCCGCTCTGGAGCAGGATATCACCAATCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATcaa                                                                                                                            <  2:366293/103‑1 (MQ=255)                                           cTGAAAATCTGTCCGCTCTGGAGCAGGATATCACCAATCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATcaa                                                                                                                            >  1:366293/1‑103 (MQ=255)                                                   cTGTCCGCTCTGGAGCAGGATATCACCAATCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTg                                                                                <  1:493204/139‑1 (MQ=255)                                                     gTCCGCTCTGGAGCAGGATATCACCAATCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGa                                                                              <  2:817677/139‑1 (MQ=255)                                                              ggAGCAGGATATCACCAATCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGtgcgtgcg                                                                     <  2:955524/139‑1 (MQ=255)                                                                    ggATATCACCAATCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGac                                                               <  1:875817/139‑1 (MQ=255)                                                                     gATATCACCAATCTGGAAAAAGCAGCGGCATTAAGTATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAgtgcg                                                                         >  2:1984943/1‑128 (MQ=255)                                                                     gATATCACCAATCTGGAAAAAGCAGCGGCATTAAGTATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAgtgcg                                                                         <  1:1984943/128‑1 (MQ=255)                                                                     gATATCACCAATCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGtgcgtgcg                                                                     >  2:308327/1‑132 (MQ=255)                                                                     gATATCACCAATCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGtgcgtgcg                                                                     <  1:308327/132‑1 (MQ=255)                                                                     gATATCACCAATCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGaca                                                              >  2:336137/1‑139 (MQ=255)                                                                       tatCACCAATCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCTAAGGtt                                                                                                               >  2:1940129/1‑88 (MQ=255)                                                                       tatCACCAATCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCCAAGGtt                                                                                                               <  1:1940129/88‑1 (MQ=255)                                                                               aTCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCAc                                                                                          >  1:1368352/1‑101 (MQ=255)                                                                               aTCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCAc                                                                                          <  2:1368352/101‑1 (MQ=255)                                                                               aTCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAgtgcg                                                                         >  1:458338/1‑118 (MQ=255)                                                                               aTCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAgtgcg                                                                         <  2:458338/118‑1 (MQ=255)                                                                                         agcagcGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTg                                          >  1:1702746/1‑139 (MQ=255)                                                                                          gcagcGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGa                                         >  1:1155684/1‑139 (MQ=255)                                                                                                 cATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTa                                  <  2:680357/139‑1 (MQ=255)                                                                                                  aTTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTAc                                 <  1:336137/139‑1 (MQ=255)                                                                                                   ttAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCTAAGGTTAGTACTTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGa                                         <  1:1400861/130‑1 (MQ=255)                                                                                                   ttAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAg                                                                                                  >  2:341930/1‑73 (MQ=255)                                                                                                   ttAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAg                                                                                                  <  1:341930/73‑1 (MQ=255)                                                                                                   ttAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGa                                         >  2:1400861/1‑130 (MQ=255)                                                                                                     aaGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGa                              <  1:1313247/139‑1 (MQ=255)                                                                                                     aaGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGa                              <  2:1702746/139‑1 (MQ=255)                                                                                                                 aaTGATTACATACCCGCGTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACgg                                                                                        <  2:1573618/69‑1 (MQ=255)                                                                                                                 aaTGATTACATACCCGCGTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACgg                                                                                        >  1:1573618/1‑69 (MQ=255)                                                                                                                     aTTACATACCCGCGTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGac                                                               >  2:336956/1‑90 (MQ=255)                                                                                                                     aTTACATACCCGCGTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGac                                                               <  1:336956/90‑1 (MQ=255)                                                                                                                           tACCCGCGTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCAc                                                                                          >  1:1051290/1‑57 (MQ=255)                                                                                                                           tACCCGCGTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCAc                                                                                          <  2:1051290/57‑1 (MQ=255)                                                                                                                           tACCCGCGTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTc        >  2:2051304/1‑139 (MQ=255)                                                                                                                           tACCCGCGTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTc        >  2:1244483/1‑139 (MQ=255)                                                                                                                             cccGCGTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTAc                                 <  1:916534/112‑1 (MQ=255)                                                                                                                             cccGCGTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTAc                                 >  2:916534/1‑112 (MQ=255)                                                                                                                             cccGCGTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAg      >  1:246620/1‑139 (MQ=255)                                                                                                                              ccGCGTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACgg                                                                                        >  1:583103/1‑56 (MQ=255)                                                                                                                              ccGCGTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACgg                                                                                        <  2:583103/56‑1 (MQ=255)                                                                                                                               cgcgTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAg                                                                                                  >  2:1226348/1‑45 (MQ=255)                                                                                                                               cgcgTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAg                                                                                                  <  1:1226348/45‑1 (MQ=255)                                                                                                                               cgcgTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGAc    <  2:1155684/139‑1 (MQ=255)                                                                                                                                gcgTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGtgcgtgcg                                                                     >  1:1290195/1‑73 (MQ=255)                                                                                                                                gcgTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGtgcgtgcg                                                                     <  2:1290195/73‑1 (MQ=255)                                                                                                                                 cgTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCa  >  1:1277195/1‑139 (MQ=255)                                                                                                                                      |                                                                                                                                      CGCAGGATCTGAAATGGGCTAATGAAGCCAGGGCTAAGCAGGCTGAAAATCTGTCCGCTCTGGAGCAGGATATCACCAATCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCA  >  NC_000913/1648908‑1649174

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