Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,521,348 A→G 11.9% intergenic (+20/+245) lysV → / ← xapR tRNA‑Lys/DNA‑binding transcriptional activator XapR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,521,3480AG11.9% 97.9 / 8.1 52intergenic (+20/+245)lysV/xapRtRNA‑Lys/DNA‑binding transcriptional activator XapR
Reads supporting (aligned to +/- strand):  ref base A (21/23);  new base G (2/4);  total (24/28)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.71e-01

ACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCAC  >  NC_000913/2521207‑2521482
                                                                                                                                             |                                                                                                                                      
ccccgcggggcgTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAg                                                                                                                                       <  2:652919/136‑1 (MQ=255)
     gaggGGGTGGGCGCTTGGAGCCCGTCTTCAGCCACTACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTGGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCTTTAGGGCGCt                                                                                                                                  <  2:1064906/143‑1 (MQ=14)
      aggGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCtt                                                                                                                                 <  1:467439/143‑1 (MQ=255)
        ggggTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCtttt                                                                                                                               <  2:711270/143‑1 (MQ=255)
        ggggTCCGCGGTTCGCTCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGCATCCTGCACGACCCACCAATGTAAAAAAGCGCCTTTAGGGCGCtttt                                                                                                                               <  1:808310/143‑1 (MQ=17)
          ggTCGGCGGGTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCTTTAGGGCGCtttttt                                                                                                                             <  1:895974/143‑1 (MQ=17)
               gcggTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCTTTAGGGCGCTTTTTTAcatt                                                                                                                        <  2:63233/143‑4 (MQ=17)
                  gTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTg                                                                                                                     <  2:358190/143‑1 (MQ=255)
                       aTCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGGTa                                                                                                                >  2:605550/1‑143 (MQ=255)
                       aTCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCACCCTAAAGGCGCTTTTTTACTATCTGCGATa                                                                                                                >  1:477101/1‑143 (MQ=255)
                        tCCCGTCATCACCCGCCACTGCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATAc                                                                                                               <  1:482217/143‑1 (MQ=255)
                            gTCATCACCCGCCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCaa                                                                                                           <  2:1061486/143‑1 (MQ=255)
                             tcatcaCCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCaaa                                                                                                          >  2:259058/1‑143 (MQ=255)
                                        cacTTCGGGTCGTTAGCTGAGTTGGTAGAGTAGTTGACTTTTAATCAATTGGTCGCAGGTTAGAATCCTGCACGACACACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTGCACGATTGGAACCt                                                                                               >  2:189834/1‑143 (MQ=11)
                                         acTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTg                                                                                              <  2:724664/143‑1 (MQ=37)
                                           ttCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTCATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAACGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTgca                                                                                            >  2:339704/1‑143 (MQ=21)
                                           ttCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTgca                                                                                            <  2:745699/143‑1 (MQ=35)
                                           ttCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTgca                                                                                            <  2:458892/143‑1 (MQ=35)
                                               ggTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGcagcag                                                                                        <  2:385129/143‑1 (MQ=32)
                                                gTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAgg                                                                                       >  2:722611/1‑143 (MQ=32)
                                                 tCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGt                                                                                      <  1:316705/143‑1 (MQ=32)
                                                  ggTTGGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGtt                                                                                     <  2:343628/142‑1 (MQ=17)
                                                     tAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTCACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttga                                                                                  <  2:700144/143‑1 (MQ=17)
                                                        cTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGACCCTGCAGCAGGTttcagtt                                                                               >  2:407406/1‑143 (MQ=17)
                                                         tCAGTGGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagttg                                                                              <  1:281537/143‑1 (MQ=17)
                                                           aGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAACGATTCGGACCTGCAGCAGGTttgagttgag                                                                            >  2:257331/1‑143 (MQ=17)
                                                           aGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgaggtgag                                                                            >  2:163029/1‑143 (MQ=17)
                                                               ggTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCCAAAGGCGCTTTTTTACTATCTGCGATACCCAAACATTCGAACCTCCAGCAGGCTTGAGGTGAgcgcg                                                                        >  2:12194/1‑142 (MQ=11)
                                                                gTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCCGGTTTGAGTTGAgcgcgg                                                                       >  1:216045/1‑141 (MQ=17)
                                                                            aCTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGa                                                           <  1:160832/143‑1 (MQ=32)
                                                                             cTTTTAATCAATTGGTGGCAGGTCCGACTCCTGCACGACCCACCAATGTACAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAg                                                          <  2:793042/143‑1 (MQ=255)
                                                                              ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCCGGTTTGAGTTGCGCGCCGCGAAACAACggcta                                                         >  2:591337/1‑140 (MQ=255)
                                                                              ttttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAGCATCGGAGc                                                         >  1:649736/1‑143 (MQ=17)
                                                                               tttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTCGAGCGCAGCGAAACAACGGAGcc                                                        >  1:84219/1‑143 (MQ=17)
                                                                                      aaTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTCAGCGGAGTGGGAGAACGGAGCCGCCcgcg                                                 >  2:414898/1‑143 (MQ=11)
                                                                                            tCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGCGCCGCTCGCGGCGACg                                           >  1:1039990/1‑143 (MQ=17)
                                                                                            tCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACg                                           >  1:125523/1‑143 (MQ=32)
                                                                                              gCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCCGGTTAGAGTAGAGCGACGACAAGCAAGGGAGCCGCGCGCGGCGACCGc                                         >  2:365312/1‑143 (MQ=11)
                                                                                                 ggTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCg                                      <  1:722611/143‑1 (MQ=32)
                                                                                                     cGAGTCCTGCACTACCCACCGATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAgg                                  <  2:438919/143‑1 (MQ=12)
                                                                                                            tGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGCACCTGCAGCCGGTTTGCGGTGAGCGAAGCGCAACAACGGAGCCGCGCGCGGCGGAGGCCCGACGGgcgaacg                           >  2:912336/1‑143 (MQ=11)
                                                                                                                  aCCCACGAACGTAACGCAGCGCCCTAATGCCGCTTTTTTACTACCTGCGATCCTCAAAGTTTCGCACCGGCAGCAGGCTTGGGCTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGagcgaagcga                     <  2:978597/143‑1 (MQ=11)
                                                                                                                  aCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGagcgaagcga                     <  1:259058/143‑1 (MQ=32)
                                                                                                                    ccaccaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAGCCGCCCTCAGGTTTGCGTTGTGCGCAGTGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGTGGCGAGt                   >  2:208055/1‑143 (MQ=11)
                                                                                                                      accaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCa                 >  1:587592/1‑143 (MQ=32)
                                                                                                                            gTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGc           <  1:903740/143‑1 (MQ=32)
                                                                                                                            gTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGc           <  1:550196/143‑1 (MQ=32)
                                                                                                                             gaaaaaaGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAGCCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCa          <  2:784715/142‑1 (MQ=17)
                                                                                                                              aaaaaaGCGCCTTTAGGGCGCTTTTTTACTATCTGCGATCCTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCAc         >  1:729321/1‑143 (MQ=11)
                                                                                                                               aaaaaGCGCCTTTAGGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCCGCAGGTTTGAGTTTTGCTCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACg        >  1:397307/1‑143 (MQ=11)
                                                                                                                                    gcgcCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACcca   >  1:496731/1‑143 (MQ=32)
                                                                                                                                     cgcCCTAAAGGCGCTTTTTTACGATCTGCGATACTCAAAGATTCGACCCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccac  <  2:12785/143‑1 (MQ=17)
                                                                                                                                             |                                                                                                                                      
ACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCAC  >  NC_000913/2521207‑2521482

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: