Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,740,474 C→A 16.7% 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,4740CA16.7% 82.6 / 10.5 36intergenic (+39/+32)hutF/PP_5037formiminoglutamate deiminase/lipocalin family lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (12/18);  new base A (3/3);  total (15/21)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.80e-01

CGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCG  >  NC_002947/5740328‑5740621
                                                                                                                                                    |                                                                                                                                                    
cgcgATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGTGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGccc                                                                                                                                                     <  3:152467/148‑1 (MQ=255)
    aTGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCcgcg                                                                                                                                                  >  6:162650/1‑147 (MQ=255)
    aTGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGaa                                                                                                                                                <  3:40765/149‑1 (MQ=255)
     tgtgATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCGAAg                                                                                                                                               <  1:63419/149‑1 (MQ=39)
           tgTGAATTGGCTATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTTCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCTAAg                                                                                                                                               <  2:228284/142‑1 (MQ=255)
           ggTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGGTGGAGCCCATAGCGTGGAAGCCCGCGCAGGGttt                                                                                                                                          >  7:193959/1‑147 (MQ=255)
           gctGAATGGGCAATGGGTGGTGCGCCAGGGGCGTCATGCACAGGAAGATGAAAGCGGGCGGGCCTTTGCCAAGCTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGCGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGAGAAg                                                                                                                                               >  1:228284/3‑143 (MQ=255)
              gAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGGGAGCCCGCGCAGGGGTggaac                                                                                                                                     >  7:241067/1‑149 (MQ=255)
                              gtgCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACt                                                                                                                     <  3:188151/149‑1 (MQ=255)
                                gcgcCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGGTTGCGCCCTTAGCGGGGGAGCCCGCGCAGGGGTTGAACGCCACCCCaaaaaaaa                                                                                                                     >  3:281431/1‑145 (MQ=255)
                                gcgcCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTTGAGCCCTTAGCGTGGCAGCCCGCGAAGGGTTGGAACGCAACCCCCAAAAAca                                                                                                                     >  1:91901/1‑146 (MQ=255)
                                    ctggggggGAAGGCCCTGTAAGTTCAAGGGGGGGGGGGGTTTGCCTTGGTGTGGCGCCTGTTTTTGGGGTGTGCGTAGCGGGGTGGGGTTGCAGCCTTTCGGGTCCAAGCCCGCTATGGTTTTAAACGCAACCCCAAAAAACTCAGTGc                                                                                                               <  8:116320/147‑1 (MQ=255)
                                      ggggCGGCATGCCCAGTATGTTTAACGCGGGGGGGGTTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGCGGTTACGGGGGTTGGGTTGCAGCCCTTCGCGTGCAAGCCCGCTATGGGTTGAAACGCAACCCCAAAAAACTCAGTGCAt                                                                                                             <  4:281431/149‑1 (MQ=255)
                                            gcATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTgcaac                                                                                                                                     >  8:72293/1‑119 (MQ=255)
                                            gcATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTgcaac                                                                                                                                     <  7:72293/119‑1 (MQ=255)
                                                   cAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCTAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGc                                                                                                <  3:70629/149‑1 (MQ=255)
                                                        aGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTCCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAAGCGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGAt                                                                                            <  1:220370/149‑1 (MQ=255)
                                                            gAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGGGGTGGAGCCCATACGCGGGAAGGCCGCGCAGGGTTTGAACGCCAACCCAAAAAAAGCAGTGGAGTTTTAAATTTTGATGATggggc                                                                                       >  3:26918/1‑147 (MQ=255)
                                                                    gcgggcgTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCa                                                                               <  5:162650/149‑1 (MQ=255)
                                                                    gcgggcgTTGGCCAAGGTGTTCCGCCAGCTTTTGGGCTGCGCGTACCGGGGTTGCGTTGCAGCCCTGCGCGTGCAAGCCCGCTAAGGTTTGAAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCa                                                                               <  4:113214/149‑1 (MQ=255)
                                                                        gcgTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGGGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACt                                                                                                                     <  5:150461/107‑1 (MQ=255)
                                                                          gTTTGCCAAGGTGTGGCGCCAGCTTTTGGGCTTAGCGTACCGGGGTGGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCaa                                                                         <  3:281810/149‑1 (MQ=255)
                                                                             tGCCAAGGGTGTGCGCCCGCTTTTGGGCTGAGCGTACCGGGGTTGCTTTGCACCCCTTCGCGTGCAAGCCCGCTATGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGc                                                                       <  5:176933/148‑1 (MQ=255)
                                                                                       ttGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTc                                                             <  1:6861/148‑1 (MQ=255)
                                                                                            ccAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGCGCAACCCCGCGACGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATAccc                                                       <  4:26918/149‑1 (MQ=255)
                                                                                                           aGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGc                                                                       <  3:113007/118‑1 (MQ=255)
                                                                                                                    gggTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGAcagcag                               <  2:91901/149‑1 (MQ=255)
                                                                                                                    gggTTGCGTTCCAGCCCTTCGCGTGCAACCCCGCTATGGGTTGAAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGAcagcag                               <  3:145409/149‑1 (MQ=255)
                                                                                                                            gttgcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCa                        >  5:157943/1‑149 (MQ=255)
                                                                                                                                 aacccTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCa                                                                          >  7:61840/3‑94 (MQ=21)
                                                                                                                                 aacccTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCa                                                                          <  8:61840/92‑1 (MQ=21)
                                                                                                                                       ttCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAgg              >  1:84715/1‑149 (MQ=255)
                                                                                                                                       ttCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAgg                              >  5:186865/1‑133 (MQ=38)
                                                                                                                                       ttCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAgg                              <  6:186865/133‑1 (MQ=38)
                                                                                                                                       ttCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAgg              >  7:136618/1‑149 (MQ=39)
                                                                                                                                                aaGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAgtg     >  2:152333/1‑149 (MQ=255)
                                                                                                                                                  gCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTc   >  4:203417/1‑149 (MQ=255)
                                                                                                                                                   cccGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCg  >  5:156431/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCG  >  NC_002947/5740328‑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: