Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 781,668 C→T 7.6% 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,6680CT 7.6% 117.3 / 3.7 53intergenic (+16/‑417)lysQ/nadAtRNA‑Lys/quinolinate synthase
Reads supporting (aligned to +/- strand):  ref base C (27/22);  new base T (2/2);  total (29/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.25e-01

ATCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACC  >  NC_000913/781526‑781803
                                                                                                                                              |                                                                                                                                       
aTCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGccc                                                                                                                                         >  2:414283/1‑143 (MQ=255)
 tCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCt                                                                                                                                        >  1:869934/1‑143 (MQ=255)
     gCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAg                                                                                                                                    <  1:299737/143‑1 (MQ=255)
         gACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGgcgc                                                                                                                                >  1:382267/1‑143 (MQ=255)
          aCCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTACAACAGCGCCCTCAACGCGCt                                                                                                                               >  2:886169/1‑143 (MQ=255)
           cccACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATGCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCtt                                                                                                                              >  2:544455/1‑143 (MQ=255)
            ccaccaTCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAACAAGCGCCCTCAAGGCGCttt                                                                                                                             >  1:501644/1‑143 (MQ=255)
            ccaccaTCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCttt                                                                                                                             <  1:753381/143‑1 (MQ=255)
            ccaacaTCATAATTGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCttt                                                                                                                             <  1:290542/143‑1 (MQ=255)
                  tCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCt                                                                                                                       >  1:609079/1‑143 (MQ=255)
                       cgcaaTTATGGCCGCCTAACCCCGGAAGGGGTCGGTAGCTCAGTTGGTAGCGCAGTTGACTTTCAATCAATCGGTCGCAGGTTCGCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAGAGGCGCTTTTTTGCTATCTg                                                                                                                  <  2:461732/139‑1 (MQ=14)
                          gATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGa                                                                                                               <  1:426739/143‑1 (MQ=255)
                             tAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATAc                                                                                                            >  2:230668/1‑143 (MQ=255)
                                  cagcaTAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCaaa                                                                                                       >  1:416458/1‑143 (MQ=255)
                                  cagcaTAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCaaa                                                                                                       <  2:336972/143‑1 (MQ=255)
                                       tatcaccGCAAGGGGCCGCTAGCTCAGTTGGTAGAGCAGTAGACTTTTAATCAATGGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTc                                                                                                  <  2:458234/138‑1 (MQ=11)
                                       tAATCCCGCGCGGGGCCGTGAGCTTAGTTTGTAGCGCAGTTGACCCTTTCTCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTc                                                                                                  <  2:628968/143‑1 (MQ=11)
                                       tAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTc                                                                                                  <  1:414283/143‑1 (MQ=255)
                                        aaTCCCGCGATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCg                                                                                                 <  2:157977/143‑1 (MQ=255)
                                           cccGCAAGGGGCCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAACGGCGCTTTTTTGCTATCTGCGATACTCACAGCTTCGGAc                                                                                              >  2:432794/1‑143 (MQ=14)
                                              gCAGGGGGTCGTGAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGGACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTg                                                                                           <  2:95037/143‑1 (MQ=255)
                                                  ggggTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCCAAGATTCGCACCTgaagc                                                                                       >  2:810986/1‑143 (MQ=12)
                                                  ggggTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGAAACCCAAAGATGCGAACGTgaagc                                                                                       >  2:114520/1‑143 (MQ=11)
                                                      tCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATATAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGt                                                                                   <  1:556788/143‑1 (MQ=17)
                                                       cGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGtt                                                                                  >  1:155682/1‑143 (MQ=32)
                                                          tAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttga                                                                               <  1:801227/143‑1 (MQ=32)
                                                                          gCAGTGGACTTTTAATCCATTGGCCGCAGGTTCGAATCCTGCACGACNCACCAATGTAAAAAAGCGCCCTTAAGGCGCTGTTTTGCGATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGaaa                                                               <  2:85154/143‑1 (MQ=11)
                                                                             gTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAaacaa                                                            <  1:49913/143‑1 (MQ=32)
                                                                               tGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCCGCAGGTTTGAGTTGAGCGCAGCGAAACAACg                                                          >  1:83273/1‑143 (MQ=17)
                                                                                     ttAATCATTTGGCGGCAGGTCGAAATCCGGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCg                                                    <  1:751526/143‑1 (MQ=11)
                                                                                             tgggccccAGGTTCGAACCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCgcggc                                            <  1:728433/136‑1 (MQ=12)
                                                                                                gTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTCTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGAc                                         >  1:870452/1‑143 (MQ=17)
                                                                                                     aGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGccc                                    >  1:139610/1‑143 (MQ=32)
                                                                                                        ttCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGaa                                 <  1:523757/143‑1 (MQ=32)
                                                                                                           gAATCCTTCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAggg                              <  2:410635/143‑1 (MQ=17)
                                                                                                           gAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTTAAAGCTCCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACACCGGAGCCGCTCGCGGCGACGGCCCGAAcgg                              >  1:81809/1‑143 (MQ=11)
                                                                                                            caTCCGGCCCGACCCCGCGATGTAACAAAGCGCCCTAACGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTGGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgc                             <  2:603371/142‑1 (MQ=11)
                                                                                                            aaTCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGCAACTCCAAGATTCGAACCTGCAGCAGGTTTGAGCTGACCACAGCGAAACAACGGCGACGCTAGCGGCGCCGGCCCCAACGgc                             >  2:162019/1‑143 (MQ=11)
                                                                                                            aaTCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACCCCAAGATCCGAACCTGCCGCCGGTTTGAGTGAAGCACAGAGAAGCCACGGAGCCCCGCGCGGCGAGGGCCCGAGGGgc                             >  2:159426/1‑143 (MQ=11)
                                                                                                            aaTCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACCCAAAGATACGATCCTGCAGCACACTTGAGTTGCGCCGAGCACTACAACAGAGCCACTACGGGCGACCGCCCGAAcgga                             >  2:343928/1‑142 (MQ=2)
                                                                                                                   cACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGagcgaa                      <  1:58500/143‑1 (MQ=32)
                                                                                                                      gACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGCCGACGGCGCGAAGGccgagcgactcg                   >  1:639505/1‑139 (MQ=255)
                                                                                                                       aCCCACCAATGTACAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGagcgaagcga                  <  1:298690/143‑1 (MQ=17)
                                                                                                                         ccaccaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACTTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACCACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGt                >  1:325284/1‑143 (MQ=17)
                                                                                                                         ccaccaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGGAGGGCGAGCGAAGCGAGt                >  1:821500/1‑143 (MQ=17)
                                                                                                                                 gTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGACGGGCGAGCGAAGCGGGTCATCCTGc        >  1:173052/1‑143 (MQ=17)
                                                                                                                                  tAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCATGCa       >  1:627405/1‑143 (MQ=17)
                                                                                                                                   aaaaaaGCGCCTTTAGGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCt                                                  >  2:292683/1‑99 (MQ=11)
                                                                                                                                   aaaaaaGCGCCTTTAGGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCt                                                  >  2:755371/1‑99 (MQ=11)
                                                                                                                                   aaaaaaGCGCCTTTAGGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCt                                                  <  1:292683/99‑1 (MQ=11)
                                                                                                                                   aaaaaaGCGCCTTTAGGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCt                                                  <  1:755371/99‑1 (MQ=11)
                                                                                                                                    aaaaaGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGACGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACg     >  2:840303/1‑143 (MQ=17)
                                                                                                                                     aaaaGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGACTCGCACCTGCAGCAGGGCTGAGTTGAGCGCAGCGAAACCACGGAGCCGCTCGCGGCGACGGCACGAAGGGAGCACGAAGCAAGTTGACCTGCACGa    >  2:35424/1‑143 (MQ=11)
                                                                                                                                       aaGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGAcc  >  1:746302/1‑143 (MQ=32)
                                                                                                                                              |                                                                                                                                       
ATCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACC  >  NC_000913/781526‑781803

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: