Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,740,484 T→C 14.7% 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,4840TC14.7% 72.5 / 6.5 34intergenic (+49/+22)hutF/PP_5037formiminoglutamate deiminase/lipocalin family lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (12/17);  new base C (2/3);  total (14/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.94e-01

ATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGGTGTGCGCGA  >  NC_002947/5740337‑5740631
                                                                                                                                                     |                                                                                                                                                   
aTGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGtt                                                                                                                                                    <  3:157541/149‑1 (MQ=255)
 tGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTg                                                                                                                                                   >  8:150704/1‑149 (MQ=255)
   gTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCGAAGGGCTgc                                                                                                                                                  <  8:173354/148‑1 (MQ=38)
      aaTGGGCAATGGGGGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTgcaac                                                                                                                                               <  7:150704/148‑1 (MQ=255)
                      tgCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTTGAGGCCCTAGCCGGGGAGCCCGCGAAGGGTTTGAACCGCACCCCaaaaaaaa                                                                                                                               >  2:276589/1‑146 (MQ=255)
                          ccAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGGGGTTGCGCCCATAGCGGGGCAAGCCCCGCACGGGTTGGAACCCCACCCCaaaaaaaacctg                                                                                                                           >  3:315022/1‑143 (MQ=255)
                                                     aaGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTGGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGAcgc                                                                                               <  3:81467/149‑1 (MQ=255)
                                                           gcgggcgTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCa                                                                                         <  8:238596/149‑1 (MQ=255)
                                                                              ttGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTc                                                                       <  1:16493/148‑1 (MQ=255)
                                                                                 cgcCAGGTTTTGGGCGGAGGGTACCGGGGTGGCGTTGCAGCCGTGGGCGTGCAAGCCGGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGTCGCCCGCCAGATCAAGCGGCTGGTGTCATa                                                                    <  1:196623/148‑1 (MQ=255)
                                                                                   ccAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTAGCGGGCAAGCCCGCGAAGGGTTGGCACGCCACCCCCAAAAACGCAGTGGATTTTTTCGATTTGCGGATCGGACGCCCGCCAGATCAAGCGGCGGGGGTCATAacc                                                                 >  2:45521/1‑149 (MQ=255)
                                                                                    cAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTgc                                                                                                                                                  >  2:316016/1‑67 (MQ=255)
                                                                                     aGCTTTTGGGCTGAGCGTACCGGGGTGGCGTTCCAGCCCGTCGCGGCCAAGCCCGCGAAGGGTGGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCct                                                               <  4:304944/149‑1 (MQ=255)
                                                                                        ttttGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCACATCAAGCGGCTGGTGTCATACCCctgct                                                            <  1:187067/149‑1 (MQ=255)
                                                                                          ttGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCctgct                                                            <  2:259385/147‑1 (MQ=255)
                                                                                                  aGCGTACCGGGGTGGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGATGGGTTGAAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCt                                                  <  4:315022/149‑1 (MQ=255)
                                                                                                   gCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTg                                                 <  1:45521/149‑1 (MQ=255)
                                                                                                             gttgcgttgcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGt                                       >  2:127446/1‑149 (MQ=255)
                                                                                                                    ttgcAACCCTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCAc                                 >  5:120314/1‑149 (MQ=38)
                                                                                                                       cAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGAAACCCCAAAAAACTCAGTGCAATTTTACTATTTGAAGAACGTACTCACGCAAGATAAAGCGGCTGGTGTCATACCCCTTCTGCCGCGCACTCGACAGCAGGTTCTCACGCt                              >  6:125138/1‑149 (MQ=255)
                                                                                                                         gCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGCTCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTcc                            >  8:188150/1‑149 (MQ=255)
                                                                                                                            cTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAg                         >  4:231948/1‑149 (MQ=255)
                                                                                                                            cTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGcc                                                         >  5:128440/1‑117 (MQ=38)
                                                                                                                            cTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGcc                                                         <  6:128440/117‑1 (MQ=38)
                                                                                                                             ttCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAgg                        >  3:6168/1‑149 (MQ=255)
                                                                                                                              tCGCGGGCTTGCACGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCt                                                  <  5:104851/122‑1 (MQ=255)
                                                                                                                              tCGCGGGCTTGCACGCGAAGGGTTGCAAAGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCt                                                  >  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)
                                                                                                                                          cacgcgAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCg            <  6:120314/147‑1 (MQ=255)
                                                                                                                                            cgcgAAGGGTTGCAACGCAACCCCATTATACTCTGTGCATTTTTACTATTTTCTGATCGGTCGCCCGCCAGATCTAGCGGCTGGTGTCATACCCATGCTGCCGCGCCCTTGTCATCAGGTTCTGACGCTCACAGGCCTGCAGTGTCGgt          <  5:125138/149‑1 (MQ=255)
                                                                                                                                                  gggTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGGTGTgcgc    <  7:11144/149‑1 (MQ=255)
                                                                                                                                                    gTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGGTGTGCGCGa  <  1:276589/149‑1 (MQ=255)
                                                                                                                                                     ttGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGATAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGGTGTGCGCGa  >  6:301125/1‑148 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
ATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGGTGTGCGCGA  >  NC_002947/5740337‑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: