Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,235,755 C→T S496N (AGC→AAC)  ycgB ← PF04293 family protein YcgB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,235,7550CT97.1% 97.4 / ‑4.4 34S496N (AGC→AAC) ycgBPF04293 family protein YcgB
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (16/17);  total (16/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.66e-01

ACCACAAACTGTCGGCGTTGTGCAGTGAAGGCGCGCACGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGT‑TTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGTTCCAGATCTGAATATTCGGCTCCAGA‑TTACTTAAGTTATATTGCGACGATAACCGGTTACGGATCTCCCGATAACCTTCTTCATTATGAATAGC  >  NC_000913/1235515‑1235999
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                     
aCCACAAACTGTCGGCGTTGTGCAGTGAAGGCGCGCACGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGT‑TTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGTTGCCGTCtt                                                                                                                                                                                                                                              <  1:501882/250‑1 (MQ=255)
                                    aCGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGT‑TTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCCACTCGATGTTTCCGTCTTCGTTTTGCTTTTCGAGCATCCCATCCAATccccccc                                                                                                                                                                                                          >  2:249993/1‑247 (MQ=255)
                                                        gTGCGTCGCTATGGGCCTGAGAGTGGCGAGATTTGTCACCGGATGCAGAAAAATCCGCCGGGTGATTTAGCCACTCAGGGTCGATAATCTCTTGCTTTTAAAGAGCAAACCCCTCAAACGAGGTGTTTTTTGTTGT‑TTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGTTGCCGTCTTCGTTTTGCTGTTCGAGCAGCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGa                                                                                                                                                                                     <  1:416076/251‑1 (MQ=255)
                                                                   tGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGT‑TTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGTTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTgcgc                                                                                                                                                                          >  1:374651/1‑251 (MQ=255)
                                                                                        tttgCACCGGATGCAGAAAAATCTGCCGTGTGTGTCAGTCATTCAGGGGCGATAATCCATTGCGTTTAAAGTGCAAACTCCTCAAACGAGGGTTTTTTTGTTGT‑TTTTACATATTGCCCATTTTTGGCGGGCAACGTTCCAGCAACTCGATGTTGCCGTCTTCTTTTTGCTGTTCGAGCATTACATCAAATCCCCACAGGCGATGCACATGCTTAAGTACTTCTTTGCGCCCCAGATCCAGCGGTgcgca                                                                                                                                                      <  2:419362/247‑2 (MQ=255)
                                                                                                     cAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGT‑TTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGTTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAatat                                                                                                                                        <  1:522583/251‑1 (MQ=255)
                                                                                                       aaaaaatttgtccgaggagttAGCCATTGAGGGGCGATATTCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTGTTGT‑GTGTTTTACAGATTTCCCATTCTTGGTGGGCATCTTTCCAGCAACTCGATGTTGCCTGCGACGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCATATGCTTCAGGACTTCTTTGCGCCCACGATCCAGCTGTGCGCGATTATGTTGACTATaa                                                                                                                                      <  2:341015/232‑1 (MQ=255)
                                                                                                                             gCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGT‑TTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGTTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGActtctt                                                                                                                                                                               >  1:236593/1‑188 (MQ=255)
                                                                                                                             gCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGT‑TTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGTTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGActtctt                                                                                                                                                                               <  2:236593/188‑1 (MQ=255)
                                                                                                                               cATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGT‑TTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGTTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACATGCGATGCACATGCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGAACGGTcgccg                                                                                                              <  1:145530/251‑1 (MQ=255)
                                                                                                                                           aTAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGATGGGTTTTTTGTTGT‑TTGTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGTTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTGCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACg                                                                                                  <  2:359931/251‑1 (MQ=255)
                                                                                                                                           aTAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGT‑TTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGTTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACg                                                                                                  <  1:34842/251‑1 (MQ=255)
                                                                                                                                            tAATCACTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGT‑TTTTACAGATTTCCAATTCTTGGCGGGCAACGTTCCAGCAACTCGATTTTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGCTGCACATTCTTCAGGCCTTCTTTGCGCCCCCTATCCAACGGGTCGCGTTTATGTGGAATATAACACAGCGTCAGCGAACAGACGCCCCGCAAATCCaatt                                                                                                 >  2:414925/1‑248 (MQ=255)
                                                                                                                                                    tCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGT‑TTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGTTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGTTcc                                                                                              >  1:542324/1‑246 (MQ=255)
                                                                                                                                                               agcaaacacattaaaCGAGGGGTATTTTGTTGT‑TTTTCCAGATTTCCCATTCTTGGCGGGCAAAGTTCCAGCAACTCGTTGTTGCCGTCTTCGTTTTGCTGTTCGCGCATCACATCAAACACACACAGGCGATGCACCTGCTTCCGGTCTTCTTTGCGCACCCGTACCAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTCCGCGAACGGCCGCCGCGCAAATCCACGATCCAGCTCTgcttcttcgg                                                                              >  2:81361/13‑242 (MQ=255)
                                                                                                                                                                   aaCCCCTCAAACGAGGGGTTTTTTGTTGT‑TTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCCGCAACTCGATGTTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGTTCCAGATCTGAATATTCGGCTcc                                                                          <  2:374651/251‑1 (MQ=255)
                                                                                                                                                                   aaCCCCTCAAACGAGGGGTTTTTTGTTGT‑TTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGTTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGTTCCAGATCTGAATATTCGGCTcc                                                                          >  1:449248/1‑251 (MQ=255)
                                                                                                                                                                          cAAACGATGGGTTTGTGGTTTT‑TTTTACAGAGTTCCCATTCTTGGCGGGAAACGTTCCAGCAACTCGATGTTGCCGTCTTCGTTTTGCGGTTCGAGCATCTCATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGCTTCAGCGAACGGTCGCCGCGCAAATCCACGTTCCAGATCTTAATATTCGGCTCCAGA‑TTAc                                                                  <  2:427601/251‑1 (MQ=255)
                                                                                                                                                                                   ggTTTTTTGTTGT‑TTTTACAGATTTCCCATTCTTGGCGGGCAAAGTTCCAGCAACTCGATGTTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGTTCCAGATCTGAATATTCGGCTCCAGA‑TTACTTAAGTtat                                                         >  1:376265/1‑251 (MQ=255)
                                                                                                                                                                                            ttgt‑tTTTACAGATTGCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGTTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCTCAGGCGATGCACATGCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGTTCCAGATCTGAATATTCGGCTCCAGA‑TTACTTAAGTTATATTGcgac                                                 <  2:109929/250‑1 (MQ=255)
                                                                                                                                                                                              gt‑tTTTACAAATGTCCCTGTCTTGGCGGGCAACCGTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTGTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGTGCCAGATCTGAATATTCGGCTCCAGA‑TTACTTAAGTTATATTGcgacga                                               <  2:186412/250‑1 (MQ=255)
                                                                                                                                                                                                     aCAGATTTCCCATTCTTGGCGTGCAACTTTCCAGCAATTCGATGTAGCCGTCTTCGTTTTGCTGATCGAGCATCACATCAAAGCCCCCCAGGCGATGCACATGCTTCAGGACTTCGTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTTAGCGAACGGTCGCCGCGCAAATCCACGTTCCAGATCTGAATATTCGGCTACAGA‑TTACTTAAGTTATATTGCGACGATAACCgg                                        <  2:46125/251‑1 (MQ=255)
                                                                                                                                                                                                             cccATTCTTGGCGGGCAACGGTCCAGCATCTCGCTGTTGCCGTCCGCGTTTTGCTGTGCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTGCGCTCCCGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGTTCCAGATCTTAATATGCGGCTCCAGA‑TTACTTAAGTTATATTGCGACGATAACCGGTTACGGAt                                <  2:273784/251‑1 (MQ=255)
                                                                                                                                                                                                              ccagtCTTGGCTGGCAACGTACCAGCACCTCGATGTTGCCGTCTTCGTTTTGCTGTTCGAGCATCACAGCAAACCCCCACAGGCGATGCACATGCTGCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGTTCCAGATCTGAATATTCGGCTCCAGA‑TTACTTAAGTTATATTGCGACGATAACCGGTTACGGAtc                               <  2:54733/247‑1 (MQ=255)
                                                                                                                                                                                                              ccATTCTTGGCGGGCAACGTTCCAGCAACTCGATGTTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACCTGATTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGCTTATTTGGAATATAAACCAGCGTCAGCGAACGGTACCCACGCAAATCCAAGTTCCAGATCTGAATATTCGGCTC‑AGACTTACTTAAGTTATAATGCGAAGATAACCGGATACgagtc                               >  2:115309/1‑247 (MQ=255)
                                                                                                                                                                                                               cATTCTTGGCGGGCAACGTTCCAGCAACTCGATGTTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGTTCCAGATCTGAATATTCGGCTCCAGA‑TTACTTAAGTTATATTGCGACGATAACCGGTTACGGAtct                              >  1:458383/1‑251 (MQ=255)
                                                                                                                                                                                                               cATTCTTGGCGGGCAACGTTCCAGCAACTCGATGTTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGTTCCAGATCTGAATATTCGGCTCCAGA‑TTACTTAAGTTATATTGCGACGATAACCGGTTACGGAtct                              >  1:84983/1‑251 (MQ=255)
                                                                                                                                                                                                                     tGGCGGGCAACGTTCCAGCAACTCGATGTTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGTTCCAGATCTGAATATTCGGCTCCAGA‑TTACTTAAGTTATATTGCGACGATAACCGGTTACGGATCTccc                           <  1:518502/248‑1 (MQ=255)
                                                                                                                                                                                                                          ggCAACGTTCCAGCAACTCGATGTTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGTTCCAGATCTGAATATTCGGCTCCAGA‑TTACTTAAGTTATATTGCGACGATAACCGGTTACGGATCTCCCGATAACct                   >  2:138555/1‑251 (MQ=255)
                                                                                                                                                                                                                            cAACGTTCCAGCAACTCGATGTTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGTTCCAGATCTGAATATTCGGCTCCAGA‑TTACTTAAGTTATATTGCGACGATAACCGGTTACGGATCTCCCGATAACcttc                 >  1:72167/1‑251 (MQ=255)
                                                                                                                                                                                                                                 ttCCAGCAACTCGATGTTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGTTCCAGATCTGAATATTCGGCTCCAGA‑TTACTTAAGTTATATTGCGACGATAACCGGTTACGGATCTCCCGATAACCTTCTTCAt            >  2:484561/1‑251 (MQ=255)
                                                                                                                                                                                                                                  taaagaaaCTGGATGTTGTCGTATTCGTTTTGCTGTTGGAGCATCACATCAAATCCCAACAGGCGATGCACATGATTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGAGATTATGTGGAATATAACGCAGCGTCAGAGAACGGTCGCCGCGCAAATCCACGTTCCAGATCTGAATATTCGGCTCAAGA‑TTACTTAAGTTATATTGCGACGATAACCGGTTACGGATCTCCCGATAACcttattaatt           >  2:251957/7‑244 (MQ=255)
                                                                                                                                                                                                                                     aGCAACTCGATGTTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGTTCCAGATCTGAATATTCGGCTCCAGA‑TTACTTAAGTTATATTGCGACGATAACCGGTTACGGATCTCCCGATAACCTTCTTCATTATg        >  1:374039/1‑251 (MQ=255)
                                                                                                                                                                                                                                            cGATGTTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGTTCCAGATCTGAATATTCGGCTCCAGA‑TTACTTAAGTTATATTGCGACGATAACCGGTTACGGATCTCCCGATAACCTTCTTCATTATGAATAGc  <  1:432495/250‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                     
ACCACAAACTGTCGGCGTTGTGCAGTGAAGGCGCGCACGATCAGGTGTACGAAACAGTGCGTCGCTATGGGCATGAGAGTGGCGAGATTTGGCACCGGATGCAGAAAAATCTGCCGGGTGATTTAGCCATTCAGGGGCGATAATCCCTTCCGTTTAAAGAGCAAACCCCTCAAACGAGGGGTTTTTTGTTGT‑TTTTACAGATTTCCCATTCTTGGCGGGCAACGTTCCAGCAACTCGATGCTGCCGTCTTCGTTTTGCTGTTCGAGCATCACATCAAATCCCCACAGGCGATGCACATGCTTCAGGACTTCTTTGCGCCCCCGATCCAGCGGTGCGCGATTATGTGGAATATAACGCAGCGTCAGCGAACGGTCGCCGCGCAAATCCACGTTCCAGATCTGAATATTCGGCTCCAGA‑TTACTTAAGTTATATTGCGACGATAACCGGTTACGGATCTCCCGATAACCTTCTTCATTATGAATAGC  >  NC_000913/1235515‑1235999

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