BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,002,928	G→A	T89I (ACA→ATA)	yigG ←	inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,002,928	0	G	A	94.7%	51.8 / ‑3.1	19	T89I (ACA→ATA)	yigG	inner membrane protein
Reads supporting (aligned to +/- strand): ref base G (1/0); new base A (7/11); total (8/11)
Fisher's exact test for biased strand distribution p-value = 4.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.64e-01

TGCTTTATATTCTGGTGTTGCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAACTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTGATGATAGAGATGATATATACTGCTAATGTACCAAAAACATAAGTTTTTATATAGATGAAACCACTATCACGGAGTCGCTGGCAATTCATGTTGATGACGAGATAATGGAGTACGATGGTAGAGACTATAACAAGAAAGCCTGCTTCTCCATCGTTAAAAAAGATAATAAGAAAGGCAAAAATGAAATTTATTAAAATAAATGAAAATATATAACGACGTCTGGAAATCTT  >  NC_000913/4002683‑4003160

tgctTTATATTCTGGTGTTGCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAACTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTatgatg                                                                                                                                                                                                                                     >  1:38461/1‑251 (MQ=255)

     tataTTCTGGTGTTGCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAACTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTgatgat                                                                                                                                                                                                                                   >  1:12872/1‑248 (MQ=255)

             ggTGTTGCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAACTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTATGATGATAGAGATGata                                                                                                                                                                                                                         <  1:330374/250‑1 (MQ=255)

                         cAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATCTATTCCTTACTCATAAATACTCCATAACTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTATGATGATAGAGATGATATATACTGCTAATg                                                                                                                                                                                                            >  1:250259/1‑251 (MQ=255)

                                          aTGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAACTACTCCATAACTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTATGATGATAGAGATGATATATACTGCTAATGTACCAAAAACATAAGtt                                                                                                                                                                                           >  2:142444/1‑251 (MQ=255)

                                            gTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAACTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTATGATGATAGAGATGATATATACTGCTAATGTACCAAAAACATAAGttt                                                                                                                                                                                          <  2:259130/250‑1 (MQ=255)

                                               gATCATAAAAAGTAAAGCGTTATTTCCATTTTTCACTTAATGAACCAACGTTCATATATTCCTTACGCATAAATACGCCATAACGATTGTTTTGATGAATCAGTAGGTGCAAGCAGTAGCATACTGAAAGTGGAAAAATAACAAATCTAAATAATCATCTAACCATTGCCTGAACAGACAAAATCGTCGGGTATCACTATGATGATTGAGATGATATATGCTGTAAATGTAGCAAAAACATATTTTTTta                                                                                                                                                                                       <  2:182609/250‑1 (MQ=255)

                                                      aaaaTTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAACTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTATGATGATAGAGATGATATATACTGCTAATGTACCAAAAACATAAGTTTTTATATAGATg                                                                                                                                                                               >  1:201687/1‑251 (MQ=255)

                                                                      ttCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAACTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTATGATGATAGAGATGATATATACTGCTAATGTACCAAAAACATAAGTTTTTATATAGATGAAACCACTATCACgg                                                                                                                                                                <  2:205079/250‑1 (MQ=255)

                                                                                         aaCCATCGTTCATATATTCCTTACTCATAAATACTCCATAACTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTATGATGATAGAGATGATATATACTGCTAATGTACCAAAAACATAAGTTTTTATATAGATGAAACCACTATCACGGAGTCGCTGGCAATTCATGtt                                                                                                                                            <  1:318341/251‑1 (MQ=255)

                                                                                             atcgttactatatTCCTTACACCTAAATACTCCATAACCATTGTTTTGATGACTCAGTCGATGCAAGCATTAGCATACTGAAAGTGGCAAAATAACAAATCAAAAAAATCATAGAACCATTGCCTGAACAGGAAAAATCTTCGGCTATCATTATGATGATAGAGATGATATATACTGCTAATGTACCAAAAACATAAGTTTTTCTATAGATGAAACCACTATCACGGAGTCGCAGGCAATTCATGTtgat                                                                                                                                         >  2:257743/9‑250 (MQ=255)

                                                                                                                                           tgatgaATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATGATGATGATAGAGATGATATATACTGCTAATGTACCAAAAACATAAGTTTTTATATAGATGAAACCACTATCACGGAGTCGCTGGCAATTCATGTTGATGACGAGATAATGGAGTACGATGGTAGAGACTATAACAAGAAAGCCt                                                                                           <  1:107667/250‑1 (MQ=255)

                                                                                                                                                aaTCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTATGATGATAGAGATGATATATACTGCTAATGTACCAAAAACATAAGTTTTTATATAGATGAAACCACTATCACGGAGTCGCTGGCAATTCATGTTGATGACGAGATAATGGAGTACGATGGTAGAGACTATAACAAGAAAGCCTGCtt                                                                                       >  1:284762/1‑249 (MQ=255)

                                                                                                                                                aaTCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTATGATGATAGAGATGATATATACTGCTAATGTACCAAAAACATAAGTTTTTATATAGATGAAACCACTATCACGGAGTCGCTGGCAATTCATGTTGATGACGAGATAATGGAGTACGATGGTAGAGACTATAACAAGAAAGCCTGCtt                                                                                       <  2:284762/249‑1 (MQ=255)

                                                                                                                                                           gCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTATGATGATAGAGATGATATATACTGCTAATGTACCAAAAACATAAGTTTTTATATAGATGAAACCACTATCACGGAGTCGCTGGCAATTCATGTTGATGACGAGATAATGGAGTACGAAGGTAGAGACTATAACAAGAAAGCCTGCTTCTCCATCGTTaa                                                                           <  2:294075/250‑1 (MQ=255)

                                                                                                                                                               gCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTATGATGATAGAGATGATATATACTGCTAATGTACCAAAAACATAAGTTTTTATATAGATGAAACCTCTATCACGGAGTCGCTGGCAATTCATGTTGATGACGAGATAATGGAGTACGATGGTAGAGACTATAACAAGAAAGCCTGCTTCTCCATCGTTAAAAAAg                                                                      <  2:235392/251‑1 (MQ=255)

                                                                                                                                                                      aagaCTGAAAGGGGAAAAATAACAAAGGAAAAAAATCATAGATACAGTGCATGAACAGGAAAAAACTTCGGCAATCATTATGAGGATAGAGATGATATATTCTGCTAATGTACCAAAAAAATAAGGTGTGATATAGATGAAACCTCGAGTACGGAGTCGATGGCAATTCATGTTGATGACGAGATAATGGAGTACGATGGTAGAGACTATAACAAGAAAGCCTGCTTCTTGAGCGTTAAAATAGAtaataa                                                               <  2:83384/248‑1 (MQ=255)

                                                                                                                                                                              aaGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTATGATGATAGAGATGATATATACTGCTAATGTACCAAAAACATAAGTTTTTATATAGATGAAACCACTATCACGGAGTCGCTGGCAATTCATGTTGATGACGAGATAATGGAGTACGATGGTAGAGACTATAACAAGAAAGCCTGCTTCTCCATCGTTAAAAAAGATAATAAGAAAGGc                                                        >  1:215897/1‑250 (MQ=255)

                                                                                                                                                                                                       aaTCATCGAATCATTGCCTGAACAGGCAAAATCTTCTTCTAAAATTAAGATGATAGAGATGAGATAAACTGCTAATGTTGCAACAACATAAATTTTTATATAGAGGAAACCACTAGCACGGAGACGCTAGCAATTCATGTTTATGACGATAAAATTTAGTACGATAGTAGAGACGATAAAAAGAAAGCCTGTTAATCCATCGTTAAAAAAGATAATAAGAAAGGCAAAAATGAAATTTATTAAAATAAATg                              <  2:359798/251‑1 (MQ=255)

                                                                                                                                                                                                                                   aaaTCTTCGGCTATCATTATGATGCTAGAGATGATATATACTGCTAATGTACCAAAAACATAAGTTTTTATATAGCTGAAACCACTATCACGGAGTCGCTGGCAATTCATGTTGATGACGAGATAATGGAGTACGATGGTAGAGACTATAACAAGAAAGCCTGCTTCTCCATCGTTAAACAAGATAATAAGAAAGGCAAAAATGAAATTTATTAAAATAAATGAAAATATATAACGACGTCTGGAAATCtt  >  2:208446/1‑251 (MQ=255)

TGCTTTATATTCTGGTGTTGCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAACTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTGATGATAGAGATGATATATACTGCTAATGTACCAAAAACATAAGTTTTTATATAGATGAAACCACTATCACGGAGTCGCTGGCAATTCATGTTGATGACGAGATAATGGAGTACGATGGTAGAGACTATAACAAGAAAGCCTGCTTCTCCATCGTTAAAAAAGATAATAAGAAAGGCAAAAATGAAATTTATTAAAATAAATGAAAATATATAACGACGTCTGGAAATCTT  >  NC_000913/4002683‑4003160

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