Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 594,946 C→A 8.0% Q166H (CAG→CAT) ‡ cusR ← DNA‑binding transcriptional activator CusR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913594,9460CA 8.0% 124.4 / 2.8 51Q166H (CAG→CAT) ‡cusRDNA‑binding transcriptional activator CusR
Reads supporting (aligned to +/- strand):  ref base C (27/20);  new base A (2/2);  total (29/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.63e-01

GAATTAGCTTCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCT‑TTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCC  >  NC_000913/594810‑595087
                                                                                                                                        |                                                                                                                                              
gAATTAGCTTCGGCTCAAAGTCGTCGTCCATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATgg                                                                                                                                          <  1:482712/143‑1 (MQ=255)
     aGCTTCGGCTCATAGTCGTTGTCGATTGTGGCGCGCAGGCGCGTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATCCCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCATCAGGGCgaag                                                                                                                                     <  2:841687/143‑1 (MQ=255)
     aGCTTCGGCTCAAAGTCGTTGTCGAGTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATGGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCATCAGGGCgaag                                                                                                                                     <  1:1029962/143‑1 (MQ=255)
      cctccGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCgaagg                                                                                                                                    <  1:1048238/139‑1 (MQ=255)
           ggCTCAAAGTCGTTGTCGATTTTTCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTTGTATCTCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGgaagaa                                                                                                                               <  1:261487/143‑1 (MQ=255)
           ggCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGgaagaa                                                                                                                               <  2:652467/143‑1 (MQ=255)
               caacgccGTTGGCGCTGTTGCCGCGCAGCCGCCTCACCGCCACATCAATAGCATTGGTATCGCTGTGAAAATTCATGTCCCATACCTGCCAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTcc                                                                                                                           <  1:550826/137‑1 (MQ=255)
                aaaGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGGTGGCGAAGGAAGAACTCca                                                                                                                          >  1:46946/1‑143 (MQ=255)
                          tCGATTTTGCCGAGCAGACGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGAGACAGCACTACACCCTGATGGCGAAGGAAGAACTCCAGCAGAGCaaa                                                                                                                >  2:1067854/1‑143 (MQ=255)
                            gctttcGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCAAGAAGCTAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAAct                                                                                                              <  1:402127/141‑1 (MQ=255)
                            gATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAGGGAAGAACTCCAGCAGAGTacacc                                                                                                              >  1:265882/1‑142 (MQ=255)
                            gATTTTGCCGCGCAGCCGCTTAACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGCAGAACTGCAGCAGAGTAAAca                                                                                                              >  2:749617/1‑142 (MQ=255)
                             aTTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGCGGCAGCACTTCGCCCTGATGGCGCAGGAAGAACTCCAGCAGAGCAACctc                                                                                                             >  1:193608/1‑143 (MQ=255)
                                  gCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTGT‑TTAc                                                                                                       >  2:959103/1‑143 (MQ=255)
                                   ccGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCT‑TTACt                                                                                                      >  1:988127/1‑143 (MQ=255)
                                          cgcGCTGCACCGCCACCTCAATAGCATTGGGATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCT‑TTACTGGTCaaa                                                                                               <  2:265882/141‑1 (MQ=255)
                                           ccGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTACACTCT‑TTACTGGTCAAAg                                                                                              >  1:859570/1‑143 (MQ=255)
                                              cTTCGCCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCT‑TTACTGGTCAAAGTGa                                                                                           <  1:938086/143‑1 (MQ=255)
                                                       cacaTCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCT‑TTACTGGTCAAAGTGATGCGCGTGc                                                                                  >  2:973180/1‑143 (MQ=255)
                                                         caTTAATAGCATGGGTATCGCTGTCAATATTCTTGTCCCATACCTGCGAGGCGATAAGCGGGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCT‑TTACTGGTCAAAGTGATGCGCGTGCCg                                                                                <  2:631852/143‑1 (MQ=255)
                                                             aaTAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATAGCGAAGGAAGAACTCCAGCAGAGTAAACAAA‑CAACTGGTGATAGGGATGCGCGTGCCGCTGc                                                                            >  1:538948/1‑143 (MQ=255)
                                                                    ttGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGCGTAAACTCT‑TTACTGGTCAAAGTGATGCGCGTGCGGCTGCGGGTGaa                                                                     >  2:347783/1‑142 (MQ=255)
                                                                    ttGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGCCGAAGGAAGAACTCCAGCAGAGCAAACTCT‑TCACTGGTCAAAGTGATGCGCGTCCCGCTGCggcccac                                                                     >  1:421586/1‑138 (MQ=255)
                                                                      ggTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCT‑TTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACtt                                                                   >  2:1009845/1‑143 (MQ=255)
                                                                         gtCGCTGTCAAAACTCATGTCCTCTACCTGCGAGGCGATAAGCGAGCGAGGCACCACTTGATCCTTATGGCGAAGGAAGAACTTCAGCAGAGTAAACTCT‑TTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGc                                                                <  2:930152/142‑1 (MQ=255)
                                                                         agcgttgtgAAAGTGCACGTCCCATACCTGCGAGGCGATAAGCGGGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCT‑TTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGc                                                                <  2:916056/138‑1 (MQ=255)
                                                                            gCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCT‑TTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGc                                                             >  2:823535/1‑143 (MQ=255)
                                                                             cTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCCGCAGAGTAAACTCT‑TTACTGGTCAACGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCt                                                            >  1:73950/1‑143 (MQ=255)
                                                                               gTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCT‑TTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGa                                                          <  1:53246/143‑1 (MQ=255)
                                                                                  aaaaTTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCT‑TTACTGGTCCAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGTGCCTGACga                                                       >  2:712197/1‑143 (MQ=255)
                                                                                  aaaaTTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCT‑TTAATGGTCAAAGTGATGCGCGTGCCGCTGCGGGGGACTTGGCGGCGGACga                                                       >  2:810560/1‑143 (MQ=255)
                                                                                   aaaTTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCT‑TTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCGGACgag                                                      >  2:338675/1‑143 (MQ=255)
                                                                                   aaaTTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCG‑TTACTGGTCAAAGTAATGCGCGTGCCGCTGCGGGGGACTTTGCGGCTGACgag                                                      >  2:540138/1‑143 (MQ=255)
                                                                                    aaTTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCT‑TTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACgaga                                                     <  1:507123/143‑1 (MQ=255)
                                                                                     aTTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAG‑ACACGATGCGCCACCTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGAt                                                    >  1:1064885/1‑143 (MQ=255)
                                                                                       tCATGTCCCGTACCTGCGAGGCGATAAGCGAGCGGGGCAGCACCTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCT‑TTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCg                                                  <  1:715424/143‑1 (MQ=255)
                                                                                            tCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCGGAGGAAACTCT‑TTACTGGTCGAAGTGAGGCGCGCGCCGCTGCGGGTGACTTCGCGGCCGACGAGCTCGAccct                                             >  1:665386/1‑141 (MQ=255)
                                                                                                   ccTCCGAGTCGAGACGCGAGCGGGTCAGCACTCCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAAGTCT‑TGACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCg                                      <  2:147845/143‑1 (MQ=255)
                                                                                                     tGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGACGAACTCCAGCAGGGCAAACTCT‑TTACTGGTCAAAGTGATGCGCGCGCCGCTGCAGGTCACGATCCGTCTGCCGAGCTCGACCACCACCTCGGc                                    >  1:829823/1‑143 (MQ=255)
                                                                                                         aGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCT‑TTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGCCTTCGCGGCTGGCGAGATCGACCACCAAATCGGCAAcc                                >  2:365627/1‑143 (MQ=255)
                                                                                                            cGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCT‑TTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGa                             >  2:721805/1‑143 (MQ=255)
                                                                                                             gATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCT‑TTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGaa                            >  2:138932/1‑143 (MQ=255)
                                                                                                               tAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCT‑TTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCACATCGGCAACCTGAAAc                          >  1:887109/1‑143 (MQ=255)
                                                                                                                  gcgagcgGTGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCT‑TTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGa                       <  1:106382/143‑1 (MQ=255)
                                                                                                                  gcgagcgGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCT‑TTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGa                       <  1:1068689/143‑1 (MQ=255)
                                                                                                                    gggcgGGGCAGCACTTCGCCCTGAGGGCGAAGGAAGAACTCCAGCAGAGTAAACTCT‑TTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACt                     <  2:1014853/141‑1 (MQ=255)
                                                                                                                     agcgGGGCAGCACTTCGCCCTGCTGGCGAAGGAAGAACTCCAGCAGAGTAAACTCT‑TTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTCTGCGGCGGCCGAGATGGACAATCCAATCCGCAACCTGaaagcaactt                    >  2:213320/1‑136 (MQ=255)
                                                                                                                        ggggCAGCACTGCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCT‑TTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAAGCGGCAACCTGAAACTGACTTTCg                 <  1:263910/143‑1 (MQ=255)
                                                                                                                             agcaCTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCT‑TTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAAt            <  2:46946/143‑1 (MQ=255)
                                                                                                                                ccTTCGCCATCAGGGCGAAGGAAGAACTCCAGCAGAGTAAACTCT‑TTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGCGACATACCACCATATCGGCACCCTGCAACTGACTTTCGATAATCAc         >  2:754720/2‑143 (MQ=255)
                                                                                                                                 cTTCGCCATCAGGGCGAAGGAAGAACTCCAGCAGAGTAAACTCT‑TTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGCCGCGATCGCCCATCGACTCGGCTGCCTCAAACTGGCTTTCGGTAGTCAcc        >  1:299063/1‑143 (MQ=255)
                                                                                                                                       ccTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCT‑TTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCCACCGGAAACTGACTTGCGATAATCACCgccgcc  >  2:750280/1‑143 (MQ=255)
                                                                                                                                        |                                                                                                                                              
GAATTAGCTTCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCT‑TTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCC  >  NC_000913/594810‑595087

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