Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,740,457 G→A 18.5% intergenic (+22/+49) 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,4570GA18.5% 48.5 / 6.3 27intergenic (+22/+49)hutF/PP_5037formiminoglutamate deiminase/lipocalin family lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (7/15);  new base A (3/2);  total (10/17)
Fisher's exact test for biased strand distribution p-value = 3.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.69e-02

CGGTGGCGATCGCCAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCAC  >  NC_002947/5740309‑5740600
                                                                                                                                                    |                                                                                                                                               
cGGTGGCGATCGCCAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTTCAg                                                                                                                                                 >  7:91677/1‑149 (MQ=255)
           gCCAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGt                                                                                                                                       <  7:240654/148‑1 (MQ=255)
                 tGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTTCAGCCCTTCGCCTGCAAAcc                                                                                                                                >  2:29867/1‑149 (MQ=255)
                        tgtgATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAg                                                                                                                         <  1:264308/149‑1 (MQ=255)
                         gtgATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAgg                                                                                                                        >  1:22292/1‑149 (MQ=255)
                          tgATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAggg                                                                                                                       <  3:86837/149‑1 (MQ=255)
                            aTGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCGaa                                                                                                                          <  3:14799/144‑1 (MQ=39)
                            aTGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGACCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGAGTACAGGGTTTGCGTTGCAACCCTTCGAGGGCTTGCACGCGaa                                                                                                                          >  4:14799/1‑144 (MQ=255)
                             tGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTg                                                                                                                    <  1:188850/149‑1 (MQ=255)
                                tGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTgcaa                                                                                                                 <  2:240160/149‑1 (MQ=255)
                                       gCAATGGGTGGTGCGCCAGCGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTTGAGCCCATACCCTGGAAGCGCGCCCCGGGGTGGAACGCCAcc                                                                                                          >  4:117549/1‑149 (MQ=255)
                                         aaTGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTgc                                                                                                                   >  5:219598/1‑138 (MQ=255)
                                         aaTGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTgc                                                                                                                   <  6:219598/138‑1 (MQ=255)
                                              gtggtgCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCCTTTCAGCCCATAGCGTGCAAGCCCGCGAAGGGTTGGAACGCCACCCCaaaa                                                                                                    >  8:261403/1‑148 (MQ=255)
                                                  tgCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTCCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTc                                                                                               <  1:29867/149‑1 (MQ=255)
                                                                     ccAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGAAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTg                                                                            <  2:22292/149‑1 (MQ=255)
                                                                                aatgcgggcggGCGTTTGCCATGGTGTTGCGCCAGGTTTTGGGCTTAGCTTACGGGGGTTGCGTTCCAGCCGTTCGCTTCCAAGCCCGCTAAGGGTTGAAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGAc                                                                  <  3:323067/145‑1 (MQ=255)
                                                                                aaaGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCTAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGAcg                                                                 <  7:261403/149‑1 (MQ=255)
                                                                                  aGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCGAAGGGCTg                                                                                                                    <  5:69301/96‑1 (MQ=21)
                                                                                  aGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCGAAGGGCTg                                                                                                                    >  6:69301/1‑96 (MQ=21)
                                                                                                   cAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTTGAGCCCATAGCGGGCCAGCCCGCGCACGGTTTGAAACCAACCCCCAAAAACAAAATGGAATTTTTATATTTTGTGATAGGACGCCCCCCAGAAAAAGCGGcg                                              >  8:63876/1‑148 (MQ=255)
                                                                                                     aGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTgg                                            >  5:268436/1‑149 (MQ=255)
                                                                                                       gtgtGGCGCCAGGTTTTGGGCTGAGCGTACGGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCTAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGgtg                                          <  2:213682/149‑1 (MQ=255)
                                                                                                          ttGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTc                                        <  3:63922/148‑1 (MQ=255)
                                                                                                             cgcCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCCGCCCTTCGCGTGCAAGCCCGCTATGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATAc                                    <  6:313398/149‑1 (MQ=255)
                                                                                                             cgcCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATAc                                    <  1:101499/149‑1 (MQ=255)
                                                                                                                                      ggggTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGAcagca           <  7:63876/149‑1 (MQ=255)
                                                                                                                                               ttgcAACCCTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCAc  >  2:201889/1‑149 (MQ=38)
                                                                                                                                                    |                                                                                                                                               
CGGTGGCGATCGCCAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCAC  >  NC_002947/5740309‑5740600

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: