Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,529,959	(T)7→6	coding (38/483 nt)	ydcD →	uncharacterized protein YdcD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,529,953	0	T	.	100.0%	223.4 / NA	50	coding (32/483 nt)	ydcD	uncharacterized protein YdcD
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (12/38);  total (12/38)

TGCCAAATGTAGCGATGCATGTATTGAAAAATATAACACTGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTAT  >  NC_000913/1529815‑1530082                                                                                                                                           |                                                                                                                                  tGCCAAATGTAGCGATGCATGTATTGAAAAATATAACACTGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCAtt                                                                                                                                  >  1:919349/1‑140 (MQ=255)  gCCAAATGTAGCGATGCATGTATTGAAAAATATAACACTGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCAttt                                                                                                                                 <  2:1184142/140‑1 (MQ=255)  gCCAAATGTAGCGATGCATGTATTGAAAAATATAACACTGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCAttt                                                                                                                                 <  1:573901/140‑1 (MQ=255)         gTAGCGATGCATGTATTGAAAAATATAACACTGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAGAGACAAACTTGTCAttttttaatt                                                                                                                          <  1:423052/140‑5 (MQ=255)              gATGCATGTATTGAAAAATATAACACTGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCA‑TTTTTTAATTACCCt                                                                                                                    <  1:1631871/140‑1 (MQ=255)               aTGCATGTATTGAAAAATATAACACTGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCA‑TTTTTTAATTACCCtt                                                                                                                   >  2:927064/1‑140 (MQ=255)               aTGCATGTATTGAAAAATATAACACTGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCA‑TTTTTTAATTACCCtt                                                                                                                   <  2:677300/140‑1 (MQ=255)                  cATGTATTGAAAAATATAACACTGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCA‑TTTTTTAATTACCCTTTGt                                                                                                                >  1:989936/1‑140 (MQ=255)                  cATGTATTGAAAAATATAACACTGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCA‑TTTTTTAATTACCCTTTGt                                                                                                                >  1:861716/1‑140 (MQ=255)                        ttGAAAAATATAACACTGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCA‑TTTTTTAATTACCCTTTGTTTAttt                                                                                                          <  1:330585/140‑1 (MQ=255)                            aaaaTATAACACTGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCA‑TTTTTTAATTACCCTTTGTTTATttgttg                                                                                                      >  2:1383197/1‑140 (MQ=255)                            aaaaTATAACACTGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCA‑TTTTTTAATTACCCTTTGTTTATttgttg                                                                                                      >  1:281465/1‑140 (MQ=255)                              aaTATAACACTGGTCAACTCTTCGGGTTATATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCA‑TTTTTTAATTACCCTTTGTTTATTTgttgtt                                                                                                    <  2:1330230/140‑1 (MQ=255)                              aaTATAACACTGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCA‑TTTTTTAATTACCCTTTGTTTATTTgttgtt                                                                                                    <  2:498821/140‑1 (MQ=255)                                   aaCACTGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCA‑TTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCt                                                                                               <  2:1194012/140‑1 (MQ=255)                                   aaCACTGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCA‑TTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCt                                                                                               <  2:456125/140‑1 (MQ=255)                                     cacTGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCA‑TTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTg                                                                                             <  1:288090/140‑1 (MQ=255)                                     cacTGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCA‑TTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTg                                                                                             <  2:938146/140‑1 (MQ=255)                                          ggCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCA‑TTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGttt                                                                                        >  2:475032/1‑140 (MQ=255)                                          ggCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCA‑TTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGttt                                                                                        <  1:1076988/140‑1 (MQ=255)                                             aaCTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCA‑TTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGAt                                                                                     <  1:1590693/140‑1 (MQ=255)                                               ctctTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCA‑TTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATaa                                                                                   <  1:54410/140‑1 (MQ=255)                                               ctctTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCA‑TTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATaa                                                                                   >  2:45737/1‑140 (MQ=255)                                               ctctTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCA‑TTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATaa                                                                                   >  2:1368941/1‑140 (MQ=255)                                                    cGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCA‑TTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAAtt                                                                              >  1:1127043/1‑140 (MQ=255)                                                        ttAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCA‑TTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCaaa                                                                          <  2:821252/140‑1 (MQ=255)                                                            aTGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCA‑TTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCt                                                                      >  1:70649/1‑140 (MQ=255)                                                            aTGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCA‑TTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCt                                                                      >  2:1511154/1‑140 (MQ=255)                                                             tGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCA‑TTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTg                                                                     <  2:281465/140‑1 (MQ=255)                                                               gTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCA‑TTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGAc                                                                   <  1:518508/140‑1 (MQ=255)                                                                tataAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCA‑TTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACa                                                                  <  2:1265273/140‑1 (MQ=255)                                                                    aaaGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCA‑TTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGt                                                              <  1:1383197/140‑1 (MQ=255)                                                                    aaaGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCA‑TTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGt                                                              <  2:1127043/140‑1 (MQ=255)                                                                           aTAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCA‑TTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTAt                                                       <  2:989936/140‑1 (MQ=255)                                                                           aTAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCA‑TTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTAt                                                       <  2:919349/140‑1 (MQ=255)                                                                           aTAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCA‑TTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTAt                                                       >  1:1087464/1‑140 (MQ=255)                                                                           aTAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCA‑TTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTAt                                                       <  2:1017923/140‑1 (MQ=255)                                                                             aaTCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCA‑TTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATgg                                                     <  1:1251297/140‑1 (MQ=255)                                                                             aaTCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCA‑TTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATgg                                                     <  2:182951/140‑1 (MQ=255)                                                                             aaTCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCA‑TTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATgg                                                     <  1:927064/140‑1 (MQ=255)                                                                                ccAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCA‑TTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACa                                                  <  1:45737/140‑1 (MQ=255)                                                                                  aaTAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCA‑TTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATa                                                <  1:475032/140‑1 (MQ=255)                                                                                  aaTAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCA‑TTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATa                                                <  2:1065994/140‑1 (MQ=255)                                                                                   aTAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCA‑TTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATaa                                               <  2:962586/140‑1 (MQ=255)                                                                                     aaaaGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCA‑TTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATaaaa                                             <  2:1262822/140‑1 (MQ=255)                                                                                                 aGGGAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCA‑TTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCg                                 <  1:835064/140‑1 (MQ=255)                                                                                                  gggAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCA‑TTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGt                                <  2:861716/140‑1 (MQ=255)                                                                                                  gggAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCA‑TTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGt                                <  1:271403/140‑1 (MQ=255)                                                                                                   ggAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCA‑TTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGtt                               <  2:70649/140‑1 (MQ=255)                                                                                                   ggAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCA‑TTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGtt                               <  1:1368941/140‑1 (MQ=255)                                                                                                    gAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCA‑TTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTg                              <  2:408313/140‑1 (MQ=255)                                                                                                                  cAAATTATTAGAGGCAAACTTGTCA‑TTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCt                <  1:1511154/140‑1 (MQ=255)                                                                                                                     attattAGAGGCAAACTTGTCA‑TTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGaa             <  2:1461259/140‑1 (MQ=255)                                                                                                                                cAAACTTGTCA‑TTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAtattat  <  2:1087464/140‑1 (MQ=255)                                                                                                                                           |                                                                                                                                  TGCCAAATGTAGCGATGCATGTATTGAAAAATATAACACTGGGCAACTCTTCGGGTTAGATGGTATAAAAGCAGATAATCCAATAAAAGCAAAGCAAGGGAGTTCAGATGCTTCAAATTATTAGAGGCAAACTTGTCATTTTTTTAATTACCCTTTGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTAT  >  NC_000913/1529815‑1530082

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