Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,740,474 C→A 26.8% intergenic (+39/+32) 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,4740CA26.8% 56.1 / 22.6 37intergenic (+39/+32)hutF/PP_5037formiminoglutamate deiminase/lipocalin family lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (13/14);  new base A (4/6);  total (17/20)
Fisher's exact test for biased strand distribution p-value = 7.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.76e-02

GCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGAT‑GAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGTTGCAGCCCTTCGC‑GTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCG  >  NC_002947/5740329‑5740621
                                                                                                                                                    |                                                                                                                                                   
gcgATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATAGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGTTGCAACCCTTCCC‑GTGCAAGcccc                                                                                                                                                   >  1:259385/1‑148 (MQ=255)
        aTGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGAT‑GAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGTTGCAGCCCTTCGC‑GTGCAAGCCCGCGAAGGGtt                                                                                                                                          <  3:157541/149‑1 (MQ=255)
        aTGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGAT‑GAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGTTGCAGCCCTTCGC‑GGGCTTGCACGCGAAgggc                                                                                                                                           <  2:240189/148‑2 (MQ=255)
        aTGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGAT‑GAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGTTGCAACCCTTCGC‑GGGCTTGCACGCGAAgggc                                                                                                                                           <  2:16289/148‑2 (MQ=38)
         tGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGAT‑GAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGTTGCAGCCCTTCGC‑GTGCAAGCCCGCGAAGGGTTg                                                                                                                                         >  8:150704/1‑149 (MQ=255)
           gTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGAT‑GAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGTTGCAACCCTTCGC‑GGGCTTGCACGCGAAGGGCTgc                                                                                                                                        <  8:173354/148‑1 (MQ=38)
              aaTGGGCAATGGGGGGTGCGCCAGGGGCGGCATGCCCAGGAAGAT‑GAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGTTGCAGCCCTTCGC‑GTGCAAGCCCGCGAAGGGTTgcaac                                                                                                                                     <  7:150704/148‑1 (MQ=255)
                   gCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGAT‑GAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGTTGCAGCCCTTCGC‑GTGCAAGCCcgc                                                                                                                                                  <  3:309237/130‑1 (MQ=255)
                   gCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGAT‑GAAAGCGGGCGGGCGTTTGCAAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGTTGCAGCCCTTCGC‑GTGCAAGCCcgc                                                                                                                                                  >  4:309237/1‑130 (MQ=255)
                              tgCGCCAGGGGCGGCATGCCCAGGAAGAT‑GAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGTTTGAGGCCCTAGC‑CGGGGAGCCCGCGAAGGGTTTGAACCGCACCCCaaaaaaaa                                                                                                                     >  2:276589/1‑146 (MQ=255)
                                  ccAGGGGCGGCATGCCCAGGAAGAT‑GAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGGGGTTGC‑GCCCATAGCGGGGCAAGCCCCGCACGGGTTGGAACCCCACCCCaaaaaaaacctg                                                                                                                 >  3:315022/1‑143 (MQ=255)
                                                  cAGGAAGAT‑GAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGTTGCAACCCTTCGC‑GGGCTTGCACGCGAAg                                                                                                                                              >  8:36291/1‑103 (MQ=37)
                                                  cAGGAAGAT‑GAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGTTGCAACCCTTCGC‑GGGCTTGCACGCGAAg                                                                                                                                              <  7:36291/103‑1 (MQ=37)
                                                              aaGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTG‑GCGTTGCAGCCCTTCGC‑GTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGAcgc                                                                                     <  3:81467/149‑1 (MQ=255)
                                                                    gcgggcgTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGTTGCAGCCCTTCGC‑GTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCa                                                                               <  8:238596/149‑1 (MQ=255)
                                                                        gcgTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGTTGCAGCCCTTCGC‑GTGCAAGCCCGCGAAg                                                                                                                                              <  5:280405/82‑1 (MQ=255)
                                                                        gcgTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGTTGCAGCCCTTCGC‑GTGCAAGCCCGCGAAg                                                                                                                                              >  6:280405/1‑82 (MQ=255)
                                                                                       ttGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGTTGCAGCCCTTCGC‑GTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTc                                                             <  1:16493/148‑1 (MQ=255)
                                                                                          cgcCAGGTTTTGGGCGGAGGGTACCGGGGTG‑GCGTTGCAGCCGTGGGC‑GTGCAAGCCGGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGTCGCCCGCCAGATCAAGCGGCTGGTGTCATa                                                          <  1:196623/148‑1 (MQ=255)
                                                                                            ccAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGTTGCAGCCCTTAGC‑GGGCAAGCCCGCGAAGGGTTGGCACGCCACCCCCAAAAACGCAGTGGATTTTTTCGATTTGCGGATCGGACGCCCGCCAGATCAAGCGGCGGGGGTCATAacc                                                       >  2:45521/1‑149 (MQ=255)
                                                                                             cAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGTTGCAGCCCTTCGC‑GTGCAAGCCCGCGAAGGGTTgc                                                                                                                                        >  2:316016/1‑67 (MQ=255)
                                                                                              aGCTTTTGGGCTGAGCGTACCGGGGTG‑GCGTTCCAGCCCGTCGC‑GGCCAAGCCCGCGAAGGGTGGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCct                                                     <  4:304944/149‑1 (MQ=255)
                                                                                                 ttttGGGCTGAGCGTACCGGGGTT‑GCGTTGCAGCCCTTCGC‑GTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCACATCAAGCGGCTGGTGTCATACCCctgct                                                  <  1:187067/149‑1 (MQ=255)
                                                                                                   ttGGGCTGAGCGTACCGGGGTT‑GCGTTGCAGCCCTTCGC‑GTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCctgct                                                  <  2:259385/147‑1 (MQ=255)
                                                                                                           aGCGTACCGGGGTG‑GCGTTGCAGCCCTTCGC‑GTGCAAGCCCGCGATGGGTTGAAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCt                                        <  4:315022/149‑1 (MQ=255)
                                                                                                            gCGTACCGGGGTT‑GCGTTGCAGCCCTTCGC‑GTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTg                                       <  1:45521/149‑1 (MQ=255)
                                                                                                                      gtt‑gcgttgcAGCCCTTCGC‑GTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGt                             >  2:127446/1‑149 (MQ=255)
                                                                                                                             ttgcAACCCTTCGC‑GGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCAc                       >  5:120314/1‑149 (MQ=38)
                                                                                                                                cAGCCCTTCGC‑GTGCAAGCCCGCGAAGGGTTGCAACGAAACCCCAAAAAACTCAGTGCAATTTTACTATTTGAAGAACGTACTCACGCAAGATAAAGCGGCTGGTGTCATACCCCTTCTGCCGCGCACTCGACAGCAGGTTCTCACGCt                    >  6:125138/1‑149 (MQ=255)
                                                                                                                                  gCCCTTCGC‑GTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGCTCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTcc                  >  8:188150/1‑149 (MQ=255)
                                                                                                                                     cTTCGC‑GTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAg               >  4:231948/1‑149 (MQ=255)
                                                                                                                                     cTTCGC‑GGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGcc                                               <  6:128440/117‑1 (MQ=38)
                                                                                                                                     cTTCGC‑GGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGcc                                               >  5:128440/1‑117 (MQ=38)
                                                                                                                                      ttCGC‑GTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAgg              >  3:6168/1‑149 (MQ=255)
                                                                                                                                       tCGC‑GGGCTTGCACGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCt                                        <  5:104851/122‑1 (MQ=255)
                                                                                                                                       tCGC‑GGGCTTGCACGCGAAGGGTTGCAAAGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCt                                        >  6:104851/1‑122 (MQ=39)
                                                                                                                                             tGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCa        <  1:127446/149‑1 (MQ=255)
                                                                                                                                                 aGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAgtgt    >  7:165564/1‑149 (MQ=255)
                                                                                                                                                   cccGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCg  <  3:93705/149‑1 (MQ=255)
                                                                                                                                                   cccGCCAAGGGTTCCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCg  <  5:283882/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGAT‑GAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTT‑GCGTTGCAGCCCTTCGC‑GTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCG  >  NC_002947/5740329‑5740621

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