Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,060,843 (G)6→7 coding (397/711 nt) yihL → putative transcriptional regulator YihL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,060,8371.G96.9% 103.9 / ‑3.1 32coding (391/711 nt)yihLputative transcriptional regulator YihL
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (17/14);  total (18/14)
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.68e-01

GCAAGGGCGAGAGCCTTCCTGGGGGTTTACCGAGAAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTACCAAA  >  NC_000913/4060719‑4060963
                                                                                                                       |                                                                                                                               
gCAAGGGCGAGAGCCTTCCTGGGGGTTTACCGAGAAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGGCGCTGGa                                                                                                                   <  2:353801/133‑1 (MQ=255)
gCAAGGGCGAGAGCCTTCCTGGGGGTTTACCGAGAAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGGCGCTGGa                                                                                                                   <  1:394199/133‑1 (MQ=255)
    gggCGAGAGCCTTCCTGGGGGTTTACCGAGAAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGGCGCTGGAAGGa                                                                                                               >  2:483186/1‑133 (MQ=255)
      gCGAGAGCCTTCCTGGGGGTTTACCGAGAAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGTATCACCGCCTGCGGGGGCGCTTGAAGGACa                                                                                                             >  1:161998/1‑133 (MQ=255)
        gagaGCCTTCCTGGGGGTTTACCGAGAAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGGCGCTGGAAGGACATa                                                                                                           >  1:55367/1‑133 (MQ=255)
        gagaGCCTTCCTGGGGGTTTACCGAGAAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGGCGCTGGAAGGACATa                                                                                                           <  1:483186/133‑1 (MQ=255)
                      gggTTTACCGAGAAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCAc                                                                                             <  2:127605/133‑1 (MQ=255)
                       ggTTTACCGAGAAAACCCGTACCACCGAAATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACg                                                                                            <  2:181541/133‑1 (MQ=255)
                       ggTTTACCGAGAAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACg                                                                                            <  1:293985/133‑1 (MQ=255)
                       ggTTTACCGAGAAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACg                                                                                            <  1:400432/133‑1 (MQ=255)
                       ggTTTACCGAGAAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACg                                                                                            <  1:278403/133‑1 (MQ=255)
                        gTTTACCGAGAAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGa                                                                                           <  1:26421/133‑1 (MQ=255)
                         tttACCGAGAAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGaa                                                                                          <  2:69159/133‑1 (MQ=255)
                          ttACCGAGAAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGaaa                                                                                         >  1:495367/1‑133 (MQ=255)
                           tACCGAGAAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAAc                                                                                        >  1:420105/1‑133 (MQ=255)
                            aCCGAGAAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACa                                                                                       <  1:448168/133‑1 (MQ=255)
                                                         cGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCtt                                                          <  2:55367/133‑1 (MQ=255)
                                                           agaCGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCt                                                                  >  1:42532/1‑123 (MQ=255)
                                                           agaCGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCt                                                                  <  2:42532/123‑1 (MQ=255)
                                                             aCGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATAAATTAATCATGACGTAGCTCCGTGCTTTAtt                                                      >  1:238990/1‑133 (MQ=255)
                                                                 tCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGaac                                                  <  1:344650/133‑1 (MQ=255)
                                                                  cgcgCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGGCGCTGGGAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGaaca                                                 >  2:123765/1‑133 (MQ=255)
                                                                                 aGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTc                                  >  2:8262/1‑133 (MQ=255)
                                                                                     aCGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCAttt                              >  2:285220/1‑133 (MQ=255)
                                                                                      cGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCAtttt                             >  2:432529/1‑133 (MQ=255)
                                                                                        ccATCAACTGAACTCTATCGCATCACCGGCTGGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGACTTTATTGAACAACTTGAAAACCACTCATTTTCt                           >  2:116813/1‑133 (MQ=255)
                                                                                                  aaCTCTATCGCATCACCGGCTGGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTggg                 >  2:414130/1‑133 (MQ=255)
                                                                                                         tCGCATCACCGGCTGGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTg          >  2:33860/1‑133 (MQ=255)
                                                                                                         tCGCATCACCGGCTGGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTg          >  2:74526/1‑133 (MQ=255)
                                                                                                            cATCACCGGCTGGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTa       >  1:384248/1‑133 (MQ=255)
                                                                                                                 ccGGCTGGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTACCaaa  >  1:93201/1‑133 (MQ=255)
                                                                                                                 ccGGCTGGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTACCaaa  >  1:128494/1‑133 (MQ=255)
                                                                                                                       |                                                                                                                               
GCAAGGGCGAGAGCCTTCCTGGGGGTTTACCGAGAAAAACCGTACCAGCGATATTCCCGAGACGCTCGCGCCACTGATTGCAGTGACGCCATCAACTGAACTCTATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTACCAAA  >  NC_000913/4060719‑4060963

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