Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 600,325 A→G T538A (ACC→GCC)  cusA → copper/silver export system RND permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913600,3250AG97.3% 119.2 / ‑3.9 37T538A (ACC→GCC) cusAcopper/silver export system RND permease
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (17/19);  total (18/19)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

GCAGTAACCCGCTCAATCGCTTTTTGATTCGTGTTTATCATCCGCTGTTGCTGAAAGTACTGCACTGGCCGAAAACCACGCTGCTGGTGGCGGCGCTTTCGGTGCTGACGGTTCTCTGGCCGCTCAATAAAGTTGGCGGGGAATTTTTACC  >  NC_000913/600251‑600401
                                                                          |                                                                            
gCAGTAACCCGCTCAATCGCTTTTTGATTCGTGTTTATCATCCGCTGTTGCTGAAAGTACTGCACTGGCCGAAAGCCACgc                                                                        <  1:130741/81‑1 (MQ=255)
gCAGTAACCCGCTCAATCGCTTTTTGATTCGTGTTTATCATCCGCTGTTGCTGAAAGTACTGCACTGGCCGAAAGCCACgc                                                                        <  2:130740/81‑1 (MQ=255)
  aGTAACCCGCTCAATCGCTTTTTGATTCGTGTTTATCATCCGCTGTTGCTGAAAGTACTGCACTGGCCGAAAGCCACgctg                                                                      <  1:130742/81‑1 (MQ=255)
  aGTAACCCGCTCAATCGCTTTTTGATTCGTGTTTATCATCCGCTGTTGCTGAAAGTACTGCACTGGCCGAAAGCCACgctg                                                                      >  1:130790/1‑81 (MQ=255)
     aaCCCGCTCAATCGCTTTTTGATTCGTGTTTATCATCCGCTGTTGCTGAAAGTACTGCACTGGCCGAAAGCCACGctgctg                                                                   >  2:130788/1‑81 (MQ=255)
          gCTCAATCGCTTTTTGATTCGTGTTTATCATCCGCTGTTGCTGAAAGTACTGCACTGGCCGAAAGCCACGCTGCTGGTggc                                                              >  2:130787/1‑81 (MQ=255)
              aaTCGTCTTTTGATTCGTGTTTAACATCCGCTGTTGCTGAAAGTACTGCACTGGCCGAAAGCCACGCTGCTGGTGgcggcg                                                          >  2:130794/1‑81 (MQ=255)
              aaTCGCTTTTTGATTCGTGTTTATCATCCGCTGTTGCTGAAAGTACTGCACTGGCCGAAAGCCACGCTGCTGGTggcggc                                                           <  1:130743/80‑1 (MQ=255)
                  gCTTTTTGATTCGTGTTTATCATCCGCTGTTGCTGAAAGTACTGCACTGGCCGAAAGCCACGCTGCTGGTGGCGGCGCttt                                                      <  1:130744/81‑1 (MQ=255)
                  gCTTTTTGATTCGTGTTTATCATCCGCTGTTGCTGAAAGTACTGCACTGGCCGAAAGCCACGCTGCTGGTGGCGGCGCtt                                                       <  1:130745/80‑1 (MQ=255)
                   cTTTTTGATTCGTGTTTATCATCCGCTGTTGCTGAAAGTACTGCACTGGCCGAAAGCCACGCTGCTGGTGGCGGCGCTTTc                                                     >  2:130762/1‑81 (MQ=255)
                   cTTTTTGATTCGTGTTTATCATCCGCTGTTGCTGAAAGTACTGCACTGGCCGAAAGCCACGCTGCTGGTGGCGGCGCTTTc                                                     <  2:130746/81‑1 (MQ=255)
                          aTTCGTGTTTATCATCCGCTGTTGCTGAAAGTACTGCACTGGCCGAAAGCCACGCTGCTGGTGGCGGCGCTTTCGGTGCt                                               <  1:130747/80‑1 (MQ=255)
                             cGTGTTTATCATCCGCTGTTGCTGAAAGTACTGCACTGGCCGAAAGCCACGCTGCTGGTGGCGGCGCTTTCGGTGCTGAc                                            <  2:130748/80‑1 (MQ=255)
                              gtgtTTATCATCCGCTGTTGCTGAAAGTACTGCACTGGCCGAAAGCCACGCTGCTGGTGGCGGCGCTTTCGGTGCTGACgg                                          >  2:130798/1‑81 (MQ=255)
                                 tttATCATCCGCTGTTGCTGAAAGTACTGCACTGGCCGAAAGCCACGCTGCTGGTGGCGGCGCTTTCGGTGCTGACGGtt                                        <  2:130749/80‑1 (MQ=255)
                                   tATCATCCGCTGTTGCTGAAAGTACTGCACTGGCCGAAAGCCACGCTGCTGGTGGCGGCGCTTTCGGTGCTGACGGTtctc                                     >  1:130791/1‑81 (MQ=255)
                                   tATCATCCGCTGTTGCTGAAAGTACTGCACTGGCCGAAAGCCACGCTGCTGGTGGCGGCGCTTTCGGTGCTGACGGTtct                                      <  2:130750/80‑1 (MQ=255)
                                   tATCATCCGCTGTTGCTGAAAGTACTGCACTGGCCGAAAGCCACGCTGCTGGTGGCGGCGCTTTCGGTGCTGACGGTtct                                      <  2:130751/80‑1 (MQ=255)
                                     tcatcCGCTGTTGCTGAAAGTACTGCACTGGCCGAAAGCCACGCTGCTGGTGGCGGCGCTTTCGGTGCTGACGGTTCTCTg                                   <  1:130752/81‑1 (MQ=255)
                                     tcatcCGCTGTTGCTGAAAGTACTGCACTGGCCGAAAGCCACGCTGCTGGTGGCGGCGCTTTCGGTGCTGACGGTTCTCTg                                   >  2:130780/1‑81 (MQ=255)
                                          cGCTGTTGCTGAAAGTACTGCACTGGCCGAAAGCCACGCTGCTGGTGGCGGCGCTTTCGGTGCTGACGGTTCTCTGGCCg                               <  2:130753/80‑1 (MQ=255)
                                           gCTGTTGCTGAAAGTACTGCACTGGCCGAAAGCCACGCTGCTGGTGGCGGCGCTTTCGGTGCTGACGGTTCTCTGGCCGCt                             >  1:130797/1‑81 (MQ=255)
                                             tGTTGCTGAAAGTACTGCACTGGCCGAAAGCCACGCTGCTGGTGGCGGCGCTTTCGGTGCTGACGGTTCTCTGGCCGCTc                            >  2:130805/1‑80 (MQ=255)
                                              gTTGCTGAAAGTACTGCACTGGCCGAAAGCCACGCTGCTGGTGGCGGCGCTTTCGGTGCTGACGGTTCTCTGGCCGCTCaa                          >  1:130792/1‑81 (MQ=255)
                                                 gCTGAAAGTACTGCACTGGCCGAAAGCCACGCTGCTGGTGGCGGCGCTTTCGGTGCTGACGGTTCTCTGGCCGCTCAATa                        >  1:130793/1‑80 (MQ=255)
                                                 gCTGAAAGTACTGCACTGGCCGAAAGCCACGCTGCTGGTGGCGGCGCTTTCGGTGCTGACGGTTCTCTGGCCGCTCAATa                        <  2:130754/80‑1 (MQ=255)
                                                    gAAAGTACTGCACTGGCCGAAAGCCACGCTGCTGGTGGCGGCGCTTTCGGTGCTGACGGTTCTCTGGCCGCTCAATAAAGt                    >  1:130786/1‑81 (MQ=255)
                                                       aGTACTGCACTGGCCGAAAGCCACGCTGCTGGTGGCGGCGCTTTCGGTGCTGACGGTTCTCTGGCCGCTCAATAAAGTTgg                 >  1:130773/1‑81 (MQ=255)
                                                           cTGCACTGGCCGAAAGCCACGCTGCTGGTGGCGGCGCTTTCGGTGCTGACGGTTCTCTGGCCGCTCAATAAAGTTGGCggg             <  2:130756/81‑1 (MQ=255)
                                                           cTGCACTGGCCGAAAGCCACGCTGCTGGTGGCGGCGCTTTCGGTGCTGACGGTTCTCTGGCCGCTCAATAAAGTTGGCgg              <  2:130755/80‑1 (MQ=255)
                                                               aCTGGCCGAAAGCCACGCTGCTGGTGGCGGCGCTTTCGGTGCTGACGGTTCTCTGGCCGCTCAATAAAGTTGGCGGGGAAt         >  2:130800/1‑81 (MQ=255)
                                                                 tGGCCGAAAGCCACGCTGCTGGTGGCGGCGCTTTCGGTGCTGACGGTTCTCTGGCCGCTCAATAAAGTTGGCGGGGAAttt       >  2:130796/1‑81 (MQ=255)
                                                                     cGAAAGCCACGCTGCTGGTGGCGGCGCTTTCGGTGCTGACGGTTCTCTGGCCGCTCAATAAAGTTGGCGGGGAATTTTTAc   <  1:130757/81‑1 (MQ=255)
                                                                     cGAAAGCCACGCTGCTGGTGGCGGCGCTTTCGGTGCTGACGGTTCTCTGGCCGCTCAATAAAGTTGGCGGGGAATTTTTAc   <  2:130758/81‑1 (MQ=255)
                                                                     cGAAAGCCACGCTGCTGGTGGCGGCGCTTTCGGTGCTGACGGTTCTCTGGCCGCTCAATAAAGTTGGCGGGGAATTTTTAc   >  2:130801/1‑81 (MQ=255)
                                                                      gAAAACCACGCTGCTGGTGGCGGCGCTTTCGGTGATGACGGTTCTCTGGCGGCTCAATCAAGTTGACAGGGAATTTTTaac  >  2:1001961/1‑79 (MQ=255)
                                                                          |                                                                            
GCAGTAACCCGCTCAATCGCTTTTTGATTCGTGTTTATCATCCGCTGTTGCTGAAAGTACTGCACTGGCCGAAAACCACGCTGCTGGTGGCGGCGCTTTCGGTGCTGACGGTTCTCTGGCCGCTCAATAAAGTTGGCGGGGAATTTTTACC  >  NC_000913/600251‑600401

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