Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 781,668 C→T 12.0% intergenic (+16/‑417) lysQ → / → nadA tRNA‑Lys/quinolinate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913781,6680CT12.0% 89.2 / 8.9 50intergenic (+16/‑417)lysQ/nadAtRNA‑Lys/quinolinate synthase
Reads supporting (aligned to +/- strand):  ref base C (28/16);  new base T (3/3);  total (31/19)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.58e-02

CTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAAT  >  NC_000913/781529‑781810
                                                                                                                                           |                                                                                                                                              
cTGCACGACCCAGCAGCCTGGAGGATCACGGCAGCATAACCCCGCTAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTGTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTaa                                                                                                                                             <  2:534979/143‑1 (MQ=17)
      gACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGCGGTCGTTAGCTCAGTTGGTAGAGCAGTCGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGgcgc                                                                                                                                       >  2:478471/1‑143 (MQ=255)
        cccACCATCCTGAATGATTAAGGCAGCAAAAACCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCTTTAGGGCGCtt                                                                                                                                     <  1:397170/143‑1 (MQ=18)
          caccaTCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCcttt                                                                                                                                   >  1:616081/1‑143 (MQ=255)
             caTCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTg                                                                                                                                >  2:854893/1‑143 (MQ=255)
              aTCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGCAAAAAAGCGCCTTTAGGGCGCttttttac                                                                                                                               <  1:634356/143‑3 (MQ=17)
                ccTTAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGCAAAAAAGCGCCTTTAGGGCGCttttttacat                                                                                                                             <  1:1015003/143‑5 (MQ=17)
                    aaTGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTACTCAATTGGTCGCAGGTTCGAATCCTGCCCGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTCACTCTCTg                                                                                                                         >  2:720190/1‑143 (MQ=17)
                    aaTGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGGCCCACCAATGTAAACAAGCGCCCTAAAGGCGCTATTTTGCTATCTg                                                                                                                         >  2:825944/1‑143 (MQ=255)
                      tGATTAAGGCAGCAGAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCg                                                                                                                       <  1:996454/143‑1 (MQ=255)
                             ggCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACCCa                                                                                                                >  2:319706/1‑143 (MQ=255)
                              gcggcaTAACCGCGGAAGCGGTCGTTCGCTCAGTTGGTAGAGCAGTCGACTTTTCATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCaa                                                                                                               <  2:638056/140‑1 (MQ=11)
                                    tAACCCCGTAATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTc                                                                                                         <  2:452731/143‑1 (MQ=255)
                                      aTCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGa                                                                                                       <  1:179938/143‑1 (MQ=40)
                                       tCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCGTGCACGACCCACCAATGTAAAAAAGCGCCCTAAGGGCGCTTTTTTGCTATCTGCGGTACCCAACGATTCGaa                                                                                                      >  1:506871/1‑143 (MQ=14)
                                          cGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTGAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCt                                                                                                   >  2:764062/1‑143 (MQ=37)
                                           gCACGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTg                                                                                                  <  2:751320/143‑1 (MQ=21)
                                              gggggTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTgcag                                                                                               <  1:20884/142‑1 (MQ=18)
                                              atgggTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTgcag                                                                                               <  2:973269/141‑1 (MQ=17)
                                               ggggTCGTTAGCTCAATTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGACTCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTGTTTGCTATCTGCGATACTCAAAGCTACGAACCTgaagc                                                                                              >  2:902700/1‑143 (MQ=11)
                                                     ggtggCTCAGTTGGTAGAGGAGTTGACTTTGAATCAGGGGGGCGCAGGTTCGAATCGTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCTGCAGGttc                                                                                        <  2:67817/139‑2 (MQ=11)
                                                        agtgcggTTGGTAGATCTGTTGACTTGTAATCGATTGTTCGCAGGTTCGAATCCTGCACGGCGCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgag                                                                                     <  2:809487/137‑1 (MQ=11)
                                                              gTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAGAAAGCGCCCTAAAGGCGCTTCTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAgc                                                                               >  2:949960/1‑143 (MQ=17)
                                                                       gCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGaaa                                                                      <  1:894304/143‑1 (MQ=32)
                                                                        cAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAaac                                                                     >  1:865767/1‑143 (MQ=32)
                                                                          gTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAaacaa                                                                   >  2:326649/1‑143 (MQ=32)
                                                                                ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGc                                                             >  2:112891/1‑143 (MQ=32)
                                                                                 tttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCCCCAATGTAACAAAGCGCCCTACCCGCGCTTTTTTGCTATCTGCGATCCTCAAAGATACGAACATGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGAAGcc                                                            >  1:225393/1‑143 (MQ=11)
                                                                                   tattCAATTGTTCGCAGGTACGAATCGTAGACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTGTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGc                                                          <  1:117259/140‑1 (MQ=11)
                                                                                     aTCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTc                                                        >  1:464577/1‑143 (MQ=32)
                                                                                         aTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCgcgg                                                    >  1:548323/1‑143 (MQ=32)
                                                                                          ttGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCgcggc                                                   >  1:694776/1‑143 (MQ=32)
                                                                                                        gAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTACGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCACGCGGCGACGGCCCGAAagg                                     >  2:322106/1‑143 (MQ=12)
                                                                                                         aagccTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTGTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgc                                    <  1:1013458/140‑1 (MQ=17)
                                                                                                                cACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTGGAGTTGAACGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGagcgaa                             >  2:525764/1‑143 (MQ=17)
                                                                                                                   gACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACAGAGCCGCTAACGACGACGGCCCGAAGGGCGagagaagcg                          >  1:462236/1‑143 (MQ=11)
                                                                                                                    aCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTCGAATTGAACGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGagcgaagcga                         <  1:902585/143‑1 (MQ=12)
                                                                                                                       caccaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGGAACAACGGAGCCGCTCGCGGCGACGGCCCGGAGGGCGAGCGACGCGAGTc                      >  2:542127/1‑143 (MQ=12)
                                                                                                                        accaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGCGTCa                     >  1:385679/1‑143 (MQ=17)
                                                                                                                        accaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCa                     >  2:252245/1‑143 (MQ=32)
                                                                                                                             tGTAAAAAAGCACCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTg                >  2:908406/1‑143 (MQ=17)
                                                                                                                              gTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCCGCGAAACAACGGCGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGAGACGCATCCTGc               >  1:301760/1‑143 (MQ=11)
                                                                                                                                aaaaaaGCGCCTTTAGGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCCTCCTGCAc             >  2:859799/1‑143 (MQ=11)
                                                                                                                                aaaaaaGCGCCTTTAGGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCAc             >  2:338047/1‑143 (MQ=255)
                                                                                                                                aaaaaaGCGCCTTTAGGGCGCTTTTTTGCTATCTGCGAAACGCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCAc             >  2:647309/1‑143 (MQ=11)
                                                                                                                                 aaaaaGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACg            >  1:143950/1‑143 (MQ=32)
                                                                                                                                    aaGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGGCATCCTGCACGAcc         >  2:908716/1‑143 (MQ=17)
                                                                                                                                      gcgcCCTAAGGGCGCTTTTTTGCTATCTGCGATACTCCAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGACGCAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACcca       <  2:357066/143‑1 (MQ=255)
                                                                                                                                      gcgcCCTAAGGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACcca       <  2:244911/143‑1 (MQ=17)
                                                                                                                                        gcCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTGCAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGACGGGCGAGCGCAGCGAGTCATCCGGCACGAcccccc     >  2:681155/1‑143 (MQ=11)
                                                                                                                                         cccTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccacca    >  2:351949/1‑143 (MQ=32)
                                                                                                                                          ccTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCaa   >  1:703898/1‑143 (MQ=32)
                                                                                                                                          ccTAAAGGCGCTTTTGTGCTATCTGCGATTCTCTATGATACGAACCTGGAGCCGGTTGGAGTCGAGCGCAGCGAAACAACGGAGCCCCTCGCGGCCACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCaa   <  2:18391/143‑1 (MQ=11)
                                                                                                                                           cTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAAt  <  1:650336/143‑1 (MQ=32)
                                                                                                                                           |                                                                                                                                              
CTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAAT  >  NC_000913/781529‑781810

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: