Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,026,764	C→T	R80C (CGC→TGC)	fre →	flavin reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,026,764	0	C	T	100.0%	126.1 / NA	36	R80C (CGC→TGC)	fre	flavin reductase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (18/18);  total (18/18)

TCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATT‑CCCCACGGAGAAGCGTGGCTGCGCGATGATGAAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTA  >  NC_000913/4026628‑4026876                                                                                                                                         |                                                                                                                  tCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACTGCa                                                                                                                <  2:240921/140‑1 (MQ=255) tCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACTGCa                                                                                                                <  1:573057/140‑1 (MQ=255) tCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACTGCa                                                                                                                <  2:369176/140‑1 (MQ=255)          gATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACTGCATCCTCAAAg                                                                                                       >  2:342724/1‑140 (MQ=255)          gATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACTGCATCCTCAAAg                                                                                                       <  1:342724/140‑1 (MQ=255)               tAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACTGCATCCTCAAAGatcat                                                                                                  >  2:452360/1‑140 (MQ=255)               tAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACTGCATCCTCAAAGatcat                                                                                                  >  1:249499/1‑140 (MQ=255)                  tGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACTGCATCCTCAAAGATCATCaa                                                                                               >  1:225309/1‑140 (MQ=255)                  tGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACTGCATCCTCAAAGATCATCaa                                                                                               >  1:48289/1‑140 (MQ=255)                     tggatgAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACTGCATCCTCAAAGATCATCAAATc                                                                                            <  1:352850/140‑1 (MQ=255)                        atgAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACTGCATCCTCAAAGATCATCAAATCGTg                                                                                         >  1:13234/1‑140 (MQ=255)                          gAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACTGCATCCTCAAAGATCATCAAATCGTGGt                                                                                       <  2:73707/140‑1 (MQ=255)                           aGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACTGCATCCTCAAAGATCATCAAATCGTGGTc                                                                                      >  1:282601/1‑140 (MQ=255)                           aGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACTGCATCCTCAAAGATCATCAAATCGTGGTc                                                                                      >  2:566970/1‑140 (MQ=255)                               cgACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACTGCATCCTCAAAGATCATCAAATCGTGGTCGACa                                                                                  <  2:226063/140‑1 (MQ=255)                                gACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACTGCATCCTCAAAGATCATCAAATCGTGGTCGACAt                                                                                 <  2:373252/140‑1 (MQ=255)                                 aCAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACTGCATCCTCAAAGATCATCAAATCGTGGTCGACAtt                                                                                <  1:148647/140‑1 (MQ=255)                                        tCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACTGCATCCTCAAAGATCATCAAATCGTGGTCGACATT‑CCCCACg                                                                        >  1:439794/1‑140 (MQ=255)                                                      cTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACTGCATCCTCAAAGATCATCAAATCGTGGTCGACATT‑CCCCACGGAGAAGCGTGGCTg                                                          <  2:282601/140‑1 (MQ=255)                                                        tCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACTGCATCCTCAAAGATCATCAAATCGTGGTCGACATT‑CCCCACGGAGAAGCGTGGCTgcg                                                        <  2:249499/140‑1 (MQ=255)                                                         cGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACTGCATCCTCAAAGATCATCAAATCGTGGTCGACATT‑CCCCACGGAGAAGCGTGGCTgcgc                                                       >  1:612590/1‑140 (MQ=255)                                                          gACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACTGCATCCTCAAAGATCATCAAATCGTGGTCGACATT‑CCCCACGGAGAAGCGTGGCTGcgcg                                                      <  2:13234/140‑1 (MQ=255)                                                          gACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACTGCATCCTCAAAGATCATCAAATCGTGGTCGACATT‑CCCCACGGAGAAGCGTGGCTGcgcg                                                      <  1:117792/140‑1 (MQ=255)                                                           aCGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACTGCATCCTCAAAGATCATCAAATCGTGGTCGACATT‑CCCCACGGAGAAGCGTGGCTgcg                                                        <  2:122864/137‑1 (MQ=255)                                                           aCGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACTGCATCCTCAAAGATCATCAAATCGTGGTCGACATT‑CCCCACGGAGAAGCGTGGCTgcg                                                        >  1:122864/1‑137 (MQ=255)                                                                   tGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACTGCATCCTCAAAGATCATCAAATCGTGGTCGACATT‑CCCCACGGAGAAGCGTGGCTGCGCGATGATGAag                                             >  1:574913/1‑140 (MQ=255)                                                                        aaGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACTGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCCACGGAGAAGCGTGGCTGCGCGATGATGAAGAGCg                                         <  2:48289/140‑1 (MQ=255)                                                                                         cATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACTGCATCCTCAAAGATCATCAAATCGTGGTCGACATT‑CCCCACGGAGAAGCGTGGCTGCGCGATGATGAAGAGCGTCCgat                                   >  2:356178/1‑128 (MQ=255)                                                                                         cATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACTGCATCCTCAAAGATCATCAAATCGTGGTCGACATT‑CCCCACGGAGAAGCGTGGCTGCGCGATGATGAAGAGCGTCCgat                                   <  1:356178/128‑1 (MQ=255)                                                                                               ggCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACTGCATCCTCAAAGATCATCAAATCGTGGTCGACATT‑CCCCACGGAGAAGCGTGGCTGCGCGATGATGAAGAGCGTCCGATGATTTTGATTgcgggcgg                 >  2:312909/1‑140 (MQ=255)                                                                                                   cTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACTGCATCCTCAAAGATCATCAAATCGTGGTCGACATT‑CCCCACGGAGAAGCGTGACTGCGCGATGATGAAGAGCGTCCGATGATTTTGATTGCGGGCGGCAcc             >  1:199869/1‑140 (MQ=255)                                                                                                     tCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACTGCATCCTCAAAGATCATCAAATCGTGGTCGACATT‑CCCCACGGAGAAGCGTGGCTGCGCGATGATGAAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCgg           >  2:140246/1‑140 (MQ=255)                                                                                                       tGAAATCAACCTTTACGCGAAAGCAGTCATGGACTGCATCCTCAAAGATCATCAAATCGTGGTCGACATT‑CCCCACGGAGAAGCGTGGCTGCGCGATGATGAAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGt         <  1:605974/140‑1 (MQ=255)                                                                                                       tGAAATCAACCTTTACGCGAAAGCAGTCATGGACTGCATCCTCAAAGATCATCAAATCGTGGTCGACATT‑CCCCACGGAGAAGCGTGGCTGCGCGATGATGAAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGt         <  1:452360/140‑1 (MQ=255)                                                                                                              aaCCTTTACGCGAAAGCAGTCATGGACTGCATCCTCAAAGATCATCAAATCGTGGTCGACATT‑CCCCACGGAGAAGCGTGGCTGCGCGATGATGAAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTa  >  2:320698/1‑140 (MQ=255)                                                                                                              aaCCTTTACGCGAAAGCAGTCATGGACTGCATCCTCAAAGATCATCAAATCGTGGTCGACATT‑CCCCACGGAGAAGCGTGGCTGCGCGATGATGAAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTa  >  1:182737/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                  TCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATT‑CCCCACGGAGAAGCGTGGCTGCGCGATGATGAAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTA  >  NC_000913/4026628‑4026876

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: ‑