Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 597,758 (G)5→6 coding (280/1224 nt) cusB → copper/silver efflux system, membrane fusion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913597,7531.G97.3% 129.1 / ‑3.9 37L92R (CTG→CGG) cusBcopper/silver efflux system, membrane fusion protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  major base G (21/15);  minor base A (1/0);  total (22/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.52e-01

TCTTATTCTGGTACGACCCAATGTATCCCAATACGCGGTTCGATAAACCAGGTAAATCGCCGTTTATGGATATGGATCTGGTGCCGAAATATGCCGATGAAGAGAGTTCTGCGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCTGGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCTGACTTTTGCCCAGAGTTTCCCGGCGAATGTCAGTTACAACGAGTATCAGTATGCCATTGTGCAGGCTCGCGCTGCCGGGTTTATCGACAAGGTGTATCCGCTTACCGTG  >  NC_000913/597603‑597898
                                                                                                                                                       |                                                                                                                                                 
tCTTATTCTGGTACGACCCAATGTATCCCAATACGCGGTTCGATAAACCAGGTAAATCGCCGTTTATGGATATGGATCTGGTGCCGAAATATGCCGATGAAGAGAGTTCTGCGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCt                                                                                                                                                    >  2:649068/1‑151 (MQ=255)
   tATTCTGGTACGACCCAATGTATCCCAATACGCGGTTCGATAAACCAGGTAAATCGCCGTTTATGGATATGGATCTGGTGCCGAAATATGCCGATGAAGAGAGTTCTGCGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCt                                                                                                                                                    <  1:51722/148‑1 (MQ=255)
   tATTCTGGTACGACCCAATGTATCCCAATACGCGGTTCGATAAACCAGGTAAATCGCCGTTTATGGATATGGATCTGGTGCCGAAATATGCCGATGAAGAGAGTTCTGCGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCTgg                                                                                                                                                 <  2:384168/150‑1 (MQ=255)
    aTTCTGGTACGACCCAATGTATCCCAATACGCGGTTCGATAAACCAGGTAAATCGCCGTTTATGGATATGGATCTGGTGCCGAAATATGCCGATGAAGAGAGTTCTGCGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCTgggt                                                                                                                                               <  1:589069/151‑2 (MQ=255)
      tCTGGTACGACCCAATGTATCCCAATACGCGGTTCGATAAACCAGGTAAATCGCCGTTTATGGATATGGATCTGGTGCCGAAATATGCCGATGAAGAGAGTTCTGCGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCTgggggg                                                                                                                                             >  2:39746/1‑150 (MQ=255)
         ggTACGACCCAATGTATCCCAATACGCGGTTCGATAAACCAGGTAAATCGCCGTTTATGGATATGGATCTGGTGCCGAAATATGCCGATGAAGAGAGTTCTGCGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCTggggggt                                                                                                                                            <  2:733896/149‑3 (MQ=255)
                         tCCCAATACGCGGTTCGATAAACCAGGTAAATCGCCGTTTATGGATATGGATCTGGTGCCGAAATATGCCGATGAAGAGAGTTCTGCGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCTGGGGGGTGAAAACGGCTACCGTCAc                                                                                                                           >  1:373679/1‑151 (MQ=255)
                          cccAATACGCGGTTCGATAAACCAGGTAAATCGCCGTTTATGGATATGGATCTGGTGCCGAAATATGCCGATGAAGAGAGTTCTGCGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCTGGGGGGTGAAAACGGCTACCGTCAc                                                                                                                           <  1:482141/150‑1 (MQ=255)
                           ccAATACGCGGTTCGATAAACCAGGTAAATCGCCGTTTATGGATATGGATCTGGTGCCGAAATATGCCGATGAAGAGAGTTCTGCGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCTGGGGGGTGAAAACGGCTACCGTCAcgc                                                                                                                         >  2:100497/1‑151 (MQ=255)
                                   cgGTTCGATAAACCAGGTAAATCGCCGTTTATGGATATGGATCTGGTGCCGAAATATGCCGATGAAGAGAGTTCTGCGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCTGGGGGGTGAAAACGGCTACCGTCACGCGCgg                                                                                                                     <  1:260012/147‑1 (MQ=255)
                                   cgGTTCGATAAACCAGGTAAATCGCCGTTTATGGATATGGATCTGGTGCCGAAATATGCCGATGAAGAGAGTTCTGCGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCTGGGGGGTGAAAACGGCTACCGTCACGCGCgg                                                                                                                     >  2:260045/1‑147 (MQ=255)
                                                      aaTCGCCGTTTATGGATATGGATCTGGTGCCGAAATATGCCGATGAAGAGAGTTCTGCGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCTGGGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCTGACTTTTGCCCAGAGtt                                                                                              >  2:176788/1‑151 (MQ=255)
                                                       aTCGCCGTTTATGGATATGGATCTGGTGCCGAAATATGCCGATGAAGAGAGTTCTGCGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCTGGGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCTGACTTTTGCCCAGAGttt                                                                                             >  2:523826/1‑151 (MQ=255)
                                                            cGTTTATGGATATGGATCTGGTGCCGAAATATGCCGATGAAGAGAGTTCTGCGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCTGGGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCTGACTTTTGCCCAGAGTTTCCCgg                                                                                        >  1:323611/1‑151 (MQ=255)
                                                               ttATGGATATGGATCTGGTGCCGAAATATGCCGATGAAGAGAGTTCTGCGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCTGGGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCTGACTTTTGCCCAGAGTTTCCCGGCg                                                                                      >  1:229511/1‑150 (MQ=255)
                                                                 atggatatggatCTGGTGCCGAAATATGCCGATGAAGAGAGTTCTGCGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCTGGGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCTGACTTTTGCCCAGAGTTTCCCGGCGAAt                                                                                   <  2:373726/151‑1 (MQ=255)
                                                                           atCTGGTGCCGAAATATGCCGATGAAGAGAGTTCTGCGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCTGGGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCTGACTTTTGCCCAGAGTTTCCCGGCGAATGTCAGTTAc                                                                          >  1:328156/1‑150 (MQ=255)
                                                                           atCTGGTGCCGAAATATGCCGATGAAGAGAGTTCTGCGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCTGGGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCTGACTTTTGCCCAGAGTTTCCCGGCGAATGTCAGTTAc                                                                          >  2:135505/1‑150 (MQ=255)
                                                                            tCTGGTGCCGAAATATGCCGATGAAGAGAGTTCTGCGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCTGGGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCTGACTTTTGCCCAGAGTTTCCCGGCGAATGTCAGTTACa                                                                         >  2:733393/1‑150 (MQ=255)
                                                                               ggTGCCGAAATATGCCGATGAAGAGAGTTCTGCGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCTGGGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCTGACTTTTGCCCAGAGTTTCCCGGCGAATGTCAGTTACAAc                                                                       <  2:704043/149‑1 (MQ=255)
                                                                               ggTGCCGAAATATGCCGATGAAGAGAGTTCTGCGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCTGGGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCTGACTTTTGCCCAGAGTTTCCCGGCGAATGTCAGTTACAACg                                                                      <  2:768429/150‑1 (MQ=255)
                                                                                 tGCCGAAATATGCCGATGAAGAGAGTTCTGCGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCTGGGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCTGACTTTTGCCCAGAGTTTCCCGGCGAATGTCAGTTACAACGAGt                                                                   >  1:181631/1‑151 (MQ=255)
                                                                                   ccGAAATATGCCGATGAAGAGAGTTCTGCGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCTGGGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCTGACTTTTGCCCAGAGTTTCCCGGCGAATGTCAGTTACAACGAGTAt                                                                 >  1:279163/1‑151 (MQ=255)
                                                                                            gCCGATGAAGAGAGTTCTGCGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCTGGGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCTGACTTTTGCCCAGAGTTTCCCGGCGAATGTCAGTTACAACGAGTATCAGTATGcc                                                        >  2:631710/1‑151 (MQ=255)
                                                                                            gCCGATGAAGAGAGTTCTGCGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCTGGGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCTGACTTTTGCCCAGAGTTTCCCGGCGAATGTCAGTTACAACGAGTATCAGTATGc                                                         <  1:649010/150‑1 (MQ=255)
                                                                                              cGATGAAGAGAGTTCTGCGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCTGGGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCTGACTTTTGCCCAGAGTTTCCCGGCGAATGTCAGTTACAACGAGTATCAGTATGCCAt                                                      >  2:535683/1‑151 (MQ=255)
                                                                                                  gAAGAGAGTTCTGCGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCTGGGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCTGACTTTTGCCCAGAGTTTCCCGGCGAATGTCAGTTACAACGAGTATCAGTATGCCATtgtg                                                  >  2:480816/1‑151 (MQ=255)
                                                                                                  gAAGAGAGTTCTGCGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCTGGGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCTGACTTTTGCCCAGAGTTTCCCGGCGAATGTCAGTTACAACGAGTATCAGTATGCCATtgtg                                                  <  1:523774/151‑1 (MQ=255)
                                                                                                   aaGAGAGTTCTGCGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCTGGGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCTGACTTTTGCCCAGAGTTTCCCGGCGAATGTCAGTTACAACGAGTATCAGTATGCCATTGTGc                                                 >  1:36160/1‑151 (MQ=255)
                                                                                                     gagagTTCTGCGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCTGGGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCTGACTTTTGCCCAGAGTTTCCCGGCGAATGTCAGTTACAACGAGTATCAGTATGCCATTGTGCAg                                               <  1:631653/151‑1 (MQ=255)
                                                                                                      agagTTCTGCGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCTGGGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCTGACTTTTGCCCAGAGTTTCCCGGCGAATGTCAGTTACAACGAGTATCAGTATGCCATTGTGCAgg                                              >  2:725002/1‑151 (MQ=255)
                                                                                                      agagTTCTGCGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCTGGGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCTGACTTTTGCCCAGAGTTTCCCGGCGAATGTCAGTTACAACGAGTATCAGTATGCCATTGTGCAg                                               <  2:704719/150‑1 (MQ=255)
                                                                                                       gagTTCTGCGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCTGGGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCTGACTTTTGCCCAGAGTTTCCCGGCGAATGTCAGTTACAACGAGTATCAGTATGCCATTGTGCAGGc                                             <  2:506831/151‑1 (MQ=255)
                                                                                                         gTTCTGCGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCTAGGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCTGACTTTTGCCCAGAGTTTCCCGGCGAATGTCAGTTACAACGAGTATCAGTATGCCATTGTGCAGGCTc                                           >  1:683693/1‑151 (MQ=255)
                                                                                                               cGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCTGGGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCTGACTTTTGCCCAGAGTTTCCCGGCGAATGTCAGTTACAACGAGTATCAGTATGCCATTGTGCAGGCTCgcgc                                       >  2:549385/1‑149 (MQ=255)
                                                                                                                    ggTGTGCGCATTGACCCGACTCAGACGCAGAATCTGGGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCTGACTTTTGCCCAGAGTTTCCCGGCGAATGTCAGTTACAACGAGTATCAGTATGCCATTGTGCAGGCTCGCGCTGcc                                   <  2:527774/148‑1 (MQ=255)
                                                                                                                       gtgCGCATTGACCCGACTCAGACGCAGAATCTGGGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCTGACTTTTGCCCAGAGTTTCCCGGCGAATGTCAGTTACAACGAGTATCAGTATGCCATTGTGCAGGCTCGCGCTGCCGGGttt                             >  1:535165/1‑151 (MQ=255)
                                                                                                                        tgCGCATTGACCCGACTCAGACGCAGAATCTGGGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCTGACTTTTGCCCAGAGTTTCCCGGCGAATGTCAGTTACAACGAGTATCAGTATGCCATTGTGCAGGCTCGCGCTGCCGGGttt                             <  1:337726/150‑1 (MQ=255)
                                                                                                                         gcgcATTGACCCGACTCAGACGCAGAATCTGGGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCTGACTTTTGCCCAGAGTTTCCCGGCGAATGTCAGTTACAACGAGTATCAGTATGCCATTGTGCAGGCTCGCGCTGCCGGGTTTAt                           <  1:479442/151‑1 (MQ=255)
                                                                                                                          cgcATTGACCCGACTCAGACGCAGAATCTGGGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCTGACTTTTGCCCAGAGTTTCCCGGCGAATGTCAGTTACAACGAGTATCAGTATGCCATTGTGCAGGCTCGCGCTGCCGGGTTTAt                           <  2:631811/150‑1 (MQ=255)
                                                                                                                           gcATTGACCCGACTCAGACGCAGAATCTGGGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCTGACTTTTGCCCAGAGTTTCCCGGCGAATGTCAGTTACAACGAGTATCAGTATGCCATTGTGCAGGCTCGCGCTGCCGGGTTTATCg                         <  2:406037/151‑1 (MQ=255)
                                                                                                                                    cGACTCAGACGCAGAATCTGGGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCTGACTTTTGCCCAGAGTTTCCCGGCGAATGTCGGTTACAACGAGTATCAGTATGCCATTGTGCAGGCTCGCGCTGCCGGGTTTATCGACAAGgtgt                >  1:181422/1‑151 (MQ=255)
                                                                                                                                                  aaTCTGGGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCTGACTTTTGCCCAGAGTTTCCCGGCGAATGTCAGTTACAACGAGTATCAGTATGCCATTGTGCAGGCTCGCGCTGCCGGGTTTATCGACAAGGGGTATCCGCTTACCGTg  >  1:503997/1‑151 (MQ=255)
                                                                                                                                                       |                                                                                                                                                 
TCTTATTCTGGTACGACCCAATGTATCCCAATACGCGGTTCGATAAACCAGGTAAATCGCCGTTTATGGATATGGATCTGGTGCCGAAATATGCCGATGAAGAGAGTTCTGCGTCTGGTGTGCGCATTGACCCGACTCAGACGCAGAATCTGGGGGTGAAAACGGCTACCGTCACGCGCGGACCGCTGACTTTTGCCCAGAGTTTCCCGGCGAATGTCAGTTACAACGAGTATCAGTATGCCATTGTGCAGGCTCGCGCTGCCGGGTTTATCGACAAGGTGTATCCGCTTACCGTG  >  NC_000913/597603‑597898

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: