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%	71.0 / NA	23	*73* (TAG→TAA)	ydfC	uncharacterized protein, Qin prophage
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (11/12);  total (11/12)

AGCGCAGGATCTGAAATGGGCTAATGAAGCCAGGGCTAAGCAGGCTGAAAATCTGTCCGCTCTGGAGCAGGATATCACCAATCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATA  >  NC_000913/1648906‑1649153                                                                                                                                        |                                                                                                                 aGCGCAGGATCTGAAATGGGCTAATGAAGCCAGGGCTAAGCAGGCTGAAAATCTGTCCGCTCTGGAGCAGGATATCACCAATCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTtaatag                                                                                                               <  2:430652/139‑1 (MQ=255)                aTGGGCTAATGAAGCCAGGGCTAAGCAGGCTGAAAATCTGTCCGCTCTGGAGCAGGATATCACCAATCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCt                                                                                                <  1:1274986/139‑1 (MQ=255)                    gCTAATGAAGCCAGGGCTAAGCAGGCTGAAAATCTGTCCGCTCTGGAGCAGGATATCACCAATCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCt                                                                                                >  1:447844/1‑135 (MQ=255)                    gCTAATGAAGCCAGGGCTAAGCAGGCTGAAAATCTGTCCGCTCTGGAGCAGGATATCACCAATCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCt                                                                                                <  2:447844/135‑1 (MQ=255)                    gCTAATGAAGCCAGGGCTAAGCAGGCTGAAAATCTGTCCGCTCTGGAGCAGGATATCACCAATCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCTAAgg                                                                                            >  2:504114/1‑139 (MQ=255)                      tAATGAAGCCAGGGCTAAGCAGGCTGAAAATCTGTCCGCTCTGGAGCAGGATATCACCAATCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCTAAGGtt                                                                                          <  2:731464/139‑1 (MQ=255)                        aTGAAGCCAGGGCTAAGCAGGCTGAAAATCTGTCCGCTCTGGAGCAGGATATCACCAATCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATcaa                                                                                                       >  2:803863/1‑124 (MQ=255)                        aTGAAGCCAGGGCTAAGCAGGCTGAAAATCTGTCCGCTCTGGAGCAGGATATCACCAATCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATcaa                                                                                                       <  1:803863/124‑1 (MQ=255)                        aTGAAGCCAGGGCTAAGCAGGCTGAAAATCTGTCCGCTCTGGAGCAGGATATCACCAATCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCTa                                                                                               >  1:1066599/1‑132 (MQ=255)                        aTGAAGCCAGGGCTAAGCAGGCTGAAAATCTGTCCGCTCTGGAGCAGGATATCACCAATCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCTa                                                                                               <  2:1066599/132‑1 (MQ=255)                              ccAGGGCTAAGCAGGCTGAAAATCTGTCCGCTCTGGAGCAGGATATCACCAATCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTa                                                                                  >  1:1024179/1‑139 (MQ=255)                                   gCTAAGCAGGCTGAAAATCTGTCCGCTCTGGAGCAGGATATCACCAATCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAg                                                                             <  1:504114/139‑1 (MQ=255)                                       aGCAGGCTGAAAATCTGTCCGCTCTGGAGCAGGATATCACCAATCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTtct                                                                         >  2:877149/1‑139 (MQ=255)                                        gCAGGCTGAAAATCTGTCCGCTCTGGAGCAGGATATCACCAATCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTtctc                                                                        >  2:398938/1‑139 (MQ=255)                                            gCTGAAAATCTGTCCGCTCTGGAGCAGGATATCACCAATCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCTAAGGtt                                                                                          >  2:1066166/1‑117 (MQ=255)                                            gCTGAAAATCTGTCCGCTCTGGAGCAGGATATCACCAATCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCTAAGGtt                                                                                          <  1:1066166/117‑1 (MQ=255)                                            gCTGAAAATCTGTCCGCTCTGGAGCAGGATATCACCAATCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACg                                                                    >  1:100947/1‑139 (MQ=255)                                             cTGAAAATCTGTCCGCTCTGGAGCAGGATATCACCAATCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACgg                                                                   <  1:877149/139‑1 (MQ=255)                                                   aTCTGTCCGCTCTGGAGCAGGATATCACCAATCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAgg                                                             <  1:398938/139‑1 (MQ=255)                                                                                                cGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTcc                                                                       >  1:868696/1‑84 (MQ=255)                                                                                                cGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTcc                                                                       <  2:868696/84‑1 (MQ=255)                                                                                                              gCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGata  <  1:399311/139‑1 (MQ=255)                                                                                                                             tACCCGCGTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTa             <  1:673481/113‑1 (MQ=255)                                                                                                                             tACCCGCGTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTa             >  2:673481/1‑113 (MQ=255)                                                                                                                                        |                                                                                                                 AGCGCAGGATCTGAAATGGGCTAATGAAGCCAGGGCTAAGCAGGCTGAAAATCTGTCCGCTCTGGAGCAGGATATCACCAATCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATA  >  NC_000913/1648906‑1649153

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