Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 596,030 T→G 10.2% I403S (ATT→AGT)  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,0300TG10.2% 138.9 / 12.9 59I403S (ATT→AGT) cusCcopper/silver efflux system, outer membrane component
Reads supporting (aligned to +/- strand):  ref base T (25/28);  new base G (2/4);  total (27/32)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.71e-01

GCGCTGGCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACTTTGCAACGGGCGCGGGCATTATATCAGCACGGCGCAGTAAGTTATCTGGAAGTGCTGGATGCCGAG  >  AP009048/595960‑596100
                                                                      |                                                                      
gcgcTGGCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAAt                                                                        >  1:6550127/1‑71 (MQ=255)
 cgcTGGCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAAtt                                                                       <  1:3874655/71‑1 (MQ=255)
 cgcTGGCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAAt                                                                        <  1:3598764/70‑1 (MQ=255)
 cgcTGGCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAAt                                                                        >  1:3448405/1‑70 (MQ=255)
 cgcTGGCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAAt                                                                        <  1:2527866/70‑1 (MQ=255)
 cgcTGGCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAAt                                                                        >  1:63098/1‑70 (MQ=255)
    tGGCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACt                                                                    >  1:4009883/1‑71 (MQ=255)
    tGGCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACt                                                                    >  1:3514905/1‑71 (MQ=255)
    tGGCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACt                                                                    >  1:3358192/1‑71 (MQ=255)
    cggCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACt                                                                    <  1:1901776/70‑1 (MQ=255)
      gCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACttt                                                                  >  1:2072698/1‑71 (MQ=255)
      gCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACttt                                                                  >  1:6352642/1‑71 (MQ=255)
      gCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACtt                                                                   <  1:5299000/70‑1 (MQ=255)
      gCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACtt                                                                   <  1:2284103/70‑1 (MQ=255)
      gCATTACGTCAAAGCCTGAACGATCAAACTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACttt                                                                  >  1:2576518/1‑71 (MQ=255)
       cATTACGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACTTTg                                                                 >  1:3980598/1‑71 (MQ=255)
         ttACGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACTTTGCa                                                               <  1:3864640/71‑1 (MQ=255)
           aCGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACTTTGCaa                                                              >  1:3958934/1‑70 (MQ=255)
               cAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACTTTGCAACGGgc                                                         >  1:5547358/1‑71 (MQ=255)
                aaaGCCTGAACGATCAAATTAGCGCCAAGCAGCGTTATCTGGCATCGCTGCAAATTACTTTGCAACGGgcg                                                        >  1:6409482/1‑71 (MQ=255)
                   gCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACTTTGCAACGGGCGCgg                                                     <  1:3318244/71‑1 (MQ=255)
                   gCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACTTTGCAACGGGCGCgg                                                     <  1:6552179/71‑1 (MQ=255)
                   gCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACTTTGCAACGGGCGCgg                                                     <  1:7173435/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:7609151/56‑1 (MQ=25)
                           gATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGTTGCAAAGTAATTTGCAGCg                                                            <  1:874019/56‑1 (MQ=25)
                           gATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGTTGCAAAGTAATTTGCAGCg                                                            <  1:3196670/56‑1 (MQ=25)
                           gATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGTTGCAAAGTAATTTGCAGCg                                                            <  1:3306117/56‑1 (MQ=25)
                           gATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACt                                                                    <  1:4879591/48‑1 (MQ=255)
                           gATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACt                                                                    <  1:2518660/48‑1 (MQ=255)
                            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)
                               aaaTTAGCGCCCACCAGCGTTATCTGGCATCGCTGCAAATTACTTTg                                                                 >  1:4170647/1‑47 (MQ=255)
                                aaTTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACttt                                                                  >  1:7533301/1‑45 (MQ=255)
                                aaTTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACttt                                                                  >  1:3406035/1‑45 (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:5482397/67‑1 (MQ=255)
                                                         gCATCGCTGCAAATTACTTTGCAACGGGCGCGGGCATTATATCAGCACGGCGCAGTAAGTTATCTgg                   <  1:5188235/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)
                                                                      |                                                                      
GCGCTGGCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACTTTGCAACGGGCGCGGGCATTATATCAGCACGGCGCAGTAAGTTATCTGGAAGTGCTGGATGCCGAG  >  AP009048/595960‑596100

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: