D66 :: BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	2,798,664	T→A	E602V (GAA→GTA)	SAUSA300_RS14355 ←	surface anchored protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	2,798,664	0	T	A	98.3%	167.0 / NA	59	E602V (GAA→GTA)	SAUSA300_RS14355	surface anchored protein
Reads supporting (aligned to +/- strand): ref base T (0/0); major base A (9/49); minor base C (0/1); total (9/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTATTTTTCAACTTATTGCGTCGTGATAACCAAATCATTAGTACACCCATTGCACCAACAATTACAGATATCGGCAACCAATGTTCTTTTATCGTTTCCCCCGCTTTAGGCAAAGATACATTACCATCAGCATTTAATAATCCACTTAACAATCCATTACCTTTACCAAGTGTTAACGTCTTTTCTGGCTTT  >  NC_007793/2798567‑2798759

tttATTTTTCAACTTATTGCGTCGTGATAACCAAATCATTAGTACACCCATTGCACCAACAATTACAGATATCGGCAACCAATGTTCTTTTATCGTTAcc                                                                                               <  1:3427043/100‑1 (MQ=255)

    tttttCAACTTATTGCGTCGTGATAACCAAATCATTAGTACACCCATTGCACCAACAATTACAGATATCGGCAACCAATGTTCTTTTATCGTTAccccct                                                                                           <  1:4287121/100‑2 (MQ=255)

    tttttCAACTTATTGCGTCGTGATAACCAAATCATTAGTACACCCATTGCACCAACAATTACAGATATCGGCAACCAATGTTCTTTTATCGTTAccccat                                                                                           <  1:359224/100‑3 (MQ=255)

    tttttCAACTTATTGCGTCGTGATAACCAAATCATTAGTACACCCATTGCACCAACAATTACAGATATCGGCAACCAATGTTCTTTTATCGTTACCCCCg                                                                                           <  1:409481/100‑1 (MQ=255)

      tttCAACTTATTGCGTCGTGATAACCAAATCATTAGTACACCCATTGCACCAACAATTACAGATATCGGCAACCAATGTTCTTTTATCGTTAcccccaag                                                                                         <  1:52339/100‑4 (MQ=255)

      tttCAACTTATTGCGTCGTGATAACCAAATCATTAGTACACCCATTGCACCAACAATTACAGATATCGGCAACCAATGTTCTTTTATCGTTACCCCCGcc                                                                                         <  1:4313664/100‑2 (MQ=255)

      tttCAACTTATTGCGTCGTGATAACCAAATCATTAGTACACCCATTGCACCAACAATTACAGATATCGGCAACCAATGTTCTTTTATCGTTACCCCCGCt                                                                                         <  1:5124931/100‑1 (MQ=255)

        tCAACTTATTGCGTCGTGATAACCAAATCATTAGTACACCCATTGCACCAACAATTACAGATATCGGCAACCAATGTTCTTTTATCGTTACCCCCGCttt                                                                                       <  1:5130170/100‑1 (MQ=255)

        tCAACTTATTGCGTCGTGATAACCAAATCATTAGTACACCCATTGCACCAACAATTACAGATATCGGCAACCAATGTTCTTTTATCGTTACCCCCGCttt                                                                                       <  1:5064951/100‑1 (MQ=255)

        tCAACTTATTGCGTCGTGATAACCAAATCATTAGTACACCCATTGCACCAACAATTACAGATATCGGCAACCAATGTTCTTTTATCGTTACCCCCGCttt                                                                                       <  1:3814146/100‑1 (MQ=255)

        tCAACTTATTGCGTCGTGATAACCAAATCATTAGTACACCCATTGCACCAACAATTACAGATATCGGCAACCAATGTTCTTTTATCGTTACCCCCGCttt                                                                                       <  1:717376/100‑1 (MQ=255)

        tCAACTTATTGCGTCGTGATAACCAAATCATTAGTACACCCATTGCACCAACAATTACAGATATCGGCAACCAATGTTCTTTTATCGTTACCCCCGCttt                                                                                       <  1:3549583/100‑1 (MQ=255)

        tCAACTTATTGCGTCGTGATAACCAAATCATTAGTACACCCATTGCACCAACAATTACAGATATCGGCAACCAATGTTCTTTTATCGTTACCCCCGCttt                                                                                       <  1:2702320/100‑1 (MQ=255)

         cAACTTATTGCGTCGCGATAACCAAATCATTAGTACACCCATTGCACCAACAATTACAGATATCGGCAACCAATGTTCTTTTATCGTTACCCCCGCTTTa                                                                                      >  1:5588230/1‑100 (MQ=255)

            tctATTGCGTCGTGATAACCAAATCATTAGTACACCCATTGCACCAACAATTACAGATATCGGCAACCAATGTTCTTTTATCGTTACCCCCGCTTTAGGc                                                                                   <  1:630465/98‑1 (MQ=255)

            cTTATTGCGTCGTGATAACCAAATCATTAGTACACCCATTGCACCAACAATTACAGATATCGGCAACCAATGTTCTTTTATCGTTACCCCCGCTTTAGGc                                                                                   <  1:5675614/100‑1 (MQ=255)

            cTTATTGCGTCGTGATAACCAAATCATTAGTACACCCATTGCACCAACAATTACAGATATCGGCAACCAATGTTCTTTTATCGTTACCCCCGCTTTAGGc                                                                                   <  1:2302952/100‑1 (MQ=255)

            cTTATTGCGTCGTGATAACCAAATCATTAGTACACCCATTGCACCAACAATTACAGATATCGGCAACCAATGTTCTTTTATCGTTACCCCCGCTTTAGGc                                                                                   <  1:62387/100‑1 (MQ=255)

             ttATTGCGTCGTGATAACCAAATCATTAGTACACCCATTGCACCAACAATTACAGATATCGGCAACCAATGTTCTTTTATCGTTACCCCCTCTTTAGGCa                                                                                  <  1:541671/100‑1 (MQ=255)

             ttATTGCGTCGTGATAACCAAATCATTAGTACACCCATTGCACCAACAATTACAGATATCGGCAACCAATGTTCTTTTATCGTTACCCCCGCTTTAGGCa                                                                                  <  1:3795191/100‑1 (MQ=255)

             ttATTGCGTCGTGATAACCAAATCATTAGTACACCCATTGCACCAACAATTACAGATATCGGCAACCAATGTTCTTTTATCGTTACCCCCGCTTTAGGCa                                                                                  <  1:846135/100‑1 (MQ=255)

             ttATTGCGTCGTGATAACCAAATCATTAGTACACCCATTGCACCAACAATTACAGATATCGGCAACCAATGTTCTTTTATCGTTACCCCCGCTTTAGGCa                                                                                  <  1:8608/100‑1 (MQ=255)

             ttATTGCGTCGTGATAACCAAATCATTAGTACACCCATTGCACCAACAATTACAGATATCGGCAACCAATGTTCTTTTATCGTTACCCCCGCTTTAGGCa                                                                                  <  1:4175857/100‑1 (MQ=255)

             ttATTGCGTCGTGATAACCAAATCATTAGTACACCCATTGCACCAACAATTACAGATATCGGCAACCAATGTTCTTTTATCGTTACCCCCGCTTTAGGCa                                                                                  <  1:5314249/100‑1 (MQ=255)

             ttATTGCGTCGTGATAACCAAATCATTAGTACACCCATTGCACCAACAATTACAGATATCGGCAACCAATGTTCTTTTATCGTTACCCCCGCTTTAGGCa                                                                                  <  1:4577981/100‑1 (MQ=255)

              tATTGCGTCGTGATAACCAACTCATTAGTACACCCATTGCACCAACAATTACAGATATCGGCAACCAATGTTCTTTTATCGTTACCCCCGCTTTAGGCaa                                                                                 <  1:2943607/100‑1 (MQ=255)

              tATTGCGTCGTGATAACCAAATCATTAGTACACCCATTGCACCAACAATTACAGATATCGGCAACCAATGTTCTTTTATCGTTCCCCCCGCTTTAGGCaa                                                                                 <  1:2050298/100‑1 (MQ=255)

              tATTGCGTCGTGATAACCAAATCATTAGTACACCCATTGCACCAACAATTACAGATATCGGCAACCAATGTTCTTTTATCGTTACCCCCGCTTTAGGCaa                                                                                 <  1:4769320/100‑1 (MQ=255)

              tATTGCGTCGTGATAACCAAATCATTAGTACACCCATTGCACCAACAATTACAGATATCGGCAACCAATGTTCTTTTATCGTTACCCCCGCTTTAGGCaa                                                                                 <  1:994367/100‑1 (MQ=255)

              tATTGCGTCGTGATAACCAAATCATTAGTACACCCATTGCACCAACAATTACAGATATCGGCAACCAATGTTCTTTTATCGTTACCCCCGCTTTAGGCaa                                                                                 <  1:4469658/100‑1 (MQ=255)

              tATTGCGTCGTGATAACCAAATCATTAGTACACCCATTGCACCAACAATTACAGATATCGGCAACCAATGTTCTTTTATCGTTACCCCCGCTTTAGGCaa                                                                                 <  1:902191/100‑1 (MQ=255)

              tATTGCGTCGTGATAACCAAATCATTAGTACACCCATTGCACCAACAATTACAGATATCGGCAACCAATGTTCTTTTATCGTTACCCCCGCTTTAGGCaa                                                                                 <  1:1980390/100‑1 (MQ=255)

              tATTGCGTCGTGATAACCAAATCATTAGTACACCCATTGCACCAACAACTACAGATATCGGCAACCAATGTTCTTTTATCGTTACCCCCGCTTTAGGCaa                                                                                 <  1:2903017/100‑1 (MQ=255)

              tATTGCGTCGTGATAACCAAATCATTAGTACACCCACTGCACCAACAATTACAGATATCGGCAACCAATGTTCTTTTATCGTTACCCCCGCTTAAGGCaa                                                                                 <  1:3067924/100‑1 (MQ=255)

               aTTGCGTCGTGATAACCAAATCATTAGTACACCCATTGCACCAACAATTACAGATATCGGCAACCAATGTTCTTTTATCGTTACCCCCGCTTTAGGCaaa                                                                                <  1:4570300/100‑1 (MQ=255)

               aTTGCGTCGTGATAACCAAATCATTAGTACACCCATTGCACCAACAATTACAGATATCGGCAACCAATGTTCTTTTATCGTTACCCCCGCTTTAGGCaaa                                                                                >  1:3561416/1‑100 (MQ=255)

                             aCCAAATCATTAGTACACCCATTGCACCAACAATTACAGATATCGGCAACCAATGTTCTTTTATCGTTACCCCCGCTTTAGGCAAAGATACATTACCATc                                                                  >  1:3809807/1‑100 (MQ=255)

                             aCCAAATCATTAGTACACCCATTGCACCAACAATTACAGATATCGGCAACCAATGTTCTTTTATCGTTACCCCCGCTTTAGGCAAAGATACATTACCATc                                                                  <  1:5093205/100‑1 (MQ=255)

                             aCCAAATCATTAGTACACCCATTGCACCAACAATTACAGATATCGGCAACCAATGTTCTTTTATCGTTACCCCCGCTTTAGGCAAAGATACATTACCATc                                                                  <  1:5380302/100‑1 (MQ=255)

                                aaaTCATTAGTACACCCATTGCACCAACAATTACAGATATCGGCAACCAATGTTCTTTTATCGTTACCCCCGCTTTAGGCAAAGATACATTACCATCAGc                                                               <  1:582266/100‑1 (MQ=255)

                                          tACACCCATTGCACCAACAATTACAGATATCGGCAACCAATGTTCTTTTATCGTTACCCCCGCTTTAGGCAAAGATACATTACCATCAGCATTTaataat                                                     >  1:904106/1‑100 (MQ=255)

                                            cacCCATTGCACCAACAATTACAGATATCGGCAACCAATGTTCTTTTATCGTTACCCCCGCTTTAGGCAAAGATACATTACCATCAGCATTTAATAATcc                                                   <  1:162078/100‑1 (MQ=255)

                                               ccATTGCACCAACAATTACAGATATCGGCAACCAATGTTCTTTTATCGTTACCCCCGCTTTAGGCAAAGATACATTACCATCAGCATTTAATAATCCagc                                                <  1:831209/100‑3 (MQ=255)

                                                cATTGCACCAACAATTACAGATATCGGCAACCAATGTTCTTTTATCGTTACCCCCGCTTTAGGCAAAGATACATTACCATCAGCATTTAATAATCCACtt                                               <  1:4565102/100‑1 (MQ=255)

                                                cATTGCACCAACAATTACAGATATCGGCAACCAATGTTCTTTTATCGTTACCCCCGCTTTAGGCAAAGATACATTACCATCAGCATTTAATAATCCACtt                                               <  1:3053109/100‑1 (MQ=255)

                                                 aTTGCACCAACAATTACAGATATCGGCAACCAATGTTCTTTTATCGTTACCCCCGCTTTAGGCAAAGATACATTACCATCAGCATTTAATAATCCACTTa                                              >  1:5587275/1‑100 (MQ=255)

                                                     cACCAACAATTACAGATATCGGCAACCAATGTTCTTTTATCGTTACCCCCGCTTTAGGCAAAGATACATTACCATCAGCATTTAATAATCCACTTAACaa                                          >  1:3881959/1‑100 (MQ=255)

                                                                aCAGATATCGGCAACCAATGTTCTTTTATCGTTACCCCCGCTTTAGGCAAAGATACATTACCATCAGCATTTAATAATCCACTTAACAATCCATTACCtt                               >  1:4459168/1‑100 (MQ=255)

                                                                     tatCGGCAACCAATGTTCTTTTATCGTTACCCCCGCTTTAGGCAAAGATACATTACCATCAGCATTTAATAATCCACTTAACAATCCATTACCTTTTCCa                          <  1:5103365/100‑1 (MQ=255)

                                                                     tatCGGCAACCAATGTTCTTTTATCGTTACCCCCGCTTTAGGCAAAGATACATTACCATCAGCATTTAATAATCCACTTAACAATCCATTACCTTTACCa                          <  1:3906508/100‑1 (MQ=255)

                                                                       acGGCAACCAATGTTCTTTTATCGTTACCCCCGCTTTAGGCAAAGATACATTACCATCAGCATTTAATAATCCACTTAACAATCCATTACCTTTACCAAg                        >  1:370400/2‑100 (MQ=255)

                                                                        cGGCAACCAATGTTCTTTTATCGTTACCCCCGCTTTAGGCAAAGATACATTACCATCAGCATTTAATAATCCACTTAACAATCCATTACCTTTACCAAgt                       <  1:3846083/100‑1 (MQ=255)

                                                                                          tATCGTTACCCCCGCTTTAGGCAAAGATACATTACCATCAGCATTTAATAATCCACTTAACAATCCATTACCTTTACCAAGTGTTAACGTCTTTTCTGGc     <  1:183982/100‑1 (MQ=255)

                                                                                           aTCGTTACCCCCGCTTTAGGCAAAGATACATTACCATCAGCATTTAATAATCCACTTAACAATCCATTACCTTTACCAAGTGTTAACGTCTTTTCTGGca    <  1:4598188/100‑2 (MQ=255)

                                                                                           aTCGTTACCCCCGCTTTAGGCAAAGATACATTACCATCAGCATTTAATAATCCACTTAACAATCCATTACCTTTACCAAGTGTTAACGTCTTTTCTGGCt    <  1:727667/100‑1 (MQ=255)

                                                                                            tCGTTACCCCCGCTTTAGGCAAAGATACATTACCATCAGCATTTAATAATCCACTTAACAATCCATTACCTTTACCAAGTGTTAACGTCTTTTCTGGcaa   <  1:3130213/100‑3 (MQ=255)

                                                                                             cGTTACCCCCGCTTTAGGCAAAGATACATTACCATCAGCATTTAATAATCCACTTAACAATCCATTACCTTTACCAAGTGTTAACGTCTTTTCTGGCttt  <  1:5342468/100‑1 (MQ=255)

                                                                                             cGTTACCCCCGCTTTAGGCAAAGATACATTACCATCAGCATTTAATAATCCACTTAACAATCCATTACCTTTACCAAGTGTTAACGTCTTTTCTGGCttt  <  1:5227615/100‑1 (MQ=255)

                                                                                             cGTTACCCCCGCTTTAGGCAAAGATACATTACCATCAGCATTTAATAATCCACTTAACAATCCATTACCTTTACCAAGTGTTAACGTCTTTTCTGGCttt  >  1:1971322/1‑100 (MQ=255)

TTTATTTTTCAACTTATTGCGTCGTGATAACCAAATCATTAGTACACCCATTGCACCAACAATTACAGATATCGGCAACCAATGTTCTTTTATCGTTTCCCCCGCTTTAGGCAAAGATACATTACCATCAGCATTTAATAATCCACTTAACAATCCATTACCTTTACCAAGTGTTAACGTCTTTTCTGGCTTT  >  NC_007793/2798567‑2798759

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 17 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 37 ≤` `ATCG/ATCG` `< 39 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`