Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 596,040 A→G 11.3% Q406Q (CAA→CAG cusC → copper/silver efflux system, outer membrane component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP009048596,0400AG11.3% 138.4 / 9.5 53Q406Q (CAA→CAGcusCcopper/silver efflux system, outer membrane component
Reads supporting (aligned to +/- strand):  ref base A (21/26);  new base G (2/4);  total (23/30)
Fisher's exact test for biased strand distribution p-value = 6.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.41e-01

ACGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACTTTGCAACGGGCGCGGGCATTATATCAGCACGGCGCAGTAAGTTATCTGGAAGTGCTGGATGCCGAGCGTTCTT  >  AP009048/595971‑596107
                                                                     |                                                                   
aCGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACTTTGCaa                                                                     >  1:3958934/1‑70 (MQ=255)
    cAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACTTTGCAACGGgc                                                                >  1:5547358/1‑71 (MQ=255)
     aaaGCCTGAACGATCAAATTAGCGCCAAGCAGCGTTATCTGGCATCGCTGCAAATTACTTTGCAACGGgcg                                                               >  1:6409482/1‑71 (MQ=255)
        gCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACTTTGCAACGGGCGCgg                                                            <  1:6552179/71‑1 (MQ=255)
        gCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACTTTGCAACGGGCGCgg                                                            <  1:7173435/71‑1 (MQ=255)
        gCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACTTTGCAACGGGCGCgg                                                            <  1:3318244/71‑1 (MQ=255)
          cTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGTTGCAAAGTAATTTGCAGCga                                                                  >  1:3535287/1‑62 (MQ=25)
          cTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGTTGCAAAGTAATTTGCAGCga                                                                  >  1:6915771/1‑62 (MQ=25)
             aaCGATCAAATTAGCGCCCAGCAGCGTTATCTTGCATCGCTGCAAATTACTTTGCAACGGGCGCGGGCAtt                                                       <  1:6683464/71‑1 (MQ=255)
             aaCGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACTTTGCAACGGGCGCGGGCAtt                                                       <  1:6499853/71‑1 (MQ=255)
              aCGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACTTTGCAACGGGCGCGGGCATta                                                      <  1:6392013/71‑1 (MQ=255)
              aCGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACTTTGCAACGGGCGCGGGCATta                                                      <  1:6865748/71‑1 (MQ=255)
               cGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACTTTGCAACGGGCGCGGGCATtat                                                     <  1:5202587/71‑1 (MQ=255)
                gATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGTTGCAAAGTAATTTGCAGCg                                                                   <  1:874019/56‑1 (MQ=25)
                gATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGTTGCAAAGTAATTTGCAGCg                                                                   <  1:3306117/56‑1 (MQ=25)
                gATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGTTGCAAAGTAATTTGCAGCg                                                                   <  1:3196670/56‑1 (MQ=25)
                gATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGTTGCAAAGTAATTTGCAGCg                                                                   <  1:7609151/56‑1 (MQ=25)
                 aTCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACTTTGCAACGGGCGCGGGCATTatat                                                   <  1:1304688/71‑1 (MQ=255)
                  tCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACTTTGCAACGGGCGCGGGCATTATATc                                                  <  1:3960053/71‑1 (MQ=255)
                   cAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACTTTGCAACGGGCGCGGGCATTatat                                                   >  1:4151418/1‑69 (MQ=255)
                    aaaTTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACTTTGCAACGGGCGCGGGCATTATATCAg                                                >  1:6068348/1‑71 (MQ=255)
                               cagcagCGTTATCTGGCATCGCTGCAAATTACTTTGCAACGGGCGCGGGCATTatat                                                   <  1:7731789/57‑1 (MQ=255)
                                  cagcGTTATCTGGCATCGCTGCAAATTACTTTGCAACGGGCGCGGGCATTatat                                                   <  1:3659883/54‑1 (MQ=255)
                                  cagcGTTATCTGGCATCGCTGCAAATTACTTTGCAACGGGCGCGGGCATTATATCAGCACGGCGCAGTAAg                                  >  1:2951240/1‑71 (MQ=255)
                                    gcGTTATCTGGCATCGCTGCAAATTACTTTGCAACGGGCGCGGGCATTATATCAGCACGGCGCAGTAAGtt                                <  1:4474330/71‑1 (MQ=255)
                                         aTCTGGCATCGCTGCAAATTACTTTGCAACGGGCGCGGGCATTATATCAGCACGGCGCAGTAAGTTATCt                            <  1:1023197/70‑1 (MQ=255)
                                         aTCTGGCATCGCTGCAAATTACTTTGCAACGGGCGCGGGCATTATATCAGCACGGCGCAGTAAGTTATCTg                           >  1:2620666/1‑71 (MQ=255)
                                         aTCTGGCATCGCTGCAAATTACTTTGCAACGGGCGCGGGCATTATATCAGCACGGCGCAGTAAGTTATCTg                           >  1:1567690/1‑71 (MQ=255)
                                         aTCTGGCATCGCTGCAAATTACTTTGCAACGGGCGCGGGCATTATATCAGCAGGCGCAGTAAGTTATCTg                           >  1:5135076/1‑70 (MQ=255)
                                          tCTGGCATCGCTGCAAATTACTTTGCAACGGGCGCGGGCATTATATCAGCACGGCGCAGTAAGTTATCTg                           <  1:2436064/70‑1 (MQ=255)
                                          tCTGGCATCGCTGCAAATTACTTTGCAACGGGCGCGGGCATTATATCAGCACGGCGCAGTAAGTTATCTg                           >  1:5587244/1‑70 (MQ=255)
                                          tCTGGCATCGCTGCAAATTACTTTGCAACGGGCGCGGGCATTATATCAGCACGGCGCAGTAAGTTATCTg                           <  1:4065360/70‑1 (MQ=255)
                                            tGGCATCGCTGCAAATTACTTTGCAACGGGCGCGGGCATTATATCAGCACGgcg                                         >  1:4735622/1‑54 (MQ=255)
                                              gCATCGCTGCAAATTACTTTGCAACGGGCGCGGGCATTATATCAGCACGgcg                                         >  1:194468/1‑52 (MQ=255)
                                              gCATCGCTGCAAATTACTTTGCAACGGGCGCGGGCATTATATCAGCACGGCGCAGTAAGTTATCTgg                          <  1:5188235/67‑1 (MQ=255)
                                              gCATCGCTGCAAATTACTTTGCAACGGGCGCGGGCATTATATCAGCACGGCGCAGTAAGTTATCTgg                          <  1:5482397/67‑1 (MQ=255)
                                                aTCGCTGCAAATTACTTTGCAACGGGCGCGGGCATTATATCAGCACGGCGCAGTAAGTTATCTGGAAGTg                     >  1:1875117/1‑70 (MQ=255)
                                                aTCGCTGCAAATCACTTTGCAACGGGCGCGGGCATTATATCAGCACGGCGCAGTAAGTTATCTGGAAGTGc                    <  1:5760480/71‑1 (MQ=255)
                                                 tCGCTGCAAATTACTTTGCAACGGGCGCGGGCATTATATCAGCACGGCGCAGTAAGTTATCTgg                          <  1:1342302/64‑1 (MQ=255)
                                                 tCGCTGCAAATTACTTTGCAACGGGCGCGGGCATTATATCAGCACGGCGCAGTAAGTTATCTGGAAGTGCt                   >  1:1823231/1‑71 (MQ=255)
                                                       cAAATTACTTTGCAACGGGCGCGGGCATTATATCAGCACGGCGCAGTAAGTTATCTGGAAGTGCTGGAt               <  1:3941568/69‑1 (MQ=255)
                                                        aaaTTACTTTGCAACGGGCGCGGGCATTATATCAGCACGGCGCAGTAAGt                                 >  1:4042549/1‑50 (MQ=255)
                                                        aaaTTACTTTGCAACGGGCGCGGGCATTATATCAGCACGGCGCAGTAAGTTATCTGGAAGTGCTGGATGcc            <  1:5316797/71‑1 (MQ=255)
                                                           ttACTTTGCAACGGGCGCGGGCATTATATCAGCACGGCGCAGTAAGTTATCTGGAAGTGCTGGATGc             <  1:3237690/67‑1 (MQ=255)
                                                           ttACTTTGCAACGGGCGCGGGCATTATATCAGCACGGCGCAGTAAGTTATCTGGAAGTGCTGGATGCCGAg         >  1:4797560/1‑71 (MQ=255)
                                                            tACTTTGCAAGGGGCGCGGGCATTATATCAGCACGGCGCAGTAAGt                                 >  1:4294887/1‑46 (MQ=255)
                                                            tACTTTGCAACGGGCGCGGGCATTATATCAGCACGGCGCAGTAAGTTATCTGGaa                        >  1:5499457/1‑55 (MQ=255)
                                                            tACTTTGCAACGGGCGCGCGCATTATATCAGCACGGCGCAGTAAGt                                 >  1:6543682/1‑46 (MQ=255)
                                                             aCTTTGCAACGGGCGCGGGCATTATATCAGCACGGCGCAGTAAGTTATCt                            >  1:196434/1‑50 (MQ=255)
                                                             aCTTTGCAACGGGCGCGGGCATTATATCAGCACGGCGCAGTAAGTTATCt                            >  1:1363338/1‑50 (MQ=255)
                                                              cTTTGCAACGGGCGCGGGCATTATATCAGCACGGCGCAGTAAGTTATCTGGAAGTGCTGGATGCCGAGCg       <  1:2146212/70‑1 (MQ=255)
                                                              cTTTGCAACGGGCGCGGGCATTATATCAGCACGGCGCAGTAAGTTATCTGGAAGTGCTGGATGCCGAGCg       <  1:1263163/70‑1 (MQ=255)
                                                                   cAACGGGCGCGGGCATTATATCAGCACGGCGCAGTAAGTTATCTGGAAGTGCTGGATGCCGAGCGTTCtt  <  1:5583730/70‑1 (MQ=255)
                                                                   cAACGGGCGCGGGCATTATATCAGCACGGCGCAGTAAGTTATCTGGAAGTGCTGGATGCCGAGCGTTCtt  >  1:2487493/1‑70 (MQ=255)
                                                                     |                                                                   
ACGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACTTTGCAACGGGCGCGGGCATTATATCAGCACGGCGCAGTAAGTTATCTGGAAGTGCTGGATGCCGAGCGTTCTT  >  AP009048/595971‑596107

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