Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,486,973	C→T	Y42Y (TAC→TAT)	unknown →	Flp pilus assembly membrane protein TadE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,486,973	0	C	T	100.0%	87.2 / NA	34	Y42Y (TAC→TAT)	unknown	Flp pilus assembly membrane protein TadE
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (15/19);  total (15/19)

GATATTATTGTGAAGGAAAACATTAATACATCTCCTGATGATTTAATTAATGCTCTCAAAAAAACTGAACGCGTATCTAATAATATCGCCAGATAAAAAAAGGTTCTAT‑AATGAAAAAATTCCTATCGAATAAGAAGGGAGCATCAACCGTTGAGTTTGCCTTGACGGTTGTATTTTATTTCTTTGTCGTATTCCTAATTCTAGAATTCTGTCGTATATCTATTACAACTGCGTACTGGGACTTGGCTATCACTGAAAGTGCAAGAATTGCTAAAAATAGAACTGCAGAGGGTAATGATTATGCGGTAGAATTTGAGAAAGCTCTAAAACAACAGCTTGTTTATCAAGAATCTTCAACTTTAGGTTATCTCGCCCGTTTGGATAAAAATGGCGGTTACAAGATTGATGTTAAATATGTTGATTGTGGTAGTGAGAGTTCTTGTGTTAAGTCTTTATTGGATGGTAGGTTTA  >  1/1486738‑1487208                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                            gATATTATTGTGAAGGAAAACATTAACACATCCCCTGATGATTTAGTTAATGCCCTCAAGAAAACTGAACGCGTGTCTAGTAATATCACCAGAT‑‑AAAAAGGTTCTGTAAATGGAAAAATTCCTATCAAATAAAAAGGGAGCGTCAGCCGTTGAGTTTGCCTTGACGGTTGTATTCTATTTCTTTATCGTATTCTTAATTCTAGAATTCTGCCGTATATCTATTACCACTGCATATTGGGatttagccataa                                                                                                                                                                                                                             >  1:17258/1‑240 (MQ=255)             aaGGAAAACATTAACACATCCCCTGATGATTTAGTTAATGCCCTCAAGAAAACTGAACGCGTGTCTAGTAATATCACCAGAT‑‑AAAAAGGTTCTGTAAATGGAAAAATTCCTATCAAATAAAAAGGGAGCGTCAGCCGTTGAGTTTGCCTTGACGGTTGTATTCTATTTCTTTATCGTATTCTTAATTCTAGAATTCTGCCGTATATCTATTACCACTGCATATTGGGATTTAGCCATAACTGAAAGTGCaa                                                                                                                                                                                                                 >  1:447683/1‑251 (MQ=255)                gAAAACATTAACACATCCCCTGATGATTTAGTTAATGCCCTCAAGAAAACTGAACGCGTGTCTAGTAATATCACCAGAT‑‑AAAAAGGTTCTGTAAATGGAAAAATTCCTATCAAATAAAAAGGGAGCGTCAGCCGTTGAGTTTGCCTTGACGGTTGTATTCTATTTCTTTATCGTATTCTTAATTCTAGAATTCTGCCGTATATCTATTACCACTGCATATTGGGATTTAGCCATAACTGAAAGTGCAAg                                                                                                                                                                                                                <  1:887285/249‑1 (MQ=255)                 aaaaCATTAACACATCCCCTGATGATTTAGTTAATGCCCTCAAGAAAACTGAACGCGTGTCTAGTAATATCACCAGAT‑‑AAAAAGGTTCTGTAAATGGAAAAATTCCTATCAAATAAAAAGGGAGCGTCAGCCGTTGAGTTTGCCTTGACGGTTGTATTCTATTTCTTTATCGTATTCTTAATTCTAGAATTCTGCCGTATATCTATTACCACTGCATATTGGGATTTAGCCATAACTGAAAGTGCAAgga                                                                                                                                                                                                              <  1:698914/250‑3 (MQ=255)                 aaaaCATTAACACATCCCCTGATGATTTAGTTAATGCCCTCAAGAAAACTGAACGCGTGTCTAGTAATATCACCAGAT‑‑AAAAAGGTTCTGTAAATGGAAAAATTCCTATCAAATAAAAAGGGAGCGTCAGCCGTTGAGTTTGCCTTGACGGTTGTATTCTATTTCTTTATCGTATTCTTAATTCTAGAATTCTGCCGTATATCTATTACCACTGCATATTGGGATTTAGCCATAACTGAAAGTGCAAgga                                                                                                                                                                                                              <  1:586392/250‑3 (MQ=255)                           caCATCCCCTGATGATTTAGTTAATGCCCTCAAGAAAACTGAACGCGTGTCTAGTAATATCACCAGAT‑‑AAAAAGGTTCTGTAAATGGAAAAATTCCTATCAAATAAAAAGGGAGCGTCAGCCGTTGAGTTTGCCTTGACGGTTGTATTCTATTTCTTTATCGTATTCTTAATTCTAGAATTCTGCCGTATATCTATTACCACTGCATATTGGGATTTAGCCATAACTGAAAGTGCAAGGATCGCTaaaaac                                                                                                                                                                                                   <  1:494047/250‑2 (MQ=255)                               tccccTGATGATTTAGTTAATGCCCTCAAGAAAACTGAACGCGTGTCTAGTAATATCACCAGAT‑‑AAAAAGGTTCTGTAAATGGAAAAATTCCTATCAAATAAAAAGGGAGCGTCAGCCGTTGAGTTTGCCTTGACGGTTGTATTCTATTTCTTTATCGTATTCTTAATTCTAGAATTCTGCCGTATATCTATTACCACTGCATATTGGGATTTAGCCATAACTGAAAGTGCAAGGATCGCTAAAAACAGaa                                                                                                                                                                                               >  1:35572/4‑251 (MQ=255)                                     gatgatTTAGTTAATGCCCTCAAGAAAACTGAACGCGTGTCTAGTAATATCACCAGAT‑‑AAAAAGGTTCTGTAAATGGAAAAATTCCTATCAAATAAAAAGGGAGCGTCAGCCGTTGAGTTTGCCTTGACGGTTGTATTCTATTTCTTTATCGTATTCTTAATTCTAGAATTCTGCCGTATATCTATTACCACTGCATATTGGGATTTAGCCATAACTGAAAGTGCAAgga                                                                                                                                                                                                              <  1:462452/230‑3 (MQ=255)                                                taatGCCCTCAAGAAAACTGAACGCGTGTCTAGTAATATCACCAGAT‑‑AAAAAGGTTCTGTAAATGGAAAAATTCCTATCAAATAAAAAGGGAGCGTCAGCCGTTGAGTTTGCCTTGACGGTTGTATTCTATTTCTTTATCGTATTCTTAATTCTAGAATTCTGCCGTATATCTATTACCACTGCATATTGGGATTTAGCCATAACTGAAAGTGCAAGGATCGCTAAAAACAGAACGGCGGAAGGTAACGAt                                                                                                                                                                              <  1:482743/251‑1 (MQ=255)                                                                    aaCGCGTGTCTAGTAATATCACCAGAT‑‑AAAAAGGTTCTGTAAATGGAAAAATTCCTATCAAATAAAAAGGGAGCGTCAGCCGTTGAGTTTGCCTTGACGGTTGTATTCTATTTCTTTATCGTATTCTTAATTCTAGAATTCTGCCGTATATCTATTACCACTGCATATTGGGATTTAGCCATAACTGAAAGTGCAAGGATCGCTAAAAACAGAACGGCGGAAGGTAACGATTATGCTGCAGAATTTgaaaa                                                                                                                                                          <  1:423956/251‑4 (MQ=255)                                                                           gtCTAGTAATATCACCAGAT‑‑AAAAAGGTTCTGTAAATGGAAAAATTCCTATCAAATAAAAAGGGAGCGTCAGCCGTTGAGTTTGCCTTGACGGTTGTATTCTATTTCTTTATCGTATTCTTAATTCTAGAATTCTGCCGTATATCTATTACCACTGCATATTGGGatttagcca                                                                                                                                                                                                                                <  1:857768/173‑9 (MQ=255)                                                                                gtaataTCACCAGAT‑‑AAAAAGGTTCTGTAAATGGAAAAATTCCTATCAAATAAAAAGGGAGCGTCAGCCGTTGAGTTTGCCTTGACGGTTGTATTCTATTTCTTTATCGTATTCTTAATTCTAGAATTCTGCCGTATATCTATTACCACTGCATATTGGGATTTAGCCATAACTGAAAGTGCAAgga                                                                                                                                                                                                              >  1:632849/2‑185 (MQ=255)                                                                                gtaataTCACCAGAT‑‑AAAAAGGTTCTGTAAATGGAAAAATTCCTATCAAATAAAAAGGGAGCGTCAGCCGTTGAGTTTGCCTTGACGGTTGTATTCTATTTCTTTATCGTATTCTTAATTCTAGAATTCTGCCGTATATCTATTACCACTGCATATTGGGATTTAGCCATAACTGAAAGTGCAAGGATCGCTAAAAACAGAACGGCGGAAGGTAACGATTATGCTGCAGAATTTGAAAAAGCTCTAAaaca                                                                                                                                              >  1:179337/2‑251 (MQ=255)                                                                                             agataaaaaGGTTCTGTAAATGGAAAAATTCCTATCAAATAAAAAGGGAGCGTCAGCCGTTGAGTTTGCCTTGACGGTTGTATTCTATTTCTTTATCGTATTCTTAATTCTAGAATTCTGCCGTATATCTATTACCACTGCATATTGGGATTTAGCCATAACTGAAAGTg                                                                                                                                                                                                                    >  1:594033/5‑170 (MQ=255)                                                                                              gataaaaaGGTTCTGTAAATGGAAAAATTCCTATCAAATAAAAAGGGAGCGTCAGCCGTTGAGTTTGCCTTGACGGTTGTATTCTATTTCTTTATCGTATTCTTAATTCTAGAATTCTGCCGTATATCTATTACCACTGCATATTGGGATTTAGCCATAACTGAAAGTGCAAGGATCGCTAAAAACAGAACGGCGGAAGGTAACGATTATGCTGCAGAATTTgaaaa                                                                                                                                                          >  1:570457/4‑224 (MQ=255)                                                                                              gataaaaaGGTTCTGTAAATGGAAAAATTCCTATCAAATAAAAAGGGAGCGTCAGCCGTTGAGTTTGCCTTGACGGTTGTATTCTATTTCTTTATCGTATTCTTAATTCTAGAATTCTGCCGTATATCTATTACCACTGCATATTGGGATTTAGCCATAACTGAAAGTGCAAGGATCGCTAAAAACAGAACGGCGGAAGGTAACGATTATGCTGCAGAATTTGAAAAAGCTCTAAAACAACAACTTATTTa                                                                                                                                  >  1:596994/4‑251 (MQ=255)                                                                                               ataaaaaGGTTCTGTAAATGGAAAAATTCCTATCAAATAAAAAGGGAGCGTCAGCCGTTGAGTTTGCCTTGACGGTTGTATTCTATTTCTTTATCGTATTCTTAATTCTAGAATTCTGCCGTATATCTATTACCACTGCATATTGGGATTTAGCCATAACTGAAAGTGCAAGGATCGCTAAAAACAGAACGGCGGAAGGTAACGATTATGCTGCAGAATTTgaaaa                                                                                                                                                          >  1:663824/3‑223 (MQ=255)                                                                                               ataaaaaGGTTCTGTAAATGGAAAAATTCCTATCAAATAAAAAGGGAGCGTCAGCCGTTGAGTTTGCCTTGACGGTTGTATTCTATTTCTTTATCGTATTCTTAATTCTAGAATTCTGCCGTATATCTATTACCACTGCATATTGGGATTTAGCCATAACTGAAAGTGCAAGGATCGCTAAAAACAGAACGGCGGAAGGTAACGATTATGCTGCAGAATTTGAAAAAGCTCTAAAACAACAACttatttat                                                                                                                                 >  1:138808/3‑251 (MQ=255)                                                                                               ataaaaaGGTTCTGTAAATGGAAAAATTCCTATCAAATAAAAAGGGAGCGTCAGCCGTTGAGTTTGCCTTGACGGTTGTATTCTATTTCTTTATCGTATTCTTAATTCTAGAATTCTGCCGTATATCTATTACCACTGCATATTGGGATTTAGCCATAACTGAAAGTGCAAGGATCGCTAAAAACAGAACGGCGGAAGGTAACGATTATGCTGCAGAATTTGAAAAAGCTCTAAAACAACAACttatttat                                                                                                                                 >  1:304702/3‑251 (MQ=255)                                                                                                                    aaaaaTTCCTATCAAATAAAAAGGGAGCGTCAGCCGTTGAGTTTGCCTTGACGGTTGTATTCTATTTCTTTATCGTATTCTTAATTCTAGAATTCTGCCGTATATCTATTACCACTGCATATTGGGATTTAGCCATAACTGAAAGTGCAAGGATCGCTAAAAACAGAACGGCGGAAGGTAACGATTATGCTGCAGAATTTGAAAAAGCTCTAAAACAACAACTTATTTATCAAGAAACCTCAACTATAggc                                                                                                            <  1:335046/251‑2 (MQ=255)                                                                                                                                  aaTAAAAAGGGAGCGTCAGCCGTTGAGTTTGCCTTGACGGTTGTATTCTATTTCTTTATCGTATTCTTAATTCTAGAATTCTGCCGTATATCTATTACCACTGCATATTGGGATTTAGCCATAACTGAAAGTGCAAGGATCGCTAAAAACAGAACGGCGGAAGGTAACGATTATGCTGCAGAATTTGAAAAAGCTCTAAAACAACAACTTATTTATCAAGAAACCTCAACtat                                                                                                                <  1:371139/233‑3 (MQ=255)                                                                                                                                  aaTAAAAAGGGAGCGTCAGCCGTTGAGTTTGCCTTGACGGTTGTATTCTATTTCTTTATCGTATTCTTAATTCTAGAATTCTGCCGTATATCTATTACCACTGCATATTGGGATTTAGCCATAACTGAAAGTGCAAGGATCGCTAAAAACAGAACGGCGGAAGGTAACGATTATGCTGCAGAATTTGAAAAAGCTCTAAAACAACAACTTATTTATCAAGAAACCTCAACTATAGGCTATCTTGCCCGttt                                                                                              >  1:111173/1‑251 (MQ=255)                                                                                                                                      aaaagGGAGCGTCAGCCGTTGAGTTTGCCTTGACGGTTGTATTCTATTTCTTTATCGTATTCTTAATTCTAGAATTCTGCCGTATATCTATTACCACTGCATATTGGGATTTAGCCATAACTGAAAGTGCAAgga                                                                                                                                                                                                              <  1:446240/133‑3 (MQ=255)                                                                                                                                                             tttGCCTTGACGGTTGTATTCTATTTCTTTATCGTATTCTTAATTCTAGAATTCTGCCGTATATCTATTACCACTGCATATTGGGATTTAGCCATAACTGAAAGTGCAAGGATCGCTAAAAACAGAACGGCGGAAGGTAACGATTATGCTGCAGAATTTGAAAAAGCTCTAAAACAACAACTTATTTATCAAGAAACCTCAACTATAGGCTATCTTGCCCGTTTGGATAAAAATGGCGGTTatgaaattg                                                                    <  1:274106/250‑9 (MQ=255)                                                                                                                                                                                         ttatCGTATTCTTAATTCTAGAATTCTGCCGTATATCTATTACCACTGCATATTGGGATTTAGCCATAACTGAAAGTGCAAGGATCGCTAAAAACAGAACGGCGGAAGGTAACGATTATGCTGCAGAATTTGAAAAAGCTCTAAAACAACAACTTATTTATCAAGAAACCTCAACTATAGGCTATCTTGCCCGTTTGGATAAAAATGGCGGTTATGAAATTGATGTTAAATATGTTGATTGTGGTAGTgag                                       <  1:690141/248‑1 (MQ=255)                                                                                                                                                                                         ttatCGTATTCTTAATTCTAGAATTCTGCCGTATATCTATTACCACTGCATATTGGGATTTAGCCATAACTGAAAGTGCAAGGATCGCTAAAAACAGAACGGCGGAAGGTAACGATTATGCTGCAGAATTTGAAAAAGCTCTAAAACAACAACTTATTTATCAAGAAACCTCAACTATAGGCTATCTTGCCCGTTTGGATAAAAATGGCGGTTATGAAATTGATGTTAAATATGTTGATTGTGGTAGTgag                                       <  1:726941/248‑1 (MQ=255)                                                                                                                                                                                           atCGTATTCTTAATTCTAGAATTCTGCCGTATATCTATTACCACTGCATATTGGGATTTAGCCATAACTGAAAGTGCAAGGATCGCTAAAAACAGAACGGCGGAAGGTAACGATTATGCTGCAGAATTTGAAAAAGCTCTAAAACAACAACTTATTTATCAAGAAACCTCAACTATAGGCTATCTTGCCCGTTTGGATAAAAATGGCGGTTATGAAATTGATGTTAAATATGTTGATTGTGGt                                             <  1:292311/242‑1 (MQ=255)                                                                                                                                                                                                 ttcttAATTCTAGAATTCTGCCGTATATCTATTACCACTGCATATTGGGATTTAGCCATAACTGAAAGTGCAAGGATCGCTAAAAACAGAACGGCGGAAGGTAACGATTATGCTGCAGAATTTGAAAAAGCTCTAAAACAACAACTTATTTATCAAGAAACCTCAACTATAGGCTATCTTGCCCGTTTGGATAAAAATGGCGGTTATGAAATTGATGTTAAATATGTTGATTGTGGTAGTGAGAGCTCTtg                               <  1:620439/247‑1 (MQ=255)                                                                                                                                                                                                          cTAGAATTCTGCCGTATATCTATTACCACTGCATATTGGGATTTAGCCATAACTGAAAGTGCAAGGATCGCTAAAAACAGAACGGCGGAAGGTAACGATTATGCTGCAGAATTTGAAAAAGCTCTAAAACAACAACTTATTTATCAAGAAACCTCAACTATAGGCTATCTTGCCCGTTTGGATAAAAATGGCGGTTATGAAATTGATGTTAAATATGTTGATTGTGGTAGTGAGAGCTCTTGTATTaaatc                      <  1:360761/251‑4 (MQ=255)                                                                                                                                                                                                                           atCTATTACCACTGCATATTGGGATTTAGCCATAACTGAAAGTGCAAGGATCGCTAAAAACAGAACGGCGGAAGGTAACGATTATGCTGCAGAATTTGAAAAAGCTCTAAAACAACAACTTATTTATCAAGAAACCTCAACTATAGGCTATCTTGCCCGTTTGGATAAAAATGGCGGTTATGAAATTGATGTTAAATATGTTGATTGTGGTAGTGAGAGCTCTTGTATTAAATCTTTAtt                <  1:768307/240‑1 (MQ=255)                                                                                                                                                                                                                           atCTATTACCACTGCATATTGGGATTTAGCCATAACTGAAAGTGCAAGGATCGCTAAAAACAGAACGGCGGAAGGTAACGATTATGCTGCAGAATTTGAAAAAGCTCTAAAACAACAACTTATTTATCAAGAAACCTCAACTATAGGCTATCTTGCCCGTTTGGATAAAAATGGCGGTTATGAAATTGATGTTAAATATGTTGATTGTGGTAGTGAGAGCTCTTGTATTAAATCTTTATTAGATGGTaagt     >  1:696224/1‑248 (MQ=255)                                                                                                                                                                                                                            tCTATTACCACTGCATATTGGGATTTAGCCATAACTGAAAGTGCAAGGATCGCTAAAAACAGAACGGCGGAAGGTAACGATTATGCTGCAGAATTTGAAAAAGCTCTAAAACAACAACTTATTTATCAAGAAACCTCAACTATAGGCTATCTTGCCCGTTTGGATAAAAATGGCGGTTATGAAATTGATGTTAAATATGTTGATTGTGGTAGTGAGAGCTCTTGTATTaaatc                      <  1:140087/233‑4 (MQ=255)                                                                                                                                                                                                                             cTATTACCACTGCATATTGGGATTTAGCCATAACTGAAAGTGCAAGGATCGCTAAAAACAGAACGGCGGAAGGTAACGATTATGCTGCAGAATTTGAAAAAGCTCTAAAACAACAACTTATTTATCAAGAAACCTCAACTATAGGCTATCTTGCCCGTTTGGATAAAAATGGCGGTTATGAAAATGATGTTAAATATGTTGATtg                                                 >  1:764995/1‑205 (MQ=255)                                                                                                                                                                                                                              tATTACCACTGCATATTGGGATTTAGCCATAACTGAAAGTGCAAGGATCGCTAAAAACAGAACGGCGGAAGGTAACGATTATGCTGCAGAATTTGAAAAAGCTCTAAAACAACAACTTATTTATCAAGAAACCTCAACTATAGGCTATCTTGCCCGTTTGGATAAAAATGGCGGTTATGAAATTGATGTTAAATATGTTGATTGTGGTAGTGAGAGCTCTTGTATTAAATCTTTATTAGATGGTAAGTTTa  >  1:731563/1‑251 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                            GATATTATTGTGAAGGAAAACATTAATACATCTCCTGATGATTTAATTAATGCTCTCAAAAAAACTGAACGCGTATCTAATAATATCGCCAGATAAAAAAAGGTTCTAT‑AATGAAAAAATTCCTATCGAATAAGAAGGGAGCATCAACCGTTGAGTTTGCCTTGACGGTTGTATTTTATTTCTTTGTCGTATTCCTAATTCTAGAATTCTGTCGTATATCTATTACAACTGCGTACTGGGACTTGGCTATCACTGAAAGTGCAAGAATTGCTAAAAATAGAACTGCAGAGGGTAATGATTATGCGGTAGAATTTGAGAAAGCTCTAAAACAACAGCTTGTTTATCAAGAATCTTCAACTTTAGGTTATCTCGCCCGTTTGGATAAAAATGGCGGTTACAAGATTGATGTTAAATATGTTGATTGTGGTAGTGAGAGTTCTTGTGTTAAGTCTTTATTGGATGGTAGGTTTA  >  1/1486738‑1487208

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