BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,937,812	A→C	22.1%	L429R (CTT→CGT) ‡	PP_2570 ←	transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,937,812	0	A	C	22.1%	26.9 / 5.8	22	L429R (CTT→CGT) ‡	PP_2570	transposase
Reads supporting (aligned to +/- strand): ref base A (8/9); new base C (2/3); total (10/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATATCACTCACACCTGAAGGTCGAGTTTGCCTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGCTTCGAACTCCTCAAG  >  NC_002947/2937670‑2937861

atatCACTCACACCTGAAGGTCGAGTTTGCCTGTGCCGGCCCCTTCGCTGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGc                                            >  1:239520/1‑150 (MQ=255)

  atCACTCACACCTGAAGGTCGAGTTTGCCTGTGCCGGCCCCTTCGCTGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGccc                                          <  2:239520/150‑1 (MQ=255)

  atCACTCACACCTGAAGGTCGAGTTTGCCTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGccc                                          <  1:928681/150‑1 (MQ=255)

      ctcacACCTGAAGGTCGAGTTTGCCTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGAc                                      <  1:1153540/150‑1 (MQ=255)

         acacCTGAAGGTCGAGTTTGCCTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTAAGAGTCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGTCCGGACCTg                                   >  1:720493/1‑150 (MQ=17)

          cacCTGAAGGTCGAGTTTGCCTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGc                                  <  2:220448/150‑1 (MQ=255)

            ccTGAAGGTCGAGTTTGCCTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATCCCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCg                                <  2:493993/150‑1 (MQ=255)

            ccTGAAGGTCGAGTTTGCCTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCAACAGTTGATACCCCAGTTTCCGGGCTTTGCGACGCAAAGCCCGGACCTGCCg                                <  2:688473/150‑1 (MQ=255)

             cTGAAGGTCGAGTTTGCCTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCgg                               >  1:80192/1‑150 (MQ=255)

             cTGAAGGTCGAGTTTGACTGTGCCGGCCCCTTCGCGGGTGAACCCGCTAATCTCAAACCAAATGAAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTGTGCGAAGCAAAGCCCGGACCTGCCgg                               >  2:500142/1‑150 (MQ=17)

                 aGGGCGAGTTTGCCTGTGACGGAAACTTCGCGGTTTTACTCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCtt                           <  2:11703/150‑1 (MQ=11)

                   gTCGAGTTTGCCTGTGCCGTCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCTTCGCAAAGCCCg                                         <  1:204581/134‑1 (MQ=17)

                   gTCGAGTTTGCCTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCTTCGCAAAGCCCg                                         >  2:204581/1‑134 (MQ=255)

                         tttGCCTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTtggcttctg                   <  1:949080/150‑1 (MQ=255)

                            gCCTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCCGCCGGTCTtggcttctggct                >  1:264552/1‑150 (MQ=255)

                              cTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCAGCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCtt                           <  2:832113/137‑1 (MQ=255)

                              cTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCTTCGCAAAGCCCgg                                        >  1:561039/1‑124 (MQ=255)

                              cTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCTTCGCAAAGCCCgg                                        <  2:561039/124‑1 (MQ=255)

                              cTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCtt                           >  1:832113/1‑137 (MQ=39)

                                     ggCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGCTTCGAActcc       <  1:601007/150‑1 (MQ=32)

                                         ccTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGCTTCGAACTCCTCaa   >  2:854874/1‑150 (MQ=31)

                                          cTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGCTTCGAACTCCTCAAg  >  2:521572/1‑150 (MQ=11)

ATATCACTCACACCTGAAGGTCGAGTTTGCCTGTGCCGGCCCCTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGCTTCGAACTCCTCAAG  >  NC_002947/2937670‑2937861

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