Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	855,490	T→A	A153A (GCA→GCT)	unknown ←	Type II/IV secretion system protein TadC, associated with Flp pilus assembly

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	855,490	0	T	A	100.0%	109.5 / NA	33	A153A (GCA→GCT)	unknown	Type II/IV secretion system protein TadC, associated with Flp pilus assembly
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (13/20);  total (13/20)

TAGTTGATGATAAATAGAGGAACCAAAATTAAGGCTTTGTTGCATTACAGTACAAAACATTCGGATTTCTGCTGTTGGTAAAGACACAGAGAGATCCTGCAATGCCTGCGACATACCGGTAATCTCGGATTTACGAATAATACGTAACATAACATAAGTCAGATCAGGGTTTAATATTTTAAAGTCAGTTGCAACGTGCTTTAGCGCAGCTTCAATACCCATTCCCGTTTGTAAATTAACTGCTACTAAATCAATGAAACCTGCTAAATCATTCATGATTTTTTTAATTTTGGATTTTAAAATAGAACCAATTAAAATTCCTGGAATAACAATAGTTACTATTAATATAGTTAGAAAGCCAAGTAACAGAGATGAGCTGTCAGTAGTTTTATCAATAAGTGAATAAAG  >  1/855250‑855657                                                                                                                                                                                                                                                 |                                                                                                                                                                        tAGTTGATGATAAATAGAGGAACCAAAATTAAGGCTTTGTTGCATTACAGTACAAAACATTCGGATTTCTGCTGTTGGTAAAGACACAGAGAGATCCTGCAATGCCTGCGACATACCGGTAATCTCGGATTTACGAATAATACGTAACATAACATAAGTCAGATCAGGGTTTAATATTTTAAAGTCAGTTGCAACGTGCTTTAGCGCAGCTTCAATACCCATTCCCGTTTGTAAATTAACAGCTACTAAAt                                                                                                                                                               <  1:510932/251‑1 (MQ=255)     tgatgaTAAATAGAGGAACCAAAATTAAGGCTTTGTTGCATTACAGTACAAAACATTCGGATTTCTGCTGTTGGTAAAGACACAGAGAGATCCTGCAATGCCTGCGACATACCGGTAATCTCGGATTTACGAATAATACGTAACATAACATAAGTCAGATCAGGGTTTAATATTTTAAAGTCAGTTGCAACGTGCTTTAGCGCAGCTTCAATACCCATTCCCGTTTGTAAATTAACAGCTACTAAATCAAt                                                                                                                                                           <  1:699177/251‑1 (MQ=255)                                                                             ggTAAAGACACAGAGAGATCCTGCAATGCCTGCGACATACCGGTAATCTCGGATTTACGAATAATACGTAACATAACATAAGTCAGATCAGGGTTTAATATTTTAAAGTCAGTTGCAACGTGCTTTAGCGCAGCTTCAATACCCATTCCCGTTTGTAAATTAACAGctacta                                                                                                                                                                  <  1:350690/172‑1 (MQ=255)                                                                                   gACACAGAGAGATCCTGCAATGCCTGCGACATACCGGTAATCTCGGATTTACGAATAATACGTAACATAACATAAGTCAGATCAGGGTTTAATATTTTAAAGTCAGTTGCAACGTGCTTTAGCGCAGCTTCAATACCCATTCCCGTTTGTAAATTAACAGCTACTAAATCAATGAAACCTGCTAAATCATTCATGATTTTTTTAATTTTGGATTTTAAAAt                                                                                                           <  1:665696/221‑1 (MQ=255)                                                                                         gagagaTCCTGCAATGCCTGCGACATACCGGTAATCTCGGATTTACGAATAATACGTAACATAACATAAGTCAGATCAGGGTTTAATATTTTAAAGTCAGTTGCAACGTGCTTTAGCGCAGCTTCAATACCCATTCCCGTTTGTAAATTAACAGCTACTAAATCAATGAAACCTGCTAAATCATTCATGATTTTTTTAATTTTGGATTTTAAAATAGAACCAATTAAAATTCCTGGAATAACAATAGTTa                                                                        <  1:795950/250‑1 (MQ=255)                                                                                           gagaTCCTGCAATGCCTGCGACATACCGGTAATCTCGGATTTACGAATAATACGTAACATAACATAAGTCAGATCAGGGTTTAATATTTTAAAGTCAGTTGCAACGTGCTTTAGCGCAGCTTCAATACCCATTCCCGTTTGTAAATTAACAGCTACTAAATCAATGAAACCTGCTAAAt                                                                                                                                             <  1:251250/179‑1 (MQ=255)                                                                                             gaTCCTGCAATGCCTGCGACATACCGGTAATCTCGGATTTACGAATAATACGTAACATAACATAAGTCAGATCAGGGTTTAATATTTTAAAGTCAGTTGCAACGTGCTTTAGCGCAGCTTCAATACCCATTCCCGTTTGTAAATTAACAGCTACTAAATCAATGAAACCTGCTAAATCATTCATGATTTTTTTAATTTTGGATTTTAAAATAGAACCAATTAAAATTCCTGGAATAACAATAGTTACTAtt                                                                   >  1:44947/1‑251 (MQ=255)                                                                                                 cTGCAATGCCTGCGACATACCGGTAATCTCGGATTTACGAATAATACGTAACATAACATAAGGCAGATCAGGGTTTAATATTTTAAAGTCAGTTGCAACGTGCTTTAGCGCAGCTTCAATACCCATTCCCGTTTGTAAATTAACAGCTACTAAATCAATGAAACCTGCTAAATCATTCATGATTTTTTTAATTTTGGATTTTAAAATAGAACCAATTAAAATTCCTGGAATAACAATAGTTACTATTAata                                                               <  1:612008/251‑1 (MQ=255)                                                                                                                       gTAATCTCGGATTTACGAATAATACGTAACATAACATAAGTCAGATCAGGGTTTAATATTTTAAAGTCAGTTGCAACGTGCTTTAGCGCAGCTTCAATACCCATTCCCGTTTGTAAATTAACAGCTACTAAATCAATGAAACCTGCTAAATCATTCATGATTTTTTTAATTTTGGATTTTAAAATAGAACCAATTAAAATTCCTGGAATAACAATAGTTACTATTAatat                                                              >  1:70082/1‑230 (MQ=255)                                                                                                                             tcGGATTTACGAATAATACGTAACATAACATAAGTCAGATCAGGGTTTAATATTTTAAAGTCAGTTGCAACGTGCTTTAGCGCAGCTTCAATACCCATTCCCGTTTGTAAATTAACAGCTACTAAATCAATGAAACCTGCTAAATCATTCATGATTTTTTTAATTTTGGATGTTAAAATAGAACCAATTAAAATTCCTGGAATAACAATAGTTACTATTAATATAGTTAGAAAGCCAAGTAACAGAGATGa                                   <  1:694927/251‑1 (MQ=255)                                                                                                                                      cGAATAATACGTAACATAACATAAGTCAGATCAGGGTTTAATATTTTAAAGTCAGTTGCAACGTGCTTTAGCGCAGCTTCAATACCCATTCCCGTTTGTAAATTAACAGCTACTAAATCAATGAAACCTGCTAAATCATTCATGATTTTTTTAATTTTGGATTTTAAAATAGAACCAATTAAAATTCCTGGAATAACAATAGTTACTATTAATATAGTTAGAAAGCCAAGTAACAGAGATGAGCTGTCagt                          <  1:205908/251‑1 (MQ=255)                                                                                                                                         ataataCGTAACATAACATAAGTCAGATCAGGGTTTAATATTTTAAAGTCAGTTGCAACGTGCTTTAGCGCAGCTTCAATACCCATTCCCGTTTGTAAATTAACAGCTACTaaa                                                                                                                                                                >  1:635977/1‑114 (MQ=255)                                                                                                                                            ataCGTAACATAACATAAGTCAGATCAGGGTTTAATATTTTAAAGTCAGTTGCAACGTGCTTTAGCGCAGCTTCAATACCCATTCCCGTTTGTAAATTAACAGCTACTAAATCAATGAAACCTGCTAAAtc                                                                                                                                            <  1:786064/131‑1 (MQ=255)                                                                                                                                            ataCGTAACATAACATAAGTCAGATCAGGGTTTAATATTTTAAAGTCAGTTGCAACGTGCTTTAGCGCAGCTTCAATACCCATTCCCGTTTGTAAATTAACAGCTACTAAATCAATGAAACCTGCTAAAt                                                                                                                                             <  1:584965/130‑1 (MQ=255)                                                                                                                                            ataCGTAACATAACATAAGTCAGATCAGGGTTTAATATTTTAAAGTCAGTTGCAACGTGCTTTAGCGCAGCTTCAATACCCATTCCCGTTTGTAAATTAACAGCTACTAAATCAATGAAACCTGCTAAATCATTCATGATTTTTTTAATTTTGGATTTTAAAATAGAACCAATTAAAATTCCTGGAATAACAATAGTTACTAtt                                                                   <  1:662688/204‑1 (MQ=255)                                                                                                                                              aCGTAACATAACATAAGTCAGATCAGGGTTTAATATTTTAAAGTCAGTTGCAACGTGCTTTAGCGCAGCTTCAATACCCATTCCCGTTTGTAAATTAACAGCTACTAAATCAATGAAACCTGCTAAATCATTCATGATTTTTTTAATTTTGGATTTTAAAATAGAACCAATTAAAATTCCTGGAATAACAATAGTTACTATTAATATAGTTAGAAAGCCAAGTAACAGAGATGAGCTGTCagt                          >  1:113262/1‑243 (MQ=255)                                                                                                                                                   acataacataAGTCAGATCAGGGTTTAATATTTTAAAGTCAGTTGCAACGTGCTTTAGCGCAGCTTCAATACCCATTCCCGTTTGTAAATTAACAGCTACTAAATCAATGAAACCTGCTAAATCATTCATGATTTTTTTAATTTTGGATTTTAAAATAGAACCAATTAAAATTCCTGGAATAACAATAGTTACTATTAatat                                                              <  1:578923/202‑1 (MQ=255)                                                                                                                                                   acataacataAGTCAGATCAGGGTTTAATATTTTAAAGTCAGTTGCAACGTGCTTTAGCGCAGCTTCAATACCCATTCCCGTTTGTAAATTAACAGCTACTAAATCAATGAAACCTGCTAAATCATTCATGATTTTTTTAATTTTGGATTTTAAAATAGAACCAATTAAAATTCCTGGAATAACAATAGTTACTATTAatat                                                              >  1:254665/1‑202 (MQ=255)                                                                                                                                                     ataacataaGTCAGATCAGGGTTTAATATTTTAAAGTCAGTTGCAACGTGCTTTAGCGCAGCTTCAATACCCATTCCCGTTTGTAAATTAACAGCTACTAAATCAATGAAACCTGCTAAATCATTCATGATTTTTTTAATTTTGGATTTTAAAATAGAACCAATTAAAATTCCTGGAATAACAATAGTTACTATTAATATAGTTAGAAAGCCAAGTAACAGAGATGAGCTGTCAGTAGTTTTATCAATAAg           <  1:527838/251‑1 (MQ=255)                                                                                                                                                      taacataaGTCAGATCAGGGTTTAATATTTTAAAGTCAGTTGCAACGTGCTTTAGCGCAGCTTCAATACCCATTCCCGTTTGTAAATTAACAGCTACTAAATCAATGAAACCTGCTAAAtc                                                                                                                                            <  1:675086/121‑1 (MQ=255)                                                                                                                                                      taacataaGTCAGATCAGGGTTTAATATTTTAAAGTCAGTTGCAACGTGCTTTAGCGCAGCTTCAATACCCATTCCCGTTTGTAAATTAACAGCTACTAAATCAATGAAACCTGCTAAATCATTCATGATTTTTTTAATTTTGGATTTTAAAATAGAACCAATTAAAATTCCTGGAATAACAATAg                                                                           <  1:320163/186‑1 (MQ=255)                                                                                                                                                        acataaGTCAGATCAGGGTTTAATATTTTAAAGTCAGTTGCAACGTGCTTTAGCGCAGCTTCAATACCCATTCCCGTTTGTAAATTAACAGCTACTAAATCAATGAAACCTGCTAAATCATTCATGATTTTTTTAATTTTGGATTTTAAAATAGAACCAATTAAAATTCCTGGAATAACAATAGTTACTATTAata                                                               >  1:341250/1‑196 (MQ=255)                                                                                                                                                               tCAGATCAGGGTTTAATATTTTAAAGTCAGTTGCAACGTGCTTTAGCGCAGCTTCAATACCCATTCCCGTTTGTAAATTAACAGCTACTAAATCAATGAAACCTGCTAAATCATTCATGATTTTTTTAATTTTGGATTTTAAAATAGAACCAATTAAAATTCCTGGAATAACAATAGTTACTATTAATATAGTTAGAAAGCCAAGTAACAGAGATGAGCTGTCAGTAGTTTTATCAATAAGTGAATAAAg  <  1:432260/250‑1 (MQ=255)                                                                                                                                                                       gggTTTAATATTTTAAAGTCAGTTGCAACGTGCTTTAGCGCAGCTTCAATACCCATTCCCGTTTGTAAATTAACAGCTACTAAATCAATGAAACCTGCTAAATCATTCATGATTTTTTTAATTTTGGATTTTAAAATAGAACCAATTAAAATTCCTGGAATAACAATAGTTACTATTAatat                                                              >  1:34736/1‑182 (MQ=255)                                                                                                                                                                                atTTTAAAGTCAGTTGCAACGTGCTTTAGCGCAGCTTCAATACCCATTCCCGTTTGTAAATTAACAGCTACTAAATCAATGAAACCTGCTaaa                                                                                                                                              <  1:229579/93‑1 (MQ=255)                                                                                                                                                                                atTTTAAAGTCAGTTGCAACGTGCTTTAGCGCAGCTTCAATACCCATTCCCGTTTGTAAATTAACAGCTACTAAATCAATGAAACCTGCTAAAt                                                                                                                                             >  1:457162/1‑94 (MQ=255)                                                                                                                                                                                atTTTAAAGTCAGTTGCAACGTGCTTTAGCGCAGCTTCAATACCCATTCCCGTTTGTAAATTAACAGCTACTAAATCAATGAAACCTGCTAAATCATTCATGATTTTTTTAATTTTGGATTTTAAAAt                                                                                                           >  1:883958/1‑128 (MQ=255)                                                                                                                                                                                atTTTAAAGTCAGTTGCAACGTGCTTTAGCGCAGCTTCAATACCCATTCCCGTTTGTAAATTAACAGCTACTAAATCAATGAAACCTGCTAAATCATTCATGATTTTTTTAATTTTGGATTTTAAAATAGAACCAAtt                                                                                                 <  1:259247/138‑1 (MQ=255)                                                                                                                                                                                 ttttAAAGTCAGTTGCAACGTGCTTTAGCGCAGCTTCAATACCCATTCCCGTTTGTAAATTAACAGCTACTAAATCAATGAAACCTGCTAAATCATTCATGATTTTTTTAATTTTGGATTTTAAAATAGAACCAATTAAAAt                                                                                            >  1:31055/1‑142 (MQ=255)                                                                                                                                                                                                  aCGTGCTTTAGCGCAGCTTCAATACCCATTCCCGTTTGTAAATTAACAGCTACTAAATCAATGAAACCTGCTAAATCATTCATGATTTTTTTAATTTTGGATTTTAAAATAGAACCAATTAAAATTCCTGGAATAACAATAGTTACTATTAATATAg                                                            <  1:163166/157‑1 (MQ=255)                                                                                                                                                                                                                 gCTTCAATACCCATTCCCGTTTGTAAATTAACAGCTACTAAATCAATGAAACCTGCTAAATCATTCATGAtt                                                                                                                                  >  1:370823/1‑72 (MQ=255)                                                                                                                                                                                                                         aCCCATTCCCGTTTGTAAATTAACAGCTACTAAATCAATGAAACCTGCTAAAt                                                                                                                                             >  1:282858/1‑53 (MQ=255)                                                                                                                                                                                                                                       gTAAATTAACAGCTACTAAATCAATGAAACCTGCTAAATCATTCATGATTTTTTTAATTTTGGATTTTAAAATAGAACCAATTAAAATTCCTGGAATAACAATAGTTACTATTAAtata                                                             >  1:424962/1‑119 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                        TAGTTGATGATAAATAGAGGAACCAAAATTAAGGCTTTGTTGCATTACAGTACAAAACATTCGGATTTCTGCTGTTGGTAAAGACACAGAGAGATCCTGCAATGCCTGCGACATACCGGTAATCTCGGATTTACGAATAATACGTAACATAACATAAGTCAGATCAGGGTTTAATATTTTAAAGTCAGTTGCAACGTGCTTTAGCGCAGCTTCAATACCCATTCCCGTTTGTAAATTAACTGCTACTAAATCAATGAAACCTGCTAAATCATTCATGATTTTTTTAATTTTGGATTTTAAAATAGAACCAATTAAAATTCCTGGAATAACAATAGTTACTATTAATATAGTTAGAAAGCCAAGTAACAGAGATGAGCTGTCAGTAGTTTTATCAATAAGTGAATAAAG  >  1/855250‑855657

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