Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,515,726	T→C	Y63C (TAT→TGT)	dam ←	DNA adenine methyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,515,726	0	T	C	100.0%	123.1 / NA	35	Y63C (TAT→TGT)	dam	DNA adenine methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (15/20);  total (15/20)

CCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCA  >  NC_000913/3515592‑3515857                                                                                                                                       |                                                                                                                                    cccGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACagag                                                                                                                                 >  1:576446/1‑139 (MQ=255)            ggaTCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGt                                                                                                                      <  2:342116/139‑1 (MQ=255)            ggaTCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGt                                                                                                                      <  2:46347/139‑1 (MQ=255)                  tGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTAt                                                                                                                <  2:28278/139‑1 (MQ=255)                   ggCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTAGTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCa                                                                                                                         >  2:600780/1‑129 (MQ=255)                   ggCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTAGTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCa                                                                                                                         <  1:600780/129‑1 (MQ=255)                     cTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGa                                                                                                             <  1:486908/139‑1 (MQ=255)                      ttttGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGat                                                                                                            <  1:63607/139‑1 (MQ=255)                        ttgttgAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGatat                                                                                                          >  1:619166/1‑139 (MQ=255)                         tgttgAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATc                                                                                                         <  1:176600/139‑1 (MQ=255)                           ttgAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCgg                                                                                                       <  2:619166/139‑1 (MQ=255)                                             tGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGc                                                                                                      >  1:378342/1‑122 (MQ=255)                                             tGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGc                                                                                                      <  2:378342/122‑1 (MQ=255)                                                tAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGATACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGt                                                                                  <  2:464113/139‑1 (MQ=255)                                                tAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGt                                                                                  <  2:121528/139‑1 (MQ=255)                                                tAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGt                                                                                  <  1:312890/139‑1 (MQ=255)                                                          cGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGt                                                                                                                      <  1:329763/93‑1 (MQ=255)                                                          cGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGt                                                                                                                      >  2:329763/1‑93 (MQ=255)                                                          cGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTgag                                                                        >  2:391101/1‑139 (MQ=255)                                                                  tttGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGat                                                                                                            >  1:547163/1‑95 (MQ=255)                                                                  tttGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGat                                                                                                            <  2:547163/95‑1 (MQ=255)                                                                                          cgTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCaa                                        >  2:706631/1‑139 (MQ=255)                                                                                            tGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAAcc                                      >  2:240768/1‑139 (MQ=255)                                                                                                               tcaGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGc                                                                                                      >  2:275315/1‑56 (MQ=255)                                                                                                               tcaGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGc                                                                                                      <  1:275315/56‑1 (MQ=255)                                                                                                               tcaGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCAcc                                                      >  1:226426/1‑104 (MQ=255)                                                                                                               tcaGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCAcc                                                      <  2:226426/104‑1 (MQ=255)                                                                                                               tcaGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCa                   >  1:697291/1‑139 (MQ=255)                                                                                                               tcaGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCa                   >  2:572060/1‑139 (MQ=255)                                                                                                                 aGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCgg                                                          >  1:725869/1‑98 (MQ=255)                                                                                                                 aGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCgg                                                          <  2:725869/98‑1 (MQ=255)                                                                                                                   tACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATg               <  1:572060/139‑1 (MQ=255)                                                                                                                           ttCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTaa       <  2:524406/139‑1 (MQ=255)                                                                                                                             cacaATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAata     >  1:594214/1‑139 (MQ=255)                                                                                                                                aaTGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATatca  <  1:683213/139‑1 (MQ=255)                                                                                                                                       |                                                                                                                                    CCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCA  >  NC_000913/3515592‑3515857

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