Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,740,470	A→T	19.6%	intergenic (+35/+36)	hutF → / ← PP_5037	formiminoglutamate deiminase/lipocalin family lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,740,470	0	A	T	19.6%	64.0 / 10.6	37	intergenic (+35/+36)	hutF/PP_5037	formiminoglutamate deiminase/lipocalin family lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (8/20);  new base T (3/3);  total (14/23)
Fisher's exact test for biased strand distribution p-value = 3.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.84e-01

CGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGC‑GTTGCAGCCC‑TTCGCGTGCAAGCCCGCGAAGGGTTGCAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTG  >  NC_002947/5740328‑5740618                                                                                                                                                 |                                                                                                                                                      cgcgATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGTGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGC‑GTTGCAGCCC‑TTCGCGTGCAAGccc                                                                                                                                                  <  3:152467/148‑1 (MQ=255)    gATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGC‑GTTGCAGCCC‑TTCGCGTGCAAGc                                                                                                                                                    <  1:144456/143‑1 (MQ=255)    gATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGC‑GTTGCAGCCC‑TTCGCGTGCAAGc                                                                                                                                                    >  2:144456/1‑143 (MQ=255)     aTGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGC‑GTTGCAGCCC‑TTCGCGTGCAAGCCcgcg                                                                                                                                               >  6:162650/1‑147 (MQ=255)     aTGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGC‑GTTGCAGCCC‑TTCGCGTGCAAGCCCGCGaa                                                                                                                                             <  3:40765/149‑1 (MQ=255)      tgtgATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGC‑GTTGCAACCC‑TTCGCGGGCTTGCACGCGAAg                                                                                                                                            <  1:63419/149‑1 (MQ=39)            tgTGAATTGGCTATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTTCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGC‑GTTGCAGCCC‑TTCGCGTGCAAGCCCGCTAAg                                                                                                                                            <  2:228284/142‑1 (MQ=255)            ggTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGC‑GGTGGAGCCC‑ATAGCGTGGAAGCCCGCGCAGGGttt                                                                                                                                       >  7:193959/1‑147 (MQ=255)            gctGAATGGGCAATGGGTGGTGCGCCAGGGGCGTCATGCACAGGAAGATGAAAGCGGGCGGGCCTTTGCCAAGCTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGCGTTGC‑GTTGCAGCCC‑TTCGCGTGCAAGCCCGAGAAg                                                                                                                                            >  1:228284/3‑143 (MQ=255)               gAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGC‑GTTGCAGCCC‑TTCGCGTGGGAGCCCGCGCAGGGGTggaa‑c                                                                                                                                  >  7:241067/1‑149 (MQ=255)                               gtgCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGC‑GTTGCAGCCC‑TTCGCGTGCAAGCCCGCGAAGGGTTGCAA‑CGCAACCCCAAAAAACt                                                                                                                  <  3:188151/149‑1 (MQ=255)                                 gcgcCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGGTTGC‑GCCC‑TTAGCGGGGGAGCCCGCGCAGGGGTTGAA‑CGCCACCCCaaaaaaaa                                                                                                                  >  3:281431/1‑145 (MQ=255)                                 gcgcCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGC‑GTTTGAGCCC‑TTAGCGTGGCAGCCCGCGAAGGGTTGGAA‑CGCAACCCCCAAAAAca                                                                                                                  >  1:91901/1‑146 (MQ=255)                                     ctggggggGAAGGCCCTGTAAGTTCAAGGGGGGGGGGGGTTTGCCTTGGTGTGGCGCCTGTTTTTGGGGTGTGCGTAGCGGGGTGGG‑GTTGCAGCCT‑TTCGGGTCCAAGCCCGCTATGGTTTTAAA‑CGCAACCCCAAAAAACTCAGTGc                                                                                                            <  8:116320/147‑1 (MQ=255)                                       ggggCGGCATGCCCAGTATGTTTAACGCGGGGGGGGTTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGCGGTTACGGGGGTTGG‑GTTGCAGCCC‑TTCGCGTGCAAGCCCGCTATGGGTTGAAA‑CGCAACCCCAAAAAACTCAGTGCAt                                                                                                          <  4:281431/149‑1 (MQ=255)                                             gcATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGC‑GTTGCAGCCC‑TTCGCGTGCAAGCCCGCGAAGGGTTgcaa‑c                                                                                                                                  <  7:72293/119‑1 (MQ=255)                                             gcATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGC‑GTTGCAGCCC‑TTCGCGTGCAAGCCCGCGAAGGGTTgcaa‑c                                                                                                                                  >  8:72293/1‑119 (MQ=255)                                                    cAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGC‑GTTGCAGCCC‑TTCGCGTGCAAGCCCGCTAAGGGTTGCAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGc                                                                                             <  3:70629/149‑1 (MQ=255)                                                         aGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGC‑GTTCCAGCCC‑TTCGCGTGCAAGCCCGCGAAGGGTTGCAAGCGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGAt                                                                                         <  1:220370/149‑1 (MQ=255)                                                             gAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGG‑GGTGGAGCCCATACGCG‑GGAAGGCCGCGCAGGGTTTGAA‑CGCCAACCCAAAAAAAGCAGTGGAGTTTTAAATTTTGATGATggggc                                                                                    >  3:26918/1‑147 (MQ=255)                                                                     gcgggcgTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGC‑GTTGCAGCCC‑TTCGCGTGCAAGCCCGCGAAGGGTTGCAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCa                                                                            <  5:162650/149‑1 (MQ=255)                                                                     gcgggcgTTGGCCAAGGTGTTCCGCCAGCTTTTGGGCTGCGCGTACCGGGGTTGC‑GTTGCAGCCC‑TGCGCGTGCAAGCCCGCTAAGGTTTGAAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCa                                                                            <  4:113214/149‑1 (MQ=255)                                                                         gcgTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGC‑GTTGCAGCCC‑TTCGCGGGCAAGCCCGCGAAGGGTTGCAA‑CGCAACCCCAAAAAACt                                                                                                                  <  5:150461/107‑1 (MQ=255)                                                                           gTTTGCCAAGGTGTGGCGCCAGCTTTTGGGCTTAGCGTACCGGGGTGGC‑GTTGCAGCCC‑TTCGCGTGCAAGCCCGCGAAGGGTTGCAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCaa                                                                      <  3:281810/149‑1 (MQ=255)                                                                              tGCCAAGGGTGTGCGCCCGCTTTTGGGCTGAGCGTACCGGGGTTGC‑TTTGCACCCC‑TTCGCGTGCAAGCCCGCTATGGGTTGCAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGc                                                                    <  5:176933/148‑1 (MQ=255)                                                                                        ttGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGC‑GTTGCAGCCC‑TTCGCGTGCAAGCCCGCGAAGGGTTGCAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTc                                                          <  1:6861/148‑1 (MQ=255)                                                                                             ccAGCTTTTGGGCTGAGCGTACCGGGGTTGC‑GTTGCAGCCC‑TTCGCGCGCAACCCCGCGACGGGTTGCAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATAccc                                                    <  4:26918/149‑1 (MQ=255)                                                                                                            aGCGTACCGGGGTTGC‑GTTGCAGCCC‑TTCGCGTGCAAGCCCGCGAAGGGTTGCAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGc                                                                    <  3:113007/118‑1 (MQ=255)                                                                                                                     gggTTGC‑GTTGCAGCCC‑TTCGCGTGCAAGCCCGCGAAGGGTTGCAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGAcagcag                            <  2:91901/149‑1 (MQ=255)                                                                                                                     gggTTGC‑GTTCCAGCCC‑TTCGCGTGCAACCCCGCTATGGGTTGAAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGAcagcag                            <  3:145409/149‑1 (MQ=255)                                                                                                                             gttgcAGCCC‑TTCGCGTGCAAGCCCGCGAAGGGTTGCAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCa                     >  5:157943/1‑149 (MQ=255)                                                                                                                                  aaccc‑TTCGCGGGCTTGCACGCGAAGGGCTGCAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCa                                                                       >  7:61840/3‑94 (MQ=21)                                                                                                                                  aaccc‑TTCGCGGGCTTGCACGCGAAGGGCTGCAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCa                                                                       <  8:61840/92‑1 (MQ=21)                                                                                                                                        ttCGCGTGCAAGCCCGCGAAGGGTTGCAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAgg           >  1:84715/1‑149 (MQ=255)                                                                                                                                        ttCGCGGGCTTGCACGCGAAGGGCTGCAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAgg                           >  5:186865/1‑133 (MQ=38)                                                                                                                                        ttCGCGGGCTTGCACGCGAAGGGCTGCAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAgg                           <  6:186865/133‑1 (MQ=38)                                                                                                                                        ttCGCGGGCTTGCACGCGAAGGGCTGCAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAgg           >  7:136618/1‑149 (MQ=39)                                                                                                                                                 aaGCCCGCGAAGGGTTGCAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAgtg  >  2:152333/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                      CGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGC‑GTTGCAGCCC‑TTCGCGTGCAAGCCCGCGAAGGGTTGCAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTG  >  NC_002947/5740328‑5740618

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