Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,614,260	G→A	*260* (TAG→TAA)	yjjV →	putative DNase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,614,260	0	G	A	100.0%	123.2 / NA	39	*260* (TAG→TAA)	yjjV	putative DNase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (21/18);  total (21/18)

GGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGT  >  NC_000913/4614132‑4614394                                                                                                                                 |                                                                                                                                       gggCAGCCTAACTGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAg                                                                                                                              >  1:180003/1‑139 (MQ=255) gggCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAg                                                                                                                              <  2:1112251/139‑1 (MQ=255)                     gCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGc                                                                                                          >  1:246700/1‑139 (MQ=255)                         gCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCgg                                                                                                                                   >  2:456415/1‑110 (MQ=255)                         gCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCgg                                                                                                                                   <  1:456415/110‑1 (MQ=255)                         gCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGtt                                                                                                      >  2:500867/1‑139 (MQ=255)                             cccGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGgcg                                                                                                  <  2:180003/139‑1 (MQ=255)                                                      gTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAAc                                                                         <  2:127762/139‑1 (MQ=255)                                                                            gagaTTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGt                                                   <  1:1303320/139‑1 (MQ=255)                                                                                   cgcAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGAt                                                                                                                                 <  1:1111061/54‑1 (MQ=255)                                                                                   cgcAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGAt                                                                                                                                 >  2:1111061/1‑54 (MQ=255)                                                                                   cgcAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCGGGCAGTTGGCGCACAATGCCTGATGCg                                                                                  >  2:602572/1‑101 (MQ=255)                                                                                   cgcAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCg                                                                                  <  1:519094/101‑1 (MQ=255)                                                                                   cgcAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCg                                                                                  <  1:602572/101‑1 (MQ=255)                                                                                   cgcAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCg                                                                                  >  2:519094/1‑101 (MQ=255)                                                                                    gcAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCg                                           <  1:670326/139‑1 (MQ=255)                                                                                    gcAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCg                                           <  2:227238/139‑1 (MQ=255)                                                                                         cggtGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCg                                      <  2:872240/136‑1 (MQ=255)                                                                                                     acacGTATACGTTGTTTAATGTACCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAAc                                                                         >  2:922002/1‑92 (MQ=17)                                                                                                     acacGTATACGTTGTTTAACGTACCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAAc                                                                         <  1:922002/92‑1 (MQ=255)                                                                                                      cacGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAg                         >  1:55769/1‑139 (MQ=255)                                                                                                        cGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTc                                                                                                                        <  1:685488/42‑1 (MQ=255)                                                                                                        cGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTc                                                                                                                        >  2:685488/1‑42 (MQ=255)                                                                                                        cGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATc                                                               >  2:1254378/1‑99 (MQ=255)                                                                                                        cGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATc                                                               <  1:1254378/99‑1 (MQ=255)                                                                                                        cGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCg                                           <  1:277599/119‑1 (MQ=255)                                                                                                        cGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCg                                           >  2:277599/1‑119 (MQ=255)                                                                                                        cGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGc                       >  1:85761/1‑139 (MQ=255)                                                                                                        cGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGc                       >  1:1272106/1‑139 (MQ=255)                                                                                                           ataCGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGtt                    >  1:163670/1‑139 (MQ=255)                                                                                                             aCGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGTACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCa                  >  2:1195373/1‑139 (MQ=255)                                                                                                               gttgttTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCAcg                <  2:1272106/139‑1 (MQ=255)                                                                                                                  gttTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGc                                                                                   <  2:1303123/69‑1 (MQ=255)                                                                                                                  gttTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGc                                                                                   >  1:1303123/1‑69 (MQ=255)                                                                                                                           gCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGcc                                            >  1:1000712/1‑99 (MQ=255)                                                                                                                           gCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGcc                                            <  2:1000712/99‑1 (MQ=255)                                                                                                                           gCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCa    >  2:1113372/1‑139 (MQ=255)                                                                                                                             cGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGt  >  2:341131/1‑139 (MQ=255)                                                                                                                             cGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGt  >  2:603184/1‑139 (MQ=255)                                                                                                                             cGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGAGCACAATGCATGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGt  >  2:809342/1‑139 (MQ=255)                                                                                                                                 |                                                                                                                                       GGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGT  >  NC_000913/4614132‑4614394

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