Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,740,457	G→A	13.8%	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	13.8%	100.4 / 7.9	46	intergenic (+22/+49)	hutF/PP_5037	formiminoglutamate deiminase/lipocalin family lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (11/27);  new base A (2/4);  total (14/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.22e-02

TGGCGATCGCCAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACC‑‑‑GGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCC  >  NC_002947/5740312‑5740605                                                                                                                                                     |                                                                                                                                                     tGGCGATCGCCAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACC‑‑‑GGGGTTGCGTTGCAGccc                                                                                                                                                   <  8:123031/149‑1 (MQ=255)    cGATCGCCAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACC‑‑‑GGGGTTGCGTTGCAGCCCTTc                                                                                                                                                <  5:246201/149‑1 (MQ=255)    cGATCGCCAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACC‑‑‑GGGGTTGCGTTGCAGCCCTTc                                                                                                                                                <  4:163074/149‑1 (MQ=255)       tCGCCAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACC‑‑‑GGGGTTGCGTTGCAGCCCTTcgcg                                                                                                                                             >  6:368397/1‑149 (MQ=255)        cGCCAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACC‑‑‑GGGGTTGCGTTGCAACCCTTcgcgg                                                                                                                                            >  7:335409/1‑148 (MQ=255)         gCCAGGGGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACC‑‑‑GGGGTTGCGTTGCAGCCCTTCGCGt                                                                                                                                            <  2:329741/148‑1 (MQ=255)           cAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACC‑‑‑GGGGTTGCGTTGCAGCCCTTCGCGTGc                                                                                                                                          >  6:276600/1‑148 (MQ=255)                gcgcGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACT‑‑‑GGGGTTGCGTTGCAGCCCTTCGCGTGCAAGccc                                                                                                                                    <  6:333224/149‑1 (MQ=255)                      tgtgATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACC‑‑‑GGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCGAAg                                                                                                                              <  8:335409/149‑1 (MQ=39)                      tgtgATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTGGCGCCAGCTTTTGGGCTGAGCGTACC‑‑‑GGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCGAAg                                                                                                                              <  6:111264/149‑1 (MQ=39)                      tggcaTGGTGAATGGGCAATGTGTGGTGCGCCAGGGGCTGCATGCCCAGGAAGATGAAAGCTGGCGGGCGTTTTCCAAGGTGTTGCGCCAGCTTTTTGGCTGAGCGTACC‑‑‑GGTGTTGCGTTTCAGCCCTTCGCGTGCAAGCCCGCGAAg                                                                                                                              <  5:6033/145‑1 (MQ=255)                       gtgATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACC‑‑‑GGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCGAAg                                                                                                                              <  4:184827/148‑1 (MQ=39)                        tgATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACC‑‑‑GGGGTTGCGTTGCCGCCCTTCGCGTGGAACCCCGCGCAggg                                                                                                                            >  8:202293/1‑149 (MQ=255)                        tgATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACC‑‑‑GGGGTTGCGTTGCAGCCCTTCGCGGGCTTGCACGCGAAgg                                                                                                                             <  4:257382/148‑1 (MQ=255)                         gATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCCTGGGGGTTGCGTTGGACCACATAGCGTGCAACCCCGCCAAg                                                                                                                              >  1:126985/1‑149 (MQ=255)                           tGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACC‑‑‑GGGGTTGCGTTTCAGCCCTTAGCGTGCAAGCCCGCGAAGGGTTg                                                                                                                         >  6:357752/1‑149 (MQ=255)                              tGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACC‑‑‑GGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTgcaa                                                                                                                      <  8:271147/149‑1 (MQ=255)                                           ggtggtGCGCCAGGGGCGGCATGCCCAGGAAGATGAACGCGGGCGGGCGTTGGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACC‑‑‑GGGGTGGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCaaa                                                                                                          <  7:202293/148‑1 (MQ=255)                                                tgCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACC‑‑‑GGGGTTGCGTTGCAGCCCATAGCGTGCAACCCCGCGAAGGGTTGCAACGCAACCCCaaaaaaaa                                                                                                     >  8:191405/1‑146 (MQ=255)                                                 gcgcCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACC‑‑‑GGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCTATGGGTTGCAACGCAACCCCAAAAAACTCa                                                                                                   <  2:251511/149‑1 (MQ=255)                                                     cAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACC‑‑‑GGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTgcaacgcaa                                                                                                                 <  2:332442/131‑1 (MQ=255)                                                        gggCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACC‑‑‑GGGGTTGCGTTGCAGCCCTTc                                                                                                                                                >  5:172692/1‑97 (MQ=255)                                                        gggCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACC‑‑‑GGGGTTGCGTTGCAGCCCTTc                                                                                                                                                <  6:172692/97‑1 (MQ=255)                                                         ggcggcATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACC‑‑‑GGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCAttt                                                                                           <  5:276600/149‑1 (MQ=255)                                                                  cccATGAAGATGAATGGGGGGGGGCGTTTGCCAAGTTGTTGCGCCTTCTTTTGGGCTGAGGTTACC‑‑‑GGGGTTGGGTTGCTCCCGTTCGCTTGCAAGCCGGCGATGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTAtt                                                                                   <  2:126985/148‑1 (MQ=255)                                                                   ccAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACC‑‑‑GGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAg                                                                                                                              >  5:354188/1‑104 (MQ=255)                                                                   ccAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACC‑‑‑GGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAg                                                                                                                              <  6:354188/104‑1 (MQ=255)                                                                         aGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACC‑‑‑GGGGTTGCGTTGCAGCCCTTCGCGTGc                                                                                                                                          <  3:44432/86‑1 (MQ=255)                                                                         aGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACC‑‑‑GGGGTTGCGTTGCAGCCCTTCGCGTGc                                                                                                                                          >  4:44432/1‑86 (MQ=255)                                                                                                     gtgtTGCTCCAGCTTTTGGGGTGAGCGTACG‑‑‑GGGGTTGCGTTGCAGGCCTTCGCGTGCATGGCGGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGgtg                                               <  8:99302/149‑1 (MQ=255)                                                                                                       gtTGCGCCAGCTTTTGGGCTGAGCGTACC‑‑‑GGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGTTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTc                                             <  7:191405/149‑1 (MQ=255)                                                                                                             ccAGCTTTTGGGCTGAGCTTACC‑‑‑GGGGTTGCGTTGCAGCCCTTCGCGTTCAAGCCCGATAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATAccc                                       <  7:136240/149‑1 (MQ=255)                                                                                                                       ggCTGAGCGTACC‑‑‑GTGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGGGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCcgc                             <  2:294726/149‑1 (MQ=255)                                                                                                                       ggCTGAGCGTACC‑‑‑GGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCcgc                             <  1:235186/149‑1 (MQ=255)                                                                                                                        gCTGAGCGTACC‑‑‑GGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTc                                             <  5:110495/132‑1 (MQ=255)                                                                                                                         cTGAGCGTACC‑‑‑GGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCgcgc                           <  5:368397/149‑1 (MQ=255)                                                                                                                               ttACC‑‑‑GGGTTTCCGTTGCAGCCGTTCGCGGGAAAGCCCGCTAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGa                     <  7:320705/148‑1 (MQ=255)                                                                                                                                        gggTTGCTTTGCAGCCCTTGGCGTACAAGCCCGCTATGGTTTGAAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGAcagcag               <  3:269557/149‑1 (MQ=255)                                                                                                                                         ggTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAgg              <  7:64448/149‑1 (MQ=255)                                                                                                                                             gcgttgcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTg                                                  <  4:70231/109‑1 (MQ=255)                                                                                                                                             gcgttgcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTg                                                  >  3:70231/1‑109 (MQ=255)                                                                                                                                               gttgcAACCCTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGAc                    <  7:4122/137‑1 (MQ=37)                                                                                                                                               gttgcAACCCTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGAc                    >  8:4122/1‑137 (MQ=37)                                                                                                                                                 tgcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGAcagcag               <  3:79263/140‑1 (MQ=255)                                                                                                                                                 tgcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGAcagcag               >  4:79263/1‑140 (MQ=255)                                                                                                                                                    aGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTc   >  3:28999/1‑149 (MQ=255)                                                                                                                                                     gCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAAGTCAGTGGATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTcc  >  1:95042/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                     TGGCGATCGCCAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACC‑‑‑GGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCC  >  NC_002947/5740312‑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: ‑