BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,606,118	G→A	14.6%	noncoding (48/87 nt)	leuV ←	tRNA‑Leu

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,606,118	0	G	A	14.6%	56.3 / 2.3	36	noncoding (48/87 nt)	leuV	tRNA‑Leu
Reads supporting (aligned to +/- strand): ref base G (11/17); new base A (3/2); total (15/21)
Fisher's exact test for biased strand distribution p-value = 6.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAG‑ATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCT  >  NC_000913/4605974‑4606251

gCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAG‑ATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGacaa                                                                                                                                   >  2:276936/1‑148 (MQ=255)

  tACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAG‑ATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTa                                                                                                                                 >  1:598792/1‑149 (MQ=255)

      tgcGTGCGTTGCACTTATGTCTGGCATACTGCGGATGATAATGTATGTTTATCAGGGATGAAAAATCCTGGAAGAAACCAGAAGGGAATTATAT‑TTAATTTGGGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAacac                                                                                                                             <  2:1049516/147‑1 (MQ=16)

      aTCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAG‑ATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAacac                                                                                                                             >  1:276711/1‑149 (MQ=255)

             cTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAG‑ATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGc                                                                                                                      <  2:526310/149‑1 (MQ=255)

               ggCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAG‑ATAATTGGTGCGAGGGGGGGGACTTGAACACCCACGACCAAAAGTAAAAAAACACCTGaatcaa                                                                                                                    >  1:1165954/1‑145 (MQ=255)

                gCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAG‑ATAATTGGTGCGAGGGGGGTGACTTGAACCCCCACGTCCGGAAGGGCACAAGCAACTGCACCTAg                                                                                                                   >  1:889854/1‑149 (MQ=38)

                   cTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAG‑ATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAgcgc                                                                                                                <  1:735806/149‑1 (MQ=255)

                   attATTGATGGAATAATTCTCATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGACTAAACAAGAAAGGAAGTATAT‑ATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAgcgc                                                                                                                <  2:634307/148‑1 (MQ=25)

                      aTTGCTGGCATCATGCTGATGATAATGAATTTGAAACAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTTAAG‑ATAATTGGTGAGAGGGGTGGGACTTCAACCCCCACGTCCGTAAGGACAAAAACACATTAAGCTAGCGCGta                                                                                                             >  1:806440/1‑148 (MQ=37)

                        tGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAG‑ATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTa                                                                                                           <  2:675795/149‑1 (MQ=255)

                        tGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAG‑ATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAAAACCTGAAGCTAGCGAGTCTa                                                                                                           >  1:668646/1‑149 (MQ=255)

                         catttcatcatGCAGATTAAAATTAATTTTAATCAGCCACCAAAAAACAAGGAACAAACATTAATTGAATTAAAC‑ATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTAc                                                                                                          <  2:514085/144‑1 (MQ=9)

                          tctttatcatGATGATTATTATTAATTTTAATCAGGCATCAAAAATCACGGAAGTAACATTATATGAAGTACAG‑AGTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTAcc                                                                                                         <  2:284697/144‑1 (MQ=9)

                           tGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAG‑ATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCa                                                                                                        >  1:145003/1‑149 (MQ=255)

                            ggCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAG‑ATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCaa                                                                                                       >  1:400035/1‑149 (MQ=255)

                             gCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAG‑ATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCTCGTCTaagcaa                                                                                                      >  1:685170/1‑144 (MQ=255)

                             gCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAG‑ATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTCAAGATAGCGCTTCTCCCaaa                                                                                                      >  1:697562/1‑148 (MQ=37)

                                  atGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAG‑ATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                                                                                 <  2:864203/149‑1 (MQ=255)

                                         tttaaaTGAATTTTAATCAGGGATGAAAAAACAAGGATGAAACATTATATGAATTATAG‑ATTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGTACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTc                                                                                          <  2:400035/145‑1 (MQ=17)

                                         tgatAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAT‑ATAATTGGGGCGCGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTc                                                                                          <  1:617935/149‑1 (MQ=255)

                                           tttaTGATTTTTATTCAGGGATGAAAAATCAAGGAAGTAACTAGATTTGATGTATAG‑ATTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGTACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGc                                                                                        <  2:806440/146‑1 (MQ=14)

                                            aaaTGAATTTGAATCAGGGATGAAAAATAAAGGCAGAAACAAGACACGAAATAAAC‑AAAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGca                                                                                       <  2:224441/148‑1 (MQ=18)

                                              aTGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAACGGAAGTAAAG‑ATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcaca                                                                                     <  2:65442/149‑1 (MQ=255)

                                               tGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAACTAAAG‑ATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGACCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATACCGCCACCTTCGcaaac                                                                                    >  2:125628/1‑146 (MQ=255)

                                                  aTTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAG‑ATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcacacaca                                                                                 >  1:178926/1‑145 (MQ=255)

                                                      taaTCAGGTATGAAAAATCAAGGAAGAAACAATAAAGGAAGTAAAT‑ATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCtt                                                                             <  2:294965/148‑1 (MQ=255)

                                                               aTGAAAAATCAATGAAGAAAAAATAAAGGAAGTAAAG‑ATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACttttttt                                                                    <  2:218157/149‑1 (MQ=255)

                                                                 gAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAG‑ATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCAACCTTCTCACAGTGAGCGTACTTTTATTGa                                                                  >  2:463223/1‑149 (MQ=21)

                                                                   aaaaTCAAGGAAGAAACAAGAAAGGAAGTAAAG‑ATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCtt                                                                             >  2:626489/1‑136 (MQ=255)

                                                                   aaaaTCAAGGAAGAAACAAGAAAGGAAGTAAAG‑ATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCtt                                                                             <  1:626489/136‑1 (MQ=255)

                                                                             cacaaaaaaGCAATGAAGTAAAGCATATTTGCT‑CGTGGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGtt                                                      <  1:1019893/146‑1 (MQ=18)

                                                                               caaaCAAGAAAGGAAGTAAAG‑AAAATTGGGGCGAGGGGGGGGACTTTAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTg                                                    <  2:913453/148‑1 (MQ=255)

                                                                                           gAAGTAAAT‑ATAATTGGTGCTAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAAATCCGCCACCTTCGCACAGGCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAgggggg                                        <  2:974803/149‑1 (MQ=255)

                                                                                                                                aCCCCCACTTGCTTAAGGAGACTTACACATGCTGCTATCTCGGATACCCATTCCGCCACCTACGCTCAGTCATCTTACTATTTTTGATTTCGCCTAGATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAacac    <  1:488738/149‑1 (MQ=18)

                                                                                                                                  ccccACGGCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCATCGCACATTCATCTTACTTTTTTTGATTTCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACAccc  <  2:913446/149‑2 (MQ=255)

GCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAG‑ATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCT  >  NC_000913/4605974‑4606251

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: `‑`