Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 601,304 T→C 7.8% L864P (CTA→CCA)  cusA → copper/silver efflux system, membrane component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913601,3040TC 7.8% 148.2 / 3.3 51L864P (CTA→CCA) cusAcopper/silver efflux system, membrane component
Reads supporting (aligned to +/- strand):  ref base T (22/25);  new base C (2/2);  total (24/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGAGAATGCGCGCCCGACGAGCTGGATTTATATCGATGCCCGCGATCGTGACATGGTGTCGGTGGTTCACGATTTGCAAAAAGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCA  >  NC_000913/601158‑601447
                                                                                                                                                  |                                                                                                                                               
cGAGAATGCGCGCCCGACGAGCTGGATTTATATCGATGCGCGGGATCGTGACATGGTGTCGGTGGTTCACGATTTGCAAAAAGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTAc                                                                                                                                               <  1:30474/149‑1 (MQ=255)
  agaATGCGCGCCCGACGAGCTGGATTTATATCGATGCCCGCGATCGTGACATGGTGTCGGTGGTTCACGATTTGCAAAAAGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTg                                                                                                                                             >  1:288201/1‑149 (MQ=255)
  agaATGCCCGCCCGCCGAGCGGGATGTATATCGAGGCCCGCGATCGTGACAGCGGGGCGGGGGGGCACGATTTGCAAAAAGCGATAGCTGAAAAAGTGCAGTTAAACCCGGGCACCAGCGTCGCATTCTCCGGGCAGTTCGAGCTACTg                                                                                                                                             <  3:13566/149‑1 (MQ=255)
        cgcgcCCGACGAGCTGGATTTATATCGATGCCCGCGATCGTGACATGGTGTCGGTGGTTCACGAATTGCAAAAAGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAgcgc                                                                                                                                       <  4:267259/149‑1 (MQ=255)
           gcCCGACGAGCTGGATTTATATCGATGCCCGCGATCGTGACATGGTGTCGGTGGTTCACGATTTGCAAAAAGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGcc                                                                                                                                    <  8:59430/149‑1 (MQ=255)
              cgacgaGCTGGATTTATATCGATGCCCGCGATCGTGACATGGTGTCGGTGGTTCACGATTTGCAAAAAGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAAc                                                                                                                                 <  5:43278/149‑1 (MQ=255)
               gacgaGCTGGATTTATATCGATGCCCGCGATCGTGACATGGTGTCGGTGGTTCACGATTTGCAAAAAGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAAcc                                                                                                                                <  1:167180/149‑1 (MQ=255)
                     cTGGATTTATATCGATGCCCGCGATCGTGACATGGTGTCGGTGGTTCACGATTTGCAAAAAGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGc                                                                                                                          >  2:145949/1‑149 (MQ=255)
                       ctactaatatCGATGCCCGCGATCGTGACATGGTGTCGGTGGTTCACGATTTGCAAAAAGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCtt                                                                                                                        <  4:28996/143‑1 (MQ=255)
                        gATTTATATCGATGCCCGCGATCGTGACATGGTGTCGGTGGTTCACGATTTGCAAAAAGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCtt                                                                                                                        <  4:93966/148‑1 (MQ=255)
                          tttATATCGATGCCCGCGATCGTGACATGGTGTCGGTGGTTCACGATTTGCAAAAAGCGATAGCTGAAAAAGTGCAGTTAAACCCGGGCACCAGCGGGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTaaa                                                                                                                     <  2:330979/149‑1 (MQ=255)
                             atatCGATGCCCGCGATCGTGACATGGTGTCGGTGGTTCACGATTTGCAAAAAGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACt                                                                                                                   >  3:8206/1‑148 (MQ=255)
                              tatCGATGCCCGCGATCGTGACATGGTGTCGGTGGTTCACGATTTGCAAAAAGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCa                                                                                                                 <  2:288201/149‑1 (MQ=255)
                                tCGATGCCCGCGATCGTGACATGGTGTCGGTGGTTCACGATTTGCAAAAAGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATg                                                                                                               <  1:131951/149‑1 (MQ=255)
                                 cGATGCCCGCGATCGTGACATGGTGTCGGTGGTTCACGATTTGCAAAAAGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCCCCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATgg                                                                                                              >  2:210062/1‑149 (MQ=255)
                                  gATGCCCGCGATCGTGACATGGTGTCGGTGGTTCACGATTTGCAAAAAGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGt                                                                                                             <  4:8206/149‑1 (MQ=255)
                                      cccGCGATCGTGACATGGTGTCGGTGGTTCACGATTTGCAAAAAGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGcc                                                                                                          >  6:224600/1‑148 (MQ=255)
                                       ccGCGATCGTGACATGGTGTCGGTGGTTCACGATTTGCAAAAAGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGa                                                                                                        >  7:145329/1‑149 (MQ=255)
                                                      ggTGTCGGTGGTTCACGATTTGCAAAAAGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCCACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTtg                                                                                                <  4:39552/142‑1 (MQ=255)
                                                      ggTGTCGGTGGTTCACGATTTGCAAAAAGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCCACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTtg                                                                                                <  4:35108/142‑1 (MQ=255)
                                                      ggTGTCGGTGGTTCACGATTTGCAAAAAGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCCACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTtg                                                                                                >  3:39552/1‑142 (MQ=255)
                                                      ggTGTCGGTGGTTCACGATTTGCAAAAAGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCCACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTtg                                                                                                >  3:35108/1‑142 (MQ=255)
                                                         gtCGGTGGTTCACGATTTGCAAAAAGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCt                                                                                      >  6:74354/1‑149 (MQ=255)
                                                         gtCGGTGGTTCACGATTTGCAAAAAGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCt                                                                                      >  1:150009/1‑149 (MQ=255)
                                                             gtggtTCACGATTTGCAAAAAGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTAATGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATAATCTTCGt                                                                                  >  6:59140/1‑149 (MQ=255)
                                                              tggtTCACGATTTGCAAAAAGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTg                                                                                 <  1:145949/149‑1 (MQ=255)
                                                                gtTCACGATTTGCAAAAAGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCt                                                                               <  6:263440/149‑1 (MQ=255)
                                                                  tCACGATTTGCAAAAAGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCtgt                                                                             >  7:15082/1‑149 (MQ=255)
                                                                    aCGATTTGCAAAAAGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTGTTGGTGCCGATGACGTTGATGATTATCTTCGTGCtgtt                                                                            <  8:12539/148‑1 (MQ=255)
                                                                     cGATTTGCAAAAAGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCtgttg                                                                           >  8:264/1‑148 (MQ=255)
                                                                               aaaGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGtt                                                                <  8:199352/149‑1 (MQ=255)
                                                                                        gCTGAAACAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTGCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTTCTGTTGTATCTGGCGTTCCGTCGGGt                                                       <  1:291075/149‑1 (MQ=255)
                                                                                               aaGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGGGGGCGaa                                                <  5:224600/149‑1 (MQ=255)
                                                                                                 gTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTcc                                                              >  1:54740/1‑133 (MQ=255)
                                                                                                 gTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTcc                                                              <  2:54740/133‑1 (MQ=255)
                                                                                                 gTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGc                                              >  2:356656/1‑149 (MQ=255)
                                                                                                   gCAGTTAAAACAGTGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTGAAAATAATGGTGCCGATGACGTTGATGATTAGCTGCGTGGGTTTGTATCTGGCGTTGCTTCGTGTGGGCGAAGCGt                                            >  4:4608/1‑149 (MQ=255)
                                                                                                    cAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGtt                                           <  3:233569/149‑1 (MQ=255)
                                                                                                    cAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGtt                                           >  2:72406/1‑149 (MQ=255)
                                                                                                    cAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGt                                            <  8:15082/148‑1 (MQ=255)
                                                                                                     aGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTg                                          >  2:7028/1‑149 (MQ=255)
                                                                                                       ttAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCt                                        <  8:273624/149‑1 (MQ=255)
                                                                                                                           ggCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGt                    >  7:238301/1‑149 (MQ=255)
                                                                                                                             cATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGttt                  <  5:74354/149‑1 (MQ=255)
                                                                                                                             cATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGttt                  <  7:264/149‑1 (MQ=255)
                                                                                                                              aTTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTAt                                                                        <  6:197634/94‑1 (MQ=255)
                                                                                                                              aTTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTAt                                                                        >  5:197634/1‑94 (MQ=255)
                                                                                                                                         aGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGggc      >  3:193099/1‑149 (MQ=255)
                                                                                                                                           ttCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGggcgg    >  5:8867/1‑149 (MQ=255)
                                                                                                                                            tCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGggcgg    >  1:156156/1‑148 (MQ=255)
                                                                                                                                             cGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCGGTTGTATCTGGCGTTGCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCGGGTGGGCGGCa  >  6:188740/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
CGAGAATGCGCGCCCGACGAGCTGGATTTATATCGATGCCCGCGATCGTGACATGGTGTCGGTGGTTCACGATTTGCAAAAAGCGATAGCTGAAAAAGTGCAGTTAAAACCGGGCACCAGCGTGGCATTCTCCGGGCAGTTCGAGCTACTGGAGCGCGCCAACCATAAGCTTAAACTCATGGTGCCGATGACGTTGATGATTATCTTCGTGCTGTTGTATCTGGCGTTCCGTCGGGTGGGCGAAGCGTTGCTGATTATCAGCAGCGTACCGTTTGCGCTGGTGGGCGGCA  >  NC_000913/601158‑601447

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑