Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,074,728	C→T	A20A (GCC→GCT)	yihW →	putative DNA‑binding transcriptional regulator YihW

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,074,728	0	C	T	100.0%	139.5 / NA	41	A20A (GCC→GCT)	yihW	putative DNA‑binding transcriptional regulator YihW
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (20/21);  total (20/21)

TGACTGTGATCAAACCCGATCTTTTTTGTCACTTTTTGTATAAAATGCCAGGGTGATGGTTTTTCGAGGAATTTTCATGAGCCTTACCGAACTAACCGGTAACCCGCGGCACGACCAACTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCC  >  NC_000913/4074593‑4074863                                                                                                                                        |                                                                                                                                        tGACTGTGATCAAACCCGATCTTTTTTGTCGCTTTTTGTATAAAATGCCAGGGTGATGGTTTTTCGAGGAATTTTCATGAGCCTTACCGAACTAACCGGTAACCCGCGGCACGACCAACTCCTCATGCTGATCGCTGAGc                                                                                                                                     <  1:194463/140‑1 (MQ=255)  gACTGTGATCAAACCCGTTCTTTTTTGTCGCATTTCGTACAAAATGCCCGGGTGATGGTTTTTCGAGGAATTTTCATGAGCCTTACCGAACTAACCGGTAACCCGCGGCACGACCAACTCCTCATGCTGATCGCTGAGCg                                                                                                                                    <  2:220000/140‑1 (MQ=255)   aCTGTGATCAAACCCGATCTTTTTTGTCGCTTTTTGTATAAAATGCCAGGGTGATGGTTTTTCGAGGAATTTTCATGAGCCTTACCGAACTAACCGGTAACCCGCGGCACGACCAACTCCTCATGCTGATCGCTGAGCGt                                                                                                                                   <  1:137931/140‑1 (MQ=255)      gtgATCAACCCCGATCTCTTTTTTCGCTTTTTGTACAAAATGCCAGGGAGATGGTTTTTCGAGGAATTTTCATGAGCCTTACCGAACTAACCGGTAACCCGCGGCACGACCAACTCCTCATGCTGATCGCTGAGCGTggg                                                                                                                                <  2:375383/140‑1 (MQ=255)       tgATCAAACCCGATCTTTTTTGTCGCTTTTTGTATAAAATGCCAGGGTGATGGTTTTTCGAGGAATTTTCATGAGCCTTACCGAACTAACCGGTAACCCGCGGCACGACCAACTCCTCATGCTGATCGCTGAGCGTGGGt                                                                                                                               >  2:581289/1‑140 (MQ=255)         aTCAAACCCGATCTTTTTTGTCGCTTTTTGTATAAAATGCCAGGGTGATGGTTTTTCGAGGAATTTTCATGAGCCTTACCGAACTAACCGGTAACCCGCGGCCCGACCAAATCCTCAAGCTGATAGCTGAGCGTGGGtat                                                                                                                             >  1:287417/1‑140 (MQ=255)         aTCAAACCCGATCTTTTTTGTCGCTTTTTGTATAAAATGCCAGGGTGATGGTTTTTCGAGGAATTTTCATGAGCCTTACCGAACTAACCGGTAACCCGCGGCACGACCAACTCCTCATGCTGATCGCTGAGCGTGGGtat                                                                                                                             >  1:822362/1‑140 (MQ=255)         aTCAAACCCGATCTTTTTTGTCGCTTTTTGTATAAAATGCCAGGGTGATGGTTTTTCGAGGAATTTTCATGAGCCTTACCGAACTAACCGGGAACCCGCGGCACGACCGACTCCTCATGCTGGTCGCTGAGCGTGGGtat                                                                                                                             >  2:225247/1‑140 (MQ=255)         aTCAAACCCGATCTTATTTGTCCCTTTTTCTATAAAATGCCAGGGTGATGGTTTTTCGAGGAATTTTCATGAGCCTTACCGAACTAACCGGTAACCCGCGGCACGACCAACTCCTCATGCTGATCGCTGAGCGTGGGtat                                                                                                                             <  1:243126/140‑1 (MQ=255)              aCCCGATCTTTTTTGTCGCTTTTTGTATAAAATGCCAGGGTGATGGTTTTTCGAGGAATTTTCATGAGCCTTACCGAACTAACCGGTAACCCGCGGCACGACCAACTCCTCATGCTGATCGCTGAGCGTGGGTatat                                                                                                                           >  1:206219/1‑137 (MQ=255)              aCCCGATCTTTTTTGTCGCTTTTTGTATAAAATGCCAGGGTGATGGTTTTTCGAGGAATTTTCATGAGCCTTACCGAACTAACCGGTAACCCGCGGCACGACCAACTCCTCATGCTGATCGCTGAGCGTGGGTatat                                                                                                                           <  2:206219/137‑1 (MQ=255)                   aTCTTTTTTGTCGCTTTTTGTATAAAATGCCAGGGTGATGGTTTTTCGAGGAATTTTCATGAGCCTTACCGAACTAACCGGTAACCCGCGGCACGACCAACTCCTCATGCTGATCGCTGAGCGTg                                                                                                                                  >  1:25120/1‑125 (MQ=255)                   aTCTTTTTTGTCGCTTTTTGTATAAAATGCCAGGGTGATGGTTTTTCGAGGAATTTTCATGAGCCTTACCGAACTAACCGGTAACCCGCGGCACGACCAACTCCTCATGCTGATCGCTGAGCGTg                                                                                                                                  <  2:25120/125‑1 (MQ=255)                   aTCTTTTTTGTCGCTTTTTGTATAAAATGCCAGGGTGATGGTTTTTCGAGGAATTTTCATGAGCCTTACCGAACTAACCGGTAACCCGCGGCACGACCAACTCCTCATGCTGATCGCTGAGCGTGGGTATATGAATATtg                                                                                                                   >  2:708876/1‑140 (MQ=255)                              cgcTTTTTGTATAAAATGCCAGGGTGATGGTTTTTCGAGGAATTTTCATGAGCCTTACCGAACTAACCGGTAACCCGCGGCACGACCAACTCCTCATGCTGATCGCTGAGCGTGGGTATATGAATATTGATGAGCTGGCa                                                                                                        >  2:186123/3‑140 (MQ=255)                                       tataAAATGCCAGGGTGATGGTTTTTCGAGGAATTTTCATGAGCCTTACCGAACTAACCGGTAACCCGCGGCACGACCAACTCCTCATGCTGATCGCTGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATctgctg                                                                                               >  2:803408/1‑140 (MQ=255)                                         taAAATGCCAGGGTGATGGTTTTTCGAGGAATTTTCATGAGCCTTACCGAACTAACCGGTAACCCGCGGCACGACCAACTCCTCATGCTGATCGCTGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGa                                                                                             <  1:581289/140‑1 (MQ=255)                                         taAAATGCCAGGGTGATGGTTTTTCGAGGAATTTTCATGAGCCTTACCGAACTAACCGGTAACCCGCGGCACGACCAACTCCTCATGCTGATCGCTGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGa                                                                                             <  1:186123/140‑1 (MQ=255)                                          aaaaTGCCAGGGTGATGGTTTTTCGAGGAATTTTCATGAGCCTTACCGAACTAACCGGTAACCCGCGGCACGACCAACTCCTCATGCTGATCGCTGAGCGTGGGTATATGAATTTTGATGAGCTGGCAAATCTGCTGGAt                                                                                            <  1:799493/140‑1 (MQ=255)                                              tGCCAGGGTGATGGTTTTTCGAGGAATTTTCATGAGCCTTACCGAACTAACCGGTAACCCGCGGCACGACCAACTCCTCATGCTGATCGCTGAGCGTGGGTATATGAATATTGATGAGCTGGCaa                                                                                                       >  1:522610/1‑125 (MQ=255)                                              tGCCAGGGTGATGGTTTTTCGAGGAATTTTCATGAGCCTTACCGAACTAACCGGTAACCCGCGGCACGACCAACTCCTCATGCTGATCGCTGAGCGTGGGTATATGAATATTGATGAGCTGGCaa                                                                                                       <  2:522610/125‑1 (MQ=255)                                                    ggTGATGGTTTTTCGAGGAATTTTCATGAGCCTTACCGAACTAACCGGTAACCCGCGGCACGACCAACTCCTCATGCTGATCGCTGAGCGTGGGTATATGAATACTGACGAGCTGGCGAATCTGCTGGATGTCTCCACGc                                                                                  >  2:352017/1‑140 (MQ=255)                                                        aTGGTTTTTCGAGGAATTTTCATGAGCCTTACCGAACTAACCGGTAACCCGCGGCACGACCAACTCCTCATGCTGATCGCTGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGAc                                                                              <  2:694686/140‑1 (MQ=255)                                                         tGGTTTTTCGAGGAATTTTCATGAGCCTTACCGAACTAACCGGTAACCCGCGGCACGACCAACTCCTCATGCTGATCGCTGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACg                                                                             >  2:824653/1‑140 (MQ=255)                                                         tGGTTTTTCGAGGAATTTTCATGAGCCTTACCGAACTAACCGGTAACCCGCGGCACGACCAACTCCTCATGCTGATCGCTGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACg                                                                             <  1:803408/140‑1 (MQ=255)                                                         tGGTTTTTCGAGGAATTTTCATGAGCCTTACCGAACTAACCGGTAACCCGCGGCACGACCAACTCCTCATGCTGATCGCTGAGCGTGGGTATATGAATATCGATTAACAGGCAAATCGGCTGCATGTCTCCACGCAGACg                                                                             >  2:352510/1‑140 (MQ=255)                                                              tttCGAGGAATTTTCATGAGCCTTACCGAACTAACCGGTAACCCGCGGCACGACCAACTCCTCATGCTGATCGCTGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTccg                                                                        >  1:659087/1‑140 (MQ=255)                                                                tCGAGGAATTTTCATGAGCCTTTCCGAACTAACCGTTAACCCGCGGCAGGACCAACTCCTCATGCTGATCGCTGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTccgcc                                                                      <  2:287417/140‑1 (MQ=255)                                                                  gAGGAATTTTCATGAGCCTTACCGAACTAACCGGTAACCCGCGGCACGACCAACTCCTCATGCTGATCGCTGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCgg                                                                    >  2:12907/1‑140 (MQ=255)                                                                    ggAATTTTCATGAGCCTTACCGAACTAACCGGTAACCCGCGGCACGACCAACTCCTCATGCTGATCGCTGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGa                                                                  <  2:659087/140‑1 (MQ=255)                                                                      aaTTTTCATGAGCCTTACCGAACTAACCGGTAACCCGCGGCACGACCAACTCCTCATGCTGATCGCTGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGata                                                                <  2:668574/140‑1 (MQ=255)                                                                                       ccGAACTAACCGGTAACCCGCGGCACGACCAACTCCTCATGCTGATCGCTGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAg                                               <  1:824653/140‑1 (MQ=255)                                                                                       ccGAACTAACCGGTAACCCGCGGCACGACCAACTCCTCATGCTGATCGCTGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGAAAATTAAGCGAg                                               <  1:12907/140‑1 (MQ=255)                                                                                            cTAACCGGTAACCCGCGGCACGACCAACTCCTCATGCTGATCGCTGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGCCTCCACGCAGACCGTCGGCCGGGATATTCGAAAATTAAGCGAGCAAgg                                          >  2:320580/1‑140 (MQ=255)                                                                                             tAACCGGTAACCCGCGGCACGAGCATCTCGTCATGCTGATCGCTGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGc                                         <  1:320580/140‑1 (MQ=255)                                                                                                      aCCCGCGGCACGACCAACTCCTCATGCTGATCGCTGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTAcg                                <  2:822362/140‑1 (MQ=255)                                                                                                                      aCTCCTCATGCTGATCGCTGAGAGTCGGTTTATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGgcg                <  1:225247/140‑1 (MQ=255)                                                                                                                             aTGCTGATCGCTGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCg         >  1:774589/1‑140 (MQ=255)                                                                                                                                  gATCGCTGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCtt    >  2:664698/1‑140 (MQ=255)                                                                                                                                    tCGCTGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTcc  >  1:785703/1‑140 (MQ=255)                                                                                                                                    tCGCTGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTcc  >  1:499179/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                        TGACTGTGATCAAACCCGATCTTTTTTGTCACTTTTTGTATAAAATGCCAGGGTGATGGTTTTTCGAGGAATTTTCATGAGCCTTACCGAACTAACCGGTAACCCGCGGCACGACCAACTCCTCATGCTGATCGCCGAGCGTGGGTATATGAATATTGATGAGCTGGCAAATCTGCTGGATGTCTCCACGCAGACGGTCCGCCGGGATATTCGTAAATTAAGCGAGCAAGGCCTGATTACGCGCCATCACGGTGGCGCGGGTCGGGCTTCC  >  NC_000913/4074593‑4074863

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