Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 594,950 T→G 11.2% H165P (CAT→CCT)  cusR ← DNA‑binding transcriptional activator CusR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913594,9500TG11.2% 67.4 / 2.6 39H165P (CAT→CCT) cusRDNA‑binding transcriptional activator CusR
Reads supporting (aligned to +/- strand):  ref base T (20/13);  new base G (2/3);  total (23/16)
Fisher's exact test for biased strand distribution p-value = 6.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.19e-01

GAATTAGCTTCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCC  >  NC_000913/594810‑595089
                                                                                                                                            |                                                                                                                                           
gAATTAGCTTCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCACATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATgg                                                                                                                                           >  1:172105/1‑143 (MQ=255)
   ttAGCTTCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTGCGCCCTGATGGCga                                                                                                                                        >  1:582160/1‑143 (MQ=255)
      gCTTCGGCTCAAAGTCGTTGTCGATTTTGCCGCTCCGCCGCTTCACTGCTACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCATCAGGGCgaagt                                                                                                                                     <  2:514662/143‑2 (MQ=255)
      gCTTCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCATCAGGGCgaagt                                                                                                                                     <  2:339806/143‑2 (MQ=255)
      gCTGCGGCGCAAAGTCGGTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCATCAGGGCgaagt                                                                                                                                     <  1:340133/143‑2 (MQ=255)
         tCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCCACATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCCATTCGCCCTGACGGCgcaggagg                                                                                                                                  >  2:526559/1‑141 (MQ=255)
         tCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAACCGAGCGGGGAAGCACTTCGCCCGACTGGCgaaggaag                                                                                                                                  >  2:511894/1‑143 (MQ=255)
         tCGGCTCAAAGTCGTTGCCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCCATAGCATTGGTATAACCGCCAAACTTGATGTCCCATACCTGCGAGGCGAGAAGCGAGCGGGGCAGCACTTCGCGCTGCGGGGgaaggggg                                                                                                                                  >  2:237390/1‑140 (MQ=255)
           ggCTCAAAGTCTTGTTCGGTTTTGCCGCGCCGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGAGAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGgaagaa                                                                                                                                <  2:229346/143‑1 (MQ=255)
                acagccGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACAGCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCGGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCca                                                                                                                           <  2:352931/138‑1 (MQ=255)
                         gTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCcagcagcgcac                                                                                                                  >  2:136558/1‑138 (MQ=255)
                          tCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATAGCGACCGAAGAACTCCAGCCGAGTaag                                                                                                                 >  2:86432/1‑142 (MQ=255)
                          tCGATTTTGCCGCGCAGCCGCATCACCGCCACATCAATAGCATTGGTATCGCTGCCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACGTCGCCCTGATGGCGCAGGAAGAACTACAGCagagacaa                                                                                                                 >  2:272587/1‑139 (MQ=255)
                          tCGATTTTGCCGCGCAGCCACTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGCTAAGCGAGCGGGGCCGCACTTCGCCCTGATGGCGAAGGGAGCACTCCAGCCGAGTcaa                                                                                                                 >  2:111447/1‑143 (MQ=255)
                           cGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGGCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAAc                                                                                                                >  2:132120/1‑143 (MQ=255)
                             aatttGCCGCGCAGCCGCTTCACCGCCACCTCAATAGCATTGGTATCGCTGTCCAACTTCATGTCCCATACCTGCGAGGCGATCAGCGAGCGGGGCAGCACTTCGCCCGGATGGCGAAGGACGCCCCCCCGCAGAGCAAAccc                                                                                                              >  2:231789/3‑141 (MQ=255)
                              ttttGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATACGAATCGCTGTCAAACTTCAAGTCCCATAACTGCGAGGCGATAAGCGAGCGGGGCAGCACCTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTACActat                                                                                                             >  1:166551/1‑141 (MQ=255)
                                  gCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGCGGCGATAAGCGCGCGGGGCAGCCCTTCGCCCTGATGGCGACGGAAGAACTCCAGCAGCGCAAACTCTTTAc                                                                                                         >  1:286539/1‑143 (MQ=255)
                                  gCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCGCTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGGAAACTCTTTAc                                                                                                         >  1:136502/1‑143 (MQ=255)
                                   ccGCGCAGCCGCTTCACTGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAAGCGATAAGCGAGCGGGCGAGCACCTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACt                                                                                                        >  2:638217/1‑143 (MQ=255)
                                          gCCGCTTCCCCGCCAAATCAATCGCAGTGGTACCGCTGTCAAAATTGATGTTCCAGACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAATGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCaaa                                                                                                 <  1:70855/143‑1 (MQ=255)
                                           ccGCTTCAACGCCCCAACAACAGCCTTGGTCTCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTTATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAg                                                                                                <  1:752334/143‑1 (MQ=255)
                                                      ccACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTg                                                                                     <  1:721930/143‑1 (MQ=255)
                                                               tagaagtGGTATCGCTGTCAAAATTCATGTCCCATGCCGGGGAGGCGATAAGCGAGCGGGGCGGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCgg                                                                            <  1:332837/137‑1 (MQ=255)
                                                               tAGCATTGGTGTCGCTGCCCCCCTCCATGTTCCATACCTGCGAGGCGATAGGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCgg                                                                            <  2:374833/143‑1 (MQ=255)
                                                                   actactATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGa                                                                        <  1:503541/138‑1 (MQ=255)
                                                                   aTTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCCGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGGGCCGCTGCGGGTGa                                                                        >  1:646929/1‑143 (MQ=255)
                                                                        tATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTg                                                                   >  1:713389/1‑143 (MQ=255)
                                                                                tCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGCACTCCAGCATAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGAc                                                           >  1:575278/1‑143 (MQ=255)
                                                                                    aaTTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTGGTCTTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACgaga                                                       <  1:443598/143‑1 (MQ=255)
                                                                                    aaTTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCCCCCTGAGGGCGAGGGAAGAACTCCAGCAGAGTAAACCCTTACCCGGCCAAGGTGAGCCGCGTGCCGCCCCGGGTGCCTTGGCGCCTGACgaga                                                       >  2:126495/1‑143 (MQ=255)
                                                                                     aTTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACCCTTTAACGGTCAAAGTGATGCGCGTGCCGCCGCGGGTGACTTTGCGGCTGACGACAt                                                      >  1:292578/1‑143 (MQ=255)
                                                                                                           gCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAGCTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACGATCACATCGGCGACCTg                                >  1:752696/1‑143 (MQ=255)
                                                                                                              aTAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGaaa                             <  1:638027/143‑1 (MQ=255)
                                                                                                                         gggCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGGGACTGACttgtgc                  >  1:47783/1‑139 (MQ=255)
                                                                                                                             agcaCGTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAAt              <  1:725667/143‑1 (MQ=255)
                                                                                                                              gggggtggCCTTGGTGGGGAAGGAGGGGCTCCAGCAGAGTAAACGCGTTACTGGTCAAAGTGATGCGCGTGCCGTTGCGGGGGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATc             <  2:335871/136‑1 (MQ=255)
                                                                                                                                 cTTCGACCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACGCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACGTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCAcc          <  1:415799/143‑1 (MQ=255)
                                                                                                                                         tGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGAACGACCATCGAATCGGCAACCTGAAACTGACTGTCGATAAACAACGCCGcccc  >  2:556350/1‑143 (MQ=255)
                                                                                                                                            |                                                                                                                                           
GAATTAGCTTCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCC  >  NC_000913/594810‑595089

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