Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 781,672 A→G 15.0% intergenic (+20/‑413) lysQ → / → nadA tRNA‑Lys/quinolinate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913781,6720AG15.0% 87.3 / 10.3 53intergenic (+20/‑413)lysQ/nadAtRNA‑Lys/quinolinate synthase
Reads supporting (aligned to +/- strand):  ref base A (28/15);  new base G (4/5);  total (33/20)
Fisher's exact test for biased strand distribution p-value = 2.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.39e-01

GACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAAT  >  NC_000913/781535‑781810
                                                                                                                                         |                                                                                                                                          
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)
                                                                                                                                         |                                                                                                                                          
GACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAAT  >  NC_000913/781535‑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: