Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,740,467	T→G	28.0%	intergenic (+32/+39)	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,467	0	T	G	28.0%	56.4 / 33.6	50	intergenic (+32/+39)	hutF/PP_5037	formiminoglutamate deiminase/lipocalin family lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (13/23);  new base G (10/4);  total (23/27)
Fisher's exact test for biased strand distribution p-value = 3.14e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.72e-01

GCCAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAG‑CCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGG  >  NC_002947/5740320‑5740613                                                                                                                                                    |                                                                                                                                                    gCCAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGt                                                                                                                                                     <  5:155145/148‑1 (MQ=255)  ccAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGGGc                                                                                                                                                   >  7:139545/1‑149 (MQ=255)      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:255240/148‑1 (MQ=255)                tgATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAG‑CCCGCGAAgg                                                                                                                                     <  8:285516/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                                                                                                                                 <  3:333909/142‑1 (MQ=255)                          tGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAG‑CCCGCGAAGGGTTg                                                                                                                                 >  4:333909/1‑142 (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                 <  5:401799/149‑1 (MQ=255)                                                                                                                                   cgttgcAGCCCTTCGCGTGCAAG‑CCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTtctc                 <  1:15111/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)                                                                                                                                                    |                                                                                                                                                    GCCAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAG‑CCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGG  >  NC_002947/5740320‑5740613

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