Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,412,281	C→A	A210D (GCC→GAC)	yhdJ →	DNA adenine methyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,412,281	0	C	A	100.0%	174.2 / NA	51	A210D (GCC→GAC)	yhdJ	DNA adenine methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (18/33);  total (18/33)

ATCCTCCACAGCCTTACAATCATCAAAAAGTACCGGGTAACGTTTGGGATTTTCCGCGCGTGCGTTATTTGATGGATGAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGTATTGAGATCAACAGCGAGTACATCAAAATGGGGCTTCGACG  >  NC_000913/3412146‑3412416                                                                                                                                        |                                                                                                                                        aTCCTCCACAGCCTTACAATCATCAAAAAGTACCGGGTAACGTTTGGGATTTTCCGCGCGTGCGTTATTTGATGGATGAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGActct                                                                                                                                     <  1:127884/140‑1 (MQ=255)        acaGCCTTACAATCATCAAAAAGTACCGGGTAACGTTTGGGATTTTCCGCGCGTGCGTTATTTGATGGATGAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGACTCTTCCAAcc                                                                                                                              <  1:675652/140‑1 (MQ=255)        acaGCCTTACAATCATCAAAAAGTACCGGGTAACGTTTGGGATTTTCCGCGCGTGCGTTATTTGATGGATGAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGACTCTTCCAAcc                                                                                                                              <  2:761242/140‑1 (MQ=255)         caGCCTTACAATCATCAAAAAGTACCGGGTAACGTTTGGGATTTTCCGCGCGTGCGTTATTTGATGGATGAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGACTCTTCCAAccc                                                                                                                             <  1:519798/140‑1 (MQ=255)            ccTTACAATCATCAAAAAGTACCGGGTAACGTTTGGGATTTTCCGCGCGTGCGTTATTTGATGGATGAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGACTCTTCCAACCCAgg                                                                                                                          <  2:735539/140‑1 (MQ=255)            ccTTACAATCATCAAAAAGTACCGGGTAACGTTTGGGATTTTCCGCGCGTGCGTTATTTGATGGATGAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGACTCTTCCAAACCAgg                                                                                                                          <  1:192115/140‑1 (MQ=255)              ttACAATCATCAAAAAGTACCGGGTAACGTTTGGGATTTTCCGCGCGTGCGTTATTTGATGGATGAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGACTCTTCCAACCCAGGCg                                                                                                                        <  1:538786/140‑1 (MQ=255)                      atcaAAAAGTACCGGGTAACGTTTGGGATTTTCCGCGCGTGCGTTATTTGATGGATGAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGACTCTTCCAACCCAGGCGATATCGtt                                                                                                                >  2:541553/1‑140 (MQ=255)                      atcaAAAAGTACCGGGTAACGTTTGGGATTTTCCGCGCGTGCGTTATTTGATGGATGAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGACTCTTCCAACCCAGGCGATATCGtt                                                                                                                >  2:164943/1‑140 (MQ=255)                      atcaAAAAGTACCGGGTAACGTTTGGGATTTTCCGCGCGTGCGTTATTTGATGGATGAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGACTCTTCCAACCCAGGCGATATCGtt                                                                                                                >  1:600542/1‑140 (MQ=255)                             aGTACCGGGTAACGTTTGGGATTTTCCGCGCGTGCGTTATTTGATGGATGAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGACTCTTCCAACCCAGGCGATATCGTTCTCGAcc                                                                                                         <  2:177842/140‑1 (MQ=255)                                aCCGGGTAACGTTTGGGATTTTCCGCGCGTGCGTTATTTGATGGATGAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGACTCTTCCAACCCAGGCGATATCGTTCTCGACCCGt                                                                                                      >  1:307002/1‑140 (MQ=255)                                 tttgatgaCGTTTGGGATTTTCCGCGCGTGCGTTATTTGATGGATGAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGACTCTTCCAACCCAGGCGATATCGTTCTCGACCCGtt                                                                                                     <  2:488760/133‑1 (MQ=255)                                       aaCGTTTGGGATTTTCCGCGCGTGCGTTATTTGATGGATGAATATGAAAACCCCCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGACTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTgg                                                                                               <  2:103972/140‑1 (MQ=255)                                        aCGTTTGGGATTTTCCGCGCGTGCGTTATTTGATGGATGAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGACTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGt                                                                                              <  2:379056/140‑1 (MQ=255)                                         cGTTTGGGATTTTCCGCGCGTGCGTTATTTGATGGATGAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGACTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTa                                                                                             >  2:444785/1‑140 (MQ=255)                                         cGTTTGGGATTTTCCGCGCGTGCGTTATTTGATGGATGAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGACTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTa                                                                                             >  2:119237/1‑140 (MQ=255)                                         cGTTTGGGATTTTCCGCGCGTGCGTTATTTGATGGATGAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGACTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTa                                                                                             >  2:734240/1‑140 (MQ=255)                                                   tttCCGCGCGTGCGTTATTTGATGGATGAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGACTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTtactac                                                                                   <  1:541553/140‑1 (MQ=255)                                                     tCCGCGCGTGCGTTATTTGATGGATGAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGACTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCg                                                                                 <  2:786315/140‑1 (MQ=255)                                                     tCCGCGCGTGCGTTATTTGATGGATGAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGACTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCg                                                                                 <  2:236154/140‑1 (MQ=255)                                                     tCCGCGCGTGCGTTATTTGATGGATGAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGACTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCg                                                                                 >  2:666051/1‑140 (MQ=255)                                                         cgcgTGCGTTATTTGATGGATGAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGACTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGc                                                                             <  1:703522/140‑1 (MQ=255)                                                           cgtgcgtTATTTGATGGATGAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGACTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCg                                                                           <  2:780578/140‑1 (MQ=255)                                                           cgtgcgtTATTTGATGGATGAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGACTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCg                                                                           <  1:214985/140‑1 (MQ=255)                                                             tgcgtTATTTGATGGATGAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGACTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTa                                                                         >  2:526273/1‑140 (MQ=255)                                                                    tttGATGGATGAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGACTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCg                                                                  >  2:525817/1‑140 (MQ=255)                                                                     ttGATGGATGAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGACTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGc                                                                 <  1:525817/140‑1 (MQ=255)                                                                          ggatgAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGACTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCg                                                            <  1:714695/140‑1 (MQ=255)                                                                            atgAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGACTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGa                                                          >  1:219349/1‑140 (MQ=255)                                                                                    gAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGACTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAAtt                                                  <  1:164943/140‑1 (MQ=255)                                                                                    gAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGACTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAAtt                                                  <  1:269817/140‑1 (MQ=255)                                                                                    gAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGACTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAAtt                                                  <  1:526273/140‑1 (MQ=255)                                                                                     aaaaCCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGACTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTc                                                 >  1:438355/1‑140 (MQ=255)                                                                                      aaaCCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGACTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCa                                                >  1:398725/1‑140 (MQ=255)                                                                                       aaCCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGACTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCAt                                               <  1:444785/140‑1 (MQ=255)                                                                                       aaCCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGACTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCAt                                               <  2:42112/140‑1 (MQ=255)                                                                                       aaCCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGACTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCAt                                               <  2:34652/140‑1 (MQ=255)                                                                                       aaCCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGACTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCAt                                               >  2:607776/1‑140 (MQ=255)                                                                                           accCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGACTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGt                                           >  1:31871/1‑140 (MQ=255)                                                                                            cccGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGACTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGTa                                          <  1:541743/140‑1 (MQ=255)                                                                                            cccGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGACTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGTa                                          <  1:248668/140‑1 (MQ=255)                                                                                             ccGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGACTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGTAt                                         >  1:484973/1‑140 (MQ=255)                                                                                              cGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGACTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGTAtt                                        <  2:398725/140‑1 (MQ=255)                                                                                              cGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGACTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGTAtt                                        <  2:729450/140‑1 (MQ=255)                                                                                              cGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGACTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGTAtt                                        <  2:275292/140‑1 (MQ=255)                                                                                              cGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGACTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGTAtt                                        <  2:639979/140‑1 (MQ=255)                                                                                                        aCCGGAAGCCTTACTGAAACGCATTATTCTCGACTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGTATTGAGATCAACa                              <  2:488475/140‑1 (MQ=255)                                                                                                                    aCTGAAACGCATTATTCTCGACTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGTATTGAGATCAACAGCGAGTACATCa                  >  1:50576/1‑140 (MQ=255)                                                                                                                                  ttCTCGACTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGTATTGAGATCAACAGCGAGTACATCAAAATGGGGCTTCGa    <  2:587013/140‑1 (MQ=255)                                                                                                                                    ctcGACTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGTATTGAGATCAACAGCGGGTACATCAAAATGGGGCTTCGACg  >  2:376802/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                        ATCCTCCACAGCCTTACAATCATCAAAAAGTACCGGGTAACGTTTGGGATTTTCCGCGCGTGCGTTATTTGATGGATGAATATGAAAACCACCCGACGCAAAAACCGGAAGCCTTACTGAAACGCATTATTCTCGCCTCTTCCAACCCAGGCGATATCGTTCTCGACCCGTTTGCTGGTAGCTTTACTACCGGTGCCGTAGCCATCGCCAGCGGACGAAAATTCATTGGTATTGAGATCAACAGCGAGTACATCAAAATGGGGCTTCGACG  >  NC_000913/3412146‑3412416

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑