Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,581,621	G→A	A476A (GCC→GCT)	hsdM ←	DNA methyltransferase M

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,581,621	0	G	A	100.0%	170.2 / NA	53	A476A (GCC→GCT)	hsdM	DNA methyltransferase M
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (27/26);  total (27/26)

CAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCGGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTT  >  NC_000913/4581492‑4581742                                                                                                                                  |                                                                                                                          cAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATa                                                                                                                  <  1:689490/139‑1 (MQ=255)          cTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGTCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTa                                                                                                         <  2:847054/139‑1 (MQ=255)               aaTCGGCCTCATCGCTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCttt                                                                                                    <  1:395980/139‑1 (MQ=255)                                            cGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATGCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGAtt                                                                       >  2:755331/1‑139 (MQ=255)                                              cATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCttt                                                                                                    >  2:1288952/1‑108 (MQ=255)                                              cATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCttt                                                                                                    <  1:1288952/108‑1 (MQ=255)                                               aTCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTtg                                                                    <  1:600072/139‑1 (MQ=255)                                                  aGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTtgcgg                                                                 <  1:1082194/139‑1 (MQ=255)                                                  aGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTtgcgg                                                                 <  1:508517/139‑1 (MQ=255)                                                  aGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTtgcgg                                                                 >  2:1213210/1‑139 (MQ=255)                                                    cgcATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTtgcggtg                                                               >  1:647935/1‑139 (MQ=255)                                                       aTCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTtgcggtg                                                               >  2:523819/1‑136 (MQ=255)                                                       aTCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTtgcggtg                                                               <  1:523819/136‑1 (MQ=255)                                                           aGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATcc                                                        >  1:479247/1‑139 (MQ=255)                                                             ttCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATc                                                                                                                      >  1:839761/1‑75 (MQ=255)                                                             ttCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATc                                                                                                                      <  2:839761/75‑1 (MQ=255)                                                              tCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCtt                                                                                                     <  1:985715/91‑1 (MQ=255)                                                              tCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCtt                                                                                                     >  2:985715/1‑91 (MQ=255)                                                                  aCAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAAt                                                                                                                   >  2:80361/1‑73 (MQ=255)                                                                  aCAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAAt                                                                                                                   <  1:80361/73‑1 (MQ=255)                                                                     gcgcCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCt                                              >  1:1119327/1‑139 (MQ=255)                                                                     gcgcCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCt                                              >  2:1023323/1‑139 (MQ=255)                                                                      cgcCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCAt                                                                                                                       >  1:36637/1‑65 (MQ=255)                                                                      cgcCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCAt                                                                                                                       <  2:36637/65‑1 (MQ=255)                                                                      cgcCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATa                                                                                                                  >  1:853613/1‑70 (MQ=255)                                                                      cgcCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCAc                                                  >  2:1155942/1‑134 (MQ=255)                                                                      cgcCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCAc                                                  <  1:1155942/134‑1 (MQ=255)                                                                      cgcCTGTACCAGTTCACCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATa                                                                                                                  <  2:853613/70‑1 (MQ=255)                                                                                  ttCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAg                                 >  1:500676/1‑139 (MQ=255)                                                                                   tCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGcc                                                                                               <  2:1170331/76‑1 (MQ=255)                                                                                   tCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGcc                                                                                               >  1:1170331/1‑76 (MQ=255)                                                                                   tCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGc                                <  1:1213210/139‑1 (MQ=255)                                                                                           cGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATa                                                                                                                  >  2:1260274/1‑49 (MQ=255)                                                                                           cGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATa                                                                                                                  <  1:1260274/49‑1 (MQ=255)                                                                                                 tGCCGCTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGa                  >  2:1017009/1‑139 (MQ=255)                                                                                                      cTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATa                                                                                                                  >  1:24610/1‑38 (MQ=38)                                                                                                      cTAATACATCCGGCTCCGGCAGGCTGTCAGCATCAATa                                                                                                                  <  2:24610/38‑1 (MQ=38)                                                                                                             aTCCGGCTCCGGCAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGgtg      <  1:70211/139‑1 (MQ=255)                                                                                                               ccggctccggcAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGc                             <  1:1107890/114‑1 (MQ=255)                                                                                                               ccggctccggcAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGc                             >  2:1107890/1‑114 (MQ=255)                                                                                                                 ggctccggcAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCAAGCGGCTGGTAGCACGATGCTGGTGGGTGtttt  >  2:1011930/1‑139 (MQ=255)                                                                                                                  gctccggcAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGc                                                                                                >  2:417601/1‑44 (MQ=255)                                                                                                                  gctccggcAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGc                                                                                                <  1:417601/44‑1 (MQ=255)                                                                                                                   ctccggcAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGcc                                                                                               >  2:264001/1‑44 (MQ=255)                                                                                                                   ctccggcAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGcc                                                                                               <  1:264001/44‑1 (MQ=255)                                                                                                                   ctccggcAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGc                                                                                                >  1:1200891/1‑43 (MQ=255)                                                                                                                   ctccggcAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGc                                                                                                <  2:1200891/43‑1 (MQ=255)                                                                                                                   ctccggcAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTg                                       <  2:1150059/100‑1 (MQ=255)                                                                                                                   ctccggcAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTg                                       >  1:1150059/1‑100 (MQ=255)                                                                                                                        gcAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGc                                                                                                <  2:376797/38‑1 (MQ=255)                                                                                                                        gcAGGCTGTCAGCATCAATACTGTCTTTATCTTTCAGc                                                                                                >  1:376797/1‑38 (MQ=255)                                                                                                                            gCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCgg        >  2:251007/1‑122 (MQ=255)                                                                                                                            gCTGTCAGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCgg        <  1:251007/122‑1 (MQ=255)                                                                                                                                  |                                                                                                                          CAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCCTGTACCAGTTCGCCCATCGCTTCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCGGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTT  >  NC_000913/4581492‑4581742

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