Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,061,752 G→T 17.4% F301L (TTC→TTA PP_3575 ← ferric siderophore receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,061,7520GT17.4% 56.7 / 3.1 23F301L (TTC→TTAPP_3575ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base G (11/8);  new base T (2/2);  total (13/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 = 1.60e-01

AGAGACGCAGGCGCTCGCCTCGAAAATCCAGGCCCAGCACCGCCATTTCGCGGTCGACGCGCTGGTGGTCCCATTCGGTATCACCGGCCTGCTTCACGCCATTGAAGCGCACGCCGAACTGCTGGCCTTCACCAAAGCGACGGCCGATGTCCACTGCGCCGCCAGCCTGGCCGGCCGAGGCGTAGCTGCCGGTGAATTCGGTGATCGGCTTGTCAGCGGCTCGCTTGGGCTGCACGTTGAT  >  NC_002947/4061637‑4061877
                                                                                                                   |                                                                                                                             
agagaCGCAGGCGCTCGCCTCGAAAATCCAGGCCCATCACCGCCATTTCGCGGTCGACGCGCTGGTGGTCCCATTCGGTATCACCGGCCTGCTTCACGCCATTGAAGCGCACGCCGAACTGCTGGCCTTCACCaaa                                                                                                           >  2:90011/1‑136 (MQ=255)
agagaCGCAGGCGCTCGCCTCGAAAATCCAGGCCCAGCACCGCCATTTCGCGGTCGACGCGCTGGTGGTCCCATTCGGTATCACCGGCCTGCTTCACGCCATTGAAGCGCACGCCGAACTGCTGGCCTTCACCaaa                                                                                                           >  2:271644/1‑136 (MQ=255)
agagaCGCAGGCGCTCGCCTCGAAAATCCAGGCACAGCACCGCCATTTCGCGGTCGACGCGCTGGTGGTCCCATTCGGTATCACCGGCCTGCTTCACGCCATTGAAGCGCACGCCGAACTGCTGGCCTTCACCaaa                                                                                                           >  1:392571/1‑136 (MQ=255)
      gCAGGCGCTCGCCTCGAAAATCCAGGCCCAGCACCGCCATTTCGCGGTCGACGCGCTGGTGGTCCCATTCGGTATCACCGGCCTGCTTCACGCCATTGAAGCGCACGCCGAACTGCTGGCCTTCACCAAAGCGACg                                                                                                     <  2:343735/136‑1 (MQ=255)
      gCAGGCGCTCGCCTCGAAAATCCAGGCCCAGCACCGCCATTTCGCGGTCGACGCGCTGGTGGTCCCATTCGGTATCACCGGCCTGCTTCACGCCATTGAAGCGCACGCCGAACTGCTGGCCTTCACCAAAGCGACg                                                                                                     >  2:297395/1‑136 (MQ=255)
        aGGCGCTCGCCTCGAAAATCCAGGCCCAGCACCGCCATTTCGCGGTCGACGCGCTGGTGGTCCCATTCGGTATCACCGGCCTGCTTCACGCCATTGAAGCGCACGCCGAACTGCTGGCCTTCACCAAAGCGACGGc                                                                                                   >  1:224770/1‑136 (MQ=255)
             cTCGCCTCGAAAATCCAGGCCCAGCACCGCCATTTCGCGGTCGACGCGCTGGTGGTCCCATTCGGTATCACCGGCCTGCTTCACGCCATTGAAGCGCACGCCGAACTGCTGGCCTTCACCAAAGCGACGGCCGATg                                                                                              >  1:229874/1‑136 (MQ=255)
                             aGGCCCAGCACCGCCATTTCGCGGTCGACGCGCTGGTGGTCCCATTCGGTATCACCGGCCTGCTTCACGCCATTGAAGCGCACGCCTAACTGCTGGCCTTCACCAAAGCGACGGCCGATGTCCACTGCGCCGCCAg                                                                              >  1:118597/1‑136 (MQ=255)
                             aGGCCCAGCACCGCCATTTCGCGGTCGACGCGCTGGTGGTCCCATTCGGTATCACCGGCCTGCTTCACGCCATTGAAGCGCACGCCTAACTGCTGGCCTTCACCAAAGCGACGGCCGATGTCCACTGCGCCGCCAg                                                                              >  1:445359/1‑136 (MQ=255)
                             aGGCCCAGCACCGCCATTTCGCGGTCGACGCGCTGGTGGTCCCATTCGGTATCACCGGCCTGCTTCACGCCATTGAAGCGCACGCCGAACTGCTGGCCTTCACCAAAGCGACGGCCGATGTCCACTGCGCCGCCAg                                                                              >  1:296685/1‑136 (MQ=255)
                               gCCCAGCACCGCCATTTCGCGGTCGACGCGCTGGTGGTCCCATTCGGTATCACCGGCCTGCTTCACGCCATTGAAGCGCACGCCGAACTGCTGGCCTTCACCAAAGCGACGGCCGATGTCCACTGCGCCGCCAGcc                                                                            >  2:242938/1‑136 (MQ=255)
                                  cAGCACCGCCATTTCGCGGTCGACGCGCTGGTGGTCCCATTCGGTATCACCGGCCTGCTTCACGCCATTGAAGCGCACGCCGAACTGCTGGCCTTCACCAAAGCGACGGCCGATGTCCACTGCGCCGCCAGCCTgg                                                                         >  1:326402/1‑136 (MQ=255)
                                             tttCGCGGTCGACGCGCTGGTGGTCCCATTCGGTATCACCGGCCTGCTTCACGCCATTGAAGCGCACGCCGAACTGCTGGCCTTCACCAAAGCGACGGCCGATGTCCACTGCGCCGCCAGCCTGGCCGGCAGAGGc                                                              >  2:402828/1‑136 (MQ=255)
                                                                  ggtCCCATTCGGTATCACCGGCCTGCTTCACGCCATTGAAGCGCACGCCGAACTGCTGGCCTTCACCAAAGCGACGGCCGATGTCCACTGCGCCGCCAGCCTGGCCGGCCGAGGCGTAGCTGCCGGTGAATTCGGt                                         <  2:229874/136‑1 (MQ=255)
                                                                                tCACCGGCCTGCTTCACGCCATTGAAGCGCACGCCTAACTGCTGGCCTTCACCAAAGCGACGGCCGATGTCCACTGCGCCGCCAGCCTGGCCGGCCGAGGCGTAGCTGCCGGTGAATTCGGTGATCGGCTTGTCAg                           <  2:118597/136‑1 (MQ=255)
                                                                                tCACCGGCCTGCTTCACGCCATTGAAGCGCACGCCTAACTGCTGGCCTTCACCAAAGCGACGGCCGATGTCCACTGCGCCGCCAGCCTGGCCGGCCGAGGCGTAGCTGCCGGTGAATTCGGTGATCGGCTTGTCAg                           <  2:445359/136‑1 (MQ=255)
                                                                                tCACCGGCCTGCTTCACGCCATTGAAGCGCACGCCGAACTGCTGGCCTTCACCAAAGCGACGGCCGATGTCCACTGCGCCGCCAGCCTGGCCGGCCGAGGCGTAGCTGCCGGTGAATTCGGTGATCGGCTTGTCAg                           <  2:296685/136‑1 (MQ=255)
                                                                                tCACCGGCCTGCTTCACGCCATTGAAGCGCACGCCGAACTGCTGGCCTTCACCAAAGCGACGGCCGATGTCCACTGCGCCGCCAGCCTGGCCGGCCGAGGCGTAGCTGCCGGTGAATTCGGTGATCGGCTTGTCAg                           <  2:446239/136‑1 (MQ=255)
                                                                                tCACCGGCCTGCTTCACGCCATTGAAGCGCACGCCGAACTGCTGGCCTTCACCAAAGCGACGGCCGATGTCCACTGCGCCGCCAGCCTGGCCGGCAGAGGCGTAGCTGCCGGTGAATTAGGTGATCGGCTTGTCAg                           <  1:402828/136‑1 (MQ=255)
                                                                                                     ttGAAGCGCACGCCGAACTGCTGGCCTTCACCAAAGCGACGGCCGATGTCCACTGCGCCGCCAGCCTGGCCGGCCGAGGCGTAGCTGCCGGTGAATTCGGTGATCGGCTTGTCAGCGGCTCGCTTGGGCTGCACGt      <  2:224770/136‑1 (MQ=255)
                                                                                                         aGCGCACGCCGAACTGCTGGCCTTCACCAAAGCGACGGCCGATGTCCACTGCGCCGCCAGCCTGGCCGGCCGAGGCGTAGCTGCCGGTGAATTTGGTGATCGGCTTGTCAGCGGCTCGCTTGGGCTGCACGTTGAt  <  2:392571/136‑1 (MQ=255)
                                                                                                                  cGAACTGCTGGCCTTCACCAAAGCGACGGCCGATGTCCACTGCGCCGCCAGCCTGGCCGGCCGAGGCGTAGCTGCCGGTGAATTCGGt                                         <  2:230174/88‑1 (MQ=255)
                                                                                                                  cGAACTGCTGGCCTTCACCAAAGCGACGGCCGATGTCCACTGCGCCGCCAGCCTGGCCGGCCGAGGCGTAGCTGCCGGTGAATTCGGt                                         >  1:230174/1‑88 (MQ=255)
                                                                                                                   |                                                                                                                             
AGAGACGCAGGCGCTCGCCTCGAAAATCCAGGCCCAGCACCGCCATTTCGCGGTCGACGCGCTGGTGGTCCCATTCGGTATCACCGGCCTGCTTCACGCCATTGAAGCGCACGCCGAACTGCTGGCCTTCACCAAAGCGACGGCCGATGTCCACTGCGCCGCCAGCCTGGCCGGCCGAGGCGTAGCTGCCGGTGAATTCGGTGATCGGCTTGTCAGCGGCTCGCTTGGGCTGCACGTTGAT  >  NC_002947/4061637‑4061877

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