Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,330,454	G→A	A72V (GCC→GTC)	yciN ←	DUF2498 protein YciN

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,330,454	0	G	A	100.0%	111.8 / NA	35	A72V (GCC→GTC)	yciN	DUF2498 protein YciN
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (15/20);  total (15/20)

GCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTAAATACCAGCACGCCATTACGCTGCGTGACGCCGGTGGCTTCAATCCCTGCG  >  NC_000913/1330328‑1330585                                                                                                                               |                                                                                                                                    gcCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGACCAGGTGTTTGa                                                                                                                         <  1:1078163/139‑1 (MQ=255) gcCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGACCAGGTGTTTGa                                                                                                                         <  1:232518/139‑1 (MQ=255)  ccGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGACCAGGTGTTTGaa                                                                                                                        <  2:729805/139‑1 (MQ=255)   cgatcgatTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGACCAGGTGTTTGAAc                                                                                                                       <  2:934448/139‑1 (MQ=255)             ttgCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTATGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGACCAGGTGTTTGAACATATTAAaca                                                                                                             <  1:745223/139‑1 (MQ=255)                        gtggCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGACCAGGTGTTTGAACATATTAAACACCGCGGTGCtt                                                                                                  >  1:954801/1‑139 (MQ=255)                                  ggTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGACCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTg                                                                                        >  2:321429/1‑139 (MQ=255)                                              ttGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGACCa                                                                                                                                  <  2:334172/85‑1 (MQ=255)                                              ttGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGACCa                                                                                                                                  >  1:334172/1‑85 (MQ=255)                                              ttGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGACCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGtt                                                                            >  2:68359/1‑139 (MQ=255)                                              ttGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGACCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGtt                                                                            >  2:100571/1‑139 (MQ=255)                                               tGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGACCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTc                                                                           <  2:954801/139‑1 (MQ=255)                                               tGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGACCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTc                                                                           <  2:843464/139‑1 (MQ=255)                                                                gcgAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGACCAGGTGTTTGACCATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGc                                                          <  1:321429/139‑1 (MQ=255)                                                                gcgAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGACCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGc                                                          <  2:141786/139‑1 (MQ=255)                                                                         aGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGACCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCg                                                           >  1:175146/1‑129 (MQ=255)                                                                         aGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGACCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCg                                                           <  2:175146/129‑1 (MQ=255)                                                                             cTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGACCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTAAATACCAg                                             <  2:1130973/139‑1 (MQ=255)                                                                               cgcgTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGACCAGGTGTTTGa                                                                                                                         <  1:920051/61‑1 (MQ=255)                                                                               cgcgTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGACCAGGTGTTTGa                                                                                                                         >  2:920051/1‑61 (MQ=255)                                                                               cgcgTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGACCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGc                                                                                             >  1:578370/1‑89 (MQ=255)                                                                               cgcgTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGACCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGc                                                                                             <  2:578370/89‑1 (MQ=255)                                                                               cgcgTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGACCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGc                                                                                      <  1:540579/96‑1 (MQ=255)                                                                               cgcgTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGACCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGc                                                                                      >  2:540579/1‑96 (MQ=255)                                                                               cgcgTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGACCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTAAATACCAGCa                                           >  2:1127070/1‑139 (MQ=255)                                                                               cgcgTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGACCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTAAATACCAGCa                                           >  1:103675/1‑139 (MQ=255)                                                                                    tttGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGACCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTAAATACCAGCACGCCa                                      <  2:84342/139‑1 (MQ=255)                                                                                      tGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGACCAGGTGTTTGAACata                                                                                                                    >  2:867418/1‑59 (MQ=255)                                                                                      tGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGACCAGGTGTTTGAACata                                                                                                                    <  1:867418/59‑1 (MQ=255)                                                                                                     gatgatATTTTTCTGAAAGCACATGGACCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTAAATACCAGCACGCCATTACGCTGCGTGACGcc                     >  1:230009/1‑139 (MQ=255)                                                                                                          tatTTTTCTGAAAGCACATGGACCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTAAATACCAGCACGCCATTACGCTGCGTGACGCCGGTgg                >  1:640615/1‑139 (MQ=255)                                                                                                                    aaaGCACATGGACCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAagta                                                               <  2:1133964/82‑1 (MQ=255)                                                                                                                    aaaGCACATGGACCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAagta                                                               >  1:1133964/1‑82 (MQ=255)                                                                                                                      aGCACATGGACCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTAAATACCAGCACGCCATTACGCTGCGTGACGCCGGTGGCTTCAATCCCTg    <  1:1127070/139‑1 (MQ=255)                                                                                                                        cacaTGGACCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTAAATACCAGCACGCCATTACGCTGCGTGACGCCGGTGGCTTCAATCCCTGCg  <  2:230009/139‑1 (MQ=255)                                                                                                                               |                                                                                                                                    GCCGATCGATTGTTGCTACGCTGGTGGCAGCGGGGTCAAAAGCCATTGATGTAAAAGACAAACGCGAGGCTAAGACCTCGCGTTTTGCTTTAATCAACCAGATGATATTTTTCTGAAAGCACATGGGCCAGGTGTTTGAACATATTAAACACCGCGGTGCTTTTGGCTGTTGGCAATCCTTGTTCATCTAAAAAGTAGTCGCCGGTAAATACCAGCACGCCATTACGCTGCGTGACGCCGGTGGCTTCAATCCCTGCG  >  NC_000913/1330328‑1330585

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: ‑