Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,648,011	T→C	L26L (CTT→CTC)	dicA →	DNA‑binding transcriptional dual regulator DicA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,648,011	0	T	C	100.0%	164.3 / NA	45	L26L (CTT→CTC)	dicA	DNA‑binding transcriptional dual regulator DicA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (20/25);  total (20/25)

AACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTTGCTAAAGCCCTGAAAATCTCCCATGTGTCTGTATCACAATGGGAACGGGGTGATAGTGAACCTACAGGGAAGAACCTTTTTGCCCTCAGTAAAGTATTGCAATGCTCACCAACATGGATTCTATTTGGCGA  >  NC_000913/1647876‑1648142                                                                                                                                        |                                                                                                                                    aactaacAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTCGCTa                                                                                                                                 >  1:22694/1‑140 (MQ=255)      acAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTCGCTAAAGcc                                                                                                                            >  2:274306/1‑140 (MQ=255)           atagatGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTCGCTAAAGCCCTGaa                                                                                                                       >  2:246151/1‑140 (MQ=255)           atagatGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTCGCTAAAGCCCTGaa                                                                                                                       >  2:122182/1‑140 (MQ=255)           atagatGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTCGCTAAAGCCCTGaa                                                                                                                       >  1:346475/1‑140 (MQ=255)              gatGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTCGCTAAAGCCCTGAAAAt                                                                                                                    >  2:47675/1‑140 (MQ=255)               atGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTCGCTAAAGCCCTGAAAAtc                                                                                                                   <  2:346475/140‑1 (MQ=255)                  ttAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTCGCTAAAGCCCTGAAAATCTcc                                                                                                                >  2:275510/1‑140 (MQ=255)                    aaCAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTCGCTAAAGCCCTGAAAATCTCCCa                                                                                                              >  2:66248/1‑140 (MQ=255)                       aaaaaCATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTCGCTAAAGCCCTGAAAATCTCCCAtgt                                                                                                           <  2:378881/140‑1 (MQ=255)                       aaaaaCATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTCGCTAAAGCCCTGAAAATCTCCCAtgt                                                                                                           <  2:298410/140‑1 (MQ=255)                       aaaaaCATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTCGCTAAAGCCCTGAAAATCTCCCAtgt                                                                                                           <  1:246151/140‑1 (MQ=255)                       aaaaaCATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTCGCTAAAGCCCTGAAAATCTCCCAtgt                                                                                                           <  1:477516/140‑1 (MQ=255)                           aCATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTCGCTAAAGCCCTGAAAATCTCCCATGTGTCt                                                                                                       <  1:52086/140‑1 (MQ=255)                               aGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTCGCTAAAGCCCTGAAAATCTCCCATGTGTCTGTAt                                                                                                   <  1:473058/140‑1 (MQ=255)                                gTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTCGCTAAAGCCCTGAAAATCTCCCATGTGTCTGTATc                                                                                                  <  2:40429/140‑1 (MQ=255)                                    aTACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTCGCTAAAGCCCTGAAAATCTCCCATGTGTCTGTATCACaa                                                                                              >  2:533185/1‑140 (MQ=255)                                    aTACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTCGCTAAAGCCCTGAAAATCTCCCATGTGTCTGTATCACaa                                                                                              <  1:270208/140‑1 (MQ=255)                                     tACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTCGCTAAAGCCCTGAAAATCTCCCATGTGTCTGTAt                                                                                                   >  1:349214/1‑134 (MQ=255)                                     tACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTCGCTAAAGCCCTGAAAATCTCCCATGTGTCTGTAt                                                                                                   <  2:349214/134‑1 (MQ=255)                                        gATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTCGCTAAAGCCCTGAAAATCTCCCATGTGTCTGTATCACAATggg                                                                                          <  1:533185/140‑1 (MQ=255)                                          tttAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTCGCTAAAGCCCTGAAAATCTCCCATGTGTCTGTATCACAATGGGaa                                                                                        <  1:66248/140‑1 (MQ=255)                                             aGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTCGCTAAAGCCCTGAAAATCTCCCATGTGTCTGTATCACAATGGGAACgg                                                                                     >  1:116826/1‑140 (MQ=255)                                                          tATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTCGCTAAAGCCCTGAAAATCTCCCATGTGTCTGTATCACAATGGGAACGGGGTGATAGTGAAc                                                                        <  2:22694/140‑1 (MQ=255)                                                            tGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTCGCTAAAGCCCTGAAAATCTCCCATGTGTCTGTATCACAATGGGaa                                                                                        <  1:154273/122‑1 (MQ=255)                                                            tGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTCGCTAAAGCCCTGAAAATCTCCCATGTGTCTGTATCACAATGGGaa                                                                                        >  2:154273/1‑122 (MQ=255)                                                            tGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTCGCTAAAGCCCTGAAAATCTCCCATGTGTCTGTATCACAATGGGAACGGGGTGATAGTGAACCt                                                                      >  1:18954/1‑140 (MQ=255)                                                             ggAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTCGCTAAAGCCCTGAAAATCTCCCATGTGTCTGTATCACAATGGGAACGGGGTGATAGTGAACCTa                                                                     <  1:274306/140‑1 (MQ=255)                                                              gAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTCGCTAAAGCCCTGAAAATCTCCCATGTGTCTGTATCACAATGGGAACGGGGTGATAGTGAACCTAc                                                                    >  2:223215/1‑140 (MQ=255)                                                                                 tCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTCGCTAAAGCCCTGAAAATCTCCCATGTGTCTGTATCACAATGGGAACGGGGTGATAGTGAACCTACAGGGAAGAACCTTTTTGcc                                                 <  1:28374/140‑1 (MQ=255)                                                                                 tCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTCGCTAAAGCCCTGAAAATCTCCCATGTGTCTGTATCACAATGGGAACGGGGTGATAGTGAACCTACAGGGAAGAACCTTTTTGcc                                                 <  1:528345/140‑1 (MQ=255)                                                                                 tCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTCGCTAAAGCCCTGAAAATCTCCCATGTGTCTGTATCACAATGGGAACGGGGTGATAGTGAACCTACAGGGAAGAACCTTTTTGcc                                                 <  1:223215/140‑1 (MQ=255)                                                                                 tCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTCGCTAAAGCCCTGAAAATCTCCCATGTGTCTGTATCACAATGGGAACGGGGTGATAGTGAACCTACAGGGAAGAACCTTTTTGcc                                                 >  2:28374/1‑140 (MQ=255)                                                                                 tCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTCGCTAAAGCCCTGAAAATCTCCCATGTGTCTGTATCACAATGGGAACGGGGTGATAGTGAACCTACAGGGAAGAACCTTTTTGcc                                                 <  2:516136/140‑1 (MQ=255)                                                                                 tCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGCTCTCTCGCTAAAGCCCTGAAAATCTCCCATGTGTGTGTATCACAATGGGAACGGGGTGATAGTGAACCTTCACGGAAGAACCTTTTTGcc                                                 <  2:401709/140‑1 (MQ=255)                                                                                           cATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTCGCTAAAGCCCTGAAAATCTCCCATGTGTCTGTATCACAATGGGAACGGGGTGATAGTGAACCTACAGGGAAGAACCTTTTTGCCCTCAGTAAAg                                       >  1:153424/1‑140 (MQ=255)                                                                                               aGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTCGCTAAAGCCCTGAAAATCTCCCATGTGTCTGTATCACAATGGGAACGGGGTGATAGTGAACCTACAGGGAAGAACCTTTTTGCCCTCAGTAAAGTAtt                                   >  2:9957/1‑140 (MQ=255)                                                                                                     cgtcgGAAAAACCTCAAACACACCCAAAGGTCTCTCGCTAAAGCCCTGAAAATCTCCCATGTGTCTGTATCACAATGGGAACGGGGTGATAGTGAACCTACAGGGAAGAACCTTTTTGCCCTCAGTAAAGTATTGCAATg                             >  1:478289/1‑140 (MQ=255)                                                                                                          gAAAAACCTCAAACACACCCAAAGGTCTCTCGCTAAAGCCCTGAAAATCTCCCATGTGTCTGTATCACAATGGGAACGGGGTGATAGTGAACCTACAGGGAAGAACCTTTTTGCCCTCAGTAAAGTATTGCAATGCTCAc                        >  1:479016/1‑140 (MQ=255)                                                                                                                ccTCAAACACACCCAAAGGTCTCTCGCTAAAGCCCTGAAAATCTCCCATGTGTCTGTATCACAATGGGAACGGGGTGATAGTGAACCTACAGGGAAGAACCTTTTTGCCCTCAGTAAAGTATTGCAATGCTCACCAACAt                  <  2:116826/140‑1 (MQ=255)                                                                                                                   cAAACACACCCAAAGGTCTCTCGCTAAAGCCCTGAAAATCTCCCATGTGTCTGTATCACAATGGGAACGGGGTGATAGTGAACCTACAGGGAAGAACCTTCTTGCCCTCAGTAAAGTATTGCAATGCTCACCAACATGGa               <  1:47675/140‑1 (MQ=255)                                                                                                                              aaaGGTCTCTCGCTAAAGCCCTGAAAATCTCCCATGTGTCTGTATCACAATGGGAACGGGGTGATAGTGAACCTACAGGGAAGAACCTTTTTGCCCTCAGTAAAGTATTGCAATGCTCACCAACATGGATTCTATTTGGc    >  1:164252/1‑140 (MQ=255)                                                                                                                               aaGGTCTCTCGCTAAAGCCCTGAAAATCTCCCATGTGTCTGTATCACAATGGGAACGGGGTGATAGTGAACCTACAGGGAAGAACCTTTTTGCCCTCAGTAAAGTATTGCAATGCTCACCAACATGGATTCTATTTGGCg   <  1:9957/140‑1 (MQ=255)                                                                                                                                aGGTCTCTCGCTAAAGCCCTGAAAATCTCCCATGTGTCTGTATCACAATGGGAACGGGGTGATAGTGAACCTACAGGGAAGAACCTTTTTGCCCTCAGTAAAGTATTGCAATGCTCACCAACATGGATTCTATTTGGCGa  <  2:153424/140‑1 (MQ=255)                                                                                                                                aGGTCTCTCGCTAAAGCCCTGAAAATCTCCCATGTGTCTGTATCACAATGGGAACGGGGTGATAGTGAACCTACAGGGAAGAACCTTTTTGCCCTCAGTAAAGTATTGCAATGCTCACCAACATGGATTCTATTTGGCGa  <  2:479016/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                                    AACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTATCGTCGGAAAAACCTCAAACACACCCAAAGGTCTCTTGCTAAAGCCCTGAAAATCTCCCATGTGTCTGTATCACAATGGGAACGGGGTGATAGTGAACCTACAGGGAAGAACCTTTTTGCCCTCAGTAAAGTATTGCAATGCTCACCAACATGGATTCTATTTGGCGA  >  NC_000913/1647876‑1648142

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