Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,582,645	C→T	D261D (GAC→GAT)	yhhZ →	putative endonuclease YhhZ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,582,645	0	C	T	100.0%	98.8 / NA	27	D261D (GAC→GAT)	yhhZ	putative endonuclease YhhZ
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (13/14);  total (13/14)

AAAAGGCTATGATCAAATACTAAATGAAAAAACAGCGGGCATAGCGATAGCAACAGCATCTATTCTTTTAACAAAGCGTTCTAATGTTGATACATATACAGAAATAAATAGTTACTTAGGCAAACTTAGAGGTCAACAAAAACTTCTTGATGGTATAGACATAATAGAAATAATATACATTAAGAGACCTTCAAAAGACTTAGCTAACTTACGAAAGGAGTTTAATAAAACTGTAAGAAAAAATTTTCTTATCAAACTTGCAAAAACCTCCGAAGCATCTGGAAGATTCAACGCCGAAGACCTTTTAAGAATGAGAAAGGGCAATGTTCCTCTAAATTATAATGTTCACCATAAACTATCTCTAGATGATGGTGGTACTAATGATTTCGAAAATT  >  NC_000913/3582486‑3582880                                                                                                                                                                |                                                                                                                                                                                                                                            aaaaGGCTATGATCAAATACTAAATGAAAAAACAGCGGGCATAGCGATAGCAACAGCATCTATTCTTTTAACAAAGCGTTCTAATGTTGATACATATACAGAAATAAATAGTTACTTAGGCAAACTTAGAGGTCAACAAAAACTTCTTGATGGTATAGATATAATAGAAATAATATACAt                                                                                                                                                                                                                         >  2:509200/1‑180 (MQ=255) aaaaGGCTATGATCAAATACTAAATGAAAAAACAGCGGGCATAGCGATAGCAACAGCATCTATTCTTTTAACAAAGCGTTCTAATGTTGATACATATACAGAAATAAATAGTTACTTAGGCAAACTTAGAGGTCAACAAAAACTTCTTGATGGTATAGATATAATAGAAATAATATACAt                                                                                                                                                                                                                         <  1:509200/180‑1 (MQ=255) aaaaGGCTATGATCAAATACTAAATGAAAAAACAGCGGGCATAGCGATAGCAACAGCATCTATTCTTTTAACAAAGCGTTCTAATGTTGATACATATACAGAAATAAATAGTTACTTAGGCAAACTTAGAGGTCAACAAAAACTTCTTGATGGTATAGATATAATAGAAATAATATACATTAAGAGACCTTCAAAAGACTTAGCTAACTTACGAAAGGAGTTTAATAAAACTGTAATAAAAAATTTTCtt                                                                                                                                                   <  2:518731/250‑1 (MQ=255)          tGATCAAATACTAAATGAAAAAACAGCGGGCATAGCGATAGCAACAGCATCTATTCTTTTAACAAAGCGTTCTAATGTTGATACATATACAGAAATAAATAGTTACTTAGGCAAACTTAGAGGTCAACAAAAACTTCTTGATGGTATAGATATAATAGAAATAATATACATTAAGAGACCTTCAAAAGACTTAGCTAACTTACGAAAGGAGTTTAATAAAACTGTAAGAAAAAATTTTCTTATCAAACTTg                                                                                                                                         >  1:348313/1‑251 (MQ=255)               aaaTACTAAATGAAAAAACAGCGGGCATAGCGATAGCAACAGCATCTATTCTTTTAACAAAGCGTTCTAATGTTGATACATATACAGAAATAAATAGTTACTTAGGCAAACTTAGAGGTCAACAAAAACTTCTTGATGGTATAGATATAATAGAAATAATATACATTAAGAGACCTTCAAAAGACTTAGCTAACTTACGAAAGGAGTTTAATAAAACTGTAAGAAAAAATTTTCTTATCAAACTTGCaaa                                                                                                                                     >  1:407238/1‑250 (MQ=255)                           aaaaaaCAGCGGGCATAGTGATAGCAACAGCATCTATTCTTTTAACAAAGCGTTCTAATGTTGATACATATACAGAAATAAATAGTTACTTAGGCAAACTTAGAGGTCAACAAAAACTTCTTGATGGTATAGATATAATAGAAATAATATACATTAAGAGACCTTCAAAAGACTTAGCTAACTTACGAAAGGAGTTTAATAAAACTGTAAGAAAAAATTTTCTTATCAAACTTGCAAAAACCTCCGAAGCa                                                                                                                        <  1:26830/251‑1 (MQ=255)                                aCAGCGGGCATAGCGATAGCAACAGCATCTATTCTTTTAACAAAGCGTTCTAATGTTGATACATATACAGAAATAAATAGTTACTTAGGCAAACTTAGAGGTCAACAAAAACTTCTTGATGGTATAGATATAATAGAAATAATATACATTAAGAGACCTTCAAAAGACTTAGCTAACTTACGAAAGGAGTTTAATAAAACTGTAAGAAAAAATTTTCTTATCAAACTTGCAAAAACCTCCGAAGCATCTgg                                                                                                                   <  2:205315/251‑1 (MQ=255)                                         aTAGCGATAGCAACAGCATCTATTCTTTTAACAAAGCGTTCTAATGTTGATACATATACAGAAATAAATAGTTACTTAGGCAAACTTAGAGGTCAACAAAAACTTCTTGATGGTATAGATATAATAGAAATAATATACATTAAGAGACCTTCAAAAGACTTAGCTAACTTACGAAAGGAGTTTAATAAAACTGTAAGAAAAAATTTTCTTATCAAACTTGCAAAAACCTCCGAAGCATCTGGAAGATTCaa                                                                                                          <  2:407238/251‑1 (MQ=255)                                           aGCGATAGCAACAGCATCTATTCTTTTAACAAAGCGTTCTAATGTTGATACATATACAGAAATAAATAGTTACTTAGGCAAACTTAGAGGTCAACAAAAACTTCTTGATGGTATAGATATAATAGAAATAATATACATTAAGAGACCTTCAAAAGACTTAGCTAACTTACGAAAGGAGTTTAATAAAACTGTAAGAAAAAATTTTCTTATCAAACTTGCAAAAACCTCCGAAGCATCTGGAAGATTCAAc                                                                                                         <  2:113857/250‑1 (MQ=255)                                               aTAGCAACAGCATCTATTCTTTTAACCAAGCGTTCTAATGTTGATACATATACAGAAATAAATAGTTACTTAGGCAAACTTAGAGGTCAACAAAAACTTCTTGATGGTATAGATATAATAGAAATAATATACATTAAGAGACCTTCAAAAGACTTAGCTAACTTACGAAAGGAGTTTAATAAAACTGTAAGAAAAAATTTTCTTATCAAACTTGCAAAAACCTCCGAAGCATCTGGAAGTTTCAACGCCg                                                                                                     >  2:349179/1‑250 (MQ=255)                                                    aaCAGCATCTATTCTTTTAACAAAGCGTTCTAATGTTGATACATATACAGAAATAAATAGTTACTTAGGCAAACTTAGAGGTCAACAAAAACTTCTTGATGGTATAGATATAATAGAAATAATATACATTAAGAGACCTTCaaa                                                                                                                                                                                                          >  2:436009/1‑144 (MQ=255)                                                    aaCAGCATCTATTCTTTTAACAAAGCGTTCTAATGTTGATACATATACAGAAATAAATAGTTACTTAGGCAAACTTAGAGGTCAACAAAAACTTCTTGATGGTATAGATATAATAGAAATAATATACATTAAGAGACCTTCaaa                                                                                                                                                                                                          <  1:436009/144‑1 (MQ=255)                                                            cTATTCTTTTAACAAAGCGTTCTAATGTTGATACATATACAGAAATAAATAGTTACTTAGGCAAACTTAGAGGTCAACAAAAACTTCTTGATGGTATAGATATAATAGAaataat                                                                                                                                                                                                                               <  1:207419/115‑1 (MQ=255)                                                            cTATTCTTTTAACAAAGCGTTCTAATGTTGATACATATACAGAAATAAATAGTTACTTAGGCAAACTTAGAGGTCAACAAAAACTTCTTGATGGTATAGATATAATAGAaataat                                                                                                                                                                                                                               >  2:207419/1‑115 (MQ=255)                                                                              gTTCTAATGTTGATACATATACAGAAATAAATAGTTACTTAGGCAAACTTAGAGGTCAACAAAAACTTCTTGATGGTATAGATATAATAGAAATAATATACATTAAGAGACCTTCAAAAGACTTAGCTAACTTACGAAAGGAGTTTAATAAAACTGTAAGAAAAAATTTTCTTATCAAACTTGCAAAAACCTCCGAAGCATCTGGAAGATTCAACGCCGAAGACCTTTTAAGAATGAGAAAGGGCAATGtt                                                                     >  1:498666/1‑251 (MQ=255)                                                                               ttCTAATGTTGATACATATACAGAAATAAATAGTTACTTAGGCAAACTTAGAGGTCAACAAAAACTTCTTGATGGTATAGATATAATAGAAATAATATACATTAAGAGACCTTCAAAAGACTTAGCTAACTTACGAAAGGAGTTTAATAAAACTGTAAGAAAAAATTTTCTTATCAAACTTGCAAAAACCTCCGAAGCATCTGGAAGATTCAACGCCGAAGACCTTTTAAGAATGAGAAAGGGCAATGTTc                                                                    >  1:163953/1‑251 (MQ=255)                                                                                          aTACATATACAGAAATAAATAGTTACTTAGGCAAACTTAGAGGTCAACAAAAACTTCTTGATGGTATAGATATAATAGAAATAATATACATTAAGAGACCTTCAAAAGACTTAGCTAACTTACGAAAGGAGTTTAATAAAACTGTAAGAAAAAATTTTCTTATCAAACTTGCAAAAACCTCCGAAGCATCTGGAAGATTCAACGCCGAAGACCTTTTAAGAATGAGAAAGGGCAATGTTCCTCTAAATtat                                                         <  1:62406/251‑1 (MQ=255)                                                                                                   cAGAAATAAATAGTTACTTAGGCAAACTTAGAGGTCAACAAAAACTTCTTGATGGTATAGATATAATAGAAATAATATACATTAAGAGACCTTCAAAAGACTTAGCTAACTTACGAAAGGAGTTTAATAAAACTGTAAGAAAAAATTTTCTTATCAAACTTGCAAAAACCTCCGAAGCATCTGGAAGATTCAACGCCGAAGACCTTTTAAGAATGAGAAAGGGCAATGTTCCTCTAAATTATAATGTTCAc                                                >  1:332275/1‑251 (MQ=255)                                                                                                     gAAATAAATAGTTACTTAGGCAAACTTAGAGGTCAACAAAAACTTCTTGATGGTATAGATATAATAGAAATAATATACATTAAGAGACCTTCAAAAGACTTAGCTAACTTACGAAAGGAGTTTAATAAAACTGTAAGAAAAAATTTTCTTATCAAACTTGCAAAAACCTCCGAAGCATCTGGAAGATTCAACGCCGAAGACCTTTTAAGAATGAGAAAGGGCAATGTTCCTCTAAATTATAATGTTCACCa                                              <  1:161015/251‑1 (MQ=255)                                                                                                                           aaCTTAGAGGTCAACAAAAACTTCTTGATGGTATAGATATAATAGAAATAATATACATTAAGAGACCTTCAAAAGACTTAGCTAACTTACGAAAGGAGTTTAATAAAACTGTAAGAAAAAATTTTCTTATCAAACTTGCAAAAACCTCCGAAGCATCTGGAAGATTCAACGCCGAAGACCTTTTAAGAATGAGAAAGGGCAATGTTCCTCTAAATTATAATGTTCACCATAAACTATCTCTAGATGAtggt                        >  1:63682/1‑251 (MQ=255)                                                                                                                           aaCTTAGAGGTCAACAAAAACTTCTTGATGGTATAGATATAATAGAAATAATATACATTAAGAGACCTTAAAAAGACTTAGCTAACTTACGAAAGGAGTTTAATAAAACTGTAAGAAAAAATTTTCTTATCAAACTTGCAAAAACCTCCGAAGCATCTGGAAGATTCAACGCCGAAGACCTTTTAAGAATGAGAAAGGGCAATGTTCCTCTAAATTATAATGTTCACCATAAACTATCTCTAGATGAtggt                        >  2:151811/1‑251 (MQ=255)                                                                                                                                   ggTCAACAAAAACTTCTTGATGGTATAGATATAATAGAAATAATATACATTAAGAGACCTTCAAAAGACTTAGCTAACTTACGAAAGGAGTTTAATAAAACTGTAAGAAAAAATTTTCTTATCAAACTTGCAAAAACCTCCGAAGCATCTGGAAGATTCAACGCCGAAGACCTTTTAAGAATGAGAAAGGGCAATGTTCCTCTAAATTATAATGTTCACCATAAACTATCTCTAGATGATGGTGGTACTaa                <  1:11428/251‑1 (MQ=255)                                                                                                                                    gTCAACAAAAACTTCTTGATGGTATAGATATAATAGAAATAATATACATTAAGAGACCTTCAAAAGACTTAGCTAACTTACGAAAGGAGTTTAATAAAACTGTAAGAAAAAATTTTCTTATCAAACTTGCAAAAACCTCCGAAGCATCTGGAAGATTCAACGCCGAAGACCTTTTAAGAATGAGAAAGGGCAATGTTCCTCTAAATTATAATGTTCACCATAAACTATCTCTAGATGATGGTGGTACTAAt               <  1:445718/251‑1 (MQ=255)                                                                                                                                      caacaaAAACTTCTTGATGGTATAGATATAATAGAAATAATATACATTAAGAGACCTTCAAAAGACTTAGCTAACTTACGAAAGGAGTTTAATAAAACTGTAAGAAAAAATTTTCTTATCAAACTTGCAAAAACCTCCGAAGCATCTGGAAGATTCAACGCCGAAGACCTTTTAAGAATGAGAAAGGGCAATGTTCCTCTAAATTATAATGTTCACCATAAACTATCTCTAGATGATGGTGGTACTAATg              >  1:463555/1‑250 (MQ=255)                                                                                                                                                    tGATGGTATAGATATAATAGAAATAATATACATTAAGAGACCTTCAAAAGACTTAGCTAACTTACGAAAGGAGTTTAATAAAACTGTAAGAAAAAATTTTCTTATCAAACTTGCAAAAACCTCCGAAGCATCTGGAAGATTCAACGCCGAAGACCTTTTAAGAATGAGAAAGGGCAATGTTCCTCTAAATTATAATGTTCACCATAAACTATCTCTAGATGATGGTGGTACTAATGATTTCGAAAAtt  <  2:6483/248‑1 (MQ=255)                                                                                                                                                     gATGGTATAGATATAATAGAAATAATATACATTAAGAGACCTTCAAAAGACTTAGCTAACTTACGAAAGGAGTTTAATAAAACTGTAAGAAAAAATTTTCTTAt                                                                                                                                                 <  2:153848/104‑1 (MQ=255)                                                                                                                                                     gATGGTATAGATATAATAGAAATAATATACATTAAGAGACCTTCAAAAGACTTAGCTAACTTACGAAAGGAGTTTAATAAAACTGTAAGAAAAAATTTTCTTAt                                                                                                                                                 >  1:153848/1‑104 (MQ=255)                                                                                                                                                                |                                                                                                                                                                                                                                            AAAAGGCTATGATCAAATACTAAATGAAAAAACAGCGGGCATAGCGATAGCAACAGCATCTATTCTTTTAACAAAGCGTTCTAATGTTGATACATATACAGAAATAAATAGTTACTTAGGCAAACTTAGAGGTCAACAAAAACTTCTTGATGGTATAGACATAATAGAAATAATATACATTAAGAGACCTTCAAAAGACTTAGCTAACTTACGAAAGGAGTTTAATAAAACTGTAAGAAAAAATTTTCTTATCAAACTTGCAAAAACCTCCGAAGCATCTGGAAGATTCAACGCCGAAGACCTTTTAAGAATGAGAAAGGGCAATGTTCCTCTAAATTATAATGTTCACCATAAACTATCTCTAGATGATGGTGGTACTAATGATTTCGAAAATT  >  NC_000913/3582486‑3582880

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