Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,740,484 T→C 15.1% intergenic (+49/+22) 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,4840TC15.1% 70.2 / 7.1 34intergenic (+49/+22)hutF/PP_5037formiminoglutamate deiminase/lipocalin family lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (8/19);  new base C (3/2);  total (13/21)
Fisher's exact test for biased strand distribution p-value = 3.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.31e-01

GGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGGTGTGCGCGA  >  NC_002947/5740339‑5740631
                                                                                                                                                   |                                                                                                                                                    
ggTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGGTGGAGCCCATAGCGTGGAAGCCCGCGCAGGGttt                                                                                                                                                    >  7:193959/1‑147 (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                                                                                    <  8:61840/92‑1 (MQ=21)
                                                                                                                      aacccTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCa                                                                                    >  7:61840/3‑94 (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)
                                                                                                                                            cgAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGgtgt        <  1:152333/149‑1 (MQ=255)
                                                                                                                                                gggTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGGTGTgcgc    <  3:203417/149‑1 (MQ=255)
                                                                                                                                                  gTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGGTGTGCGCGa  <  6:157943/149‑1 (MQ=255)
                                                                                                                                                   ttGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGGTGTGCGCGa  >  7:28938/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGGTGTGCGCGA  >  NC_002947/5740339‑5740631

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: