Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 596,023 C→T 9.1% L401L (CTG→TTG)  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,0230CT 9.1% 202.2 / 7.5 66L401L (CTG→TTG) cusCcopper/silver efflux system, outer membrane component
Reads supporting (aligned to +/- strand):  ref base C (28/32);  new base T (2/4);  total (30/36)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

GCAGATGCGCTGGCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACTTTGCAACGGGCGCGGGCATTATATCAGCACGGCGCAGTAAGTTATCTGGAAGTGCT  >  AP009048/595954‑596090
                                                                     |                                                                   
gCAGATGCGCTGGCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGc                                                                     <  1:6396504/70‑1 (MQ=255)
   gATGCGCTGGCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCa                                                                 <  1:5493046/71‑1 (MQ=255)
    aTGCGCTGGCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCt                                                                    <  1:1535598/67‑1 (MQ=255)
     tGCGCTGGCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCAAGCAGCGTTATCTGGCATCGCTGCaaa                                                               >  1:2748721/1‑71 (MQ=255)
      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:63098/1‑70 (MQ=255)
       cgcTGGCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAAt                                                              <  1:2527866/70‑1 (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:2284103/70‑1 (MQ=255)
            gCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACtt                                                         <  1:5299000/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:7173435/71‑1 (MQ=255)
                         gCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACTTTGCAACGGGCGCgg                                           <  1:6552179/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)
                             gAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCa                                                                 <  1:2625994/45‑1 (MQ=255)
                              aaCGATCAAATTAGCGCCCAGCAGCGTTATCTTGCATCGCTGCAAATTACTTTGCAACGGGCGCGGGCAtt                                      <  1:6683464/71‑1 (MQ=255)
                              aaCGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACTTTGCAACGGGCGCGGGCAtt                                      <  1:6499853/71‑1 (MQ=255)
                               aCGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACTTTGCAACGGGCGCGGGCATta                                     <  1:6865748/71‑1 (MQ=255)
                               aCGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACTTTGCAACGGGCGCGGGCATta                                     <  1:6392013/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:5587244/1‑70 (MQ=255)
                                                           tCTGGCATCGCTGCAAATTACTTTGCAACGGGCGCGGGCATTATATCAGCACGGCGCAGTAAGTTATCTg          <  1:4065360/70‑1 (MQ=255)
                                                           tCTGGCATCGCTGCAAATTACTTTGCAACGGGCGCGGGCATTATATCAGCACGGCGCAGTAAGTTATCTg          <  1:2436064/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)
                                                                     |                                                                   
GCAGATGCGCTGGCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACTTTGCAACGGGCGCGGGCATTATATCAGCACGGCGCAGTAAGTTATCTGGAAGTGCT  >  AP009048/595954‑596090

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