Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,233,540 A→G intergenic (‑307/‑218) lysC ← / → pgi aspartate kinase III/glucose‑6‑phosphate isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,233,5400AG92.9% 34.3 / ‑6.0 14intergenic (‑307/‑218)lysC/pgiaspartate kinase III/glucose‑6‑phosphate isomerase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (7/6);  total (8/6)
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

GCTACTATTTTCACCCAGAAGTGCTCCACCACTTGCGAAACGCCCGACTGCGAACGCTTCTGGTGACAACCCAGGGGATTCAGCCCCTGTAGCCGATGATGAACGTGGCCAGCCGTTCAATCACCTCGGCGATGCACCCCCTCAGGTGTTATCACAGGACTGGCTCCTCCAACACCGTTACTTGGGCAACGCGCCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACGC  >  NC_000913/4233315‑4233779
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                               
gCTACTATTTTCACCCAGAAGTGCTCCACCACTTGCGAAACGCCCGACTGCGAACGCTTCTGGTGACAACCCAGGGGATTCAGCCCCTGTAGCCGATGATGAACGTGGCCAGCCGTTCAATCACCTCGGCGATGCACCCCCTCAGGTGTTATCACAGGACTGGCTCCTCCAACACCGTTACTTGGGCAACTCGCCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACGTTTATAAATAAAGGGTGAGCGGGGc                                                                                                                                                                                                                        <  1:184781/251‑1 (MQ=255)
                aGAAGTGCTCCACCACTTGCGAAACGCCCGACTGCGAACGCTTCTGGTGACAACCCAGGGGATTCAGCCCCTGTAGCCGATGATGAACGTGGCCAGCCGTTCAATCACCTCGGCGATGCACCCCCTCAGGTGTTATCACAGGACTGGCTCCTCCAACACCGTTACTTGGGCAACGCGCCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACGTTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATggn                                                                                                                                                                                                        <  1:29275/251‑2 (MQ=255)
                                                      cGCTTCTGGTGACAACCCAGGGGATTCAGCCCCTGTAGCCGATGATGAACGTGGCCAGCCGTTCAATCACCTCGGCGATGCACCCCCTCAGGTGTTATCACAGGACTGGCTCCTCCAACTCCGTTACTTTGGCAACGCGCCTCTTCTGGCCTGCGCTAGCGCAGGTAGTGCGGTTATAAATAAAGGGTTAGTGTGGCTGTTGTCAACGATGGGGTCATGCTGATTTTTCATCCACTCCTGGCGGTCagtag                                                                                                                                                                  <  1:136584/251‑1 (MQ=255)
                                                                                     tctcgtccgatctaggaACGTGGCCAGCCTTTCAATCACCTCGGCGATGCACCCCCGCAGGTGTTATCACAGGACTGGCTCCTCCAACACCGTTACTTGGGCAACGCGCGTCGTCTGGACGGCGCTAGTGCAGGTAGTACGGTTATAAAGAAAGGGTGAGCGGGTCGGTTGTCAACGATGGGGTCATGTGGAGTTTTCATCCACTCCTCGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAg                                                                                                                                    <  1:171494/235‑1 (MQ=255)
                                                                                     ncTGTAGCCGATGATGAACGTGGCCAGCCGTTCAATCACCTCGGCGATGCACCCCCTCAGGTGTTATCACAGGACTGGCTCCTCCAACACCGTTACTTGGGCAACGCGCCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACGTTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATggg                                                                                                                                                                                                        >  1:61111/2‑182 (MQ=255)
                                                                                                                         cACCTCGGCGATGCACCCCCTCAGGTGTTATCACAGGACTGGCTCCTCCAACACCGTTACTTGGGCAACGCGCCTCTTCTGTCCTGCGCTAGCGCAGGTAGTACGTTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAAcc                                                                                               <  1:180096/251‑1 (MQ=255)
                                                                                                                                       aCCCCCTCAGGTGTTATCACAGGACTGGCTCCTCCAACACCGTTACTTGGGCAACGCGCCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACGTTTATAAATAAAGGGTGAGCGGGGCTGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAAc                                                                                                           <  1:267641/225‑1 (MQ=255)
                                                                                                                                                            ngACTGGCTCCTCCAACACCGTTACTTGGGCAACGCGCCTCTTCTGGCCTGCGCTAGCGCAGGTAGTCCGTTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCCGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCCACATTACGCTAACGGCACTAAAACCATCACCTTTTTCTTTGACTGGCGCTACAATCTTcc                                                            >  1:64531/2‑251 (MQ=255)
                                                                                                                                                                                    attGGGCAACGCGCCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACGTTTATAAATAAAGGGTGAGCGGGGCGGTTGTTAACGATGGGGTCATGCGGATTTTTCAt                                                                                                                                                                                      >  1:524/2‑105 (MQ=255)
                                                                                                                                                                                        ggCAACGCGCCTCTTCTGGCCTGCGCTAGCGCAGGCAGTACGTTTATAAATAAAGGGTGAGCGGGGCGGTTGTCACCGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCCGCTAATCAATGCTTCACTGCGCTACGGGTTTACAATCAACATTATGCTAACGGCACTAACACCATCACATTTTTCTGTGACTGGCGCTACACTCTTCCAAAGTCACAATTCTCAAAATCagaag                                  >  1:212311/1‑249 (MQ=255)
                                                                                                                                                                                                   tcttctGGCCTGCGCTAGCGCAGGTAGTAAGTTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAAAGATAGGGGCATGCGGCTTTTTTATCCACTCCTTGCGGTATGTAGTTCAGCTAATAAATGTTTACCTGTGCTAAGGGTTTGCATTCAACATTTCGCTAACGGACCTCAAACCATCATATTTTTCTGTGACTGGCGCTATAATCTTCCAAAGTCAACATTTTCAAAATCATAAGAGTATTTCTAATg                     >  1:319207/1‑251 (MQ=255)
                                                                                                                                                                                                                      gcAGGTAGTACGTTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGCAGTTCAGCTAATAAATGCTTCAc                                                                                                                                             >  1:180040/1‑112 (MQ=255)
                                                                                                                                                                                                                       cAGGTAGTACGTTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGAAGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTCGTTAAGCTAATAAATGATTAACTGCGCTAAGGGATTACACTCAACATTACGCTAACGGCACTAAAACCCTCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACGc  >  1:151485/1‑250 (MQ=255)
                                                                                                                                                                                                                          gtagtnCATTTNTNNCTCACGCGTGCGCGGGCCGNGTGTCAACGATGGGGNCATGNGGAttnttn                                                                                                                                                                                        >  2:130092/1‑64 (MQ=255)
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                               
GCTACTATTTTCACCCAGAAGTGCTCCACCACTTGCGAAACGCCCGACTGCGAACGCTTCTGGTGACAACCCAGGGGATTCAGCCCCTGTAGCCGATGATGAACGTGGCCAGCCGTTCAATCACCTCGGCGATGCACCCCCTCAGGTGTTATCACAGGACTGGCTCCTCCAACACCGTTACTTGGGCAACGCGCCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACGC  >  NC_000913/4233315‑4233779

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