Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 4,184,570 G→T 7.6% S90Y (TCT→TAT)  gspC ← general secretory pathway component, cryptic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090484,184,5700GT 7.6% 344.8 / 7.8 105S90Y (TCT→TAT) gspCgeneral secretory pathway component, cryptic
Reads supporting (aligned to +/- strand):  ref base G (39/58);  new base T (6/2);  total (45/60)
Fisher's exact test for biased strand distribution p-value = 7.14e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01

Only 100 of 120 total aligned reads displayed.

CCCCTCATTAATTAAGACAAAAGAAGTCTCATCATTCGATGTGAGTACAATGCCATTAAGCGCTATAGAGAGCTTAGGTGTTTCCTGATTAACTGCCACGGGATTTGCATCTTTTACCCCAGCATGATTTTCGTTA  >  AP009048/4184503‑4184638
                                                                   |                                                                    
ccccTCATTAATTAAGACAAAAGAAGTCTCATCATTCGATGTGAGTACAATGCCATTAAGCGCTATaga                                                                     >  1:4524401/1‑69 (MQ=255)
  ccTCATTAATTAAGACAAAAGAAGTCTCATCATTCGATGTGAGTACAATGCCATTAAGCGCTATagagag                                                                  <  1:5987243/70‑1 (MQ=255)
   cTCATTAATTAAGACAAAAGAAGTCTCATCATTCGATGTGAGTACAATGCCATTAAGCTCTCTATAGCGCt                                                                >  1:5817377/1‑71 (MQ=255)
      attaattaAGACAAAAGAAGTCTCATCATTCGATGTGAGTACAATGCCATTAAGCTCTCTATAGCGCTTa                                                              <  1:7207890/70‑1 (MQ=38)
      attaattaAGACAAAAGAAGTCTCATCATTCGATGTGAGTACAATGCCATTAAGCGCTATAGAGAGCTTAg                                                             >  1:462704/1‑71 (MQ=255)
      attaattaAGACAAAAGAAGTCTCATCATTCGATGTGAGTACAATGCCATTAAGCGCTATAGAGAGCTTAg                                                             >  1:6971277/1‑71 (MQ=255)
       ttaattaaGACAAAAGAAGTCTCATCATTCGATGTGAGTACAATGCCATTAAGCTCTCTATAGCGCTTaat                                                            >  1:5572189/1‑69 (MQ=37)
           ttaaGACAAAAGAAGTCTCATCATTCGATGTGAGTACAATGCCATTAAGCGCTATAGAGAGCTTAGGTGtt                                                        >  1:3293153/1‑71 (MQ=255)
            taaGACAAAAGAAGTCTCATCATTCGATGTGAGTACAATGCCATTAAGCGCTATaga                                                                     >  1:5927617/1‑57 (MQ=255)
            taaGACAAAAGAAGTCTCATCATTCGATGTGAGTACAATGCCATTAAGCGCTATaga                                                                     >  1:129124/1‑57 (MQ=255)
             aaGACAAAAGAAGTCTCATCATTCGATGTGAGTACAATGCCATTAAGCGCTATAGAGAGCt                                                                <  1:3372984/61‑1 (MQ=255)
             aaGACAAAAGAAGTCTCATCATTCGATGTGAGTACAATGCCATTAAGCGCTATAGAGAGCt                                                                <  1:3476300/61‑1 (MQ=255)
             aaGACAAAAGAAGTCTCATCATTCGATGTGAGTACAATGCCATTAAGCGCTATAGAGAGCTTAGGTGTTTc                                                      >  1:1247002/1‑71 (MQ=255)
             aaGACAAAAGAAGTCTCATCATTCGATGTGAGTACAATGCCATTAAGCGCTATAGAGAGCTTAGGTGTTTc                                                      >  1:4790044/1‑71 (MQ=255)
             aaGACAAAAGAAGTCTCATCATTCGATGAGAGTACAATGCCATTAAGCGCTATAGAGAGCTTAGGTGTTTc                                                      >  1:7721733/1‑71 (MQ=255)
                  aaaaGAAGTCTCATCATTCGATGTGAGTACAATGCCATTAAGCGCTATAGAGAGCTTAGGTGTTTCCTGAt                                                 >  1:2073131/1‑71 (MQ=255)
                   aaaGAAGTCTCATCATTCGATGTGAGTACAATGCCATTAAGCGCTATAGAGAGCTTAGGTGTTTCCTGAtt                                                <  1:3706105/71‑1 (MQ=255)
                   aaaGAAGTCTCATCATTCGATGTGAGTACAATGCCATTAAGCGCTATAGAGAGCTTAGGTGTTTCCTGAtt                                                <  1:5026325/71‑1 (MQ=255)
                    aagaagTCTCATCATTCGATGTGAGTACAATGCCATTAAGCTCTCTatag                                                                    >  1:1597508/1‑50 (MQ=255)
                    aagaagTCTCATCATTCGATGTGAGTACAATGCCATTAAGCTCTCTatag                                                                    >  1:7857770/1‑50 (MQ=255)
                     agaagTCTCATCATTCGATGTGAGTACAATGCCATTAAGCGCTATATAGAGCTTAGGTGTTTCCTGATTaa                                              <  1:1871918/71‑1 (MQ=255)
                     agaagTCTCATCATTCGATGTGAGTACAATGCCATTAAGCGCTATAGAGAGCTTa                                                              <  1:7916754/55‑1 (MQ=255)
                     agaagTCTCATCATTCGATGTGAGTACAATGCCATTAAGCGCTATAGAGAGCTTAGGTGTTTCCTGATTaa                                              <  1:4137830/71‑1 (MQ=255)
                     agaagTCTCATCATTCGATGTGAGTACAATGCCATTAAGCGCTATAGAGAGCTTAGGTGTTTCCTGATTaa                                              <  1:3022148/71‑1 (MQ=255)
                     agaagTCTCATCATTCGATGTGAGTACAATGCCATTAAGCGCTATAGAGAGCTTAGGTGTTTCCTGATTaa                                              <  1:2671956/71‑1 (MQ=255)
                      gaagTCTCATCATTCGATGTGAGTACAATGCCATTAAGCGCTATAGAGAGCTTAGGTGTTTCCTGATTAAc                                             <  1:2177471/71‑1 (MQ=255)
                       aagTCTCATGATTCGATGTGAGTACAATGCCATTAAGCGCTATAGAGAGCTTAGGTGTTTCCTGATTAACt                                            <  1:5583356/71‑1 (MQ=255)
                       aagTCTCATCATTCGATGTGAGTACAATGCCATTAAGCGCTATAGAGAGCTTAGGTGTTTCCTGATTAACt                                            <  1:3134777/71‑1 (MQ=255)
                       aagTCTCATCATTCGATGTGAGTACAATGCCATTAAGCGCTATAGAGAGCTTAGGTGTTTCCTGATTAACt                                            <  1:1275916/71‑1 (MQ=255)
                        agTCTCATCATTCGATGTGAGTACAATGCCATTAAGCGCTATAGAGAGCTTAGGTGTTTCCTGATTAACt                                            <  1:2316132/70‑1 (MQ=255)
                         gTCTCATCATTCGATGTGAGTACAATGCCATTAAGCGCTATAGAGAGCTTAgg                                                            <  1:6643235/53‑1 (MQ=255)
                         gTCTCATCATTCGATGTGAGTACAATGCCATTAAGCGCTATAGAGAGCTTAgg                                                            <  1:3715444/53‑1 (MQ=255)
                         gTCTCATCATTCGATGTGAGTACAATGCCATTAAGCGCTATAGAGAGCTTAgg                                                            <  1:1315421/53‑1 (MQ=255)
                         gTCTCATCATTCGATGTGAGTACAATGCCATTAAGCGCTATAGAGAGCTTAgg                                                            <  1:3495242/53‑1 (MQ=255)
                          tctcATCATTCGATGTGAGTACAATGCCATTAAGCGCTATAGAGAGCTTAgg                                                            <  1:1717440/52‑1 (MQ=255)
                          tctcATCATTCGATGTGAGTACAATGCCATTAAGCGCTATAGAGAGCTTAGGTGTTTCCTGATTAACTGcc                                         <  1:2497103/71‑1 (MQ=255)
                          tctcATCATTCGATGTGAGTACAATGCCATTAAGCGCTATAGAGAGCTTAGGTGTTTCCTGATTAACTGcc                                         <  1:6440155/71‑1 (MQ=255)
                           ctcATCATTCGATGTGAGTACAATGCCATTAAGCGCTATAGAGAGCTTAGGTGTTTCCTGATTAACTGCCa                                        >  1:7680516/1‑71 (MQ=255)
                            tcatcaTTCGATGTGAGTACAATGCCATTAAGCTCTCTATAGCGCTTaat                                                            >  1:2158074/1‑48 (MQ=21)
                            tcatcaTTCGATGTGAGTACAATGCCATTAAGCTCTCTATAGCGCTTaat                                                            >  1:5059164/1‑48 (MQ=21)
                            tcatcaTTCGATGTGAGTACAATGCCATTAAGCTCTCTATAGCGCTTaat                                                            >  1:5100275/1‑48 (MQ=21)
                            tcatcaTTCGATGTGAGTACAATGCCATTAAGCTCTCTATAGCGCTTaat                                                            >  1:2190831/1‑48 (MQ=21)
                            tcatcaTTCGATGTGAGTACAATGCCATTAAGCGCTATagaga                                                                   >  1:5120851/1‑43 (MQ=255)
                              atcatTCGATGTGAGTACAATGCCATTAAGCGCTATAGAGAGCTTAGGTGttt                                                       >  1:7896899/1‑53 (MQ=255)
                              atcatTCGATGTGAGTACAATGCCATTAAGCGCTATAGAGAGCTTAGGTGTTTCCTGAt                                                 >  1:4197705/1‑59 (MQ=255)
                              atcatTCGATGTGAGTACAATGCCATTAAGCGCTATAGAGAGCTTAGGTGTTTCCTGATTAACTGCCa                                        >  1:5409187/1‑68 (MQ=255)
                              atcatTCGATGTGAGTACAATGCCATTAAGCGCTATAGAGAGCTTAGGTGTTTCCTGATTAACTGCCAc                                       <  1:2660450/69‑1 (MQ=255)
                                catTCGATGTGAGTACAATGCCATTAAGCGCTATGGAGAGCTTAg                                                             <  1:7427049/45‑1 (MQ=255)
                                catTCGATGTGAGTACAATGCCATTAAGCGCTATAGAGAGCTTAg                                                             <  1:4185169/45‑1 (MQ=255)
                                catTCGATGTGAGTACAATGCCATTAAGCGCTATAGAGAGCTTAg                                                             <  1:4116320/45‑1 (MQ=255)
                                catTCGATGTGAGTACAATGCCATTAAGCGCTATAGAGAGCTTAg                                                             <  1:5054669/45‑1 (MQ=255)
                                 atTCGATGTGAGTACAATGCCATTAAGCGCTATAGAGAGCt                                                                >  1:4465862/1‑41 (MQ=255)
                                 atTCGATGTGAGTACAATGCCATTAAGCGCTATAGAGAGCt                                                                >  1:6532666/1‑41 (MQ=255)
                                 atTCGATGTGAGTACAATGCCATTAAGCGCTATAGAGAGCt                                                                >  1:6119892/1‑41 (MQ=255)
                                 atTCGATGTGAGTACAATGCCATTAAGCGCTATAGAGAGCt                                                                >  1:1011996/1‑41 (MQ=255)
                                  ttCGATGTGAGTACAATGCCATTAAGCGCTATAGAGAGCTTAGGTGtt                                                        <  1:2278766/48‑1 (MQ=255)
                                    cGATGTGAGTACAATGCCATTAAGCGCTATAGAGAGCTTAGGTGTTTc                                                      >  1:3961625/1‑48 (MQ=255)
                                      aTGTGAGTACAATGCCATTAAGCGCTATAGAGAGCTTAGGTGTTTCCTGATTAACTGCCACggg                                    <  1:3395715/64‑1 (MQ=255)
                                        gtgAGTACAATGCCATTAAGCGCTATAGAGAGCTTAGGTGTTTCCt                                                    <  1:111545/46‑1 (MQ=255)
                                        gtgAGTACAATGCCATTAAGCGCTATAGAGAGCTTAGGTGTTTCCt                                                    <  1:7215175/46‑1 (MQ=255)
                                        gtgAGTACAATGCCATTAAGCGCTATAGAGAGCTTAGGTGTTTCCTGAtt                                                <  1:4858926/50‑1 (MQ=255)
                                           aGTACAATGCCATTAAGCGCTATAGAGAGCTTAGGTGtt                                                        >  1:4448571/1‑39 (MQ=255)
                                           aGTACAATGCCATTAAGCGCTATAGAGAGCTTAGGTGtt                                                        >  1:6550530/1‑39 (MQ=255)
                                           aGTACAATGCCATTAAGCGCTATAGAGAGCTTAGGTGtt                                                        >  1:1761004/1‑39 (MQ=255)
                                            gTACAATGCCATTAAGCGCTATAGAGAGCTTAGGTGTTTCCTGa                                                  >  1:4830091/1‑44 (MQ=255)
                                            gTACAATGCCATTAAGCGCTATAGAGAGCTTAGGTGTTTCCTGATTAACTGCCACGGGATTTGCATCtttt                       <  1:2445984/71‑1 (MQ=255)
                                             tACAATGCCATTAAGCGCTATAGAGAGCTTAGGTGTTTCCTGa                                                  >  1:230301/1‑43 (MQ=255)
                                             tACAATGCCATTAAGCGCTATAGAGAGCTTAGGTGTTTCCTGATTAACTGCCACGGGATTTGCATCTTTTa                      <  1:5413207/71‑1 (MQ=255)
                                               cAATGCCATTAAGCGCTATAGAGAGCTTAGGTGTTTCCt                                                    >  1:1904573/1‑39 (MQ=255)
                                               cAATGCCATTAAGCGCTATAGAGAGCTTAGGTGTTTCCt                                                    >  1:1839687/1‑39 (MQ=255)
                                               cAATGCCATTAAGCGCTATAGAGAGCTTAGGTGTTTCCTGATTAACTGCCACGGGATTTGCATCTTTTAcc                    <  1:191614/71‑1 (MQ=255)
                                               cAATGCCATTAAGCGCTATAGAGAGCTTAGGTGTTTCCTGATTAACTGCCACGGGATTTGCATCTTTTAc                     <  1:7652246/70‑1 (MQ=255)
                                                aaTGCCATTAAGCGCTATAGAGAGCTTAGGTGTTTCCTGATTAACTGCCAc                                       <  1:6107562/51‑1 (MQ=255)
                                                  tGCCATTAAGCGCTATAGAGAGCTTAGGTGTTTCCTGAtt                                                >  1:4633808/1‑40 (MQ=255)
                                                  tGCCATTAAGCGCTATAGAGAGCTTAGGTGTTTCCTGAt                                                 >  1:2032273/1‑39 (MQ=255)
                                                  tGCCATTAAGCGCTATAGAGAGCTTAGGTGTTTCCTGATTAACTGCCACGGGATTTGCATCTTTTAccc                   >  1:820227/1‑69 (MQ=255)
                                                  tGCCATTAAGCGCTATAGAGAGCTTAGGTGTTTCCTGATTAACTGCCACGGGATTTGCATCTTTTACCCCa                 <  1:7540516/71‑1 (MQ=255)
                                                   gCCATTAAGCGCTATAGAGAGCTTAGGTGTTTCCTGATTaa                                              >  1:631734/1‑41 (MQ=255)
                                                     cATTAAGCGCTATAGAGAGCTTAGGTGTTTCCTGATTAACTGc                                          <  1:2328816/43‑1 (MQ=255)
                                                      aTTAAGCGCTATAGAGAGCTTAGGTGTTTCCTGATTAACTGCCACGGGATTTGCATCTTTTACCCCAGCAt             >  1:2512653/1‑71 (MQ=255)
                                                       ttAAGCGCTATAGAGAGCTTAGGTGTTTCCTGATTAACTGCCACGGGATTTGCATCTTTTAc                     >  1:4114897/1‑62 (MQ=255)
                                                       ttAAGCGCTATAGAGAGCTTAGGTGTTTCCTGATTAACTGCCACGGGATTTGCATCTTTTACCCCAGCATg            <  1:6832301/71‑1 (MQ=255)
                                                       ttAAGCGCTATAGAGAGCTTAGGTGTTTCCTGATTAACTGCCACGGGATTTGCATCTTTTACCCCAGCATg            <  1:3383280/71‑1 (MQ=255)
                                                         aaGCGCTATAGAGAGCTTAGGTGTTTCCTGATTAACTGc                                          >  1:1763152/1‑39 (MQ=255)
                                                         aaGCGCTATAGAGAGCTTAGGTGTTTCCTGATTAACTGc                                          >  1:2479446/1‑39 (MQ=255)
                                                          aGCGCTATAGAGAGCTTAGGTGTTTCCTGATTAACTGc                                          <  1:4946938/38‑1 (MQ=255)
                                                          aGCGCTATAGAGAGCTTAGGTGTTTCCTGATTAACTGc                                          <  1:2615344/38‑1 (MQ=255)
                                                          aGCGCTATAGAGAGCTTAGGTGTTTCCTGATTAACTGc                                          <  1:4773447/38‑1 (MQ=255)
                                                           gcgcTATAGAGAGCTTAGGTGTTTCCTGATTAACTGCCACGGGATTTGCATCTTTTACCCCAGCATg            >  1:1590219/1‑67 (MQ=255)
                                                           gcgcTAAAGAGAGCTTAGGTGTTTCCTGATTAACTGCCACGGGATTTGCATCTTTTACCCCAGCATGAttt        >  1:1020781/1‑71 (MQ=255)
                                                            cgcTATAGAGAGCTTAGGTGTTTCCTGATTAACTGCCACGGGATTTGCATCTTTTACCCCAGCATGAtttt       <  1:1073382/71‑1 (MQ=255)
                                                             gcTATAGAGAGCTTAGGTGCTTCCTGATTAACTGCCACGGGATTTGCATCTTTTACCCCAGCATGATTTTc      >  1:5404167/1‑71 (MQ=255)
                                                               tataGAGAGCTTAGGTGTTTCCTGATTAACTGCCACGGGATTTGCATCTTTTACCACAGCATGATTTTCGt    >  1:1268066/1‑71 (MQ=255)
                                                                ataGAGAGCTTAGGTGTTTCCTGATTAACTGCCACGGGATTTGCATCtttt                       <  1:716995/51‑1 (MQ=255)
                                                                 taGAGAGCTTAGGTGTTTCCTGATTAACTGCCACGGGATTTGCATCTTTTACCCCAGCATGATTTTCGt    <  1:2045282/69‑1 (MQ=255)
                                                                 taGAGAGCTTAGGTGTTTCCTGATTAACTGCCACGGGATTTGCATCTTTTACCCCAGCATGATTTTCGTTa  >  1:2267066/1‑71 (MQ=255)
                                                                   gagagCTTAGGTGTTTCCTGTTTAACTGCCACGGGAtt                                 <  1:1062890/38‑1 (MQ=255)
                                                                   gagagCTTAGGTGTTTCCTGATTAACTGCCACGGGAtt                                 <  1:3593625/38‑1 (MQ=255)
                                                                   gagagCTTAGGTGTTTCCTGATTAACTGCCACGGGAtt                                 <  1:3931278/38‑1 (MQ=255)
                                                                   gagagCTTAGGTGTTTCCTGATTAACTGCCACGGGATTTGCATCtttt                       >  1:7788920/1‑48 (MQ=255)
                                                                   |                                                                    
CCCCTCATTAATTAAGACAAAAGAAGTCTCATCATTCGATGTGAGTACAATGCCATTAAGCGCTATAGAGAGCTTAGGTGTTTCCTGATTAACTGCCACGGGATTTGCATCTTTTACCCCAGCATGATTTTCGTTA  >  AP009048/4184503‑4184638

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: