Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,260,746	G→C	10.1%	G445A (GGC→GCC)	fimV →	motility protein FimV

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,260,746	0	G	C	10.1%	116.5 / 2.3	38	G445A (GGC→GCC)	fimV	motility protein FimV
Reads supporting (aligned to +/- strand):  ref base G (17/17);  new base C (2/2);  total (19/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCTGGTGCTGGCACTGCTGCTGTGGCTGCTGGCGCGCAAACGCAAGGCCCAGCAAGAAGCCGAGAAGCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAGTGGCAGAACCGGTCGCCGAGGCAGCACCCGAGCCTGAACCACAAGCAGACCC  >  NC_002947/2260605‑2260881                                                                                                                                              |                                                                                                                                        tCCTGGTGCTGGCACTGCTGCTGTGGCTGCTGGCGCGCAAACGCAAGGCCCAGCAAGAAGCCGAGAAGCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGaa                                                                                                                                  <  4:203603/149‑1 (MQ=255)       tGCTGGCACTGCTGCTGTGGCTGCTGGCGCGCAAACGCAAGGCCCAGCAAGAAGCCGAGAAGCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGACGCCGCCGAAGctctc                                                                                                                            <  4:138271/149‑5 (MQ=255)              aCTGCTGCTGTGGCTGCTGGCGCGCAAACGCAAGGCCCAGCAAGAAGCCGAGAAGCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGACGCCGCCGAAGctc                                                                                                                              >  1:146773/1‑138 (MQ=255)              aCTGCTGCTGTGGCTGCTGGCGCGCAAACGCAAGGCCCAGCAAGAAGCCGAGAAGCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGACGCCGCCGAAGctc                                                                                                                              <  2:146773/140‑3 (MQ=255)                 gctgctgTGGCTGCTGGCGCGCAAACGCAAGGCCCAGCAAGAAGCCGAGAAGCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCggcggc                                                                                                                                        >  6:177758/1‑127 (MQ=255)                 gctgctgTGGCTGCTGGCGCGCAAACGCAAGGCCCAGCAAGAAGCCGAGAAGCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCggcggc                                                                                                                                        <  5:177758/127‑1 (MQ=255)                  ctgctgTGGCTGCTGGCGCGCAAACGCAAGGCCCAGCAAGAAGCCGAGAAGCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGgt                                                                                                                   <  3:114479/147‑1 (MQ=255)                        tgGCTGCTGGCGCGCAAACGCAAGGCCCAGCAAGAAGCCGAGAAGCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTTTAAcc                                                                                                             <  7:299663/147‑1 (MQ=255)                        tgGCTGCTGGCGCGCAAACGCAAGGCCCAGCAAGAAGCCGAGAAGCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAAcc                                                                                                             >  8:299663/1‑147 (MQ=255)                            tactgGCGCGCAAACGCAAGGCCCAGCAAGAAGCCGAGAAGCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCg                                                                                                       >  3:71370/3‑149 (MQ=255)                                   ccgcAAACGCAAGGCCCAGCAAGAAGCCGAGAAGCACCTGCGCAACCCGCGGGCGCTGGAAGAAGAGCAGCGCGCAGGTTTCGACAGCGAAATCGAGAGCTTCGGCGCCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGcag                                                                                                >  7:142776/2‑147 (MQ=255)                                          cGCAATGCCCAGCAAGATGCCGAGAAGCACCTGCGGATGGCGCGTGCGCTGGAAGAAGAGCAGGGCGGAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGGCGGCTGTGGTTGc                                                                                         <  8:142776/149‑1 (MQ=255)                                           gCAAGGCCCAGCAAGAAGCCGAGAAGCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGcc                                                                                        >  7:233467/1‑149 (MQ=255)                                            cAAGGCCCAGCAAGAAGCCGAGAAGCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCgg                                                                                                                     <  5:136409/119‑1 (MQ=255)                                            cAAGGCCCAGCAAGAAGCCGAGAAGCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCgg                                                                                                                     >  6:136409/1‑119 (MQ=255)                                             aaGGCCCAGCAAGAAGCCGAGAAGCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGc                                                                                      >  3:71920/1‑149 (MQ=255)                                                  ccAGCAAGAAGCCGAGAAGCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAAcc                                                                                                             >  2:183708/1‑121 (MQ=255)                                                  ccAGCAAGAAGCCGAGAAGCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCgg                                                                                 <  8:30205/149‑1 (MQ=255)                                                          aagCCGAGAAGCACCTGCGCATGGCGCGGGCGATGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCt                                                                         <  8:53132/149‑1 (MQ=255)                                                               gagaAGCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGt                                                                    <  4:71920/149‑1 (MQ=255)                                                                    gCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTc                                                                                                        <  1:118702/108‑1 (MQ=255)                                                                    gCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTc                                                                                                        >  2:118702/1‑108 (MQ=255)                                                                             cATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAg                                                      >  2:278824/1‑149 (MQ=255)                                                                              aTGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAGt                                                     <  8:233467/149‑1 (MQ=255)                                                                                 gcgcGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAGTGGc                                                  <  4:71370/149‑1 (MQ=255)                                                                                               agaagaGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAGTGGCAGAACCGGTCGcc                                     >  4:136507/1‑148 (MQ=255)                                                                                                gaagaGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAGTGGCAGAACCGGTCGcc                                     >  5:186492/1‑147 (MQ=255)                                                                                                      cAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAGTGGCAGa                                               <  5:326105/131‑1 (MQ=255)                                                                                                      cAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAGTGGCAGa                                               >  6:326105/1‑131 (MQ=255)                                                                                                              aGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAGTGGCAGAACCGGTCGCCGAGGCAGCACCCGAGc                     <  6:186492/149‑1 (MQ=255)                                                                                                               ggTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCATCGCTGCGGAAAAAGTGGCAGAACCGGTCGCCGAGGCAGCACCCGAGcc                    <  2:309355/149‑1 (MQ=255)                                                                                                                  ttCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAGTGGCAGAACCGGTCGCCGa                                   <  6:94414/131‑1 (MQ=255)                                                                                                                  ttCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAGTGGCAGAACCGGTCGCCGa                                   >  5:94414/1‑131 (MQ=255)                                                                                                                            atatCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAGTGGCAGAACCGGTCGCCGAGGCAGCACCCGAGCCTGAACCACAAg         >  4:196166/1‑147 (MQ=255)                                                                                                                              atCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAGTGGCAGAACCGGTCGCCGAGGCAGCACCCGAGCCTGAAc               >  7:15550/1‑139 (MQ=255)                                                                                                                              atCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAGTGGCAGAACCGGTCGCCGAGGCAGCACCCGAGCCTGAAc               <  8:15550/139‑1 (MQ=255)                                                                                                                                 gagaGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAGTGGCAGAACCGGTc                                        <  4:16488/111‑1 (MQ=255)                                                                                                                                 gagaGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAGTGGCAGAACCGGTc                                        >  3:16488/1‑111 (MQ=255)                                                                                                                                 gagaGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAGTGGCAGAACCGGTCGCCGAGGCAGCACCCGAGCCTGAACCACAAGCAGAccc  >  2:151216/1‑149 (MQ=255)                                                                                                                                 gagaGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAGTGGCAGAACCGGTCGCCGAGGCAGCACCCGAGCCTGAACCACAAGCAGAccc  >  2:146451/1‑149 (MQ=255)                                                                                                                                              |                                                                                                                                        TCCTGGTGCTGGCACTGCTGCTGTGGCTGCTGGCGCGCAAACGCAAGGCCCAGCAAGAAGCCGAGAAGCACCTGCGCATGGCGCGGGCGCTGGAAGAAGAGCAGGGCGCAGGTTTCGACAGCGATATCGAGAGCTTCGGCGGCGTCGAAGCCCCCGAGCCGGGTGTAACCCTGTCGCCGGCTGTGGTTGCCGCTTCGGCTGCAGCTGCCGTCGCTGCGGAAAAAGTGGCAGAACCGGTCGCCGAGGCAGCACCCGAGCCTGAACCACAAGCAGACCC  >  NC_002947/2260605‑2260881

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