Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 303,911 C→T G91E (GGA→GAA)  ykgJ ← putative metal‑chelating domain‑containing protein YkgJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913303,9110CT95.7% 65.1 / ‑6.2 23G91E (GGA→GAA) ykgJputative metal‑chelating domain‑containing protein YkgJ
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (14/8);  total (14/9)
Fisher's exact test for biased strand distribution p-value = 3.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

GCTTTTTATATTCTGATGCAGATAGTCAGTGAGTATATCGCGCTACTTCAGGATGATGTAGATCCGAAGAACGCTACAGAAGAGAGGCATTGTTGCTGGCAAATAGAAGAAGTATCGGGTTTTGTTACCCCTGAAAAACGAAGCCCCGCTATTATCGCTGGCGGGGCAGTGCAATTAATTATTGTATGTTTATAACGGTGTCAGCCCGTATTTAGCCCTTGCACGATTGCAAGCCTCATTGACGACTCCGTTCTCACCAGACATGGCGAATTCTCTGCATGTGGACGATCGATTTTTATATATCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAAC  >  NC_000913/303664‑304145
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                           
gCTTTTTATATTCTGATGCAGATAGTCAGTGAGTATATCGCGCTACTTCAGGATGATGTAGATCCGAAGAACGCTACAGAAGAGAGGCATTGTTGCTGGCAAATAGAAGAAGTATCGGGTTTTGTTACCCCTGAAAAACGAAGCCCCGCTATTATCGCTGGCGGGGCAGTGCAATTAATTATTGTATGTTTATAACGGTGTCATCCCGTATTTAGCCCTTGCACGATTGCAAGCCTCATTGACGACTTCGt                                                                                                                                                                                                                                          <  1:165635/251‑1 (MQ=255)
                  taGATAGTCAGTGAGAATATCGCGCTACTTCAGGATGATGTAGATCCGAAGAACGCTACAGAAGAGAGGCATTGTTGCTGGAAAATAGAAGAAGTATCGGGTTTTGTTACCCCTGAAAAACGAAGCCCCGCTATTATCGCTGGAGGGGCAGGGCAATTAATTATTGTATGTTTATAACGGTGTCAGCCCGAATTTAACCCTTGCACGATTGCAATCCTCATTGACGACTTCGTTCTCACCAGACATGGCgg                                                                                                                                                                                                                        >  2:60518/2‑250 (MQ=255)
                     atagtaagtgaatcaatCGCGCTACTTCAGGATGATGTAGAACCGAAGAACGCTACAGAAGAGAGGGATTGTTGCTGGCAAATAGAAGACGTATCGGGATTTGTTACCCCTGAAAAACGAAGCCCCGCTATTATAGCGTGCGGGGCAGTGCAATTAATTCTTGTGTGTTTATAACGGTGTCAGCCCGAATTTAGCCCGTGCACGTCTGCAAGCCTCATTGACGACTTCGTTCTCCCCCGACATGGCCAAtt                                                                                                                                                                                                                     >  2:248627/16‑251 (MQ=255)
                     aTAGTCAGTGAGTATATCGCGCTACTTCAGGATGATGTAGATCCGAAGAACGCTACAGAAGAGAGGCATTGTTGCTGGCAAATAGAAGAAGTATCGGGTTTTGTTACCCCTGAAAAACGAAGCCCCGCTATTATCGCTGGCGGGGCAGTGCAATTAATTATTGTATGTTTATAACGGTGTCAGCCCGTATTTAGCCCTTGCACGATTGCAAGCCTCATTGACGACTTCGTTCTCACCAGACATGGCGAAtt                                                                                                                                                                                                                     >  2:116255/1‑251 (MQ=255)
                          cAGTGAGTATATCGCGCTACTTCAGGATGATGTAGATCCGAAGAACGCTACAGAAGAGAGGCATTGTTGCTGGCAAATAGAAGAAGTATCGGGTTTTGTTACCCCTGAAAAACGAAGCCCCGCTATTATCGCTGGCGGGGCAGTGCAATTAATTATTGTATGTTTATAACGGTGTCAGCCCGTATTTAGCCCTTGCACGATTGCAAGCCTCATTGACGACTTCGTTCTCACCAGACATGGCGAATTctct                                                                                                                                                                                                                 <  1:116255/250‑1 (MQ=255)
                                 tataTCGCGCTACTTCAGGATGATGTAGATCCGAAGAACGCTACAGAAGAGAGGCATTGTTGCTGGCAAATAGAAGCAGTATCGGGTTTTGTTACCCCTGAAAAACGAAGCCCCGCTATTATCGCTGGCGGGGCAGTGCAATTAATTATTGTATGTTTATAACGGTGTCAGCCCGTATTTAGCCCTTGCACGATTGCAAGCCTCATTGACGACTTCGTTCTCACCAGACCTGGCGAATTCTCTGCAtgtg                                                                                                                                                                                                          <  2:233536/250‑1 (MQ=255)
                                        cgcTACTTCAGGATGATGTAGATCCGAAGAACGCTACAGAAGAGAGGCATTGTTGCTGGCAAATAGAAGAAGTATCGGGTTTTGTTACCCCTGAAAAACGAAGCCCCGCTATATCGCTGGCGGGGCAGTGCAATTAATTATTGTATGTTTATAACGGTGTCAGCCCGTATTTAGCCCTTGCACGATTGCAAGCCTCATTGACGACTTCGTTCTCACCAGACATGGCGAATTCTCTGCATGTGGAcgatcga                                                                                                                                                                                                 <  1:84520/251‑1 (MQ=255)
                                          cTACTTCAGGATGATGTAGATCCGAAGAACGCTACAGAAGAGAGGCATTGTTGCTGGCAAATAGAAGAAGTATCGGGTTTTGTTACCCCTGAAAAACGAAGCCCCGCTATTATCGCTGGCGGGGCAGTGCAATTAATTCTTGTATGTTTATAACGGTGTCAGCCCGTATTTAGCCCTTGCACGATTGCAAGCCTCATTGACGACTTCGTTCTCACCAGACATGGCGAATTCTCTGAATGTGGAcgatcgat                                                                                                                                                                                                >  1:30037/1‑251 (MQ=255)
                                           tACTTCAGGATGATGTAGATCCGAAGAACGCTACAGAAGCGAGGCATTGTTGCTGGCAACTAGAAGACGTATCGGGTTTTGTTACCCCTGAAAAACGAAGCCCCGCTATTATCGCTGGCGGGGCAGTGCCATTAATTATTGTATGTTTATAAAGGTGTCAGCCCGTATTTAGCCCTTGCCCGATTGCAAGCCTCATTGACGCCTTCGTTCTCCCCAGACATGGCGCATTCTCTGCATGTGGAccatcgctt                                                                                                                                                                                               >  2:137156/1‑248 (MQ=255)
                                                                                                             aagTATCGGGTTTTGTTACCCCTGAAAAACGAAGCCCCGCTATTATCGCTGGCGGGGCAGAGCAATTAATTATTGTATGTTTATAACGGTGTCAGCCCGTATTTAGCCCTTGCACGATTGCAAGCCTCATTGACGACTTCGTTTTCACCAGACATGGCGAATTCTCTGCATGTGGACGATCGATTTTTATATATCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGttttt                                                                                                                             >  1:145620/1‑251 (MQ=255)
                                                                                                                   cGGGTTTTGATATCCCTGAAAAACGAAGCCACGCTATTATCGCTGGCGGGGCAGTGCAATTAATTATTGTATGTTTCTAACGGTGTCAGACCGTATTTAGCCCTTGCACGATTGCAAGCCTCATTGACGACTTCGTTCTCACCAGACCTGGCGAATTCTCTGCATGTGGCCGATCGCTTTTTATATATCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGAtt                                                                                                                       >  2:119914/1‑251 (MQ=255)
                                                                                                                              aCCCCTGAAAAACGAAGCCCCGCTATTATCGCTGGCGGGGCAGTGCAATTAATTATTGTATGTTTATAACGGTGTCAGCCCGTATTTAGCCCTTGCACGATTGCAAGCCTCATTGACGACTTCGTTCTCACCAGACATGGCGAATTCTCTGCATGTGGACGATCGATTTTTATATATCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGGTTTTCTGATTGGTACCGCTc                                                                                                             <  2:124852/250‑1 (MQ=255)
                                                                                                                                             aGCCCCGCTATTATCGCTGGCGGGGCAGTGCAATTAATTATTGTATGTTTATAACGGTGTCAGCCCGTATTTAGCCCTTGCACGATTGCAAGCCTCATTGACGACTTCGTGCTCACCAGACATGGCGAATTCTCTGCATGTGGACGATCGATTTTTATATATCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACAGCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAg                                                                                            <  1:93254/251‑1 (MQ=255)
                                                                                                                                                   gCTATTATCGCTGGCGGGGCCGTGCAATTAATTATTGTATGTTTATAACGGTGTCAGCCCGTATTTAGCCCTTGCACGATTGCAAGCCTCATTGACGACTTCGTTCTCACCAGACATGGCGAATTCTCTGCCTGTGGACGATCGATTTTTATATATCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGGAAAGGGGatt                                                                                      >  1:88619/1‑250 (MQ=255)
                                                                                                                                                                                           gTTTATAACGGAGTCAGCCCGTATTTAGCCCTTGCACGATTGCAAGCCTCATTGACGACTTCGTTCTCACCAGACATGGCGAATTCTCTGCATGTGGGCGATCGATTTTTATATATCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAATGGCAATAAATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATAGGTGAAAAGGGGATATTTGCTCAGTGAGCCTGGCGGGAATCGTACCGCCAGCAtg                                              >  2:22560/1‑250 (MQ=255)
                                                                                                                                                                                                        tCAGCCCGTATTTAGCCCTTGCACGATTGCAAGCCTCATTGACGCCTTCGTTCTCACCAGACATGGCGAATTCTCTGCATGTGGACGATCGATTTTTATATATCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGcc                                 >  2:115387/1‑251 (MQ=255)
                                                                                                                                                                                                            cccGTATTTAGCCCTTGCACGATTGCAAGCCTCATTGACGACTTCGTTCTCACCAGACATGGCGAATTCTCTGCATGTGGACGATCGATTTTTATATATCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAgg                             >  2:54442/1‑250 (MQ=255)
                                                                                                                                                                                                                      gCCCTTGCACGATTGCAAGCCTCATTGACGACTTCGTTCTCACCAGACATGGCGAATTCTCTGCATGTGGACGATCGATTTTTATATATCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCCGCATCGTCGGCTTCGGCCCAGTAAAAAgaga                    >  1:158866/1‑250 (MQ=255)
                                                                                                                                                                                                                          gtgcacgatggcaaccctgttggccgacTCCTTTCTCACAAGACATGGCGATTTCTCTGCATGTGTATGACCGATTTTTATATATCGTGCAGCAGGCTTTTTTGCCCGGGTTTCCTGAAAGGGCAATACATCGGGGGTTTTTCGGATTGTTACCGCTCATGCACCGTTGAAAGGGGGATTTTTGTCCAGGGAGCCGGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCgg               <  2:88619/227‑1 (MQ=255)
                                                                                                                                                                                                                             cACGATTGCAAGCCGCATTGGCGACTTAGTTCTCCACAGAAAGGGCAAAGTCTATGCATGGGGGCGAGCTTTTTTGATATATCGTGCAGCAGGCATTTTTGCCCGGGTTTCCTGCAGGGGCAATACATGTGGGTTTTTTCTGATGGGTACCTCGCTGACATCGGTGAAAAGGGGATGTTTGGTCAGTGATCCGTGTGTGAATAGTACCGCCAGCATCGTCGGTTTCGGCCCAGTAAAAAGAGACGCGGaaa            <  2:126767/251‑1 (MQ=255)
                                                                                                                                                                                                                                     cAAGCCTCATGGACGACTTCGTCCTCACAAGCCATGGCGATTTCTCGGCTTGTGGACGATCGATTTTTATATAGGGGGCAGCAGGCTTTTTTGCCCGGGGTCCCTGAAAGGGCAATACATCGGGGGTTTTTTTGATTGGTCCCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGGGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGcaca    <  1:54442/251‑1 (MQ=255)
                                                                                                                                                                                                                                      aaGCCTCATTGACGACTTCGTTCTCACCAGACATGGCGAATTCTCTGCATGTGGACGATCGATTTTTATATATCGTGCAGCAGGCATTTTTGCCAGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAACAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGAATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGcaca    >  1:160346/1‑250 (MQ=255)
                                                                                                                                                                                                                                       aGCCTCATTGACGACTTCGATCTCACCAGACATGGCGAATTCTCTGCATGTGGACGATCGATTTTTCTATATCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAACAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGGAAAAAGAGAAGCGGAAAAATGcacacc  >  2:236612/1‑249 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                           
GCTTTTTATATTCTGATGCAGATAGTCAGTGAGTATATCGCGCTACTTCAGGATGATGTAGATCCGAAGAACGCTACAGAAGAGAGGCATTGTTGCTGGCAAATAGAAGAAGTATCGGGTTTTGTTACCCCTGAAAAACGAAGCCCCGCTATTATCGCTGGCGGGGCAGTGCAATTAATTATTGTATGTTTATAACGGTGTCAGCCCGTATTTAGCCCTTGCACGATTGCAAGCCTCATTGACGACTCCGTTCTCACCAGACATGGCGAATTCTCTGCATGTGGACGATCGATTTTTATATATCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAAC  >  NC_000913/303664‑304145

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