Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	BS168_reference_genome	1,704,223	C→T	Q5* (CAA→TAA)	fliZ →	flagellar regulatory protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	BS168_reference_genome	1,704,223	0	C	T	100.0%	163.5 / NA	46	Q5* (CAA→TAA)	fliZ	flagellar regulatory protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (27/19);  total (27/19)

ATCAATGGTTATTGATGCCATTCAGGCCGGCGCTAAGGATTTCATTGTGAAACCTTTCCAGGCTGACCGTGTGCTTGAAGCAATCAACAAAACATTAAATTAAAGGGTGTACGACTGTTGAAAAAGAGTCAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAAT  >  BS168_reference_genome/1704094‑1704355                                                                                                                                  |                                                                                                                                     atcaatGGTTATTGATGCCATTCAGGCCGGCGCTAAGGATTTCATTGTGAAACCTTTCCAGGCTGACCGTGTGCTTGAAGCAATCAACAAAACATTAAATTAAAGGGTGTACGACTGTTGAAAAAGAGTTAATATTTTa                                                                                                                             >  8:77227/1‑139 (MQ=255) atcaatGGTTATTGATGCCATTCAGGCCGGCGCTAAGGATTTCATTGTGAAACCTTTCCAGGCTGACCGTGTGCTTGAAGCAATCAACAAAACATTAAATTAAAGGGTGTACGACTGTTGAAAAAGAGTTAATATTTTa                                                                                                                             <  3:53028/139‑1 (MQ=255)                 gCCATTCAGGCCGGCGCTAAGGATTTCATTGTGAAACCTTTCCAGGCTGACCGTGTGCTTGAAGCAATCAACAAAACATTAAATTAAAGGGTGTACGACTGTTGAAAAAGAGTTAATATTTTATTGTTTTTATTTGttt                                                                                                             >  4:149985/1‑139 (MQ=255)                   cATTCAGGCCGGCGCTAAGGATTTCATTGTGAAACCTTTCCAGGCTGACCGTGTGCTTGAAGCAATCAACAAAACATTAAATTAAAGGGTGTACGACTGTTGAAAAAGAGTTAATATTTTATTGTTTTTATTTGttttt                                                                                                           <  8:44558/139‑1 (MQ=255)                   cATTCAGGCCGGCGCTAAGGATTTCATTGTGAAACCTTTCCAGGCTGACCGTGTGCTTGAAGCAATCAACAAAACATTAAATTAAAGGGTGTACGACTGTTGAAAAAGAGTTAATATTTTATTGTTTTTATTTGttttt                                                                                                           <  7:233902/139‑1 (MQ=255)                   cATTCAGGCCGGCGCTAAGGATTTCATTGTGAAACCTTTCCAGGCTGACCGTGTGCTTGAAGCAATCAACAAAACATTAAATTAAAGGGTGTACGACTGTTGAAAAAGAGTTAATATTTTATTGTTTTTATTTGttttt                                                                                                           <  1:135634/139‑1 (MQ=255)                         ggCCGTCGCTAAGGATTTCATTGTGAAACCTTTCCAGTCTGACCGTGTGCTTGAAGCAATCAACAAAACATTAAATTAAAGGGTGTACGAATGTTGAAAAAGAGTTAATATTTTATTGTTTTTATTTGTTTTTTCGttt                                                                                                     <  1:26580/139‑1 (MQ=255)                             ggCGCTAAGGATTTCATTGTGAAACCTTTCCAGGCTGACCGTGTGCTTGAAGCAATCAACAAAACATTAAATTAAAGGGTGTACGACTGTTGAAAAAGAGTTAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTAtt                                                                                                 <  6:86417/139‑1 (MQ=255)                             ggCGCTAAGGATTTCATTGTGAAACCTTTCCAGGCTGACCGTGTGCTTGAAGCAATCAACAAAACATTAAATTAAAGGGTGTACGACTGTTGAAAAAGAGTTAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTAtt                                                                                                 <  7:90803/139‑1 (MQ=255)                             ggCGCTAAGGATTTCATTGTGAAACCTTTCCAGGCTGACCGTGTGCTTGAAGCAATCAACAAAACATTAAATTAAAGGGTGTACGACTGTTGAAAAAGAGTTAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTAtt                                                                                                 <  5:104404/139‑1 (MQ=255)                               cgcTAAGGATTTCATTGTGAAACCTTTCCAGGCTGACCGTGTGCTTGAAGCAATCAACAAAACATTAAATTAAAGGGTGTACGACTGTTGAAAAAGAGTTAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCa                                                                                               >  8:247161/1‑139 (MQ=255)                                           cATTGTGAAACCTTTCCAGGCTGACCGTGTGCTTGAAGCAATCAACAAAACATTAAATTAAAGGGTGTACGACTGTTGAAAAAGAGTTAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTAtt                                                                                                 <  3:19187/125‑1 (MQ=255)                                           cATTGTGAAACCTTTCCAGGCTGACCGTGTGCTTGAAGCAATCAACAAAACATTAAATTAAAGGGTGTACGACTGTTGAAAAAGAGTTAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTAtt                                                                                                 >  4:19187/1‑125 (MQ=255)                                              tgtgAAACCTTTCCAGGCTGACCGTGTGCTTGAAGCAATCAACAAAACATTAAATTAAAGGGTGTACGACTGTTGAAAAAGAGTTAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATtg                                                                                >  4:83303/1‑139 (MQ=255)                                                tgAAACCTTTCCAGGCTGACCGTGTGCTTGAAGCAATCAACAAAACATTAAATTAAAGGGTGTACGACTGTTGAAAAAGAGTTAATATTTTATTGTTTTTATTTGTTTTTTCGtt                                                                                                      >  1:302126/1‑115 (MQ=255)                                                tgAAACCTTTCCAGGCTGACCGTGTGCTTGAAGCAATCAACAAAACATTAAATTAAAGGGTGTACGACTGTTGAAAAAGAGTTAATATTTTATTGTTTTTATTTGTTTTTTCGtt                                                                                                      <  2:302126/115‑1 (MQ=255)                                                tgAAACCTTTCCAGGCTGACCGTGTGCTTGAAGCAATCAACAAAACATTAAATTAAAGGGTGTACGACTGTTGAAAAAGAGTTAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATtgct                                                                              >  7:271479/1‑139 (MQ=255)                                                  aaaCCTTTCCAGGCTGACCGTGTGCTTGAAGCAATCAACAAAACATTAAATTAAAGGGTGTACGACTGTTGAAAAAGAGTTAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTAtt                                                                                                 >  6:321514/1‑118 (MQ=255)                                                  aaaCCTTTCCAGGCTGACCGTGTGCTTGAAGCAATCAACAAAACATTAAATTAAAGGGTGTACGACTGTTGAAAAAGAGTTAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTAtt                                                                                                 <  5:321514/118‑1 (MQ=255)                                                                tGACCGTGTGCTTGAAGCAATCAACAAAACATTAAATTAAAGGGTGTACGACTGTTGAAAAAGAGTTAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGAct                                                              >  5:229191/1‑139 (MQ=255)                                                                  aCCGTGTGCTTGAAGCAATCAACAAAACATTAAATTAAAGGGTGTACGACTGTTGAAAAAGAGTTAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGActct                                                            <  8:271479/139‑1 (MQ=255)                                                                       gtgCTTGAAGCAATCAACAAAACATTAAATTAAAGGGTGTACGACTGTTGAAAAAGAGTTAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATaa                                                       >  4:224741/1‑139 (MQ=255)                                                                            tGAAGCAATCAACAAAACATTAAATTAAAGGGTGTACGACTGTTGAAAAAGAGTTAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCTGTGTACATCCGATTGCTGCTGCGGCGGCAGACTCTGATAATTCaa                                                  >  5:41794/1‑139 (MQ=255)                                                                             gAAGCAATCAACAAAACATTAAATTAAAGGGTGTACGACTGTTGAAAAAGAGTTAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCaa                                                  <  3:224741/138‑1 (MQ=255)                                                                                    tCAACAAAACATTAAATTAAAGGGTGTACGACTGTTGAAAAAGAGTTAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAAc                                          >  4:70488/1‑139 (MQ=255)                                                                                                taaattaaaGGGTGTACGACTGTTGAAAAAGAGTTAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGa                                                          >  5:232833/1‑111 (MQ=255)                                                                                                taaattaaaGGGTGTACGACTGTTGAAAAAGAGTTAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGa                                                          <  6:232833/111‑1 (MQ=255)                                                                                                taaattaaaGGGTGTACGACTGTTGAAAAAGAGTTAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCa                               >  6:195822/1‑138 (MQ=255)                                                                                                       aaGGGTGTACGACTGTTGAAAAAGAGTTAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGa                                        >  7:24159/1‑122 (MQ=255)                                                                                                       aaGGGTGTACGACTGTTGAAAAAGAGTTAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGa                                        <  8:24159/122‑1 (MQ=255)                                                                                                                    tGTTGAAAAAGAGTTAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCTACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTGc          <  6:229191/139‑1 (MQ=255)                                                                                                                    tGTTGAAAAAGAGTTAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTg           >  2:214230/1‑138 (MQ=255)                                                                                                                    tGTTGAAAAAGAGTTAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTg           >  8:210978/1‑138 (MQ=255)                                                                                                                      ttGAAAAAGAGTTAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTg           >  5:338601/1‑136 (MQ=255)                                                                                                                      ttGAAAAAGAGTTAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTg           >  6:316740/1‑136 (MQ=255)                                                                                                                      ttGAAAAAGAGTTAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTg           >  5:37984/1‑136 (MQ=255)                                                                                                                      ttGAAAAAGAGTTAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTg           >  4:145662/1‑136 (MQ=255)                                                                                                                      ttGAAAAAGAGTTAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTg           >  3:62952/1‑136 (MQ=255)                                                                                                                      ttGAAAAAGAGTTAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTg           >  3:230345/1‑136 (MQ=255)                                                                                                                      ttGAAAAAGAGTTAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTg           >  3:181565/1‑136 (MQ=255)                                                                                                                      ttGAAAAAGAGTTAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTg           >  1:170354/1‑136 (MQ=255)                                                                                                                      ttGAAAAAGAGTTAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGAAAATTCAACTGTAAACGAATGGTTACAAAAGAAAGATGaaaaaaatt           >  5:53890/1‑133 (MQ=255)                                                                                                                      ttGAAAAAGAGTTAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGAGGCAGACTCTGATATTTCATCTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTg           >  1:14559/1‑136 (MQ=255)                                                                                                                            aaGAGTTAATATTTTTTTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTGCAGatcaat  <  8:50055/139‑1 (MQ=255)                                                                                                                            aaGAGTTAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTAGATCCGATTGCTGGTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTGCAGatcaat  <  5:9610/139‑1 (MQ=255)                                                                                                                            aaGAGTTAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCTGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTGCAGatcaat  <  7:247161/139‑1 (MQ=255)                                                                                                                                  |                                                                                                                                     ATCAATGGTTATTGATGCCATTCAGGCCGGCGCTAAGGATTTCATTGTGAAACCTTTCCAGGCTGACCGTGTGCTTGAAGCAATCAACAAAACATTAAATTAAAGGGTGTACGACTGTTGAAAAAGAGTCAATATTTTATTGTTTTTATTTGTTTTTTCGTTTTATTCAGTGTACATCCGATTGCTGCTGCCGCGGCAGACTCTGATAATTCAACTGTAAACGAATGGTTTCAAAAGAAAGATGAAAAAACTGCAGATCAAT  >  BS168_reference_genome/1704094‑1704355

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑