BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,392,506	C→G	16.3%	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	16.3%	75.0 / 15.2	43	noncoding (35/76 nt)	glyX	tRNA‑Gly
Reads supporting (aligned to +/- strand): ref base C (14/22); new base G (3/4); total (17/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.61e-01

CGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAG‑CACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCT  >  NC_000913/4392359‑4392644

cgcgcgAATAGCTCAGTTGATAGAGCACGACATTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAATTCATACAAAGCAAAAACCACCAAATCGGAAATAGATAATTTAGTAGAACAAGACCTTGcaa                                                                                                                                           >  2:1892837/1‑148 (MQ=11)

cgcgGGAATAGCTCAGTTGGTAGAGCAGGACCTTGCCAAGGGCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAATCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCa                                                                                                                                           <  2:120487/150‑1 (MQ=255)

       aTAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGGGAAGGTCgt                                                                                                                                    <  1:60277/150‑2 (MQ=255)

         aGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCgggg                                                                                                                                  <  2:1958689/150‑1 (MQ=39)

         aGCTCAGTTGGTAGAGCACGACCTTGCCAAGGGCGGGGTCGCGAGATCGAGGCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAACCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGGCAAGGTCgtgc                                                                                                                                  <  1:423390/150‑4 (MQ=14)

         aGCTCACTTTGTAGAGCACGACCTTCCCAAGGTCGGGGTCGCTACTTCGAGTCGCGTTTCCCGCGCCAAAATTTGAAAAGTCCTCCAAATCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGGCAAGGTCgtgc                                                                                                                                  <  1:816998/150‑4 (MQ=14)

              aGTTTGTAGCGCACGACCTTCCCCACCGCGGTGTCGCGAGTTCGAGGCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTCCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTcgcg                                                                                                                             <  1:1825191/150‑1 (MQ=14)

                   gTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGGCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAACCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTc                                                                                                                        <  1:173586/150‑1 (MQ=255)

                      gagCACGACCTTGCCAATGGCGGGGTCGCGAGTTCGCGTCTCGGTTCCCGCTCCTAAATTTGAAAAGTGCTTCAAAGCACAGACCTCCCACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCgcg                                                                                                                     <  2:2210776/150‑3 (MQ=14)

                                     aatgTCTGGGTCGCGATTTCGTGTCTCGGTTCCAGCTCCAAAATTTGAAAAGTGCTGCAAATCACATACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTc                                                                                                      <  2:2399916/147‑1 (MQ=14)

                                      agctaggggTCGCTACATCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAATCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTcc                                                                                                     <  1:1540070/145‑1 (MQ=14)

                                              gTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGGCAAGGTcctgct                                                                                                                                 <  2:1329571/114‑6 (MQ=255)

                                              gTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGGCAAGGTcctgct                                                                                                                                 >  1:1329571/1‑109 (MQ=255)

                                                     tttCGCGTCTCGGTTACGTCTCATAAATTTGAAAATTGCTTAATATAACTGACCACCGAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGc                                                                                      <  1:1675191/149‑1 (MQ=14)

                                                       tCGAGTCTAGTTTCCCGCTCAATCATTTGTAATGTGATTCATCTCACAGACCTCTCTAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTg                                                                                    <  1:1892837/150‑1 (MQ=14)

                                                             ttcGTTTACCTCTCATAAATTTGAAAAGTTCTTAAAATAAAAGACCAACCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGc                                                                              <  2:29327/149‑1 (MQ=17)

                                                             ctcGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGGGCCAAAATATGAAAGTGCGGGAAGGc                                                                              >  1:2182106/1‑150 (MQ=17)

                                                                 tttCCGGCTCCTAAATTCGAAAATTTCTTAACATAACATAAGACGCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAg                                                                          <  1:1686206/150‑1 (MQ=14)

                                                                 gttCCCGCGCAAAAATTTGAAAATTGCTCCAACTCACAGACCACGCACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAg                                                                          <  1:1676606/149‑1 (MQ=17)

                                                                           cAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCACGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAg                                                                >  2:8730/1‑150 (MQ=255)

                                                                            aaaaTTTGAAAAGTTCGGCACAACACTTACCACCCACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGc                                                               <  1:258613/150‑1 (MQ=17)

                                                                               atgggAAAAGTGCTCCAAACAACTGACCACCCACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTTGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCggg                                                            <  1:864260/146‑1 (MQ=17)

                                                                               aTTTGAAAATTTCTCAACATCAAATACCACGCACGCGGGAATAGCTCGGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCggg                                                            <  1:2436870/150‑1 (MQ=14)

                                                                               aTTTGAAAATTCCTCCAAACAACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCggg                                                            <  1:1145991/150‑1 (MQ=17)

                                                                                   gAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGAACAGACCACCCCAGCGGGAATa                                                        >  1:1766690/1‑150 (MQ=255)

                                                                                     aaaGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGAGCTGTAAGGACCAGACCACCCCAGCGGGAATAgg                                                      >  2:1914008/1‑149 (MQ=255)

                                                                                       tgTGCTGGAACGCACAGCCCCCCCAATCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGGCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTc                                                    <  1:1480800/149‑1 (MQ=14)

                                                                                        ttgctccACACCACATACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCa                                                   <  2:1229270/149‑1 (MQ=255)

                                                                                        gTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGtc                                                                                                                   <  1:2058642/86‑1 (MQ=38)

                                                                                        gTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGtc                                                                                                                   >  2:2058642/1‑86 (MQ=38)

                                                                                        gTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATGTGATATTGATGAAAGGGCCAGACCGCCCAGGGGGGATTATCtgg                                                   >  1:525171/1‑148 (MQ=14)

                                                                                                   cacaGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGTCCCAGACCACCCACGCGGGAATAGCTCCGTTGGTAGAG‑c                                       >  1:929524/1‑150 (MQ=255)

                                                                                                       gACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATATGATATGGATATAAGGAACAGACCACCGCAGAGGGATTAGGTGAGTTGTTagagtaagt                                    >  2:833781/1‑145 (MQ=14)

                                                                                                                   gggAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAG‑CACGACCTTGCCAAGGt                       >  2:1230724/1‑150 (MQ=38)

                                                                                                                        tAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAG‑CACGACCTTGGCAAGGTCgtgc                  <  1:8730/150‑4 (MQ=255)

                                                                                                                        tAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGGCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAG‑CACGACCTTGCCAAGGTCgggg                  <  1:1966829/150‑1 (MQ=38)

                                                                                                                          gCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTTTAAGGAACAGACCACCCAAGCGGGAATAGCTCAGTTGTTAGAG‑CACGACCATGCCAAGGTAGGGGTc                >  2:1928577/1‑150 (MQ=11)

                                                                                                                          gCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTCTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAG‑CACGACATTGCCAAGGTCGGGGTc                >  2:916702/1‑150 (MQ=255)

                                                                                                                          gCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATATTAAAGTGCTATAAGGAACAGACCACCCCAGAGCGAATAGCTGAGTTTGTAGAGACAC‑AGAATGCCAAGTTCGGGGTc                >  2:405297/1‑150 (MQ=11)

                                                                                                                             cAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAG‑CACGACCTTGCCAAGGTCGGGGTcgcg             >  2:1229176/1‑150 (MQ=38)

                                                                                                                             aattttgTATAGCACTACCTTCCCTAGTTCGTTGTCGCTTTTTCTTGTGTCGTTTCCCGCTCAAAAATTTTAAAGTGCTGTAATGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAG‑CACGACCTTGCCAAGGTCGGGGTcgcg             <  2:531500/144‑1 (MQ=11)

                                                                                                                                   gTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATATGAAAGAGCTGTAAGGAACAGACCACCCAAGCGGGAATAGCTCAGTTTGTAGAG‑CAAGACCTTACCAATGTCGGGGACGCGACTTcc       >  2:1225190/1‑149 (MQ=11)

                                                                                                                                    gcgaccccGACCTTGGCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGCCCACCCAAGCGGGAATAGCTCAGTTGGTAGAG‑CACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGa      >  1:1944899/8‑150 (MQ=255)

                                                                                                                                        ccccGACCTTGGCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCCCAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAG‑CACGCCCTTGCGAAGGTCGGGGTCGCGAGTTCGAGtct  >  1:2287822/4‑150 (MQ=11)

CGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAG‑CACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCT  >  NC_000913/4392359‑4392644

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