Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,792,768 A→G 18.0% intergenic (‑217/+58) kbl ← / ← yibB glycine C‑acetyltransferase/YibB family protein, function unknown

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,792,7680AG18.0% 119.3 / 17.6 50intergenic (‑217/+58)kbl/yibBglycine C‑acetyltransferase/YibB family protein, function unknown
Reads supporting (aligned to +/- strand):  ref base A (21/20);  new base G (4/5);  total (25/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

TTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATC  >  NC_000913/3792629‑3792903
                                                                                                                                           |                                                                                                                                       
ttCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGGAAAAATGATGCTTTGGCATCAtt                                                                                                                                 <  3:546381/148‑1 (MQ=255)
                      aaCGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAg                                                                                                                    <  3:1071683/139‑1 (MQ=255)
                      aaCGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAg                                                                                                                    >  4:1071683/1‑139 (MQ=255)
                                 tGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCTTTGGCATCAtt                                                                                                                                 <  2:125992/115‑1 (MQ=255)
                                 tGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCTTTGGCATCAtt                                                                                                                                 >  1:125992/1‑115 (MQ=255)
                                      tctctcAACGCCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAGTGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGca                                                                                                                                      >  4:739232/1‑105 (MQ=255)
                                      tctctcAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCTTTGGCATCAtt                                                                                                                                 <  3:924451/110‑1 (MQ=255)
                                      tctctcAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGt                                                                                                                   >  2:1067939/1‑124 (MQ=255)
                                      tctctcAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGt                                                                                                                   <  1:1067939/124‑1 (MQ=255)
                                      tctctcAACACCAACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGca                                                                                                                                      <  3:739232/105‑1 (MQ=255)
                                         ctcAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTAt                                                                                                             <  2:93771/127‑1 (MQ=255)
                                         ctcAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTAt                                                                                                             >  1:93771/1‑127 (MQ=255)
                                             acacCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCTTTGGCATCAttt                                                                                                                                <  1:5389/104‑1 (MQ=255)
                                             acacCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCTTTGGCATCAttt                                                                                                                                >  2:5389/1‑104 (MQ=255)
                                             acacCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGcatc                                                                                                                                    <  1:211799/100‑1 (MQ=255)
                                             acacCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGcatc                                                                                                                                    >  2:211799/1‑100 (MQ=255)
                                                       aTAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGt                                                                                                                   <  3:982843/107‑1 (MQ=255)
                                                       aTAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGt                                                                                                                   >  4:982843/1‑107 (MQ=255)
                                                              aaaaaCTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTcc                                                                                                                            <  4:644204/91‑1 (MQ=255)
                                                              aaaaaCTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTcc                                                                                                                            >  3:644204/1‑91 (MQ=255)
                                                               aaaaCTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCTTTGGCATCAtt                                                                                                                                 >  3:424388/1‑85 (MQ=39)
                                                                     gTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAg                                                                                                                                        >  3:25986/1‑72 (MQ=255)
                                                                     gTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAg                                                                                                                                        <  4:25986/72‑1 (MQ=255)
                                                                                            gatgatTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAg                                                                                      <  2:475409/99‑1 (MQ=255)
                                                                                            gatgatTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAg                                                                                      >  1:475409/1‑99 (MQ=255)
                                                                                             atgatTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGaa                                                                                  >  1:946134/1‑102 (MQ=255)
                                                                                             atgatTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGaa                                                                                  <  2:946134/102‑1 (MQ=255)
                                                                                             atgatTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTg                                       >  4:837206/1‑145 (MQ=255)
                                                                                             atgatTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTg                                       <  3:837206/145‑1 (MQ=255)
                                                                                                                     taAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCt                                                                       <  2:914553/89‑1 (MQ=255)
                                                                                                                     taAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCt                                                                       >  1:914553/1‑89 (MQ=255)
                                                                                                                       aaaaGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAgtg                            >  1:304532/1‑130 (MQ=255)
                                                                                                                       aaaaGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAgtg                            <  2:304532/130‑1 (MQ=255)
                                                                                                                       aaaaGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTCGaa                                                                   >  1:1013613/1‑91 (MQ=255)
                                                                                                                       aaaaGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTCGaa                                                                   <  2:1013613/91‑1 (MQ=255)
                                                                                                                        aaaGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGaa                                                                   <  2:940836/90‑1 (MQ=255)
                                                                                                                        aaaGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGaa                                                                   >  1:940836/1‑90 (MQ=255)
                                                                                                                         aaGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCt                                                                       >  2:165574/1‑85 (MQ=255)
                                                                                                                         aaGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCt                                                                       <  1:165574/85‑1 (MQ=255)
                                                                                                                               aaaTGATGCTTTGGCATCATTTTTCCTATCTCAGTTGGTATCATAACGCTACTCACTAATGGAGGGAAACTTATTTTCTTGaa                                                                   >  2:229890/1‑83 (MQ=37)
                                                                                                                               aaaTGATGCTTTGGCATCATTTTTCCTATCTCAGTTGGTATCATAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGaa                                                                   <  1:229890/83‑1 (MQ=38)
                                                                                                                               aaaTGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTg                                       >  3:802002/1‑111 (MQ=255)
                                                                                                                               aaaTGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTg                                       <  4:802002/111‑1 (MQ=255)
                                                                                                                               aaaTGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAGATc  >  4:219191/1‑148 (MQ=255)
                                                                                                                                aaTGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCt                                                                       >  2:431772/1‑78 (MQ=255)
                                                                                                                                aaTGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCt                                                                       <  1:431772/78‑1 (MQ=255)
                                                                                                                                    atgCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCt                                                                                                    <  1:395811/45‑1 (MQ=255)
                                                                                                                                    atgCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCt                                                                                                    >  2:395811/1‑45 (MQ=255)
                                                                                                                                      gCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTa                                 >  4:192038/1‑110 (MQ=255)
                                                                                                                                      gCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTa                                 <  3:192038/110‑1 (MQ=255)
                                                                                                                                         aaaGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCt                                                                       <  2:783164/69‑1 (MQ=255)
                                                                                                                                         aaaGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCt                                                                       >  1:783164/1‑69 (MQ=255)
                                                                                                                                          aaGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATtc                                                         >  3:754307/1‑82 (MQ=255)
                                                                                                                                          aaGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATtc                                                         <  4:754307/82‑1 (MQ=255)
                                                                                                                                          aaGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAAt   >  3:324479/1‑136 (MQ=255)
                                                                                                                                          aaGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTCTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCACTAAAGTGTTACTCCTGAAGCAGCGAAACaatt   <  4:324479/136‑3 (MQ=255)
                                                                                                                                           |                                                                                                                                       
TTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTACAGTTTTGGCTGGGCATAATGATGATTGATTTTACTGCATTGATATAAAAAGAAAAAATGATGCCAAAGCATCATTTTTCCTATCTCAGTTGGTATCAAAACGCTACCCACTAATGGAGGGAAACCTATTTTCTTGAAAGGAAAATTCTCAGCGCCTGCATTTTTGCCCCTAAAGTGTTACTCCTGAAGCAGCGAAACAAATC  >  NC_000913/3792629‑3792903

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