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

CGGTGGCGATCGCCAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCC  >  NC_002947/5740309‑5740605                                                                                                                                                      |                                                                                                                                                     cGGTGGCGATCGCCAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGTTGCAg                                                                                                                                                      <  8:228598/149‑1 (MQ=255)        gATCGCCAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGTTGCAGCCCTTcg                                                                                                                                               >  6:251557/1‑149 (MQ=255)            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                                                                                                                                      >  2:337051/1‑122 (MQ=255)                                        gCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGTTGCAACCCTTCGCGGGcttgc                                                                                                                                      <  1:337051/126‑5 (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                                                                                                                            <  1:127531/84‑1 (MQ=21)                                                                                            gcgTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGTTGCAACCCTTCGCGGGCTTGCACGCGAAggg                                                                                                                            >  2:127531/1‑84 (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        >  6:187726/1‑149 (MQ=38)                                                                                                                                                gttgcAACCCTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCa        >  7:137137/1‑149 (MQ=38)                                                                                                                                                 ttgcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCAc       >  4:174085/1‑149 (MQ=255)                                                                                                                                                   gcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGAcag                  >  3:312109/1‑136 (MQ=255)                                                                                                                                                   gcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGAcag                  <  4:312109/136‑1 (MQ=255)                                                                                                                                                     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)                                                                                                                                                      |                                                                                                                                                     CGGTGGCGATCGCCAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCC  >  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: ‑