Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,740,484 T→C 11.3% 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,4840TC11.3% 101.6 / 5.7 46intergenic (+49/+22)hutF/PP_5037formiminoglutamate deiminase/lipocalin family lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (15/22);  new base C (3/2);  total (21/25)
Fisher's exact test for biased strand distribution p-value = 6.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.62e-01

ATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGC‑TTTTGGGCTGAGCGTACCGGGG‑TTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGGT  >  NC_002947/5740337‑5740623
                                                                                                                                                     |                                                                                                                                           
aTGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGC‑TTTTGGGCTGAGCGTACCGGGGTTTGGGTTG‑AGCCCCTAGCCTGGAAGCCCGCCAAAGGtt                                                                                                                                            >  2:356442/1‑149 (MQ=255)
 ttttGAATGGGGAATGGGTGTTGCGTCAGGGGCGGCATGCGCAGGTAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGC‑TTTTGGTCTGAGCGTACCGGGG‑TTGCGTTGCAACGCTTCGCGGGGTTGCACGCGAAGGGGTg                                                                                                                                           <  8:373868/146‑1 (MQ=255)
 tGGTGAATGGGCAATGGGGGGTGCGCCAGGGGCGGCAAGCCCATGAAGTTGTAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTTCTGCTGAGCGTACCGGGG‑TTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGtt                                                                                                                                            <  1:304224/149‑1 (MQ=255)
  ggTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGC‑TTTTGGGCTGAGCGTACCGGGG‑TTTCGTTTGAGCCCATACCCTGGAAGCCCCGGAAAGGTTg                                                                                                                                           >  8:261094/1‑148 (MQ=255)
    tGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGC‑TTTTGGGCTGAGCGTACCGGGG‑TTGGGGTTGAGCCCCTGGCGGGGCAGCCCGCGCAGGGGTggca                                                                                                                                        >  5:345102/1‑146 (MQ=255)
      aaTGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGC‑TTTTGGGCTGAGCGTACCGGGG‑TTGCGTTTGAGCCCATAGCGTGGAAGCCCGCGAAGGGTTtgaacg                                                                                                                                      >  2:21489/1‑149 (MQ=255)
      aaTGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGC‑TTTTGGGCTGAGCGTACCGGGG‑TTGCCTTGCAGCCCATAGCGTGCCACCCCGCGCAGGGTTgccacg                                                                                                                                      >  6:173765/1‑149 (MQ=255)
                  gtggggCGCCAGGGGCGGCATGCCCAGGAAGAGGAACGCGGGCGGGCGTTTGCCAAGGTGGTGCGCCAGC‑TTTTGGGCTGAGCGTACCGGGG‑TTGCGTGGCAGCCCTGCGCGTGCAACCCCGCGAAGGGTTGCAACGCAACCCCaaaaa                                                                                                                          <  7:93242/149‑1 (MQ=255)
                   tcgtccgcCAGGGGCGGCATGCCCAGGAAGATGAACGCGGCCGGGCGTTTGCCAAGGTGTTGCCCCAGC‑TTTTGGGCTGAGCGTACCGGGG‑TTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCaaaaaa                                                                                                                         <  5:198105/144‑1 (MQ=255)
                         gcCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGC‑TTTTGGGCTGAGCGTACCGGGG‑TTGCGTTTCAGCCCATAGCGTGGAAGGCCGCGCAGGGGTGCAACGCCACCCCACAAAACTCAGt                                                                                                                   >  6:360659/1‑149 (MQ=255)
                          ccAGGGGCCGCATGCCCACGAAGATGAAAGCGGGCGGGCGTTTGCCAAGCTGTTGCGCCAGC‑TTTTGGGCTGAGCGTACCGGGG‑TTGCGTTGCAGCCCTAACCGGGCAACGCCCCCCAGGGTTGCAAAGCCACCCCAAAAAACtagctg                                                                                                                  >  8:283920/1‑144 (MQ=255)
                                   gcATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGC‑TTTTGGGCTGAGCGTACCGGGG‑TTGCCATTCAGCCCATACCGTGGCACCCCCCCAAGGGTTGCAACCCAACCCCAAAAAACAGAGAGCAtttttt                                                                                                          >  6:183937/1‑147 (MQ=255)
                                        cccAGGAAGATGAAAGGGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGC‑TTTTGGGCTGAGCGTACCGGGG‑TTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTAtt                                                                                                     <  7:261094/148‑1 (MQ=255)
                                        cccAGGAAGATGAAAGCGGGCGGTCGTTTGCCAAGGTGTTGCGCCAGC‑TTTTGGGCTGAGCGTACCGGGG‑TTGCGTTGCAGCCCTTCTCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTAttt                                                                                                    <  5:183937/149‑1 (MQ=255)
                                        cccAGGAAGATGAAAGCGGGCGGGCTTTTGCCAAGGTGTTGCGCCTGC‑TTTTGGGCTGAGCGTACCGGGG‑TTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTAttt                                                                                                    <  5:173765/149‑1 (MQ=255)
                                        cccAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGC‑TTTTGGGCTGAGCGTACCGGGG‑TTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTAttt                                                                                                    <  1:21489/149‑1 (MQ=255)
                                         ccAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGC‑TTTTGGGCTGAGCGTACCGGGG‑TTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTg                                                                                                   <  6:377057/149‑1 (MQ=255)
                                          cAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGC‑TTTTGGGCTGAGCGTACCGGGG‑TTGCGTTTCCGCCCCTCGCCGGGCAGCCCCCGCCCGGGTGCCCCGCCACCCCAAAAAACACAGTGCCATTTTAATAttttc                                                                                                  >  7:4128/1‑147 (MQ=255)
                                          cAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGC‑TTTTGGGCTGAGCGTACCGGGG‑TTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGc                                                                                                  <  5:360659/149‑1 (MQ=255)
                                            ggAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGC‑TTTTGGGCTGAGCGTACCGGGG‑TTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTg                                                                                                <  6:374808/149‑1 (MQ=255)
                                             taaGATCAAAGCGGGGGGGCGTTTGCCATGGTGTTGCGCCTGG‑TTTTGGGCTGAGCGTACGGGGG‑TTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCTAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGa                                                                                               <  6:131083/148‑1 (MQ=255)
                                                    aaaGCGGGCGGGCGTTTGCAAAGGTGTTGCGCCAGC‑TTTTGGGCTTAGCGTACCGGGG‑TTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCTAAGGGTTGAAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGAc                                                                                         <  1:8980/148‑1 (MQ=255)
                                                          ggcgggcgTTTGCCAAGGTGTTGCGCCAGC‑TTTTGGGCTGAGCGTACCGGGG‑TTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCAAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGAcgcccgcc                                                                                  <  1:356442/149‑1 (MQ=255)
                                                          ggcgggcgTTTGCCAAGGTGTGGCGCCAGC‑TTTTGGGCTGAGCGTACCGGGG‑TGGCGTTGCAGCCCTTCGCGGCCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGAcgcccgcc                                                                                  <  3:360534/149‑1 (MQ=255)
                                                                           gtgtTGCGCCAGC‑TTGTGGGCTGCGCGTACCGGGG‑TTGGGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGgt                                                                  <  8:4128/148‑1 (MQ=255)
                                                                                gcgcCAGC‑TTTTGGGCTGCGCGTACCGGGG‑TGGCGTGGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCAt                                                             <  7:283920/148‑1 (MQ=255)
                                                                                 cgcCAGC‑TTTTGGGCTGAGCGTACCGGGG‑TTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGAAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATAc                                                           <  3:40293/149‑1 (MQ=255)
                                                                                                              g‑ttgcgttgcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCACTGGACAGCAGGt                               >  2:34171/1‑149 (MQ=255)
                                                                                                                  gcgttgcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTtct                            >  3:362277/1‑149 (MQ=255)
                                                                                                                    gttgcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCcg                                                <  8:360020/127‑1 (MQ=255)
                                                                                                                    gttgcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCcg                                                >  7:360020/1‑127 (MQ=255)
                                                                                                                        cAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCt                      >  7:283667/1‑149 (MQ=255)
                                                                                                                        cAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCt                      >  1:224756/1‑149 (MQ=255)
                                                                                                                         aacccTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTc                     >  4:211815/3‑149 (MQ=38)
                                                                                                                         aGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTc                     >  3:63730/1‑149 (MQ=255)
                                                                                                                          acccTTCGCGGGCTTGCACGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGACGCGCCCTGGACAGCAGGTTCTCACGCTcc                    >  7:367206/2‑149 (MQ=39)
                                                                                                                             cTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAg                 >  4:18890/1‑149 (MQ=39)
                                                                                                                              ttCGCGGGCTTGCACGCGAAGGGTCGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAgg                >  2:365499/1‑149 (MQ=255)
                                                                                                                                cgcgTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGc               <  7:39653/148‑1 (MQ=255)
                                                                                                                                      cAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAgt        <  4:362277/149‑1 (MQ=255)
                                                                                                                                         gCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGAca                                     >  7:317463/1‑117 (MQ=255)
                                                                                                                                         gCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGAca                                     <  8:317463/117‑1 (MQ=255)
                                                                                                                                           acgcgAAGGGCTGCAACGCAACCCCGAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTtctc                           >  2:222502/2‑125 (MQ=255)
                                                                                                                                           acgcgAAGGGCTGCAACGCAACCCCGAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTtctc                           <  1:222502/124‑1 (MQ=255)
                                                                                                                                            cgcgAAGGGTCGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGgt  <  1:365499/149‑1 (MQ=255)
                                                                                                                                            cgcgAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGgt  <  3:18890/149‑1 (MQ=255)
                                                                                                                                             gcgAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGgt  <  5:142761/148‑1 (MQ=255)
                                                                                                                                             gcgAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGgt  >  6:142761/1‑148 (MQ=255)
                                                                                                                                                     |                                                                                                                                           
ATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGC‑TTTTGGGCTGAGCGTACCGGGG‑TTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGGT  >  NC_002947/5740337‑5740623

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