Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,474,138	(A)6→7	coding (15/594 nt)	bdcR →	transcriptional repressor for divergent bdcA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,474,132	1	.	A	100.0%	139.1 / NA	43	T3T (ACT→ACA)	bdcR	transcriptional repressor for divergent bdcA
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base A (30/13);  total (30/13)

GCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACT‑AAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGG  >  NC_000913/4473995‑4474263                                                                                                                                           |                                                                                                                                    gCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACT‑a                                                                                                                                    >  2:630264/1‑139 (MQ=255)                   ccACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGtt                                                                                                                   >  2:235103/1‑139 (MQ=255)                   ccACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGtt                                                                                                                   >  1:702543/1‑139 (MQ=255)                                      tCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGtt                                                                                                >  2:269217/1‑139 (MQ=255)                                       cTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGttt                                                                                               <  1:300875/139‑1 (MQ=255)                                           ccTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCt                                                                                           >  2:672181/1‑139 (MQ=255)                                              gTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCt                                                                                        >  2:677441/1‑139 (MQ=255)                                                aaaaGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGa                                                                                      >  1:175761/1‑139 (MQ=255)                                                 aaaGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAg                                                                                     <  2:277407/139‑1 (MQ=255)                                                      gcCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGc                                                                                >  1:676119/1‑139 (MQ=255)                                                       cccATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAAtctc                                                                                                                        <  1:304076/98‑1 (MQ=255)                                                       cccATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAAtctc                                                                                                                        >  2:304076/1‑98 (MQ=255)                                                       cccATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGgt                                                                                                             >  1:149819/1‑109 (MQ=255)                                                       cccATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGgt                                                                                                             <  2:149819/109‑1 (MQ=255)                                                                                  tAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCAc                                                    >  2:330282/1‑139 (MQ=255)                                                                                  tAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCAc                                                    >  1:602500/1‑139 (MQ=255)                                                                                   aaCATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCAcc                                                   >  1:165011/1‑139 (MQ=255)                                                                                     cATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTg                                                                                       <  1:218569/101‑1 (MQ=255)                                                                                     cATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTg                                                                                       >  2:218569/1‑101 (MQ=255)                                                                                       tGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCaaaa                                               >  2:272237/1‑139 (MQ=255)                                                                                       tGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCaaaa                                               >  2:177651/1‑139 (MQ=255)                                                                                       tGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCaaaa                                               >  2:39594/1‑139 (MQ=255)                                                                                       tGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCaaaa                                               >  2:398553/1‑139 (MQ=255)                                                                                       tGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAAAAATCTCGCGTTCCAGGTCGTCCAAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCaaaa                                               >  1:574933/1‑139 (MQ=255)                                                                                          aTTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAACAGTGCTTTTTCACCAAAAAgg                                            >  2:470254/1‑139 (MQ=255)                                                                                          aTTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAgg                                            >  2:449479/1‑139 (MQ=255)                                                                                          aTTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAgg                                            >  1:643256/1‑139 (MQ=255)                                                                                          aTTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAgg                                            >  2:375629/1‑139 (MQ=255)                                                                                            tACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGtt                                          >  2:397441/1‑139 (MQ=255)                                                                                            tACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGtt                                          >  2:780303/1‑139 (MQ=255)                                                                                            tACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGtt                                          >  2:653584/1‑139 (MQ=255)                                                                                               aaaaCTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCg                                       <  2:165011/139‑1 (MQ=255)                                                                                                aaaCTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGa                                      <  1:272237/139‑1 (MQ=255)                                                                                                   cTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGc                                   >  1:204601/1‑139 (MQ=255)                                                                                                   cTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGc                                   >  2:103366/1‑139 (MQ=255)                                                                                                      ttttAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGt                                >  1:477616/1‑139 (MQ=255)                                                                                                      ttttAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGt                                <  1:447804/139‑1 (MQ=255)                                                                                                      ttttAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGt                                >  2:275167/1‑139 (MQ=255)                                                                                                                 tGTCAACAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTg                     <  1:672181/139‑1 (MQ=255)                                                                                                                     aaCAGGTTTTGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGt                 <  1:756627/139‑1 (MQ=255)                                                                                                                            tttGTAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGAt          <  1:375629/139‑1 (MQ=255)                                                                                                                                tAATGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATc      <  2:700158/139‑1 (MQ=255)                                                                                                                                 aaTGGTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCt     <  1:39594/139‑1 (MQ=255)                                                                                                                                    ggTCACTAAAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTgg  >  2:559519/1‑139 (MQ=255)                                                                                                                                           |                                                                                                                                    GCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACT‑AAAAAACAATCTCGCGTTCCAGGTCGTCCCAGACGTTTCGCTCCTGAGCAGGCAATCTCTGCGGCAAAAGTGCTTTTTCACCAAAAAGGTTTCGATGCTGTCAGTGTTGCTGAAGTTACTGATTATCTTGG  >  NC_000913/4473995‑4474263

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