Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,740,467 T→G 14.6% 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_0029475,740,4670TG14.6% 54.8 / 6.4 34intergenic (+32/+39)hutF/PP_5037formiminoglutamate deiminase/lipocalin family lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (8/21);  new base G (3/2);  total (11/23)
Fisher's exact test for biased strand distribution p-value = 3.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.68e-01

CCAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGC  >  NC_002947/5740321‑5740610
                                                                                                                                                  |                                                                                                                                               
ccAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGt                                                                                                                                                 >  1:111410/1‑147 (MQ=255)
ccAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGt                                                                                                                                                 <  2:111410/147‑1 (MQ=255)
ccAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGc                                                                                                                                               <  2:213645/149‑1 (MQ=255)
    gTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGc                                                                                                                                           <  6:331849/149‑1 (MQ=255)
                aTGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGTCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGtt                                                                                                                               <  8:197682/149‑1 (MQ=255)
                aTGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGt                                                                                                                                <  4:129968/148‑1 (MQ=255)
                  ggTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCCAAGGGTTg                                                                                                                              >  6:351068/1‑148 (MQ=255)
                    tGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTgcaa                                                                                                                           <  3:55165/149‑1 (MQ=255)
                     gAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGGGCTTGCACGCGaa                                                                                                                                    >  1:280092/1‑139 (MQ=255)
                     gAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGGGCTTGCACGCGaa                                                                                                                                    <  2:280092/139‑1 (MQ=255)
                                       gcgcCTGGGGCGGATTGCCCTGTAAGTTTAACGGGGGGGGGGGTTGGCCATGGTGTGGCGCCTGTTTTGGGGCTGTGCGTACGGGGGTTGCGTTCCAGCCCTTCGCGTGCAAGCCCGCTATGGTTTGAAACGCAACCCCAAAAAACTCa                                                                                                        <  8:223746/149‑1 (MQ=255)
                                        cgcCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAg                                                                                                       <  4:119581/149‑1 (MQ=255)
                                           cAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTTCAGCCCATAGCCTGCAAACCCGCGCAGGGTTGGAACGCCACCCCCAAAAAcacaagg                                                                                                     >  4:309340/1‑145 (MQ=255)
                                           cAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGt                                                                                                                                                 <  2:103664/104‑1 (MQ=255)
                                           cAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGt                                                                                                                                                 >  1:103664/1‑104 (MQ=255)
                                                                   gAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCCCGCCCGCGCAGGGTTGCAACGCAACCCCCAAAAACTCAGGGCATTTTTAAAAGTTGCTGATCGGAc                                                                            >  2:270011/1‑149 (MQ=255)
                                                                     aaGCGGGCGGGCGTTTGCCAAGGTGTGGCGGCAGGTTTTGGGCTGAGCGTACCGGGGTTGCGGTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTGGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGAcgc                                                                          <  4:233794/149‑1 (MQ=255)
                                                                                        aaGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTg                                                       <  7:65907/149‑1 (MQ=255)
                                                                                        aaGGGGTTGCGCCAGCTTTTGGGCTGCGCGTACGGGGGTGGCGTGGCAGCCCTTCGCGTTCAAGCCCGCGAAGGGTTGAAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCTAGCGGCTg                                                       <  5:105969/149‑1 (MQ=255)
                                                                                                    cAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCACGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCt                                                        <  5:327891/136‑1 (MQ=255)
                                                                                                    cAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGc                                                         <  1:276848/135‑1 (MQ=255)
                                                                                                          ttGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCctgct                                       <  1:270011/147‑1 (MQ=255)
                                                                                                                  aGCGTACCGGGGTTGCGTTTCAGCCCATAGCGTGGAAACCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTAATATTTTGTGATAGGACGCCCGCCAGATAAAGCCGCTGGTGTCATACCCCTGCTGCCGCGcccg                             >  6:377132/1‑148 (MQ=255)
                                                                                                                   gggTACCGGGGTTGCGTTGCAGCCCTTGGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTg                            <  3:105197/147‑1 (MQ=255)
                                                                                                                   gCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTg                            <  5:377132/149‑1 (MQ=255)
                                                                                                                     gTACCGGGGTTGAGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGa                          <  3:50732/149‑1 (MQ=255)
                                                                                                                         cGGGGTTGCGTTTCAGCCCATAGCGTGCAAGCCCGCGAAGGGTTTGAACGCAACCCCAAAAAAATAAGTGCATTTTTTATATTTTGTGGTAGGACGCCCCCCAGATCAAGCGGcgg                                                       >  4:314572/1‑114 (MQ=255)
                                                                                                                         cGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTg                                                       <  3:314572/116‑1 (MQ=255)
                                                                                                                                gcgttgcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTg                            <  8:76456/136‑1 (MQ=255)
                                                                                                                                gcgttgcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTg                            >  7:76456/1‑136 (MQ=255)
                                                                                                                                gcgttgcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTtc                >  6:272962/1‑148 (MQ=255)
                                                                                                                                  gttgcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCa             >  4:329410/1‑149 (MQ=255)
                                                                                                                                  gttgcAACCCTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCAGTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCa             >  6:169781/1‑149 (MQ=38)
                                                                                                                                      cAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGc          <  1:301731/148‑1 (MQ=255)
                                                                                                                                       aacccTTCGCGGGCTTGCACGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTc        >  8:181732/3‑149 (MQ=255)
                                                                                                                                        gCCCTTGGCGGGCAAGCCCGCTAAGGGTGGCAACGCTACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTcc       <  5:351068/149‑1 (MQ=255)
                                                                                                                                          ccTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCa     <  5:272962/149‑1 (MQ=255)
                                                                                                                                              cgcgTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGc  <  3:310613/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
CCAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGC  >  NC_002947/5740321‑5740610

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