Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,740,467	T→G	36.7%	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	36.7%	31.4 / 56.7	55	intergenic (+32/+39)	hutF/PP_5037	formiminoglutamate deiminase/lipocalin family lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (12/23);  new base G (9/11);  total (21/34)
Fisher's exact test for biased strand distribution p-value = 5.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.43e-01

GCCAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGG  >  NC_002947/5740320‑5740613                                                                                                                                                     |                                                                                                                                                   gCCAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGTCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGTTGCAGCCCTTCGCGt                                                                                                                                                    <  1:191336/148‑1 (MQ=255)  ccAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGTTGCAGCCCTTCGCGTGc                                                                                                                                                  <  5:251557/149‑1 (MQ=255)  ccAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGTTGCAGCCCTTCGCGTGc                                                                                                                                                  <  2:298678/149‑1 (MQ=255)      gTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGTTGCAGCCCTTCGCGGGcttgc                                                                                                                                              <  3:252943/149‑5 (MQ=255)             aTGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGTTGCAGCCCTTCGCGTGCAAGCCCGCGaa                                                                                                                                       <  2:318465/149‑1 (MQ=255)              tgtgATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAg                                                                                                                                      <  3:134779/149‑1 (MQ=255)              tgtgATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGTTGCAACCCTTCGCGGGCTTGCATGCGAAg                                                                                                                                      <  2:276062/149‑1 (MQ=39)               gtgATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGTGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGTTGCAACCCTTCGCGGGCTTGCACGCGAAg                                                                                                                                      <  2:6388/148‑1 (MQ=39)                  aTGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGtt                                                                                                                                  <  3:254551/149‑1 (MQ=255)                  aTGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGt                                                                                                                                   <  6:123274/148‑1 (MQ=255)                  aTGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGTTGCAACCCTTCGCGGGCTTGCACGCGAAgggc                                                                                                                                   <  2:98698/148‑2 (MQ=38)                   tGGTGAATGGGCAATGGGGGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGTTGCAACCCTTCGCGGGCTTGCACGCGAAGGGCTg                                                                                                                                 <  6:290079/149‑1 (MQ=38)                           gggCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGTTGCAGCCCTTCGCGGGCTTGCACGCGAAg                                                                                                                                      <  6:66378/136‑1 (MQ=255)                           gggCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGTTGCAGCCCTTCGCGGGCTTGCACGCGAAg                                                                                                                                      >  5:66378/1‑136 (MQ=255)                             gCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGTTGCAACCCTTCGCGGGcttgc                                                                                                                                              <  1:337051/126‑5 (MQ=255)                             gCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGTTGCAACCCTTCGCGGGcttgc                                                                                                                                              >  2:337051/1‑122 (MQ=255)                                      ggtgCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGGCGTTGGGGCCCTAGGGGGGGACGCCCCGCACGGGTTGGACCGCCACCCCcaaaa                                                                                                                >  1:125982/1‑148 (MQ=255)                                         gcgcCTGGGGGGGCTTGCCAAGTATTATTATCGGGGGCGGGCTTTTCCAATGGTGTTGGGCAAGCTTTTGGTCTGTGTTTAGGGGGTTT‑GGGTTCCAGCCCTTGGCGTGCAAGCCGGCTATGGGTTCAAACGCAACCCCAAAAAACTCa                                                                                                           <  4:139346/149‑1 (MQ=255)                                             cAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTg                                                                                                                                 >  7:12027/1‑123 (MQ=255)                                             cAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTg                                                                                                                                 <  8:12027/123‑1 (MQ=255)                                                        tGCCCAGGAATATGACAGCGGGCGGGCGTTTGCAAAGGTGTTGCGCCAGCTTTTGGGCTGCGCGTACCGGGGTT‑GCGTTGCAGCCCTTCGCGTGCAAGCCCGCTAAGGGTTGAAACGCAACCCCAAAAAACTCAGTGCATTTTTACTAt                                                                                            <  5:294683/149‑1 (MQ=255)                                                             aGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGTTGCAGCCCTTAGCGTGGAAGCCCGCGAAGGGTTTGAACGCAACCCCAAAAAAATCAGTGCATTTTTACTAttttgt                                                                                       >  5:61143/1‑146 (MQ=255)                                                                    tGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGATTTTGGGCTGAGCGTACCGGGGTT‑GCGTTGCATCCCTTCGCTTGAAAGCCCGCTAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGa                                                                                <  5:181877/149‑1 (MQ=255)                                                                     gAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAg                                                                                                                                      >  6:210747/1‑94 (MQ=255)                                                                     gAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAg                                                                                                                                      <  5:210747/94‑1 (MQ=255)                                                                           gggcgggcGTTTGCCTAGGGGTTGCTCCTGCTTTTGGGCTGAGCTTACCGGGGTT‑GCGTTGCAGCCCTTCGCGTGCACGCCGGCGATGGGTTGAAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGAcgcccgc                                                                         <  6:126466/149‑1 (MQ=255)                                                                                 gcgTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGTTGCAACCCTTCGCGGGCTTGCACGCGAAggg                                                                                                                                    >  2:127531/1‑84 (MQ=21)                                                                                 gcgTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGTTGCAACCCTTCGCGGGCTTGCACGCGAAggg                                                                                                                                    <  1:127531/84‑1 (MQ=21)                                                                                   gTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGTTTCAGCCCTTAGCGTGCAAGCCCCCGAAGGGTTGGAACGCAACCCCAAAAAACAGAGGGGATTTTTAATATTTGCTGAAAGGACGCCCCCCCGATCaa                                                                 >  6:158771/1‑149 (MQ=255)                                                                                         cAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGGTTCAGCCCATAGCGGGGAAGCCCGCGAACGGTTTGAACCCAAACCCCAAAAAAAAACGGCATTTTTTATATTTTGTGAAAGGCCCCCCGCCACAGAAAGCGGcg                                                           >  2:208472/1‑148 (MQ=255)                                                                                          aaGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTg                                                          <  1:208472/149‑1 (MQ=255)                                                                                               gtTGCGCCAGCTTTTGGGCTGAGCGTACCGNGGTT‑GCGTTGCAGCCCTTCGCGTGGAAGCCCGCGAAGGGTTGGGACGCAACCCCAAAAAACTCAGTGCATTTTTAGTATTTGGTGATGGGGCGCCCGCCCGATGAAGCGGgtggtgg                                                      >  1:27501/1‑147 (MQ=255)                                                                                                   cgcCTGCTTTTGGGGTCCGCGACCCGGGGTG‑GCGTGGCAGCCTTCCGCGTCCACGCCCGATAAGGGTTCAAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATAc                                                 <  2:325968/149‑1 (MQ=255)                                                                                                   cgcCAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATAc                                                 <  2:100901/149‑1 (MQ=255)                                                                                                     ccAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATAccc                                               <  5:158771/149‑1 (MQ=255)                                                                                                                       gTACCGGGGTT‑GCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGa                             <  8:208215/149‑1 (MQ=255)                                                                                                                                tt‑gcgttgcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGtt                    >  5:78056/1‑149 (MQ=255)                                                                                                                                     gttgcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCa                <  6:78056/149‑1 (MQ=255)                                                                                                                                     gttgcAACCCTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCa                >  7:137137/1‑149 (MQ=38)                                                                                                                                     gttgcAACCCTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCa                >  6:187726/1‑149 (MQ=38)                                                                                                                                      ttgcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCAc               >  4:174085/1‑149 (MQ=255)                                                                                                                                        gcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGAcag                          <  4:312109/136‑1 (MQ=255)                                                                                                                                        gcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGAcag                          >  3:312109/1‑136 (MQ=255)                                                                                                                                         caacccTTCGCGGGCTTGCACGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCAc               >  1:214526/4‑146 (MQ=39)                                                                                                                                         caacccTTCGCGGGCTTGCACGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCAc               <  2:214526/143‑1 (MQ=39)                                                                                                                                          aacccTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGa                                                                     <  6:165934/89‑1 (MQ=21)                                                                                                                                          aacccTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGa                                                                     >  5:165934/3‑91 (MQ=21)                                                                                                                                          aGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAgg                      >  6:206733/1‑138 (MQ=255)                                                                                                                                          aGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAgg                      <  5:206733/138‑1 (MQ=255)                                                                                                                                          aGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTc           >  6:13583/1‑149 (MQ=255)                                                                                                                                           gCCCTTCGCGTGCAAGCCCGCTATGGGTTGAAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTcc          <  1:324200/149‑1 (MQ=255)                                                                                                                                            cccTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCCACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCTGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCGcc          >  6:39569/1‑148 (MQ=255)                                                                                                                                            cccTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTg                               <  3:296459/127‑1 (MQ=255)                                                                                                                                            cccTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTg                               >  4:296459/1‑127 (MQ=255)                                                                                                                                              cTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAg       <  8:137137/149‑1 (MQ=39)                                                                                                                                                   cgTGCATCCCCGCTATGGGTTCAAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCgg  <  4:175931/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                   GCCAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGG  >  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: ‑