Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,333,655	A→T	L99Q (CTG→CAG)	basR ←	DNA‑binding transcriptional dual regulator BasR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,333,655	0	A	T	100.0%	106.8 / NA	30	L99Q (CTG→CAG)	basR	DNA‑binding transcriptional dual regulator BasR
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (18/12);  total (18/12)

TCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCATGTAACTCTTCCAGCGCAAAAGGCTTCACCAGATAGTCGTCGGCACCGACATCCAGCCCGGCGATTTTGTCGGTCAGCGTATCGCGAGCGGTGAGGATCAGTACCGG  >  NC_000913/4333580‑4333731                                                                            |                                                                             tCGCCCTGATTATGATGGCGTCGTAGCAGGGCGCGGATACGGGCATGTAACTCTTCCAGCGCAAAAGGCTTCACCTGAt                                                                           <  2:866214/79‑1 (MQ=255)   gCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCATGTAACTCTTCCAGCGCAAAAGGCTTCACCTGATAgtc                                                                       >  1:866263/1‑81 (MQ=255)    cccTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCATGTAACTCTTCCAGCGCAAAAGGCTTCACCTGATAgtcg                                                                      >  2:866269/1‑81 (MQ=255)             ttatGGCGTCGTAGCAGGGCGCGGATACGGGCATGTAACTCTTCCAGCGCAAAAGGCTTCACCTGATAGTCGTCGGCACCg                                                             <  1:866215/81‑1 (MQ=255)               atGGCGTCGTAGCAGGGCGCGGATACGGGCATGTAACTCTTCCAGCGCAAAAGGCTTCACCTGATAGTCGTCGGCACCGa                                                            >  2:866267/1‑80 (MQ=255)               atGGCGTCGTAGCAGGGCGCGGATACGGGCATGTAACTCTTCCAGCGCAAAAGGCTTCACCTGATAGTCGTCGGCACCGa                                                            >  1:866259/1‑80 (MQ=255)                tGGCGTCGTAGCAGGGCGCGGATACGGGCATGTAACTCTTCCAGCGCAAAAGGCTTCACCTGATAGTCGTCGGCACCGACa                                                          >  2:866272/1‑81 (MQ=255)                  gCGTCGTAGCAGGGCGCGGATACGGGCATGTAACTCTTCCAGCGCAAAAGGCTTCACCTGATAGTCGTCGGCACCGACATc                                                        >  2:866275/1‑81 (MQ=255)                       gtAGCAGGGCGCGGATACGGGCATGTAACTCTTCCAGCGCAAAAGGCTTCACCTGATAGTCGTCGGCACCGACATCCAGcc                                                   >  1:866280/1‑81 (MQ=255)                        tAGCAGGGCGCGGATACGGGCATGTAACTCTTCCAGCGCAAAAGGCTTCACCTGATAGTCGTCGGCACCGACATCCAGccc                                                  <  2:866216/81‑1 (MQ=255)                          gCAGGGCGCGGATACGGGCATGTAACTCTTCCAGCGCAAAAGGCTTCACCTGATAGTCGTCGGCACCGACATCCAGCCCgg                                                >  2:866268/1‑81 (MQ=255)                          gCAGGGCGCGGATACGGGCATGTAACTCTTCCAGCGCAAAAGGCTTCACCTGATAGTCGTCGGCACCGACATCCAGCCCg                                                 <  1:866218/80‑1 (MQ=255)                          gCAGGGCGCGGATACGGGCATGTAACTCTTCCAGCGCAAAAGGCTTCACCTGATAGTCGTCGGCACCGACATCCAGCCCg                                                 <  1:866217/80‑1 (MQ=255)                           cAGGGCGCGGATACGGGCATGTAACTCTTCCAGCGCAAAAGGCTTCACCTGATAGTCGTCGGCACCGACATCCAGCCCGGc                                               >  2:866279/1‑81 (MQ=255)                           cAGGGCGCGGATACGGGCATGTAACTCTTCCAGCGCAAAAGGCTTCACCTGATAGTCGTCGGCACCGACATCCAGCCCGGc                                               >  2:866270/1‑81 (MQ=255)                           cAGGGCGCGGATACGGGCATGTAACTCTTCCAGCGCAAAAGGCTTCACCTGATAGTCGTCGGCACCGACATCCAGCCCGGc                                               <  1:866219/81‑1 (MQ=255)                                cgcgGATACGGGCATGTAACTCTTCCAGCGCAAAAGGCTTCACCTGATAGTCGTCGGCACCGACATCCAGCCCGGCGAttt                                          >  2:866276/1‑81 (MQ=255)                                  cgGATACGGGCATGTAACTCTTCCAGCGCAAAAGGCTTCACCTGATAGTCGTCGGCACCGACATCCAGCCCGGCGATTTTg                                        >  1:866278/1‑81 (MQ=255)                                    gATACGGGCATGTAACTCTTCCAGCGCAAAAGGCTTCACCTGATAGTCGTCGGCACCGACATCCAGCCCGGCGATTTTGTc                                      >  2:866277/1‑81 (MQ=255)                                     aTACGGGCATGTAACTCTTCCAGCGCAAAAGGCTTCACCTGATAGTCGTCGGCACCGACATCCAGCCCGGCGATTTTGTCg                                     >  2:866273/1‑81 (MQ=255)                                      tACGGGCATGTAACTCTTCCAGCGCAAAAGGCTTCACCTGATAGTCGTCGGCACCGACATCCAGCCCGGCGATTTTGTCgg                                    >  2:866271/1‑81 (MQ=255)                                        cGGGCATGTAACTCTTCCAGCGCAAAAGGCTTCACCTGATAGTCGTCGGCACCGACATCCAGCCCGGCGATTTTGTCGGTc                                  <  2:866220/81‑1 (MQ=255)                                                   ctctTCCAGCGCAAAAGGCTTCACCTGATAGTCGTCGGCACCGACATCCAGCCCGGCGATTTTGTCGGTCAGCGTATcgcg                       <  2:866221/81‑1 (MQ=255)                                                   ctctTCCAGCGCAAAAGGCTTCACCTGATAGTCGTCGGCACCGACATCCAGCCCGGCGATTTTGTCGGTCAGCGTATcgc                        >  2:866266/1‑80 (MQ=255)                                                        ccAGCGCAAAAGGCTTCACCTGATAGTCGTCGGCACCGACATCCAGCCCGGCGATTTTGTCGGTCAGCGTATCGCGAg                     <  2:866222/78‑1 (MQ=255)                                                            cgcAAAAGGCTTCACCTGATAGTCGTCGGCACCGACATCCAGCCCGGCGATTTTGTCGGTCAGCGTATCGCGAGCGGTGAg              <  2:866223/81‑1 (MQ=255)                                                             gcAAAAGGCTTCACCTGATAGTCGTCGGCACCGACATCCAGCCCGGCGATTTTGTCGGTCAGCGTATCGCGAGCGGTGAgg             <  2:866224/81‑1 (MQ=255)                                                               aaaaGGCTTCACCTGATAGTCGTCGGCACCGACATCCAGCCCGGCGATTTTGTCGGTCAGCGTATCGCGAGCGGTGAGGAt           >  2:866288/1‑81 (MQ=255)                                                                 aaGGCTTCACCTGATAGTCGTCGGCACCGACATCCAGCCCGGCGATTTTGTCGGTCAGCGTATCGCGAGCGGTGAGGATCa         <  2:866225/81‑1 (MQ=255)                                                                        cACCTGATAGTCGTCGGCACCGACATCCAGCCCGGCGATTTTGTCGGTCAGCGTATCGCGAGCGGTGAGGATCAGTACCgg  >  2:866274/1‑81 (MQ=255)                                                                            |                                                                             TCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCATGTAACTCTTCCAGCGCAAAAGGCTTCACCAGATAGTCGTCGGCACCGACATCCAGCCCGGCGATTTTGTCGGTCAGCGTATCGCGAGCGGTGAGGATCAGTACCGG  >  NC_000913/4333580‑4333731

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