Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 596,033 C→A 9.0% T404N (ACT→AAT)  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,0330CA 9.0% 208.9 / 5.0 66T404N (ACT→AAT) cusCcopper/silver efflux system, outer membrane component
Reads supporting (aligned to +/- strand):  ref base C (31/29);  new base A (2/4);  total (33/33)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

TGGCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACTTTGCAACGGGCGCGGGCATTATATCAGCACGGCGCAGTAAGTTATCTGGAAGTGCTGGATGCCGAGCG  >  AP009048/595964‑596102
                                                                     |                                                                     
tGGCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACt                                                                      >  1:3514905/1‑71 (MQ=255)
tGGCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACt                                                                      >  1:3358192/1‑71 (MQ=255)
tGGCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACt                                                                      >  1:4009883/1‑71 (MQ=255)
cggCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACt                                                                      <  1:1901776/70‑1 (MQ=255)
  gCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACttt                                                                    >  1:6352642/1‑71 (MQ=255)
  gCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACttt                                                                    >  1:2072698/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:7173435/71‑1 (MQ=255)
               gCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACTTTGCAACGGGCGCgg                                                       <  1:3318244/71‑1 (MQ=255)
               gCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACTTTGCAACGGGCGCgg                                                       <  1:6552179/71‑1 (MQ=255)
                 cTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGTTGCAAAGTAATTTGCAGCga                                                             >  1:6915771/1‑62 (MQ=25)
                 cTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGTTGCAAAGTAATTTGCAGCga                                                             >  1:3535287/1‑62 (MQ=25)
                    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: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)
                       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:4065360/70‑1 (MQ=255)
                                                 tCTGGCATCGCTGCAAATTACTTTGCAACGGGCGCGGGCATTATATCAGCACGGCGCAGTAAGTTATCTg                      >  1:5587244/1‑70 (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)
                                                                     |                                                                     
TGGCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACTTTGCAACGGGCGCGGGCATTATATCAGCACGGCGCAGTAAGTTATCTGGAAGTGCTGGATGCCGAGCG  >  AP009048/595964‑596102

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: