Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,462,842	A→G	intergenic (‑182/+212)	nrdD ← / ← treC	anaerobic ribonucleoside‑triphosphate reductase/trehalose‑6‑phosphate hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,462,842	0	A	G	100.0%	155.4 / NA	42	intergenic (‑182/+212)	nrdD/treC	anaerobic ribonucleoside‑triphosphate reductase/trehalose‑6‑phosphate hydrolase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (15/27);  total (15/27)

AAGGCGTTTTCCGGTTTTTTCGTTTCGTTTTCCGGTGCTTTCCGCAATCGCCATCCGCTTATCCACAACGCCGCACCTTTTTCGTGGGGAGCTGTAGTAGCAATTATAGTCGATTAATACAACATATTGGGTTGGGACGCATTTTAAAGTCTATATATAGTGCTTTGCATCAAGGATGTTTGAGCTTTTTTTGATGTAGCTCAAAGTAAAAAGCAGAGCGTACGGATGACGGGCGCTACAGCGATATGTAAATTTTTTAATGAATTTGC  >  NC_000913/4462706‑4462974                                                                                                                                         |                                                                                                                                     aaGGCGTTTTCCGGTTTTTTCGTTTCGTTTTCCGGTGCTTTCCGCAATCGCCATCCGCTTATCCACAACGCCGCACCTTTTTCGTGGGGAGCTGTAGTAGCAATTATAGTCGATTAATACAACATATTGGGTTGGGgcgc                                                                                                                                   <  1:220159/140‑1 (MQ=255)  aGGCGTTTTCCGGTTTTTTCGTTTCGTTTTCCGGTGCTTTCCGCAATCGCCATCCGCTTATCCACAACGCCGCACCTTTTTCGTGGGGAGCTGTAGTAGCAATTATAGTCGATTAATACAACATATTGGGTTGGGGCGCa                                                                                                                                  >  2:256860/1‑140 (MQ=255)  aGGCGTTTTCCGGTTTTTTCGTTTCGTTTTCCGGTGCTTTCCGCAATCGCCATCCGCTTATCCACAACGCCGCACCTTTTTCGTGGGGAGCTGTAGTAGCAATTATAGTCGATTAATACAACATATTGGGTTGGGGCGCa                                                                                                                                  >  1:596081/1‑140 (MQ=255)    gCGTTTTCCGGTTTTTTCGTTTCGTTTTCCGGTGCTTTCCGCAATCGCCATCCGCTTATCCACAACGCCGCACCTTTTTCGTGGGGAGCTGTAGTAGCAATTATAGTCGATTAATACAACATATTGGGTTGGGGCGCAtt                                                                                                                                >  1:136882/1‑140 (MQ=255)       ttttCCGGTTTTTTCGTTTCGTTTTCCGGTGCTTTCCGCAATCGCCATCCGCTTATCCACAACGCCGCACCTTTTTCGTGGGGAGCTGTAGTAGCAATTATAGTCGATTAATACAACATATTGGGTTGGGGCGCATTTTa                                                                                                                             <  1:152870/140‑1 (MQ=255)       ttttCCGGTTTTTTCGTTTCGTTTTCCGGTGCTTTCCGCAATCGCCATCCGCTTATCCACAACGCCGCACCTTTTTCGTGGGGAGCTGTAGTAGCAATTATAGTCGATTAATACAACATATTGGGTTGGGGCGCATTTTa                                                                                                                             <  2:416740/140‑1 (MQ=255)       ttttCCGGTTTTTTCGTTTCGTTTTCCGGTGCTTTCCGCAATCGCCATCCGCTTATCCACAACGCCGCACCTTTTTCGTGGGGAGCTGTAGTAGCAATTATAGTCGATTAATACAACATATTGGGTTGGGGCGCATTTTa                                                                                                                             <  1:273054/140‑1 (MQ=255)        tttCCGGTTTTTTCGTTTCGTTTTCCGGTGCTTTCCGCAATCGCCATCCGCTTATCCACAACGCCGCACCTTTTTCGTGGGGAGCTGTAGTAGCAATTATAGTCGATTAATACAACATATTGGGTTGGGGCGCATTTTaa                                                                                                                            <  2:87288/140‑1 (MQ=255)        tttCCGGTTTTTTCGTTTCGTTTTCCGGTGCTTTCCGCAATCGCCATCCGCTTATCCACAACGCCGCACCTTTTTCGTGGGGAGCTGTAGTAGCAATTATAGTCGATTAATACAACATATTGGGTTGGGGCGCATTTTaa                                                                                                                            <  2:517026/140‑1 (MQ=255)        tttCCGGTTTTTTCGTTTCGTTTTCCGGTGCTTTCCGCAATCGCCATCCGCTTATCCACAACGCCGCACCTTTTTCGTGGGGAGCTGTAGTAGCAATTATAGTCGATTAATACAACATATTGGGTTGGGGCGCATTTTaa                                                                                                                            <  2:204358/140‑1 (MQ=255)        tttCCGGTTTTTTCGTTTCGTTTTCCGGTGCTTTCCGCAATCGCCATCCGCTTATCCACAACGCCGCACCTTTTTCGTGGGGAGCTGTAGTAGCAATTATAGTCGATTAATACAACATATTGGGTTGGGGCGCATTTTaa                                                                                                                            <  2:573035/140‑1 (MQ=255)          tCCGGTTTTTTCGTTTCGTTTTCCGGTGCTTTCCGCAATCGCCATCCGCTTATCCACAACGCCGCACCTTTTTCGTGGGGAGCTGTAGTAGCAATTATAGTCGATTAATACAACATATTGGGTTGGGGCGCATTTTAAAg                                                                                                                          <  2:501709/140‑1 (MQ=255)          tCCGGTTTTTTCGTTTCGTTTTCCGGTGCTTTCCGCAATCGCCATCCGCTTATCCACAACGCCGCACCTTTTTCGTGGGGAGCTGTAGTAGCAATTATAGTCGATTAATACAACATATTGGGTTGGGGCGCATTTTAAAg                                                                                                                          <  1:82164/140‑1 (MQ=255)          tCCGGTTTTTTCGTTTCGTTTTCCGGTGCTTTCCGCAATCGCCATCCGCTTATCCACAACGCCGCACCTTTTTCGTGGGGAGCTGTAGTAGCAATTATAGTCGATTAATACAACATATTGGGTTGGGGCGCATTTTAAAg                                                                                                                          <  2:596081/140‑1 (MQ=255)          tCCGGTTTTTTCGTTTCGTTTTCCGGTGCTTTCCGCAATCGCCATCCGCTTATCCACAACGCCGCACCTTTTTCGTGGGGAGCTGTAGTAGCAATTATAGTCGATTAATACAACATATTGGGTTGGGGCGCATTTTAAAg                                                                                                                          <  1:256860/140‑1 (MQ=255)          tCCGGTTTTTTCGTTTCGTTTTCCGGTGCTTTCCGCAATCGCCATCCGCTTATCCACAACGCCGCACCTTTTTCGTGGGGAGCTGTAGTAGCAATTATAGTCGATTAATACAACATATTGGGTTGGGGCGCATTTTAAAg                                                                                                                          <  1:312343/140‑1 (MQ=255)           ccGGTTTTTTCGTTTCGTTTTCCGGTGCTTTCCGCAATCGCCATCCGCTTATCCACAACGCCGCACCTTTTTCGTGGGGAGCTGTAGTAGC‑ATTATAGTCGATTAATACAACATATTGGGTTGGGGCGCATTTTAAAGTc                                                                                                                        <  2:562806/140‑1 (MQ=255)             ggTTTTTTCGTTTCGTTTTCCGGTGCTTTCCGCAATCGCCATCCGCTTATCCACAACGCCGCACCTTTTTCGTGGGGAGCTGTAGTAGCAATTATAGTCGATTAATACAACATATTGGGTTGGGGCGCATTTTAAAGTCt                                                                                                                       <  2:136882/140‑1 (MQ=255)                                    tGCTTTCCGCAATCGCCATCCGCTTATCCACAACGCCGCACCTTTTTCGTGGGGAGCTGTAGTAGCAATTATAGTCGATTAATACAACATATTGGGTTGGGGCGCATTTTAAAGTCTATATATAGTGCTTTGCATCAAgg                                                                                                >  1:644727/1‑140 (MQ=255)                                     gCTTTCCGCAATCGCCATCCGCTTATCCACAACGCCGCACCTTTTTCGTGGGGAGCTGTAGTAGCAATTATAGTTGATTAATACAACATATTGGGTTGGGGCGCATTTTAAAGTCTATATATAGTGCTTTGCATCAAGGa                                                                                               >  1:651604/1‑140 (MQ=255)                                         tCCGCAATCGCCATCCGCTTATCCACAACGCCGCACCTTTTTCGTGGGGAGCTGTAGTAGCAATTATAGTCGATTAATACAACATATTGGGTTGGGGCGCATTTTAAAGTCTATATATAGTGCTTTGCATCAAGGATGtt                                                                                           >  2:427925/1‑140 (MQ=255)                                           cGCAATCGCCATCCGCTTATCCACAACGCCGCACCTTTTTCGTGGGGAGCTGTAGTAGCAATTATAGTCGATTAATACAACATATTGGGTTGGGGCGCATTTTAAAGTCTATATATAGTGCTTTGCATCAAGGATGTTTg                                                                                         >  2:498349/1‑140 (MQ=255)                                             cAATCGCCATCCGCTTATCCACAACGCCGCACCTTTTTCGTGGGGAGCTGTAGTAGCAATTATAGTCGATTAATACAACATATTGGGTTGGGGCGCATTTTAAAGTCTATATATAGTGCTTTGCATCAAGGATGTTTGAg                                                                                       <  2:644727/140‑1 (MQ=255)                                              aaTCGCCATCCGCTTATCCACAACGCCGCACCTTTTTCGTGGGGAGCTGTAGTAGCAATTATAGTTGATTAATACAACATATTGGGTTGGGGCGCATTTTAAAGTCTATATATAGTGCTTTGCATCAAGGATGTTTGAGc                                                                                      <  2:651604/140‑1 (MQ=255)                                                                cacaACGCCGCACCTTTTTCGTGGGGAGCTGTAGTAGCAATTATAGTCGATTAATACAACATATTGGGTTGGGGCGCATTTTAAAGTCTATATATAGTGCTTTGCATCAAGGATGTTTGAGCTTTTTTTGATGTAGCTCa                                                                    <  1:210482/140‑1 (MQ=255)                                                                cacaACGCCGCACCTTTTTCGTGGGGAGCTGTAGTAGCAATTATAGTCGATTAATACAACATATTGGGTTGGGGCGCATTTTAAAGTCTATATATAGTGCTTTGCATCAAGGATGTTTGAGCTTTTTTTGATGTAGCTCa                                                                    <  1:364463/140‑1 (MQ=255)                                                                cacaACGCCGCACCTTTTTCGTGGGGAGCTGTAGTAGCAATTATAGTCGATTAATACAACATATTGGGTTGGGGCGCATTTTAAAGTCTATATATAGTGCTTTGCATCAAGGATGTTTGAGCTTTTTTTGATGTAGCTCa                                                                    <  2:570443/140‑1 (MQ=255)                                                                cacaACGCCGCACCTTTTTCGTGGGGAGCTGTAGTAGCAATTATAGTCGATTAATACAACATATTGGGTTGGGGCGCATTTTAAAGTCTATATATAGTGCTTTGCATCAAGGATGTTTGAGCTTTTTTTGATGTAGCTCa                                                                    <  1:498349/140‑1 (MQ=255)                                                                  caACGCCGCACCTTTTTCGTGGGGAGCTGTAGTAGCAATTATAGTCGATTAATACAACATATTGGGTTGGGGCGCATTTTAAAGTCTATATATAGTGCTTTGCATCAAGGATGTTTGAGCTTTTTTTGATGTAGCTCaaa                                                                  >  1:246325/1‑140 (MQ=255)                                                                   aaCGCCGCACCTTTTTCGTGGGGAGCTGTAGTAGCAATTATAGTCGATTAATACAACATATTGGGTTGGGGCGCATTTTAAAGTCTATATATAGTGCTTTGCATCAAGGATGTTTGAGCTTTTTTTGATGTAGCTCAAAg                                                                 <  2:110547/140‑1 (MQ=255)                                                                            ccTTTTTCGTGGGGAGCTGTAGTAGCAATTATAGTCGATTAATACAACATATTGGGTTGGGGCGCATTTTAAAGTCTATATATAGTGCTTTGCATCAAGGATGTTTGAGCTTTTTTTGATGTAGCTCAAAGTAAAAAGCa                                                        <  1:427925/140‑1 (MQ=255)                                                                            ccTTTTTCGTGGGGAGCTGTAGTAGCAATTATAGTCGATTAATACAACATATTGGGTTGGGGCGCATTTTAAAGTCTATATATAGTGCTTTGCATCAAGGATGTTTGAGCTTTTTTTGATGTAGCTCAAAGTAAAAAGCa                                                        <  1:635120/140‑1 (MQ=255)                                                                            ccTTTTTCGTGGGGAGCTGTAGTAGCAATTATAGTCGATTAATACAACATATTGGGTTGGGGCGCATTTTAAAGTCTATATATAGTGCTTTGCATCAAGGATGTTTGAGCTTTTTTTGATGTAGCTCAAAGTAAAAAGCa                                                        <  1:525647/140‑1 (MQ=255)                                                                                            cTGTAGTAGCAATTATAGTCGATTAATACAACATATTGGGTTGGGGCGCATTTTAAAGTCTATATATAGTGCTTTGCATCAAGGATGTTTGAGCTTTTTTTGATGTAGCTCAAAGTAAAAAGCAGAGCGTACGGATGACg                                        >  2:616286/1‑140 (MQ=255)                                                                                              gtagtaGCAATTATAGTCGATTAATACAACATATTGGGTTGGGGCGCATTTTAAAGTCTATATATAGTGCTTTGCATCAAGGATGTTTGAGCTTTTTTTGATGTAGCTCAAAGTAAAAAGCAGAGCGTACGGATGACggg                                      <  1:653737/140‑1 (MQ=255)                                                                                              gtagtaGCAATTATAGTCGATTAATACAACATATTGGGTTGGGGCGCATTTTAAAGTCTATATATAGTGCTTTGCATCAAGGATGTTTGAGCTTTTTTTGATGTAGCTCAAAGTAAAAAGCAGAGCGTACGGATGACggg                                      <  1:420901/140‑1 (MQ=255)                                                                                              gtagtaGCAATTATAGTCGATTAATACAACATATTGGGTTGGGGCGCATTTTAAAGTCTATATATAGTGCTTTGCATCAAGGATGTTTGAGCTTTTTTTGATGTAGCTCAAAGTAAAAAGCAGAGCGTACGGATGACggg                                      <  1:241494/140‑1 (MQ=255)                                                                                                         tataGTCGATTAATACAACATATTGGGTTGGGGCGCATTTTAAAGTCTATATATAGTGCTTTGCATCAAGGATGTTTGAGCTTTTTTTGATGTAGCTCAAAGTAAAAAGCAGAGCGTACGGATGACGGGCGCTACAGCGa                           >  1:279396/1‑140 (MQ=255)                                                                                                           taGTCGATTAATACAACATATTGGGTTGGGGCGCATTTTAAAGTCTATATATAGTGCTTTGCATCAAGGATGTTTGAGCTTTTTTTGATGTAGCTCAAAGTAAAAAGCAGAGCGTACGGATGACGGGCGCTACAGCGata                         >  2:427108/1‑140 (MQ=255)                                                                                                           taGTCGATTAATACAACATATTGGGTTGGGGCGCATTTTAAAGTCTATATATAGTGCTTTGCATCAAGGATGTTTGAGCTTTTTTTGATGTAGCTCAAAGTAAAAAGCAGAGCGTACGGATGACGGGCGCTACAGCGata                         >  1:12659/1‑140 (MQ=255)                                                                                                                          acaTATTGGGTTGGGGCGCATTTTAAAGACTATATATAGTGCTTTGCATCAAGGATGTTTGAGCTTTTTTTGATGTAGCTCAAAGTAAAAAGCAGAGCGTACGGATGACGGGCGCTACAGCGATATGTAAATTTTTTAAt          >  2:206105/1‑140 (MQ=255)                                                                                                                              atTGGGTTGGGGCGCATTTTAAAGTCTATATATAGTGCTTTGCATCAAGGATGTTTGAGCTTTTTTTGATGTAGCTCAAAGTAAAAAGCAGAGCGTACGGATGACGGGCGCTACAGCGATATGTAAATTTTTTAATGAAt      >  1:393488/1‑140 (MQ=255)                                                                                                                                tgggttgggGCGCATTTTAAAGTCTATATATAGTGCTTTGCATCAAGGATGTTTGAGCTTTTTTTGATGTAGCTCAAAGTAAAAAGCAGAGCGTACGGATGACGGGCGCTACAGCGATATGTAAATTTTTTAATGAAttt    >  1:111209/1‑140 (MQ=255)                                                                                                                                  ggttgggGCGCATTTTAAAGTCTATATATAGTGCTTTGCATCAAGGATGTTTGAGCTTTTTTTGATGTAGCTCAAAGTAAAAAGCAGAGCGTACGGATGACGGGCGCTACAGCGATATGTAAATTTTTTAATGAATTTGc  >  2:185902/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                     AAGGCGTTTTCCGGTTTTTTCGTTTCGTTTTCCGGTGCTTTCCGCAATCGCCATCCGCTTATCCACAACGCCGCACCTTTTTCGTGGGGAGCTGTAGTAGCAATTATAGTCGATTAATACAACATATTGGGTTGGGACGCATTTTAAAGTCTATATATAGTGCTTTGCATCAAGGATGTTTGAGCTTTTTTTGATGTAGCTCAAAGTAAAAAGCAGAGCGTACGGATGACGGGCGCTACAGCGATATGTAAATTTTTTAATGAATTTGC  >  NC_000913/4462706‑4462974

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑