Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,893,403	C→T	S213S (AGC→AGT)	yidZ →	putative DNA‑binding transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,893,403	0	C	T	100.0%	153.2 / NA	40	S213S (AGC→AGT)	yidZ	putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (20/20);  total (20/20)

TCGCTACCGTTAGCCATTGATTATGAAGTGCTGTTTAGTGATGTGCCCTGCGTCTGGTTACGCAAAGATCATCCGGCACTGCATCAAACGTGGAATCTGGACACCTTCTTACGTTATCCGCATATCAGCATTTGCTGGGAACAGAGCGATACCTGGGCGCTGGACAATGTGTTACAGGAGCTGGGACGCGAACGCACGATTGCTATGAGCCTGCCGGAATTCGAGCAGTCACTGTTTATGGCAGCGCAACCCGACAATCTGCTACTGGCGACCGCGCCGCGCTACT  >  NC_000913/3893257‑3893542                                                                                                                                                   |                                                                                                                                            tCGCTACCGTTAGCCATTGATTATGAAGTGCTGTTTAGTGATGTGCCCTGCGTCTGGTTACGCAAAGATCATCCGGCACTGCATCAAACGTGGAATCTGGACACCTTCTTACGTTATCCGCATATCAGCATTTGCTGGGAACAGAGTGAt                                                                                                                                          >  1:216702/1‑150 (MQ=255)     tACCGTTAGCCATTGATTATGAAGTGCTGTTTAGTGATGTGCCCTGCGTCTGGTTACGCAAAGATCATCCGGCACTGCATCAAACGTGGAATCTGGACACCTTCTTACGTTATCCGCATATCAGCATTTGCTGGGAACAGAGTGATACCTg                                                                                                                                     >  2:7035/1‑151 (MQ=255)      aCCGTTAGCCATTGATTATGAAGTGCTGTTTAGTGATGTGCCCTGCGTCTGGTTACGCAAAGATCATCCGGCACTGCATCAAACGTGGAATCTGGACTCCTTCTTACGTTATCCGCATATCAGCATTTGCTGGGAACAGAGTGATACCt                                                                                                                                      <  1:524360/149‑1 (MQ=255)      aCCGTTAGCCATTGATTATGAAGTGCTGTTTAGTGATGTGCCCTGCGTCTGGTTACGCAAAGATCATCCGGCACTGCATCAAACGTGGAATCTGGACACCTTCTTACGTTATCCGCATATCAGCATTTGCTGGGAACAGAGTGATACCt                                                                                                                                      <  1:10881/149‑1 (MQ=255)       ccGTTAGCCATTGATTATGAAGTGCTGTTTAGTGATGTGCCCTGCGTCTGGTTACGCAAAGATCATCCGGCACTGCATCAAACGTGGAATCTGGACACCTTCTTACGTTATCCGCATATCAGCATTTGCTGGGAACAGAGTGATACCTggg                                                                                                                                   >  2:153873/1‑151 (MQ=255)           tAGCCATTGATTATGAAGTGCTGTTTAGTGATGTGCCCTGCGTCTGGTTACGCAAAGATCATCCGGCACTGCATCAAACGTGGAATCTGGACACCTTCTTACGTTATCCGCATATCAGCATTTGCTGGGAACAGAGTGATACCTGGGCGCt                                                                                                                               >  2:431651/1‑151 (MQ=255)                   gATTATGAAGTGCTGTTTAGTGATGTGCCCTGCGTCTGGTTACGCAAAGATCATCCGGCACTGCATCAAACGTGGAATCTGGACACCTTCTTACGTTATCCGCATATCAGCATTTGCTGGGAACAGAGTGATACCTGGGCGCTGGACAAtg                                                                                                                       <  1:396276/151‑1 (MQ=255)                   gATTATGAAGTGCTGTTTAGTGATGTGCCCTGCGTCTGGTTACGCAAAGATCATCCGGCACTGCATCAAACGTGGAATCTGGACACCTTCTTACGTTATCCGCATATCAGCATTTGCTGGGAACAGAGTGATACCTGGGCGCTGGACAAt                                                                                                                        <  2:190068/150‑1 (MQ=255)                              gCTGTTTAGTGATGTGCCCTGCGTCTGGTTACGCAAAGATCATCCGGCACTGCATCAAACGTGGAATCTGGACACCTTCTTACGTTATCCGCATATCAGCATTTGCTGGGAACAGAGTGATACCTGGGCGCTGGACAATGTGTTACAGGa                                                                                                             <  1:315326/150‑1 (MQ=255)                                     aGTGATGTGCCCTGCGTCTGGTTACGCAAAGATCATCCGGCACTGCATCAAACGTGGAATCTGGACACCTTCTTACGTTATCCGCATATCAGCATTTGCTGGGAACAGAGTGATACCTGGGCGCTGGACAATGTGTTACAGGAGCTGGGAc                                                                                                     >  1:171191/1‑151 (MQ=255)                                        gATGTGCCCTGCGTCTGGTTACGCAAAGATCATCCGGCACTGCATCAAACGTGGAATCTGGACACCTTCTTACGTTATCCGCATATCAGCATTTGCTGGGAACAGAGTGATACCTGGGCGCTGGACAATGTGTTACAGGAGCTGGGAcgcg                                                                                                  >  1:255707/1‑151 (MQ=255)                                           gtgCCCTGCGTCTGGTTACGCAAAGATCATCCGGCACTGCATCAAACGTGGAATCTGGACACCTTCTTACGTTATCCGCATATCAGCATTTGCTGGGAACAGAGTGATACCTGGGCGCTGGACAATGTGTTACAGGAGCTGGGACGCGAAc                                                                                               <  2:28896/151‑1 (MQ=255)                                             gCCCTGCGTCTGGTTACGCAAAGATCATCCGGCACTGCATCAAACGTGGAATCTGGACACCTTCTTACGTTATCCGCATATCAGCATTTGCTGGGAACAGAGTGATACCTGGGCGCTGGACAATGTGTTACAGGAGCTGGGACGCGAACg                                                                                              <  1:444684/150‑1 (MQ=255)                                             gCCCTGCGTCTGGTTACGCAAAGATCATCCGGCACTGCATCAAACGTGGAATCTGGACACCTTCTTACGTTATCCGCATATCAGCATTTGCTGGGAACAGAGTGATACCTGGGCGCTGGACAATGTGTTACAGGAGCTGGGACGCGAACGc                                                                                             >  1:131047/1‑151 (MQ=255)                                             gCCCTGCGTCTGGTTACGCAAAGATCATCCGGCACTGCATCAAACGTGGAATCTGGACACCTTCTTACGTTATCCGCATATCAGCATTTGCTGGGAACAGAGTGATACCTGGGCGCTGGACAATGTGTTACAGGAGCTGGGACGCGAACGc                                                                                             >  1:237067/1‑151 (MQ=255)                                                 tGCGTCTGGTTACGCAAAGATCATCCGGCACTGCATCAAACGTGGAATCTGGACACCTTCTTACGTTATCCGCATATCAGCATTTGCTGGGAACAGAGTGATACCTGGGCGCTGGACAATGTGTTACAGGAGCTGGGACGCGAACGCACg                                                                                          >  1:430721/1‑150 (MQ=255)                                                         gTTACGCAAAGATCATCCGGCACTGCATCAAACGTGGAATCTGGACACCTTCTTACGTTATCCGCATATCAGCATTTGCTGGGAACAGAGTGATACCTGGGCGCTGGACAATGTGTTACAGGAGCTGGGACGCGAACGCACGATTGCTATg                                                                                 <  2:583595/151‑1 (MQ=255)                                                                aaaGATCATCCGGCACTGCATCAAACGTGGAATCTGGACACCTTCTTACGTTATCCGCATATCAGCATTTGCTGGGAACAGAGTGATACCTGGGCGCTGGACAATGTGTTACAGGAGCTGGGACGCGAACGCACGATTGCTATGAGCCTGc                                                                          >  1:379470/1‑151 (MQ=255)                                                                         ccGGCACTGCATCAAACGTGGAATCTGGACACCTTCTTACGTTATCCGCATATCAGCATTTGCTGGGAACAGAGTGATACCTGGGCGCTGGACAATGTGTTACAGGAGCTGGGACGCGAACGCACGATTGCTATGAGCCTGCCGGAATTc                                                                  <  1:535003/150‑1 (MQ=255)                                                                         ccGGCACTGCATCAAACGTGGAATCTGGACACCTTCTTACGTTATCCGCATATCAGCATTTGCTGGGAACAGAGTGATACCTGGGCGCTGGACAATGTGTTACAGGAGCTGGGACGCGAACGCACGATTGCTATGAGCCTGCCGGAATTc                                                                  <  1:7036/150‑1 (MQ=255)                                                                         ccGGCACTGCATCAAACGTGGAATCTGGACACCTTCTTACGTTATCCGCATATCAGCATTTGCTGGGAACAGAGTGATACCTGGGCGCTGGACAATGTGTTACAGGAGCTGGGACGCGAACGCACGATTGCTATGAGCCTGCCGGAATTCg                                                                 <  1:218643/151‑1 (MQ=255)                                                                             cACTGCATCAAACGTGGAATCTGGACACCTTCTTACGTTATCCGCATATCAGCATTTGCTGGGAACAGAGTGATACCTGGGCGCTGGACAATGTGTTACAGGAGCTGGGACGCGAACGCACGATTGCTATGAGCCTGCCGGAATTCGAGCa                                                             <  1:99201/151‑1 (MQ=255)                                                                                 gCATCAAACGTGGAATCTGGACACCTTCTTACGTTATCCGCATATCAGCATTTGCTGGGAACAGAGTGATACCTGGGCGCTGGACAATGTGTTACAGGAGCTGGGACGCGAACGCACGATTGCTATGAGCCTGCCGGAATTCGAGCAg                                                            >  2:445765/1‑148 (MQ=255)                                                                                 gCATCAAACGTGGAATCTGGACACCTTCTTACGTTATCCGCATATCAGCATTTGCTGGGAACAGAGTGATACCTGGGCGCTGGACAATGTGTTACAGGAGCTGGGACGCGAACGCACGATTGCTATGAGCCTGCCGGAATTCGAGCAGTc                                                          <  1:375975/150‑1 (MQ=255)                                                                                 gCATCAAACGTGGAATCTGGACACCTTCTTACGTTATCCGCATATCAGCATTTGCTGGGAACAGAGTGATACCTGGGCGCTGGACAATGTGTTACAGGAGCTGGGACGCGAACGCACGATTGCTATGAGCCTGCCGGAATTCGAGCAGTc                                                          >  1:554461/1‑150 (MQ=255)                                                                                  cATCAAACGTGGAATCTGGACACCTTCTTACGTTATCCGCATATCAGCATTTGCTGGGAACAGAGTGATACCTGGGCGCTGGACAATGTGTTACAGGAGCTGGGACGCGAACGCACGATTGCTATGAGCCTGCCGGAATTCGAGCAGTCAc                                                        >  2:299165/1‑151 (MQ=255)                                                                                         cGTGGAATCTGGACACCTTCTTACGTTATCCGCATATCAGCATTTGCTGGGAACAGAGTGATACCTGGGCGCTGGACAATGTGTTACAGGAGCTGGGACGCGAACGCACGATTGCTATGAGCCTGCCGGAATTCGAGCAGTCACTGTTTAt                                                 <  2:444432/151‑1 (MQ=255)                                                                                          gTGGAATCTGGACACCTTCTTACGTTATCCGCATATCAGCATTTGCTGGGAACAGAGTGATACCTGGGCGCTGGACAATGTGTTACAGGAGCTGGGACGCGAACGCACGATTGCTATGAGCCTGCCGGAATTCGAGCAGTCACTGTTTATg                                                <  2:139999/151‑1 (MQ=255)                                                                                             gAATCTGGACACCTTCTTACGTTATCCGCATATCAGCATTTGCTGGGAACAGAGTGATACCTGGGCGCTGGACAATGTGTTACAGGAGCTGGGACGCGAACGCACGATTGCTATGAGCCTGCCGGAATTCGAGCAGTCACTGTTTATGGCa                                             <  2:404825/151‑1 (MQ=255)                                                                                              aaTCTGGACACCTTCTTACGTTATCCGCATATCAGCATTTGCTGGGAACAGAGTGATACCTGGGCGCTGGACAATGTGTTACAGGAGCTGGGACGCGAACGCACGATTGCTATGAGCCTGCCGGAATTCGAGCAGTCACTGTTTATGGc                                              >  2:139856/1‑149 (MQ=255)                                                                                                 cTGGACACCTTCTTACGTTATCCGCATATCAGCATTTGCTGGGAACAGAGTGATACCTGGGCGCTGGACAATGTGTTACAGGAGCTGGGACGCGAACGCACGATTGCTATGAGCCTGCCGGAATTCGAGCAGTCACTGTTTATGGCAgcgc                                         >  2:620942/1‑151 (MQ=255)                                                                                                       acCTTCTTACGTTATCCGCATATCAGCATTTGCTGGGACCAGAGTGATACCTGGGCGCTGGACAATGTGTTACAGGAGCTGGGACGCGAACGCACGATTGCTATGAGCCTGCCGGAATTCGAGCAGTCTCTGTTTATGGCAGCGCAAccc                                    <  1:68953/150‑1 (MQ=255)                                                                                                                    aTCCGCATATCAGCATTTGCTGGGAACAGAGTGATACCTAGGCGCTGGACAATGTGTTACAGGAGCTGGGACGCGAACGCACGATTGCTATGAGCCTGCCGGAATTCGAGCAGTCACTGTTTATGGCAGCGCAACCCGACAATCTGCTAc                       <  1:619688/150‑1 (MQ=255)                                                                                                                      ccGCATATCAGCATTTGCTGGGAACAGAGTGATACCTGGGCGCTGGACAATGTGTTACAGGAGCTGGGACGCGAACGCACGATTGCTATGAGCCTGCCGGAATTCGAGCAGTCACTGTTTATGGCAGCGCAACCCGACAATCTGCTACTgg                    >  1:522450/1‑151 (MQ=255)                                                                                                                          atatCAGCATTTGCTGGGAACAGAGTGATACCTGGGCGCTGGACAATGTGTTACAGGAGCTGGGACGCGAACGCACGATTGCTATGAGCCTGCCGGAATTCGAGCAGTCACTGTTTATGGCAGCGCAACCCGACAATCTGCTACTGGCGa                 <  1:364408/150‑1 (MQ=255)                                                                                                                              cAGCATTTGCTGGGAACAGAGTGATACCTGGGCGCTGGACAATGTGTTACAGGAGCTGGGACGCGAACGCACGATTGCTATGAGCCTGCCGGAATTCGAGCAGTCACTGTTTATGGCAGCGCAACCCGACAATCTGCTACTGGCGAccgc             <  1:456007/150‑1 (MQ=255)                                                                                                                                gCATTTGCTGGGAGCAGAGTGATACCTGGGCGCTGGACAATGTGTTACAGGAGCTGGGACGCGAACGCACGATTGCTATGAGCCTGCCGGAATTCGAGCAGTCACTGTTTATGGCAGCGCAACCCGACAATCTGCTACTGGCGAccgcgcc          >  1:21648/1‑151 (MQ=255)                                                                                                                                 cATTTGCTGGGAACAGAGTGATACCTGGGCGCTGGACAATGTGTTACAGGAGCTGGGACGCGAACGCACGATTGCTATGAGCCTGCCGGAATTCGAGCAGTCACTGTTTATGGCAGCGCAACCCGACAATCTGCTACTGGCGAccgcgcc          >  1:357343/1‑150 (MQ=255)                                                                                                                                     tGCTGGGAACAGAGTGATACCTGGGCGCTGGACAATGTGTTACAGGAGCTGGGACGCGAACGCACGATTGCTATGAGCCTGCCGGAATTCGAGCAGTCACTGTTTATGGCAGCGCAACCCGACAATCTGCTACTGGCGACCGCGCCGCGCt     >  1:521423/1‑151 (MQ=255)                                                                                                                                        tGGGAACAGAGTGATACCTGGGCGCTGGACAATGTGTTACAGGAGCTGGGACGCGAACGCACGATTGCTATGAGCCTGCCGGAATTCGAGCAGTCACTGTTTATGGCAGCGCAACCCGACAATCTGCTACTGGCGACCGCGCCGCGCTACt  >  2:602667/1‑151 (MQ=255)                                                                                                                                                   |                                                                                                                                            TCGCTACCGTTAGCCATTGATTATGAAGTGCTGTTTAGTGATGTGCCCTGCGTCTGGTTACGCAAAGATCATCCGGCACTGCATCAAACGTGGAATCTGGACACCTTCTTACGTTATCCGCATATCAGCATTTGCTGGGAACAGAGCGATACCTGGGCGCTGGACAATGTGTTACAGGAGCTGGGACGCGAACGCACGATTGCTATGAGCCTGCCGGAATTCGAGCAGTCACTGTTTATGGCAGCGCAACCCGACAATCTGCTACTGGCGACCGCGCCGCGCTACT  >  NC_000913/3893257‑3893542

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