Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,487,693	2 bp→GT	coding (225‑226/603 nt)	unknown →	Flp pilus assembly surface protein TadF, ATP/GTP‑binding motif

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,487,693	0	A	G	100.0%	96.5 / NA	34	Q75Q (CAA→CAG)	unknown	Flp pilus assembly surface protein TadF, ATP/GTP‑binding motif
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (22/12);  total (22/12)
*	1	1,487,694	0	C	T	97.1%	81.2 / ‑4.5	34	L76L (CTA→TTA)	unknown	Flp pilus assembly surface protein TadF, ATP/GTP‑binding motif
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (21/12);  total (22/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

ATAAAAGGGGCGCAGTTAGTCTTGAGTTCCTTTTTATGCTTATACTTTTAGTATTTATTTTCGCTTTTTTAACAGACTTAGTTATTGTTCGTACGACCCAAGGGAAACTGGACAATGCCTCATACTCACTGGTCAATGTTCTTCGTGAACGTAATCAACTATACAATGATAATGGGAGCAAACCACTAAAAAGTACTGATTTAACTGAATATGAAAAAATGGCAAAACTTATTTTGTTTGGAGATAAAGACAGCCCAAATAAAGTCGGTGTCACTATTGAGCACTGGGAAGAGAAAAAAAGTCCTGAATCATTAACGACTCTAAGTACCTGCCAGCCTTATCGTAAGTTGGATAATAATCTCTCTTATTTGTCACCACGCTCAGAAGCAGCAAACC  >  1/1487536‑1487931                                                                                                                                                              ||                                                                                                                                                                                                                                              ataaAAAGGGCGCAGTTACGCTTGAGCTTCTTTTTATGCTTATTCTTTTAGTATTTATTTTCGCCTTTTTAACAGACCTAGTTATTGTTCGTACTACCCAGGGAAAACTGGACAATGCCTCGTATTCACTGGTCAATATTCTTCGCGAGCGTAATCAGTTATACAATGATAATGGGACCGAAAAACTAAAAAGTACTGACTTAACTGAATATGAAAAAATGGCAAAACTTATTTTGTTTGGAGATAAAGAc                                                                                                                                                   <  1:884263/251‑1 (MQ=255)                   cgcTTGAGCTTCTTTTTATGCTTATTCTTTTAGTATTTATTTTCGCCTTTTTAACAGACCTAGTTATTGTTCGTACTACCCAGGGAAAACTGGACAATGCCTCGTATTCACTGGTCAATATTCTTCGCGAGCGTAATCAGTTATACAATGATAATGGGACCGAAAAACTAAAAAGTACTGACTTAACTGAATATGAAAAAATGGCAAAACTTATTTTGTTTGGAGATAAAGACAGTCCAAATAAAGTCGgt                                                                                                                                 >  1:779644/3‑251 (MQ=255)                             tcTTTTTATGCTTATTCTTTTAGTATTTATTTTCGCCTTTTTAACAGACCTAGTTATTGTTCGTACTACCCAGGGAAAACTGGACAATGCCTCGTATTCACTGGTCAATATTCTTCGCGAGCGTAATCAGTTATACAATGATAATGGGACCGAAAAACTAAAAAGTACTGACTTAACTGAATATGAAAAAATGGCAAAACTTATTTTGTTTGGAGATAAAGAc                                                                                                                                                   <  1:459000/222‑1 (MQ=255)                              cTTTTTATGCTTATTCTTTTAGTATTTATTTTCGCCTTTTTAACAGACCTAGTTATTGTTCGTACTACCCAGGGAAAACTGGACAATGCCTCGTATTCACTGGTCAATATTCTTCGCGAGCGTAATCAGTTATACAATGATAATGGGACTGAAAAACTAAAAAGTACTGACTTAACTGAATATGAAAAAATGGCAAAACTTATTTTGTTTGGAGATAAAGACAGTCCAAATAAAGTCGGTGTCACTATTGa                                                                                                                      >  1:67448/1‑251 (MQ=255)                                ttttATGCTTATTCTTTTAGTATTTATTTTCGCCTTTTTAACAGACCTAGTTATTGTTCGTACTACCCAGGGAAAACTGGATAATGCCTCGTATTCACTGGTCAATATTCTTCGCGAGCGTAATCAGTTATACAATGATAATGGGACCGAAAAACTAAAAAGTACTGACTTAACTGAATATGAAAAAATGGCAAAACTTATTTTGTTTGGAGATAAAGACAGTCCAAATAAAGTCGGTGTCACTATTGAGc                                                                                                                    <  1:442647/251‑1 (MQ=255)                                           ttcTTTTAGTATTTATTTTCGCCTTTTTAACAGACCTAGTTATTGTTCGTACTACCCAGGGAAAACTGGACAATGCCTCGTATTCACTGGTCAATATTCTTCGCGAGCGTAATCAGTTATACAATGATAATGGGACCGAAAAACTAAAAAGTACTGACTTAACTGAATATGAAAAAATGGCAAAACTTATTTTGTTTGGAGATAAAGACAGTCCAAATAAAGTCGGTGTCACTATTGAGCACTGGGAagag                                                                                                         >  1:779739/3‑251 (MQ=255)                                           ttcTTTTAGTATTTATTTTCGCCTTTTTAACAGACCTAGTTATTGTTCGTACTACCCAGGGAAAACTGGACAATGCCTCGTATTCACTGGTCAATATTCTTCGCGAGCGTAATCAGTTATACAATGATAATGGGACCGAAAAACTAAAAAGTACTGACTTAACTGAATATGAAAAAATGGCAAAACTTATTTTGTTTGGAGATAAAGACAGTCCAAATAAAGTCGGTGTCACTATTGAGCACTGGGAagag                                                                                                         >  1:574359/3‑251 (MQ=255)                                                   gTATTTATTTTCGCCTTTTTAACAGACCTAGTTATTGTTCGTACTACCCAGGGAAAACTGGACAATGCCTCGTATTCACTGGTCAATATTCTTCGCGAGCGTAATCAGTTATACAATGATAATGGGACCGAAAAACTAAAAAGTACTGACTTAACTGAATATGAAAAAATGGCAAAACTTATTTTGTTTGGAGATAAAGACAGTCCAAATAAAGTCGGTGTCACTAt                                                                                                                         <  1:147517/227‑1 (MQ=255)                                                                cctttttAACAGACCTAGTTATTGTTCGTACTACCCAGGGAAAACTGGACAATGCCTCGTATTCACTGGTCAATATTCTTCGCGAGCGTAATCAGTTATACAATGATAATGGGACCGAAAAACTAAAAAGTACTGACTTAACTGAATATGAAAAAATGGCAAAACTTATTTTGTTTGGAGATAAAGAc                                                                                                                                                   >  1:281579/3‑188 (MQ=255)                                                                       aaCAGACCTAGTTATTGTTCGTACTACCCAGGGAAAACTGGACAATGCCTCGTATTCACTGGTCAATATTCTTCGCGAGCGTAATCAGTTATACAATGATAATGGGACCGAAAAACTAAAAAGTACTGACTTAACTGAATATGAAAAAATGGCAAAACTTATTTTGTTTGGAGATAAAGACAGTCCAAATAAAGTCGGTGTCACTATTGAGCACTGGGAAGAGAAAAAAAATCCTGaaatgttaacgaatc                                                                             <  1:651093/251‑14 (MQ=255)                                                                           gacctAGTTATTGTTCGTACTACCCAGGGAAAACTGGACAATGCCTCGTATTCACTGGTCAATATTCTTCGCGAGCGTAATCAGTTATACAATGATAATGGGACCGAAAAACTAAAAAGTACTGACTTAACTGAATATGAAAAAATGGCAAAACTTATTTTGTTTGGAGATAAAGACAGTCCAAATaa                                                                                                                                        <  1:892138/184‑1 (MQ=255)                                                                           gacctAGTTATTGTTCGTACTACCCAGGGAAAACTGGACAATGCCTCGTATTCACTGGTCAATATTCTTCGCGAGCGTAATCAGTTATACAATGATAATGGGACCGAAAAACTAAAAAGTACTGACTTAACTGAATATGAAAAAATGGCAAAACTTATTTTGTTTGGAGATAAAGACAGTCCAAATAAAGTCGGTGTCACTATTg                                                                                                                       >  1:664346/5‑205 (MQ=255)                                                                               tAGTTATTGTTCGTACTACCCAGGGAAAACTGGACAATGCCTCGTATTCACTGGTCAATATTCTTCGCGAGCGTAATCAGTTATACAATGATAATGGGACCGAAAAACTAAAAAGTACTGACTTAACTGAATATGAAAAAATGGCAAAACTTATTTTGTTTGGAGATAAAGAc                                                                                                                                                   <  1:589346/173‑1 (MQ=255)                                                                               tAGTTATTGTTCGTACTACCCAGGGAAAACTGGACAATGCCTCGTATTCACTGGTCAATATTCTTCGCGAGCGTAATCAGTTATACAATGATAATGGGACCGAAAAACTAAAAAGTACTGACTTAACTGAATATGAAAAAATGGCAAAACTTATTTTGTTTGGAGATAAAGACAGTCCAAATAAAGTCGGTGTCATTATTGAGCACTGGGAAGAGAAAAAAAATCCTGaaatgttaacgaatctaggcaca                                                                     >  1:720347/1‑230 (MQ=255)                                                                                aGTTATTGTTCGTACTACCCAGGGAAAACTGGACAATGCCTCGTATTCACTGGTCAATATTCTTCGCGAGCGTAATCAGTTATACAATGATAATGGGACCGAAAAACTAAAAAGTACTGACTTAACTGAATATGAAAAAATGGCAAAACTTATTTTGTTTGGAGATAAAGACAGTCCAAATAAAGTCGGTGTCACTATTGAGCACTGGGAAGAGAAAAAAAATCCTGaaatgttaacgaa                                                                               >  1:524724/1‑229 (MQ=255)                                                                                 gTTATTGTTCGTACTACCCAGGGAAAACTGGACAATGCCTCGTATTCACTGGTCAATATTCTTCGCGAGCGTAATCAGTTATACAATGATAATGGGACCGAAAAACTAAAAAGTACTGACTTAACTGAATATGAAAAAATGGCAAAACTTATTTTGTTTGgaga                                                                                                                                                          <  1:158211/164‑1 (MQ=255)                                                                                 gTTATTGTTCGTACTACCCAGGGAAAACTGGACAATGCCTCGTATTCACTGGTCAATATTCTTCGCGAGCGTAATCAGTTATACAATGATAATGGGACCGAAAAACTAAAAAGTACTGACTTAACTGAATATGAAAAAATGGCAAAACTTATTTTGTTTGGAGATAAAGGCAGTCCAAATAAAGTCGGTGTCACTATTGAGCACTGGGAAGAGaaaaaaaa                                                                                                 >  1:431347/1‑220 (MQ=255)                                                                                           gTACTACCCAGGGAAAACTGGACAATGCCTCGTATTCACTGGTCAATATTCTTCGCGAGCGTAATCAGTTATACAATGATAATGGGACCGAAAAACTAAAAAGTACTGACTTAACTGAATATGAAAAAATGGCAAAACTTATTTTGTTTGGAGATAAAGAc                                                                                                                                                   <  1:360430/161‑1 (MQ=255)                                                                                                aCCCAGGGAAAACTGGACAATGCCTCGTATTCACTGGTCAATATTCTTCGCGAGCGTAATCAGTTATACAATGATAATGGGACCGAAAAACTAAAAAGTACTGACTTAACTGAATATGAAAAAATGGCAAAACTTATTTTGTTTGGAGATAAAGAc                                                                                                                                                   >  1:577717/1‑156 (MQ=255)                                                                                                aCCCAGGGAAAACTGGACAATGCCTCGTATTCACTGGTCAATATTCTTCGCGAGCGTAATCAGTTATACAATGATAATGGGACCGAAAAACTAAAAAGTACTGACTTAACTGAATATGAAAAAATGGCAAAACTTATTTTGTTTGGAGATAAAGACAGTCCAAATAAAGTCGGTGTc                                                                                                                              >  1:682741/1‑177 (MQ=255)                                                                                                aCCCAGGGAAAACTGGACAATGCCTCGTATTCACTGGTCAATATTCTTCGCGAGCGTAATCAGTTATACAATGATAATGGGACCGAAAAACTAAAAAGTACTGACTTAACTGAATATGAAAAAATGGCAAAACTTATTTTGTTTGGAGATAAAGACAGTCCAAATAAAGTCGGTGTCACTATTGAGCACTGGGAAGAGAAAAAAAATCCTGAAATGTTAACGAATCTAGGCACATGTCGGCCTTATCGTaa                                                    >  1:254629/1‑251 (MQ=255)                                                                                                  ccagggaaaaCTGGACAATGCCTCGTATTCACTGGTCAATATTCTTCGCGAGCGTAATCAGTTATACAATGATAATGGGACCGAAAAACTAAAAAGTACTGACTTAACTGAATATGAAAAAATGGCAAAACTTATTTTGTTTGGAGATAAAGACAGTCCAAATAAAGTCGGTGTCACTATTGAGCACTGGGAAGAGAAAAAAAATCCTGAAATGTTAACGAATCTAGGCACATGTCGGCCTTATCGTAAGt                                                  <  1:276791/244‑1 (MQ=255)                                                                                                         aaaCTGGACAATGCCTCGTATTCACTGGTCAATATTCTTCGCGAGCGTAATCAGTTATACAATGATAATGGGACCGAAAAACTAAAAAGTACTGACTTAACTGAATATGAAAAAATGGCAAAACTTATTTTGTTTGGAGATAAAGACAGTCCAAATAAAGTCGGTGTCACTATTGAGCACTGGGAAGAGAAAAAAAATCCTGAAATGTTAACGAATCTAGGCACATGTCGGCCTTATCGTAAGTTGGatga                                           >  1:675512/1‑249 (MQ=255)                                                                                                         aaaCTGGACAATGCCTCGTATTCACTGGTCAATATTCTTCGCGAGCGTAATCAGTTATACAATGATAATGGGACCGAAAAACTAAAAAGTACTGACTTAACTGAATATGAAAAAATGGCAAAACTTATTTTGTTTGGAGATAAAGACAGTCCAAATAAAGTCGGTGTCACTATTGAGCACTGGGAAGAGAAAAAAAATCCTGAAATGTTAACGAATCTAGGCACATGTCGGCCTTATCGTAAGTTGGatga                                           >  1:18778/1‑249 (MQ=255)                                                                                                           aCTGGACAATGCCTCGTATTCACTGGTCAATATTCTTCGCGAGCGTAATCAGTTATACAATGATAATGGGACCGAAAAACTAAAAAGTACTGACTTAACTGAATATGAAAAAATGGCAAAACTTATTTTGTTTGGAGATAAAGACAGTCCAAATAAAGt                                                                                                                                     >  1:568382/1‑159 (MQ=255)                                                                                                                     gCCTCGTATTCACTGGTCAATATTCTTCGCGAGCGTAATCAGTTATACAATGATAATGGGACCGAAAAACTAAAAAGTACTGACTTAACTGAATATGAAAAAATGGCAAAACTTATTTTGTTTGGAGATAAAGACAGTCCAAATa                                                                                                                                         >  1:33869/1‑145 (MQ=255)                                                                                                                     gCCTCGTATTCACTGGTCAATATTCTTCGCGAGCGTAATCAGTTATACAATGATAATGGGACCGAAAAACTAAAAAGTACTGACTTAACTGAATATGAAAAAATGGCAAAACTTATTTTGTTTGGAGATAAAGACAGTCCAAATAAAGTCGGTGTCACTATTGAGCACTGGGAAGAGAAAAAAAATCCTGAAATGTTAACGAATCTAGGCACATGTCGGCCTTATCGTAAGTTGGATGATAATCTATCTTa                               >  1:794139/1‑251 (MQ=255)                                                                                                                         cgtattCACTGGTCAATATTCTTCGCGAGCGTAATCAGTTATACAATGATAATGGGACCGAAAAACTAAAAAGTACTGACTTAACTGAATATGAAAAAATGGCAAAACTTATtttg                                                                                                                                                                  <  1:790982/111‑1 (MQ=255)                                                                                                                         cgtattCACTGGTCAATATTCTTCGCGAGCGTAATCAGTTATACAATGATAATGGGACCGAAAAACTAAAAAGTACTGACTTAACTGAATATGAAAAAATGGCAAAACTTATTTTGTTTGGAGATAAAGACAGTCCAAATAAAGTCGGTGTCACTATTGAGCACTGGGAAGAGAAAAAAAATCCTGaaa                                                                                         >  1:126619/6‑188 (MQ=255)                                                                                                                                    gTCAATATTCTTCGCGAGCGTAATCAGTTATACAATGATAATGGGACCGAAAAACTAAAAAGTACTGACTTAACTGAATATGAAAAAATGGCAAAACTTATTTTGTTTGGAGATAAAGACAGTCCAAATAAAGTCGGTGTCACTATTGAGCACTGGGAAGAGAAAAAAAATCCTGAAATGTTAACGAATCTAGGCACATGTCGGCCTTATCGTAAGTTGGATGATAatctatct                                 <  1:670121/234‑1 (MQ=255)                                                                                                                                        atattcttcGCGAGCGTAATCAGTTATACAATGATAATGGGACCGAAAAACTAAAAAGTACTGACTTAACTGAATATGAAAAAATGGCAAAACTTATTTTGTTTGGAGATAAAGACAGTCCAAATAAAGTCGGTGTCACTATTGAGCACTGGGAAGAGAAAAAAAATCCTGAAATGTTAACGAATCTAGGCACATGTCGGCCTTATCGTAAGTTGGATGATAATCTATCTTAttt                            >  1:225663/4‑235 (MQ=255)                                                                                                                                        atattcttcGCGAGCGTAATCAGTTATACAATGATAATGGGACCGAAAAACTAAAAAGTACTGACTTAACTGAATATGAAAAAATGGCAAAACTTATTTTGTTTGGAGATAAAGACAGTCCAAATAAAGTCGGTGTCACTATTGAGCACTGGGAAGAGAAAAAAAATCCTGAAATGTTAACGAATCTAGGCACATGTCGGCCTTATCGTAAGTTGGATGATAATCTATCTTATTTGTCGCCACGCTCAGaa            >  1:727878/4‑251 (MQ=255)                                                                                                                                                    agcGTAATCAGTTATACAATGATAATGGGACCGAAAAACTAAAAAGTACTGACTTAACTGAATATGAAAAAATGGCAAAACTTATTTTGTTTGGAGATAAAGACAGTCCAAATAAAGTCGGTGTCACTATTGAGCACTGGGAAGAGAAAAAAAATCCTGAAATGTTAACGAATCTAGGCACATGTCGGCCTTATCGTAAGTTGGATGATAATCTATCTTATTTGTCGCCACGCTCAGAATCAGCTAAcc  >  1:385261/3‑249 (MQ=255)                                                                                                                                                      cGTAATCAGCTATACAATGATAATGGGACCGAAAAACTAAAAAGTACTGACTTAACTGAATATGAAAAAATGGCAAAACTTATTTTGTTTGGAGATAAAGACAGTCCAAATAAAGTCGGTGTCACTATTGAGCACTGGGAAGAGAAAAAAAATCCTGAAATGTTAACGAATCTAGGCACATGTCGGCCTTATCGTAAGTTGGATGATAAtc                                      >  1:539316/1‑211 (MQ=255)                                                                                                                                                              ||                                                                                                                                                                                                                                              ATAAAAGGGGCGCAGTTAGTCTTGAGTTCCTTTTTATGCTTATACTTTTAGTATTTATTTTCGCTTTTTTAACAGACTTAGTTATTGTTCGTACGACCCAAGGGAAACTGGACAATGCCTCATACTCACTGGTCAATGTTCTTCGTGAACGTAATCAACTATACAATGATAATGGGAGCAAACCACTAAAAAGTACTGATTTAACTGAATATGAAAAAATGGCAAAACTTATTTTGTTTGGAGATAAAGACAGCCCAAATAAAGTCGGTGTCACTATTGAGCACTGGGAAGAGAAAAAAAGTCCTGAATCATTAACGACTCTAAGTACCTGCCAGCCTTATCGTAAGTTGGATAATAATCTCTCTTATTTGTCACCACGCTCAGAAGCAGCAAACC  >  1/1487536‑1487931

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