Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,740,457	G→A	19.5%	intergenic (+22/+49)	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,457	0	G	A	19.5%	71.6 / 13.4	42	intergenic (+22/+49)	hutF/PP_5037	formiminoglutamate deiminase/lipocalin family lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (14/19);  new base A (2/6);  total (16/26)
Fisher's exact test for biased strand distribution p-value = 4.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.92e-03

CGGTGGCGATCGCCAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCG‑AAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCC  >  NC_002947/5740309‑5740605                                                                                                                                                     |                                                                                                                                                      cGGTGGCGATCGCCAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAg                                                                                                                                                       <  5:96583/149‑1 (MQ=255)    tGGCGATCGCCAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGccc                                                                                                                                                    <  1:291373/149‑1 (MQ=255)      gCGATCGCCAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCtt                                                                                                                                                  >  2:9234/1‑149 (MQ=255)         aTCGCCAGGTGCGCTATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTcgc                                                                                                                                               <  1:105235/149‑1 (MQ=255)         aTCGCCAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTcgc                                                                                                                                               <  4:59887/149‑1 (MQ=255)            gCCAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGt                                                                                                                                             <  1:9234/148‑1 (MQ=255)             ccAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGc                                                                                                                                           <  7:2594/149‑1 (MQ=255)                ggTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGcttg                                                                                                                                        <  7:36121/149‑4 (MQ=255)                 gTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTTCAGCCCATCGCGTGCaacc                                                                                                                                       >  8:66333/1‑147 (MQ=255)                     gcgATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTcgc                                                                                                                                               <  3:174930/137‑1 (MQ=255)                     gcgATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTcgc                                                                                                                                               >  4:174930/1‑137 (MQ=255)                       gATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCG‑a                                                                                                                                <  6:263849/149‑1 (MQ=39)                        aTGTGTTGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAATGCGGGCGGGGGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTGGCGTTGCAGCGCTTCGCGGGCTTGCACGCG‑aa                                                                                                                               <  8:184392/149‑1 (MQ=255)                        aTGTGATGGTGTATGGGCAATGGGTGGTGCGCCAGGGGCGGCAGGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCG‑aa                                                                                                                               <  2:112370/149‑1 (MQ=38)                              tGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCG‑AAGGGCTg                                                                                                                         <  7:137894/149‑1 (MQ=38)                                gTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGGGGTTGGGCCCATGGCGTGGGCGGCCGCG‑CAGGGTTgg                                                                                                                        >  3:23226/1‑147 (MQ=255)                                gTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTCCGTTGCAGCCCTTCGCGTGCAAGCCCGCG‑AAGGGTTgca                                                                                                                       <  1:203351/149‑1 (MQ=255)                                     tGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTTCAGCCCTTAGCGTGGAAGCCCGCG‑AAGGGTTgccacgc                                                                                                                   >  3:210640/1‑148 (MQ=255)                                                   tgCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCGTGGCGTGCAAGCCCGCG‑GAGGGTTGGACCGCCACCCCAAAAAACTc                                                                                                    >  1:9849/1‑149 (MQ=255)                                                       ccAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCG‑AAGGGCTgcaa                                                                                                                      <  6:304099/127‑1 (MQ=37)                                                       ccAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCG‑AAGGGCTgcaa                                                                                                                      >  5:304099/1‑127 (MQ=37)                                                            gcgggcATGCCCGGTATGTTGAACGCGGGGTGGCGTTTGAAATGGGTTTGCGCATGGTTTTTGGCTGTGTGTACCTGGGTGGTGTTGCTGCCTTTGGTGGCCAACCCCGTG‑ATGGGTTTCAACGCAACCCCAAAAAACTCAGTGCAttt                                                                                           <  4:258301/146‑1 (MQ=255)                                                             gcggcATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCG‑AAggg                                                                                                                            <  7:258896/115‑1 (MQ=37)                                                             gcggcATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCG‑AAggg                                                                                                                            >  8:258896/1‑115 (MQ=37)                                                                gcATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTTGAGCCCGTAGCGGGGAAGCCCGCG‑CAGGGTTGCAACGCCACCCCAAAAAACGCAGTGCCATTTTAc                                                                                       >  6:151195/1‑149 (MQ=255)                                                                gcATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCG‑TAGGGTTGCAAAGCAACCCCAAAAAACTCAGTGCATTTTTa                                                                                        >  6:295247/1‑148 (MQ=255)                                                                      ccAGGAAGATGAACGCGCGCGGGCTGTTGCCAAGGGGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCG‑AAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTg                                                                                 <  4:210640/149‑1 (MQ=255)                                                                            aGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCG‑AAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATc                                                                           <  5:151195/149‑1 (MQ=255)                                                                                            gcgTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCG‑AAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATc                                                           <  7:66333/149‑1 (MQ=255)                                                                                                    cAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCCTTTCAGCCCATAGCGGGCAATCCCCCG‑AACGGATGCAACGCAACCCCAAAAAACACACGTCAATTTTAATATTTTCTGATCGGCCCCCCCCCCCATCAAACGGcg                                                   >  8:173540/1‑148 (MQ=255)                                                                                                               gcCAGGTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCGGCC‑TAGGTTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATAcc                                        <  1:12012/149‑1 (MQ=255)                                                                                                               gcCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGGAGCCCATAGCGGGGCAGCCCGCGCACGGTTTGCAA‑GCAACCCCCAAAAAAAAAGTGCATTTTTTATATTTTGTGATAGGCCCCCCCCCACATAAAGAGGGGGGGGTGATAcc                                        >  2:78256/1‑149 (MQ=255)                                                                                                                               aGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCT‑AAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCt                        <  4:23226/149‑1 (MQ=255)                                                                                                                                 cttACGGGGTTTGCGTTGAATCCCTTCGTGTTAAAGCCCTCT‑ATTGTTTTAAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTg                       <  4:24987/146‑1 (MQ=255)                                                                                                                                   tACCGGGGTTGCGTGGCAGCCCTTCGCGTGCAAGCCCGCG‑AAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGAc                    <  2:9849/149‑1 (MQ=255)                                                                                                                                    aCCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCG‑AAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGAcgcccgc                                                                 <  2:43820/103‑1 (MQ=255)                                                                                                                                    aCCGGGGTTGCGTTGCAGCCCATCGCGTGCCAGCCCGCG‑AAGGGTTGGAACGCAACCCCAAAAAACTCAGTGGATTTTTTATATTTGCTGATCGGccgcccgc                                                                 >  1:43820/1‑103 (MQ=255)                                                                                                                                       ggggTTGCGTTGCAGCCCTTCGCGTGCATGCCCGCG‑AAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGAcagca                <  7:173540/149‑1 (MQ=255)                                                                                                                                       ggggTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCG‑AAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGAcagca                <  4:217484/149‑1 (MQ=255)                                                                                                                                                   cAGCCCTTCGCGTGCAAGCCCGCG‑AAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCt    >  6:157745/1‑149 (MQ=255)                                                                                                                                                   cAGCCCTTCGCGTGCAAGCCCGCG‑AAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTAACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCctg                                    <  7:232746/117‑1 (MQ=255)                                                                                                                                                   cAGCCCTTCGCGTGCAAGCCCGCG‑AAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTAACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCctg                                    >  8:232746/1‑117 (MQ=255)                                                                                                                                                    aGCCCTTCGCGTGCAAGCCCGCG‑AAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTc   >  1:244675/1‑149 (MQ=255)                                                                                                                                                     gCCCTTCGCGTGCAAGCCCGCG‑AAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGGGTCATACCCCTGCTGCCGCGCCCGGGCCAGCAGGGTCTCACGCAcc  >  8:60650/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                      CGGTGGCGATCGCCAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCG‑AAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCC  >  NC_002947/5740309‑5740605

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