BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,392,506	C→G	13.2%	noncoding (35/76 nt)	glyX →	tRNA‑Gly

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,392,506	0	C	G	13.2%	76.9 / 6.6	37	noncoding (35/76 nt)	glyX	tRNA‑Gly
Reads supporting (aligned to +/- strand): ref base C (15/17); new base G (3/2); total (18/19)
Fisher's exact test for biased strand distribution p-value = 6.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCG  >  NC_000913/4392367‑4392646

aagatCAGTGGTTAGCGCAAGACCTTCCCTCGGGCGGGGTCGCGATTTCGAGTTTCGTTTCCAGATCCTAAATTTGAAAACTGCTGCAAAGAACAGACCACCCACGCGGGAATAGCTGAGTTGGTAGAGCACGACCTTGCCAAGGTCggg                                                                                                                                    <  2:519744/146‑1 (MQ=14)

 aGCTCAGTTGGTAGAGCACGACCTTTCCAAGGGCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATGTGAAAAGTGCTGCAAAGCACAGACCACCCACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCgggg                                                                                                                                   <  1:550468/150‑1 (MQ=255)

 aGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGATCGAGTCTCGTTTCCCTCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGGCAAGGTCgtgc                                                                                                                                   <  1:1093717/150‑4 (MQ=255)

  gCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTTAAAAGTGATGAAAATCACCCACCACCCAAGCGGGAATAGCTCAGTTGTTATAGCACTACCTTACAAAGGTCGGGGt                                                                                                                                  >  2:2014/1‑150 (MQ=11)

   cTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGCGGCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGGCGGGGTc                                                                                                                                 <  1:1142311/150‑1 (MQ=255)

      aGTTGGTAGAGCACGACCTTGCCAAGGGCGGGGTCGCGAGTTCGTGTCTCGTTTCCCGCTCCAAAATTTGAAAATTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTcgcg                                                                                                                              <  2:569921/150‑1 (MQ=255)

            tAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCg                                                                                                                        <  1:377923/150‑1 (MQ=39)

                      ccTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAATGCCCAGACAACCCATGGCGGAATAGCTCAGTTGGTATCGCACGACCTAGCCCAGGTCGGGGTCGCGAGTTCGAGTCTCGttt                                                                                                              >  1:411192/1‑150 (MQ=11)

                         tGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTACCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTccc                                                                                                           <  1:1106276/150‑1 (MQ=255)

                             aagtgaggggTCGCTAGTTCGTGGAGCGTTTCCCGCTCAAAAATTTGAAAAGTGCTGAACATCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGGCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTc                                                                                                       <  1:775576/144‑1 (MQ=11)

                                       tCGCGAGATCGAGTCTCGTTTCCCGCTCAAAAATTTGAAAAGTGCTGCAAACCACAGACCTCCTAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCAAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGa                                                                                             <  2:888647/150‑1 (MQ=14)

                                                cGAGTCTCGTTTCCCGCTCCAAAATTTGAAAACTGCTGCATAGCACAGACCACCCAAGCGGGAATAGCTCAGTTTGTAGAGCACGCCCTTGCCAAGGTCGGGGTCGGGAGTTCGAGTCTCGGTTCCGGCTCCAAAATTTGaaaatacagt                                                                                    >  2:552157/1‑143 (MQ=14)

                                                   gTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTATAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTAGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTaat                                                                                 >  2:724606/1‑149 (MQ=255)

                                                    tctcGTATCCCGCTCCAAAATTTGAAAAGTGCTGCAAACCACAGACCACGCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAgg                                                                                <  2:446553/150‑1 (MQ=255)

                                                    gcacGTTTCCCGCTCCAAAATTTGAAAAGTTCTGCAACGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAgg                                                                                <  1:1012764/147‑1 (MQ=255)

                                                       cGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGCGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGATAGTGCTGTAAGGcac                                                                             >  2:722825/1‑150 (MQ=255)

                                                           tCCCGCTCCAAAATTTGAAAAGTGCTGCAAACCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGac                                                                         <  2:1225289/150‑1 (MQ=255)

                                                               gCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTTAAATTGATGTAATGACCAGaccacc                                                                     >  2:487243/1‑150 (MQ=17)

                                                               ccTCCAAAATTTGAAAAGTGCTTCACCTAAAAGACCACCTCAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGaccacc                                                                     <  2:156157/149‑1 (MQ=17)

                                                                   cAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAg                                                                                                                      >  2:888028/1‑97 (MQ=255)

                                                                   cAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAg                                                                                                                      <  1:888028/97‑1 (MQ=255)

                                                                    aaaaTTTGAAAAGTGCTCCAAAGAACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACTGACCACCCAAGc                                                                <  1:1263997/150‑1 (MQ=255)

                                                                        ttagAAAAGTTCTGCACACAACATACCACCAACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGa                                                            <  1:1241358/147‑1 (MQ=14)

                                                                              aaGTTCTGAACATCACAGACCTCCGACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCt                                                      <  1:505498/150‑1 (MQ=14)

                                                                                     gcAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTTAAATAGCTGTAAGGCACAGACCACCCCAGAGGGTAAAGCAGAGTTgg                                               >  1:384103/1‑150 (MQ=17)

                                                                                      cacaacacaTACCTCCAAAGCGGCAATAGCTCAGTTGGTTGCGCCCGACATTGCCAAGGTCGGGGGCGCGAGTGCGCGTCCCGTTTCCCGCCCCAAAATTGGAAAGTGCTGTAAGGCACAGACCACCCACGCGGGAATAGCTCAGTTGGt                                              <  1:816006/145‑1 (MQ=14)

                                                                                               gACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTTGAGCACGa                                     >  1:593878/1‑150 (MQ=255)

                                                                                               gACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGAGCTGTAATGCACAGACCACCCAAGCGGGAATAGCTCAGTTGTTAGAGCACGa                                     >  1:1090434/1‑150 (MQ=255)

                                                                                                            ggAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTAGAGTATCGTTTCCCGCTCCAAAATATGAAAGTTATGAAAGGAACAGACCACCCAAACGCGAATATCACAGTTGGTATAGCAATACCATTCCaaacta                        >  1:279720/1‑146 (MQ=11)

                                                                                                                  gCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAATTGATTAAATTGACAGACACCCCAAGCGGGAATAGCTCACTTTGTAGAACAAGACCTTGCCAAGGTCGGCGta                  >  2:958256/1‑149 (MQ=11)

                                                                                                                   cTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTTAAAGTGATGTAAGGCACAGACCACCCAAGAGGGAATAGCTGAGTTTGTAGAGCAAGACCTTTACAAGGTATGGGTcg                 >  2:1108424/1‑150 (MQ=11)

                                                                                                                   cTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGAACAGACCCCCCACGCGGGAATAGCTCAGTTGGTAGAGAACGACCTAGCAAAGGTCGGGGTcg                 >  2:941558/1‑150 (MQ=11)

                                                                                                                                gCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGCGtct    >  2:988956/1‑150 (MQ=38)

                                                                                                                                ccccGACCTTGGCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTTAAAGTGCTGTAATGAACACACCACCCAAGCGGGATTAGCTCAGTTGGTAGAGCACGAGATAGCCAAGGTAGGTGTCGCGAGTTCGAGtca    >  2:763093/4‑149 (MQ=11)

                                                                                                                                ccccGACCTTGGCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGtct    >  1:420492/4‑150 (MQ=255)

                                                                                                                                 cccGACCTTGGCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGTATAGCTCAGTTTGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTGGAggcgg   >  1:354867/3‑146 (MQ=11)

                                                                                                                                  agccccTTGGCAAGGTCTGGGCCGATCTTGCGCGTCTCGTTTCCCCCGCCAAATTTTCAAAGTGCGGAACCTCACAGACCACCCACGCGGGAATAGCGCAGTTGGTAGAGCACGACCTTGCCAAGGGCGGGGTCGCGAGTTCGAGTCTCg  <  1:679926/146‑1 (MQ=2)

TAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCG  >  NC_000913/4392367‑4392646

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 14 ≤` `ATCG/ATCG` `< 15 ≤` `ATCG/ATCG` `< 27 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`