Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,015,369 C→T G267G (GGG→GGA ysgA ← putative dienelactone hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,015,3690CT95.5% 56.7 / NA 23G267G (GGG→GGAysgAputative dienelactone hydrolase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (7/14);  minor base G (0/1);  total (8/15)
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.00e+00

GGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATA  >  NC_000913/4015160‑4015574
                                                                                                                                                                                                                 |                                                                                                                                                                                                             
ggATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTTCCACCATACTGc                                                                                                                                                                                                   <  1:13033/222‑1 (MQ=255)
                                cGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGTCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTTCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGc                                                                                                                                      <  1:464956/251‑1 (MQ=255)
                                     cccGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTTCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATa                                                                                                                                  >  1:38042/2‑250 (MQ=255)
                                                           aaCCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTTCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGc                                                                                                                                      >  1:32089/1‑224 (MQ=255)
                                                                         ccTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGGGCAGGCCGCGCCGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGCAATACCACCCGGCCTTTTCTCATTACAGCGACTTCTTTCCACCATACTGCTTAATCCATTCCCGCATACGCTACCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAAAGCATGCACGGCGTCCGGGTACACg                                                                                             >  1:219944/1‑251 (MQ=255)
                                                                                        ntGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTTCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGGACACGATAATCTCTGCTTTc                                                                              >  1:129567/2‑251 (MQ=255)
                                                                                            ccAGAAACCCGCGCGGCACTGGCTAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTTCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCg                                                                                                         <  1:162444/220‑1 (MQ=255)
                                                                                            ccAGAAAACCGCGCGGCACTGGCGGACTAGGCGCAGGCCGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGAGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTTCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCAt                                                                          >  1:420782/1‑251 (MQ=255)
                                                                                               gAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTTCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTan                                                                       <  1:301995/251‑2 (MQ=255)
                                                                                                     cgcgcgGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTTCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAg                                                                                                                             <  1:282/191‑1 (MQ=255)
                                                                                                            cACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTTCCACCAAACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCTTGGAATGCAGGCCCGTCGTCCGGTTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGcc                                                         <  1:88372/251‑1 (MQ=255)
                                                                                                                            gCAGGCGGCGCAGAACTTTCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTTCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTc                                           <  1:439089/250‑1 (MQ=255)
                                                                                                                             cAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTTCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCg                                                                             <  1:161777/215‑1 (MQ=255)
                                                                                                                                 nggcgCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTTCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCAtcttc                                                                                                                                                                >  1:92972/2‑128 (MQ=255)
                                                                                                                                   gcgcAGAACTTGCGTCGGGTGTAAAATCCAAACCGTGTGGTAATACCACCCGGTCGTTTCTCATTACAGCGACTTCTGGCCGCCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGGAGGAGAGCGCGGGCGATAATGAGCTTTGAATGCATGCCCGGCGTCCGGGTACGCGATTATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTCAACGctct                                   <  1:493825/251‑1 (MQ=255)
                                                                                                                                    ctcAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTTTTTTATCATTACAGAGACTTATTACCACCACACTGCTTACACCATTCCAGCATCCGCTgccagcca                                                                                                                                                                     >  1:468103/3‑120 (MQ=255)
                                                                                                                                      cAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTTCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATg                                                <  1:288650/235‑1 (MQ=255)
                                                                                                                                                    ggggTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTTCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCGCGGGCGATAATCATCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCgn                                                                            <  1:105796/193‑2 (MQ=255)
                                                                                                                                                       gTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTTCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATcn               <  1:133452/251‑2 (MQ=255)
                                                                                                                                                                    cGGGTGGTAAGACCTCCCGGGCTTTTCTCATTACAGCGACTTCTGTCCAGCATCCTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCAtata  <  1:245062/251‑1 (MQ=255)
                                                                                                                                                                          gtAATACCACCCGGTCTTTTCTCATTACAGCGCCTTCTTTCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGACAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGCATGCTGTTATCCTGACCACCata    >  1:387437/1‑243 (MQ=255)
                                                                                                                                                                                              ctcATTACAGCGACTTCTTTCCACCATACTGCTTAAAc                                                                                                                                                                                             <  1:245835/38‑1 (MQ=255)
                                                                                                                                                                                                        cGACTTCTTTCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCgg                            <  1:178881/189‑1 (MQ=255)
                                                                                                                                                                                                                 |                                                                                                                                                                                                             
GGATTGAAGCCTTGCTGAAGAAAGCGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATA  >  NC_000913/4015160‑4015574

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: