Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 623,110 C→A 23.6% intergenic (+59/+103) PP_0535 → / ← PP_0536 ferric siderophore receptor/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947623,1100CA23.6% 25.9 / 5.1 17intergenic (+59/+103)PP_0535/PP_0536ferric siderophore receptor/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (5/8);  new base A (2/2);  total (7/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.05e-01

GGATTACCTGTCGCGGGTGGGCTCGACCTCCACGTTCAACTTCTACGGGCCGTCGCGCAGCATGATGGTTGGAGCGCGGTATGACTTCTGACGGGGTCTGTTGCCGCTAGGTGCACCGCGACACAAGGCCGTGCCTG‑‑CAAAGGGCGCGGCCTTGATGCTGAATCATGAGAAATCTAGAAGTGTATCTGCACTATCGGGTATCAGTGCAATAAAGTTGTGAGAGTAAGTGTTAATCAGGCGCTTGGTTTTAATTCAAGAAGTCTCCGA  >  NC_002947/622961‑623227
                                                                                                                                                       |                                                                                                                     
ggATTACCTGTCGCGGGTGGGCTCGACCTCCACGTTCAACTTCTACGGGCCGTCGCGCAGCATGATGGTTGGAGCGCGGTATGACTTCTGACGGGGTCTGTTGCCGCTAGGTGCACCGCGACACAAGGCCGTGCCTG‑‑CAAAGGGCGCGGc                                                                                                                       >  1:467166/1‑150 (MQ=255)
        tGTCGCGGGTGGGCTCGACCTCCACGTTCAACTTCTACGGGCCGTCGCGCAGCATGATGGTTGGAGCGCGGTATGACTTCTGACGGGGTCTGTTGCCGCTAGGTGCACCGCGACACAAGGCCGTGCCTG‑‑CAAAGGGCGCGGCCTTGATGc                                                                                                               <  2:442865/150‑1 (MQ=255)
                gTGGGCTCGACCTCCACGTTCAACTTCTACGGGCCGTCGCGCAGCATGATGGTTGGAGCGCGGTATGACTTCTGACGGGGTCTGTTGCCGCTAGGTGCACCGCGACACAAGGCCGTGCCTG‑‑CAAAGGGCGCGGCCTTGATGCTGAATCAt                                                                                                       <  1:6173/150‑1 (MQ=255)
                     cTCGACCTCCACGTTCAACTTCTACGGGCCGTCGCGCAGCATGATGGTTGGAGCGCGGTATGACTTCTGACGGGGTCTGTTGCCGCTAGGTGCACCGCGACACAAGGCCGTGCCTG‑‑CAAAGGGCGCGGCCTTGATGCTGAATCATGAGaa                                                                                                  >  2:352287/1‑150 (MQ=255)
                      tCGACCTCCACGTTCAACTTCTACGGGCCGTCGCGCAGCATGATGGTTGGAGCGCGGTATGACTTCTGACGGGGTCTGTTGCCGCTAGGTGCACCGCGACACAAGGCCGTGCCTG‑‑CAAAGGGCGCGGCCTTGATGCTGAATCATGAGaaa                                                                                                 >  1:435857/1‑150 (MQ=255)
                             ccACGTTCAACTTCTACGGGCCGTCGCGCAGCATGATGGTTGGAGCGCGGTATGACTTCTGACGGGGTCTGTTGCCGCTAGGTGCACCGCGACACAAGGCCGTGCCTG‑‑CAAAGGGCGCGGCCTTGATGCTGAATCATGAGAAATCTAGaa                                                                                          <  1:350150/150‑1 (MQ=255)
                                                         cagtaTGATGGTTGGGGCGGGGTATGACTTCTGACGGGGTCTGTTGCCGCTAGGTGCACCGCGACACAAGGCCGTGCCTG‑‑CAAAGGGCGCGGCCTTGATGCTGAACCATGAGAAATCTAGAAGGGTATCTGCACTATCGGGTATCAGTGc                                                              <  2:406098/150‑1 (MQ=255)
                                                                  ggTTGGAGCGCGGTTTGACTTCTGACGGGGTCTGTTGCCGCTAGGTGCACTGCGCCTTTGCTGAACAGGATAAAAGGTCGCGGCCTTGATTCTGAATCATGAAAAATCTAGAAATGTATCTGCACTATCGGGTATCAGTGCAATAAAGtt                                                     <  1:369409/150‑1 (MQ=255)
                                                                                          aCGGGGTCTGTTGCCGCTAGGTGCACCGCGACACAAGGCCGTGCCTG‑‑CAAAGGGCGCGGACTTGATGCTGAATCATGAGAAATCTAGAAGTGTATCTGCACTATCGGGTATCAGTGCAATaaa                                                        <  1:121678/123‑1 (MQ=255)
                                                                                          aCGGGGTCTGTTGCCGCTAGGTGCACCGCGACACAAGGCCGTGCCTG‑‑CAAAGGGCGCGGACTTGATGCTGAATCATGAGAAATCTAGAAGTGTATCTGCACTATCGGGTATCAGTGCAATaaa                                                        >  2:121667/1‑123 (MQ=255)
                                                                                          aCGGGGTCTGTTGCCGCTAGGTGCACCGCGACACAAGGCCGTGCCTG‑‑CAAAGGGCGCGGACTTGATGCTGAATCATGAGAAATCTAGAAGTGTATCTGCACTATCGGGTATCAGTGCAATaaa                                                        >  2:121678/1‑123 (MQ=255)
                                                                                          aCGGGGTCTGTTGCCGCTAGGTGCACCGCGACACAAGGCCGTGCCTG‑‑CAAAGGGCGCGGACTTGATGCTGAATCATGAGAAATCTAGAAGTGTATCTGCACTATCGGGTATCAGTGCAATaaa                                                        <  1:121667/123‑1 (MQ=255)
                                                                                                  tGTTGCCGCTAGGTGCACCGCGACACAAGGCCGTGCCTG‑‑CAAAGGGCGCGGCCTTGATGCTGAATCATGAGAAATCTAGAAGTGTATCTGCACTATCGGGTATCAGTGCAATAAAGTTGTGAGAGTAAGTGTTAATCAGGCGCTTGGttt                     >  1:465955/1‑150 (MQ=255)
                                                                                                         gCTAGGTGCACCGCGACACAAGGCCGTGCCTG‑‑CAAAGGGCGCGGCCTTGATGCTGAATCATGAGAAATCTAGAAGTGTATCTGCACTATCGGGTATCAGTGCAATAAAGTTGTGAGAGTAAGTGTTAATCAGGCGCTTGGTTTTAATTCa              >  1:378224/1‑150 (MQ=255)
                                                                                                          cTAGGTGCACCGCGACACAAGGCCGTGCCTG‑‑CAAAGGGCGCGGCCTTGATGCTGAATCATGAGAAATCTAGAAGTGTATCTGCACTATCGGGTATCAGTGCAATAAAGTTGTGAGAGTAAGTGTTAATCAGGCGCTTGGTTTTAATTCaa             <  2:435857/150‑1 (MQ=255)
                                                                                                                     gcgACACAAGGCCGTGCCTG‑‑CAAAGGGCGCGGCCTTGATGCTGAATCATGAGAAATCTAGAAGTGTATCTGCACTATCGGGTATCAGTGCAATAAAGTTGTGAGAGTAAGTGTTAATCAGGCGCTTGGTTTTAATTCAAGAAGTCTCCgg  <  1:352287/150‑2 (MQ=255)
                                                                                                                                  gTGCCTG‑‑CAAAGGGCGCGGCCTTGATGCTGAATCATGAGAAATCTAGAAGTGTATCTGCACTATCGGGTATCAGTGCAATAAAGTTGTGAGAGTAAg                                          >  1:140876/1‑97 (MQ=255)
                                                                                                                                  gTGCCTG‑‑CAAAGGGCGCGGCCTTGATGCTGAATCATGAGAAATCTAGAAGTGTATCTGCACTATCGGGTATCAGTGCAATAAAGTTGTGAGAGTAAg                                          <  2:140876/97‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                     
GGATTACCTGTCGCGGGTGGGCTCGACCTCCACGTTCAACTTCTACGGGCCGTCGCGCAGCATGATGGTTGGAGCGCGGTATGACTTCTGACGGGGTCTGTTGCCGCTAGGTGCACCGCGACACAAGGCCGTGCCTG‑‑CAAAGGGCGCGGCCTTGATGCTGAATCATGAGAAATCTAGAAGTGTATCTGCACTATCGGGTATCAGTGCAATAAAGTTGTGAGAGTAAGTGTTAATCAGGCGCTTGGTTTTAATTCAAGAAGTCTCCGA  >  NC_002947/622961‑623227

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: