Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,740,470	A→T	23.0%	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	23.0%	88.3 / 25.0	52	intergenic (+35/+36)	hutF/PP_5037	formiminoglutamate deiminase/lipocalin family lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (15/25);  new base T (8/4);  total (23/29)
Fisher's exact test for biased strand distribution p-value = 1.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.49e-01

GTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAG‑CCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAG  >  NC_002947/5740325‑5740616                                                                                                                                                  |                                                                                                                                                    gTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAG‑c                                                                                                                                                 <  4:347103/149‑1 (MQ=255)    cgcgATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAG‑ccc                                                                                                                                               <  7:347572/148‑1 (MQ=255)        aTGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAG‑CCCGCGaa                                                                                                                                          <  5:119309/149‑1 (MQ=255)         tgtgATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAG‑CCCGCGAAg                                                                                                                                         <  3:26436/149‑1 (MQ=255)          gtgATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGGGCTTG‑CACGCGAAg                                                                                                                                         <  8:139545/148‑1 (MQ=255)           tgaggctGAATGGGCATTGGGTGGCGCGCCATGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAG‑CCCGCGAAgg                                                                                                                                        <  8:264309/142‑1 (MQ=255)           tgATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAG‑CCCGCGAAgg                                                                                                                                        <  8:285516/148‑1 (MQ=255)           tgATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAG‑CCCGCGAAgg                                                                                                                                        <  8:255240/148‑1 (MQ=255)           tgATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAG‑CCCGCGAAgg                                                                                                                                        <  3:20407/148‑1 (MQ=255)             aTGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAG‑CCCGCGAAGGGtt                                                                                                                                     <  3:278666/149‑1 (MQ=255)              tGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGCTTG‑CACGCGAAGGGCTg                                                                                                                                    >  4:19698/1‑149 (MQ=38)               ggTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGCTTG‑CACGCGAAGGGCTg                                                                                                                                    <  3:19698/148‑1 (MQ=38)                gTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGCTTG‑CACGCGAAGGGCTgc                                                                                                                                   <  5:179484/148‑1 (MQ=38)                     tGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAG‑CCCGCGAAGGGTTg                                                                                                                                    >  4:333909/1‑142 (MQ=255)                     tGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAG‑CCCGCGAAGGGTTg                                                                                                                                    <  3:333909/142‑1 (MQ=255)                               gtggtgCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAG‑CCCGCGAAGGGTTGCAACGCAACCCCaaaa                                                                                                                    >  5:418315/1‑148 (MQ=255)                                     cgcCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAG‑CCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAg                                                                                                             <  3:49634/149‑1 (MQ=255)                                      gcCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGCGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAG‑CCCGCTAAGGGTTGAAACGCAACCCCAAAAAACTCAGt                                                                                                            <  4:327970/149‑1 (MQ=255)                                                        aGGAAGATGAAAGCGGGCGGGCGTTGGCCAAGGTGTTGCGCCAGCTTTTGGGCGGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAG‑CCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCt                                                                                          <  4:313844/149‑1 (MQ=255)                                                                  aaGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAG‑CCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGacgc                                                                                >  6:418315‑M2/1‑145 (MQ=255)                                                                    gcgggcggGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTTGAACCCATACCGGGGAAGCCCCCCCAAGGGTTG‑AACCCAAACCCCAAAAAAAAAATGCATTTTTTATATTTTGTTATCGGAcgccc                                                                              >  2:15111/1‑149 (MQ=255)                                                                      gggggggcGTTTGCCATGGTGTTGCGCCAGCTTTTGGGCGGAGCGTACCGGGGTGGCGTTGCAGCCCTTCGCGTGCAAG‑CCCGCGATGGGTTGAAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGacgcccgc                                                                            >  4:298328‑M2/1‑141 (MQ=255)                                                                                      aGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAG‑CCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTgg                                                            <  3:289787/149‑1 (MQ=255)                                                                                             gcgcCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAG‑CCCGCGAAGGGTTgc                                                                                                                                   >  5:223800/1‑71 (MQ=255)                                                                                             gcgcCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAG‑CCCGCGAAGGGTTgc                                                                                                                                   <  6:223800/71‑1 (MQ=255)                                                                                               gcCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTTCAGCCCTTAGCGTGCAAG‑CCCGCCAAGGGTTGCAACGCCACCCCAAAAAAAGAAGTGGATTTTTTATATTTGGTGATAGGACGCCCGCCACATCAAGCGGCGGGTGTAATaac                                                   >  6:401799/1‑147 (MQ=255)                                                                                               gcCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTAGCGTGCAAG‑CCCGCGAACGGTTGCAACCCAACCCCAAAAAACTCACTGCAATTTTACTATTTGCGGATCGGACGCCCCCCAGATCAAGCGGAAGGTGTCATAcc                                                   >  7:386471/1‑149 (MQ=255)                                                                                                     ttttGGGCGGAGCGTACCGGGGTGGCGTTGCAGCCCTTCGCGTGCAAG‑CCCGCTATGGGTTGAAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCctgct                                             <  2:361478/149‑1 (MQ=255)                                                                                                          ggCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAG‑CCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCcg                                         <  1:113586/148‑1 (MQ=255)                                                                                                               aGCTTACCGGGGTTGCGTTGCAGCCGTGGGCGTGCATC‑CCCGGGAAGGTTTGAAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGcc                                     <  5:26406/147‑1 (MQ=255)                                                                                                                 cGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAG‑CCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCt                                   <  1:82444/147‑1 (MQ=255)                                                                                                                              ctttgcAGCCCTTCGCGTGCAAG‑CCCGCTAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTtctc                    <  2:127287/147‑1 (MQ=255)                                                                                                                              cgttgcAGCCCTTCGCGTGCAAG‑CCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTtctc                    <  1:15111/149‑1 (MQ=255)                                                                                                                              cgttgcAGCCCTTCGCGTGCAAG‑CCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTtctc                    <  5:401799/149‑1 (MQ=255)                                                                                                                               gttgcAGCCCTTGGCGTGCAAG‑CCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTTCTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCTCCCTGGACAGCAGGTTCTCa                   <  8:386471/149‑1 (MQ=255)                                                                                                                                     gCCCTTCGCGTGCAAG‑CCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTcc             >  8:63228/1‑149 (MQ=255)                                                                                                                                     acccTTCGCGGGCTTG‑CACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTcc             >  5:134336/2‑149 (MQ=38)                                                                                                                                     acccTTCGCGGGCTTG‑CACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTcc             >  3:264582/2‑149 (MQ=38)                                                                                                                                       ccTTCGCGTGCAAG‑CCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTcc             >  3:202247/1‑147 (MQ=255)                                                                                                                                       ccTTCGCGTGCAAG‑CCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACACCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTcc             >  3:359135/1‑147 (MQ=255)                                                                                                                                       ccTTCGCGGGCTTG‑CACGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTcc             >  8:234595/1‑147 (MQ=255)                                                                                                                                       ccTTCGCGGGCTTG‑CACGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCAGGACAGCAGGTTCTCACGCTcc             >  2:53415/1‑147 (MQ=39)                                                                                                                                        cTTCGCGGGCTTG‑CACGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAg          >  5:118184/1‑149 (MQ=255)                                                                                                                                        cTTCGCGGGCTTG‑CACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAg          >  2:304780/1‑149 (MQ=39)                                                                                                                                         ttCGCGTGCAAG‑CCCGCGAAGGGTTGCAACGCAACCACAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATCCCCCTGCTGCCGCGccc                                    >  2:67209/1‑122 (MQ=255)                                                                                                                                         ttCGCGGGCTTG‑CACGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGAcagca                           <  7:207271/131‑1 (MQ=255)                                                                                                                                         ttCGCGGGCTTG‑CACGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGAcagca                           >  8:207271/1‑131 (MQ=255)                                                                                                                                              gTGCAAG‑CCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCgg     >  6:199676/1‑148 (MQ=255)                                                                                                                                              gTGCAAG‑CCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCgg     >  6:126152/1‑148 (MQ=255)                                                                                                                                                gCAAG‑CCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGacgcccgcccgctttcatcttcctgggcatgccgcccctggcgcaccacccattgcccattcaccatcacatcgcgcacctg  <  6:7229‑M2/149‑83 (MQ=255)                                                                                                                                                gCAAG‑CCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGacgcccgcccgctttcatcttcctgggcatgccgcccctgg                                           <  5:43600‑M2/108‑42 (MQ=255)                                                                                                                                                gCAAG‑CCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGacgcccgcccgctttcatcttcctgggcatgccgcccctgg                                           >  6:43600‑M2/1‑67 (MQ=255)                                                                                                                                                  |                                                                                                                                                    GTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAG‑CCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAG  >  NC_002947/5740325‑5740616

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: ‑