D66 :: BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	1,632,345	T→G	E28A (GAA→GCA)	SAUSA300_RS08070 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	1,632,345	0	T	G	94.1%	83.0 / ‑1.7	34	E28A (GAA→GCA)	SAUSA300_RS08070	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/2); new base G (4/28); total (4/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

CAACATTTTTATACTTCTTTGTATTTTTCACATCTTGAATACTTTCTTTTTCTAACTTTGTAACTTCATCAATTTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAACCATCGCCTCATGATAAATTCACTC  >  NC_007793/1632280‑1632440

cAACATTTTTATACTTCTTTGTATTTTTCACATCTTGAATACTTTCTTTTTCTAACTTTGTAACTTCATCAATTTCTTTATCATATTTCTGACCACATGc                                                               <  1:4721570/100‑1 (MQ=255)

                                tcactcaTACTTTTTCTTTCAGAATTTAAAACTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATa                               >  1:2390683/7‑100 (MQ=255)

                                tcactcaTACTTTTTCTTTCAGAATTTAAAACTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATa                               >  1:2946612/7‑100 (MQ=255)

                                 cactcaTACTTTTTCTTTCAGAATTTAAAACTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATAc                              >  1:4568314/6‑100 (MQ=255)

                                  actcaTACTTTTTCTTTCAGAATTTAAAACTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACa                             <  1:2683281/96‑1 (MQ=255)

                                   ctcaTACTTTTTCTTTCAGAATTTAAAACTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACaa                            <  1:5736548/97‑1 (MQ=255)

                                   ctcaTACTTTTTCTTTCAGAATTTAAAACTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACaa                            >  1:2131754/4‑100 (MQ=255)

                                      aTACTTTTTCTTTCAGAATTTAAAACTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAAcc                         <  1:2887236/100‑1 (MQ=255)

                                      aTACTTTTTCTTTCAGAATTTAAAACTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAAcc                         <  1:484466/100‑1 (MQ=255)

                                      aTACTTTTTCTTTCAGAATTTAAAACTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAAcc                         <  1:1147150/100‑1 (MQ=255)

                                      aTACTTTTTCTTTCAGAATTTAAAACTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAAcc                         <  1:521497/100‑1 (MQ=255)

                                      aTACTTTTTCTTTCAGAATTTAAAACTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAAcc                         <  1:2720558/100‑1 (MQ=255)

                                      aTACTTTTTCTTTCAGAATTTAAAACTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAAcc                         <  1:2204886/100‑1 (MQ=255)

                                       tACTTTTTCTTTCAGAATTTAAAACTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAACCa                        <  1:4389652/100‑1 (MQ=255)

                                                tttCTAACTTTGTAACTTCATCAATTTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAACCATCGCCtagt               <  1:2599159/100‑4 (MQ=255)

                                                   cagaatttaaaaCTGCATCAATCTCTTTATCATATTTCTGTCCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAACCATCGCCTCATGAt            <  1:5182523/90‑1 (MQ=255)

                                                   cagaatttaaaaCTGCATCAATCTCTTTATCATATTTCTGACCACGTGCAGATAATAATATAACTAATCCTAAAATTAATACAAACCATCGCCTCATGAt            <  1:3133114/90‑1 (MQ=255)

                                                   cagaatttaaaaCTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAACCATCGCCTCATGAt            <  1:4723097/90‑1 (MQ=255)

                                                    agaatttaaaaCTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAACCATCGCCTCATGATa           <  1:3958322/91‑1 (MQ=255)

                                                    agaatttaaaaCTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAACCATCGCCTCATGATa           <  1:499003/91‑1 (MQ=255)

                                                    agaatttaaaaCTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAACCATCGCCTCATGATa           <  1:5320058/91‑1 (MQ=255)

                                                    agaatttaaaaCTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAACCATCGCCTCATGATa           <  1:5750230/91‑1 (MQ=255)

                                                    agaatttaaaaCTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAACCATCGCCTCATGATa           <  1:748042/91‑1 (MQ=255)

                                                    agaatttaaaaCTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAACCATCGCCTCATGATa           <  1:4556545/91‑1 (MQ=255)

                                                    agaatttaaaaCTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAACCATCGCCTCATGATa           <  1:3093557/91‑1 (MQ=255)

                                                     gaatttaaaaCTGCATCAATCTCTTTAGCATATTTCTGACCACATGCAGATAATAATATAACTAAGCCTAAAATTAATACAAGCCATCGCCTCATGATaa          <  1:5477808/92‑1 (MQ=255)

                                                     gaatttGAAACTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAACCATCGCCTCATGATaa          <  1:3280461/97‑1 (MQ=255)

                                                      aatttaaaaCTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAACCATCGCCTCATGATaag         <  1:3838150/93‑2 (MQ=255)

                                                      aatttaaaaCTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAACCATCGCCTCATGATaag         <  1:1675201/93‑2 (MQ=255)

                                                      aatttaaaaCTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAACCATCGCCTCATGATaaa         <  1:1627176/93‑1 (MQ=255)

                                                          taaaaCTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAACCATCGCCTCATGATAAATTCa     <  1:2054829/97‑1 (MQ=255)

                                                           aaaaCTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAACCATCGCCTCATGATAAATTCAc    <  1:4724733/98‑1 (MQ=255)

                                                           aaaaCTGCATCAATCTCTTTACCATATTTCTGCCCACATGCATATAATAATATAACTAATCCTAAAACTAATACAATCCATCGCCTCATGATAAATTCAc    <  1:1246990/98‑1 (MQ=255)

                                                             aaCTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAACCATCGCCTCATGATAAATTCACtg  <  1:4634092/100‑2 (MQ=255)

CAACATTTTTATACTTCTTTGTATTTTTCACATCTTGAATACTTTCTTTTTCTAACTTTGTAACTTCATCAATTTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAACCATCGCCTCATGATAAATTCACTC  >  NC_007793/1632280‑1632440

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