Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,740,457 G→A 20.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_0029475,740,4570GA20.5% 69.2 / 14.1 39intergenic (+22/+49)hutF/PP_5037formiminoglutamate deiminase/lipocalin family lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (11/20);  new base A (3/5);  total (14/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.07e-01

CGATCGCCAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGG‑GTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCA  >  NC_002947/5740315‑5740599
                                                                                                                                              |                                                                                                                                               
cGATCGCCAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCTGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTc                                                                                                                                           <  8:83518/149‑1 (MQ=255)
cGATCGCCAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGCGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTc                                                                                                                                           <  4:194637/149‑1 (MQ=255)
 gATCGCCAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGATAGCGGGCGGCCGTTTGCCAAGGTGTGGCGCCCGCGTTTGGGCAGAGCGTACCGGGGTTGCGTTGAAGCCCATcg                                                                                                                                          >  6:371743/1‑149 (MQ=255)
      ccAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTg                                                                                                                                      >  6:154796/1‑148 (MQ=255)
             cgcgATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGccc                                                                                                                               <  2:311631/148‑1 (MQ=255)
                 aTGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGaa                                                                                                                          <  6:68340/149‑1 (MQ=255)
                 aTGTGATGGAGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGaa                                                                                                                          <  7:333261/149‑1 (MQ=255)
                      aTGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGG‑Gtt                                                                                                                    <  2:142292/149‑1 (MQ=255)
                           gAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGGGcttgc                                                                                                                                 >  3:119387/1‑128 (MQ=255)
                           gAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGGGcttgc                                                                                                                                 <  4:119387/132‑5 (MQ=255)
                                 gCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAAccc                                                                                                                                              >  1:361991/1‑113 (MQ=255)
                                 gCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAAccc                                                                                                                                              <  2:361991/113‑1 (MQ=255)
                                          ggtgCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCGAAgg‑gct                                                                                                                    <  4:339649/129‑3 (MQ=38)
                                          ggtgCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCGAAgg‑gct                                                                                                                    >  3:339649/1‑127 (MQ=38)
                                           gtgCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCaa                                                                                                                                   <  2:22790/114‑1 (MQ=255)
                                           gtgCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCaa                                                                                                                                   >  1:22790/1‑114 (MQ=255)
                                             gcgcCCGGGGCCGCATGCCCAGGAAGATGACCCCGGGCGGGCGGTTGCCAAGGGGTGCCGCCAGCGTTTGGGCGGAGCGTACCGGGGTTGCGTTGCAGCCCTGCGCGTGCAAGCCCGCGAAGG‑GTTGCAACGCAACCCCAAAAAACTCa                                                                                             <  4:231708/149‑1 (MQ=255)
                                             gcgcCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCCTTGCAGCCCTTCGCGGGGAAGCCCGCGAACG‑GGTGGAACGCCACCCCCAAAAAca                                                                                               >  5:379090/1‑146 (MQ=255)
                                                   ggggCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGG‑GTTGCAACGCAACCCCAAAAAACTCAGTGCAt                                                                                       <  1:266875/149‑1 (MQ=255)
                                                               ccAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTTCAGCCCTTAGCGTGCAAGCCCGCGAAGG‑GTTGGCACGCAACCCCAAAAAAATCAGTGCATTTTTTCTAtttt                                                                           >  2:24704/1‑148 (MQ=255)
                                                               ccAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGG‑GTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTg                                                                           <  8:291953/149‑1 (MQ=255)
                                                               ccAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCGAAg                                                                                                                         <  6:177473/104‑1 (MQ=37)
                                                                cAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTTGAGCCCATAGCGGGGAAGCCCGCGAAAG‑GTTTGAACGCCAACCCAAAAAAAAAAGTGGAATTTTTATAttttg                                                                          >  6:142546/1‑147 (MQ=255)
                                                                 aGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTAGCGTGGCAGCCCGGGCAGG‑GTTGGCACGCCACCCCAAAAAACGCAGGGGCTTTTTACTATTTTCt                                                                         >  1:239005/1‑149 (MQ=255)
                                                                     aGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGG‑GTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATc                                                                     <  1:24704/149‑1 (MQ=255)
                                                                         gAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCGAAg                                                                                                                         >  3:5015/1‑94 (MQ=25)
                                                                         gAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCGAAg                                                                                                                         <  4:5015/94‑1 (MQ=25)
                                                                                               tggTGTGGCGCCAGCTTTTGGGCGTAGCTTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCTAAGG‑GTTTCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTgg                                           <  4:11881/148‑1 (MQ=255)
                                                                                                 gtgtTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGGGTGCAAGCCCGCGAAGG‑GTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGgtg                                         <  7:278520/149‑1 (MQ=255)
                                                                                                 gtgtTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGCGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTAttt                                                                            <  1:71643/115‑1 (MQ=255)
                                                                                                   gtTCCGCCATCTTTTGGGCTGTGCGTACCGGGGTTGCGTTGCAGCCGTTCGCGGCCAAGCCCGCGATGG‑GTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTc                                       <  8:125227/149‑1 (MQ=255)
                                                                                                       cgcCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGG‑GTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATa                                    <  5:154796/148‑1 (MQ=255)
                                                                                                              ttttGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGG‑GTTGCAACGCAACCCCAAAAAACACAGTGCATTTTTAGTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTGATACCCctgct                            >  7:265776/1‑149 (MQ=255)
                                                                                                                   ggCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGG‑GTTGAAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCcgc                       <  2:239005/149‑1 (MQ=255)
                                                                                                                    gCTGTGCGTACGGGGGTTGCGTTCCTGCCTGTCGGGTGAATGCCGGCTATGG‑TTTCAAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGcc                         <  5:25644/146‑1 (MQ=255)
                                                                                                                                   gttgcgttgcAGCCCTTCGCGTGCAAGCCCGCGAAGG‑GTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGt       >  6:231423/1‑149 (MQ=255)
                                                                                                                                    ttgcgttgcAGCCCTTCGCGTGCAAGCCCGCGAAGG‑GTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGAcagccccc       >  8:277990/1‑144 (MQ=255)
                                                                                                                                       cgttgcAGCCCTTCGCGTGCAAGCCCGCGAAGG‑GTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTtctc   <  5:231423/149‑1 (MQ=255)
                                                                                                                                        gttgcAGCCCTTCGCGTGCAAGCCCGCGAAGG‑GTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCa  <  3:217784/149‑1 (MQ=255)
                                                                                                                                        gttcccgCCCTTCGCGTGCACGCCCGCTACGG‑GTTGAAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCa  <  3:329459/143‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                               
CGATCGCCAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGG‑GTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCA  >  NC_002947/5740315‑5740599

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