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%	208.2 / NA	57	Y63C (TAT→TGT)	dam	DNA adenine methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (26/31);  total (26/31)

CGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATC  >  NC_000913/3515594‑3515859                                                                                                                                     |                                                                                                                                      cGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGAc                                                                                                                                 <  1:1083389/139‑1 (MQ=255)         cggaTCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAgg                                                                                                                         <  2:562205/139‑1 (MQ=255)          ggaTCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGt                                                                                                                        <  1:404141/139‑1 (MQ=255)          ggaTCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGt                                                                                                                        <  1:1699553/139‑1 (MQ=255)            aTCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCg                                                                                                                      <  1:723784/139‑1 (MQ=255)                tGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCGCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTAt                                                                                                                  <  2:1576305/139‑1 (MQ=255)                 ggCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAgg                                                                                                                         <  1:802020/131‑1 (MQ=255)                 ggCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAgg                                                                                                                         >  2:802020/1‑131 (MQ=255)                 ggCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTa                                                                                                                   <  1:1691127/137‑1 (MQ=255)                 ggCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTa                                                                                                                   >  2:1691127/1‑137 (MQ=255)                    ttttGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGat                                                                                                              <  2:796447/139‑1 (MQ=255)                     tttGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGata                                                                                                             >  2:241785/1‑139 (MQ=255)                           gAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCa                                                                                                       <  2:810274/139‑1 (MQ=255)                           gAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCa                                                                                                       <  1:241785/139‑1 (MQ=255)                                 ttCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGata                                                                                                 >  1:1071773/1‑139 (MQ=255)                                         aCTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACgaga                                                                                         <  1:1698336/139‑1 (MQ=255)                                             aTAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAg                                                                                     >  2:1016609/1‑139 (MQ=255)                                              tAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGt                                                                                    <  1:1069406/139‑1 (MQ=255)                                                   aaCCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGgtg                                                                               <  1:758435/139‑1 (MQ=255)                                                    aCCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGgtgt                                                                              >  2:1153453/1‑139 (MQ=255)                                                       tCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCgg                                                                                 <  2:102403/133‑1 (MQ=255)                                                       tCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCgg                                                                                 >  1:102403/1‑133 (MQ=255)                                                        cGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTgag                                                                          >  1:1600728/1‑139 (MQ=255)                                                        cGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTgag                                                                          >  2:1787459/1‑139 (MQ=255)                                                           cgcAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGt                                                                                                                        >  1:650504/1‑90 (MQ=255)                                                           cgcAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGt                                                                                                                        <  2:650504/90‑1 (MQ=255)                                                           cgcAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGc                                                                                                        >  1:991766/1‑106 (MQ=255)                                                           cgcAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGc                                                                                                        <  2:991766/106‑1 (MQ=255)                                                                tttGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCg                                                                  >  2:997320/1‑139 (MQ=255)                                                                   gTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGAcc                                                               <  1:997320/139‑1 (MQ=255)                                                                     ttCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCg                                                             <  1:1787459/139‑1 (MQ=255)                                                                       cGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGActtttc                                                           <  2:1600728/139‑6 (MQ=255)                                                                               aacaGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACaa                                                   <  2:1071773/139‑1 (MQ=255)                                                                               aacaGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACaa                                                   <  2:1248196/139‑1 (MQ=255)                                                                                          tGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAAcc                                        >  2:238057/1‑139 (MQ=255)                                                                                                 tGTACGTACTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACAtt                                 <  1:238057/139‑1 (MQ=255)                                                                                                   tacgtacTCATCAGTACGCATCTTCACAA‑GTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGc                              >  2:1404735/1‑139 (MQ=255)                                                                                                         cTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCa                                                                                                       >  2:947015/1‑61 (MQ=255)                                                                                                         cTCATCAGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCa                                                                                                       <  1:947015/61‑1 (MQ=255)                                                                                                           catcaGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCgg                                                                                 <  2:473037/81‑1 (MQ=255)                                                                                                           catcaGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCgg                                                                                 >  1:473037/1‑81 (MQ=255)                                                                                                             tcaGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGAtata                                                                                               <  1:1044929/65‑1 (MQ=255)                                                                                                             tcaGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGAtata                                                                                               >  2:1044929/1‑65 (MQ=255)                                                                                                             tcaGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCa                     >  2:1754127/1‑139 (MQ=255)                                                                                                             tcaGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCa                     >  1:851835/1‑139 (MQ=255)                                                                                                             tcaGTACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCa                     >  1:1156796/1‑139 (MQ=255)                                                                                                                gTACGCATCTTCACAA‑GTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATg                 <  1:1404735/139‑1 (MQ=255)                                                                                                                 tACGCATCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATg                 <  2:1156796/139‑1 (MQ=255)                                                                                                                      aTCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGc                                                                                                        <  1:995114/47‑1 (MQ=255)                                                                                                                      aTCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGc                                                                                                        >  2:995114/1‑47 (MQ=255)                                                                                                                      aTCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGtt            >  1:1548949/1‑139 (MQ=255)                                                                                                                      aTCTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGtt            >  1:1282248/1‑139 (MQ=255)                                                                                                                        cTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGata                                                                                                             >  2:258555/1‑40 (MQ=255)                                                                                                                        cTTCACAATGTTACAGAGACTGATCAGGTCGCTATTGata                                                                                                             <  1:258555/40‑1 (MQ=255)                                                                                                                           cacaATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGtt            <  2:1327686/134‑1 (MQ=255)                                                                                                                           cacaATGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGtt            >  1:1327686/1‑134 (MQ=255)                                                                                                                                tGTTACAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATatcatc  >  1:707351/1‑139 (MQ=255)                                                                                                                                     |                                                                                                                                      CGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATC  >  NC_000913/3515594‑3515859

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