BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,227,582	C→T	10.2%	A88T (GCA→ACA)	PP_3706 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,227,582	0	C	T	10.2%	104.3 / 3.3	39	A88T (GCA→ACA)	PP_3706	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (20/15); new base T (2/2); total (22/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-01

GTGAGCCATCATGAGCTGCATTGGCTTATCTTCAGCCGTGAAAATCTTGGATGGTGAGTTAAGCCAACGACCGTCCAAGCGTGAGTAAAGGCTTTTTAGAAAGCTACTCCACTCCGCTTCGATGTATGCACGCTCCCCAGCATCGGTCACGTTTTCTGATACGACCGGTGCGCCAGGTCTACGAAAGTACGCGGCTGTTACCTCTTCGCCTCTAACGACTCTGCCACCCAGAGAAATAGACCACGCATCGCGAGGATGACTAGCCATTGTGCACA  >  NC_002947/4227454‑4227728

gTGAGCCATCATGAGCTGCATTGGCTTATCTTCAGCCGTGAAAATCTTGGATGGTGAGTTAAGCCAACGACCGTCCAAGCGTGAGTAAAGGCTTTTTAGAAAGCTACTCCACTCCGCTTCGATGTATGCACGCTCCCCAGCATCGGTc                                                                                                                                 <  2:243901/148‑1 (MQ=255)

          atGAGCTGCATTGGCTTATCTTCAGCCGTGAAAATCTTGGATGGTGAGTTAAGCCAACGACCGTCCAAGCGTGAGTAAAGGCTTTTTAGAAAGCTACTCCACTCCGCTTCGATGTATGCACGCTCCCCAGCATCGGTCACGTTTTCTg                                                                                                                       <  4:115814/148‑1 (MQ=255)

                  cATTGGCTTATCTTCAGCAGTGAAAATCTTGGATGGTGAGTTAAGCCAACGACCGTCCAAGCGTGAGTAAAGGCTTTTTATAAAGATACTCAACACCGCTTCTATGTATGCACGCTCCCCAGCAACGGAAACGTTTTCTGATACTACCg                                                                                                              >  1:71196/1‑149 (MQ=255)

                       gCTTATCTTCAGCCGTGAAAATCTTGGATGGTGAGTTAAGCCAACGACCGTCCAAGCGTGAGTAAAGGCTTTTTAGAAAGCTACTCCACTCCGCTTCGATGTATGCACGCTCCCCAGCATCGGTCACGTTTTCTGATACGACCGGTgcg                                                                                                         >  6:240284/1‑149 (MQ=255)

                        cTTATCTTCAGCCGTGAAAATCTAGGATGGTGAGTTAAGCCAACGACAGTCCAAGCGTGAGTAAAGGATATTAAGAAAGCTACTCCACTCAGCTTCGATGTATGCACGCTCCCCCGCATCGGTCGCGTATTCTCATACGACCGGTgcg                                                                                                         >  8:222843/1‑148 (MQ=255)

                         ttATCTTCAGCCGTGAAAATCTTGGATGGTGAGTTAAGCCAACGACCGTCCAAGCGTGAGTAAAGGCTTTTTAGAAAGCTACTCCACTCCGCTTCGATGTATGCACGCTCCCCAGCATCGGTCACGTTTTCTGATACGACCg                                                                                                              <  1:144993/142‑1 (MQ=255)

                         ttATCTTCAGCCGTGAAAATCTTGGATGGTGAGTTAAGCCAACGACCGTCCAAGCGTGAGTAAAGGCTTTTTAGAAAGCTACTCCACTCCGCTTCGATGTATGCACGCTCCCCAGCATCGGTCACGTTTTCTGATACGACCg                                                                                                              >  2:144993/1‑142 (MQ=255)

                          tATCTTCAGCCGTGAAAATCTTGGATGGTGAGTTAAGCCAACGACCGTCCAAGCGTGAGTAAAGGCTTTTTAGAAAGCTACTCCACTCCGCTTCGATGTATGTACGCTCCTCAGCATCGGTCACGttt                                                                                                                           <  4:144899/128‑1 (MQ=255)

                          tATCTTCAGCCGTGAAAATCTTGGATGGTGAGTTAAGCCAACGACCGTCCAAGCGTGAGTAAAGGCTTTTTAGAAAGCTACTCCACTCCGCTTCGATGTATGTACGCTCCTCAGCATCGGTCACGttt                                                                                                                           >  3:144899/1‑128 (MQ=255)

                                     gtgtAAATCTTGGATGGTGTGTTAAGGCAACGAACGTCATATCGTGAGTAAATGCTTTTTAGATAGCTACTCCACTCCGCTTCGATGTATGCACGCTCCCCAGCATCGGTCACGTTTTCTGATACGACCGGTGCGCCAGGTCTACGaaa                                                                                           <  2:71196/149‑1 (MQ=255)

                                          aaTCTTGGATGGTGAGTTAAGCCAACGACCGTCCAAGCGTGAGTAAAGGCTTTTTAGAAAGCTACTCCACTCCGCTTCGATGTATGCACGCTCCCCAGCATCGGTCACGTTTTCTGATACGACCGGTGCGCCAGGTCTACGaaa                                                                                           >  6:245660/1‑144 (MQ=255)

                                          aaTCTTGGATGGTGAGTTAAGCCAACGACCGTCCAAGCGTGAGTAAAGGCTTTTTAGAAAGCTACTCCACTCCGCTTCGATGTATGCACGCTCCCCAGCATCGGTCACGTTTTCTGATACGACCGGTGCGCCAGGTCTACGaaa                                                                                           <  5:245660/144‑1 (MQ=255)

                                           aTCTTGGATGGTGAGTTAAGCCAACGACCGTCCAAGCGTGAGTAAAGGCTTTTTAGAAAGCTACTCCACTCCGCTTCGATGTATGCACGCTCCCCAGCATCGGTCACGTTTTCTGATACGACCGGTGCGCCAGGTCTACGAAAGTAcgc                                                                                     <  7:66211/149‑1 (MQ=255)

                                             cTTGGATGGTGAGTTAAGCCAACGACCGTCCAAGCGTGAGTAAAGGCTTTTTAGAAAGCTACTCCACTCCGCTTCGATGTATGCACGCTCCCCAGCATCGGTCACGTTTTCTGATACGACCGGTGCGCCAGGTCTACGAAAGTACGCgg                                                                                   >  5:229877/1‑149 (MQ=255)

                                                   tGGTGAGTTAAGCCAACGACCGTCCAAGCGTGAGTAAAGGCTTTTTAGAAAGCTACTCCACTCCGCTTCGATGTATGCACGCTCCCCAGCATCGGTCACGTTTTCTGATACGACCGGTGCGCCAGGTCTACGAAAGTACGCGGCTGTTa                                                                             >  3:145554/1‑149 (MQ=255)

                                                     gTGAGTTAAGCCAACGACCGTCCAAGCGTGAGTAAAGGCTTTTTAGAAAGCTACTCCACTCCGCTTCGATGTATGCACGCTCCCCAGCATCGGTCACGTTTTCTGATACGACCGGTGCGCCAGGTCTACGAAAGTACGCGGCTGTTAc                                                                            <  4:145554/148‑1 (MQ=255)

                                                         gTTAAGCCAACGACCGTCCAAGCGTGAGTAAAGGCTTTTTAGAAAGCTACTCCACTCCGCTTCGATGTATGCACGCTCCCCAGCATCGGTCACGTTTTCTGATACGACCGGTGCGCCAGGTCTACGAAAGTACGCGGCTGTTACCTCtt                                                                       >  7:97512/1‑149 (MQ=255)

                                                         gTTAAGCCAACGACCGTCCAAGCGTGAGTAAAGGCTTTTTAGAAAGCTACTCCACTCCGCTTCGATGTATGCACGCTCCCCAGCATCGGTCACGTTTTCTGATACGACCGGTGCGCCAGGTCTACGAAAGTACGCGGCTGTTACCTCtt                                                                       >  1:126151/1‑149 (MQ=255)

                                                          ttAAGCCAACGACCGTCCAAGCGTGAGTAAAGGCTTTTTAGAAAGCTACTCCACTCCGCTTCGATGTATGCACGCTCCCCAGCATCGGTCACGTTTTCTGATACGACCGGTGCGCCAGGTCTACGAAAGTACGCGGCTGTTACct                                                                          <  7:225273/145‑1 (MQ=255)

                                                          ttAAGCCAACGACCGTCCAAGCGTGAGTAAAGGCTTTTTAGAAAGCTACTCCACTCCGCTTCGATGTATGCACGCTCCCCAGCATCGGTCACGTTTTCTGATACGACCGGTGCGCCAGGTCAACGAAAGTACGCGGCTGTTACct                                                                          >  8:225273/1‑145 (MQ=255)

                                                                         tCCAAGCGTGAGTAAAGGCTTTTTAGAAAGCTACTCCACTCCGCTTCGATGTATGCACGCTCCCCAGCATCGGTCACGTTTTCTGATACGACCGGTGCGCCAGGTCTACGAAAGTACGCGGCTGTTACCTCTTCGCCTCTAACGActct                                                       <  8:97512/149‑1 (MQ=255)

                                                                               cGTGAGTAAAGGCTTTTTAGAAAGCTACTCCACTCCGCTTCGATGTATGCACGCTCCCCAGCATCGGTCACGTTTTCTGATACGACCGGTGCGCCAGGTCTACGAAAGTACGCGGCTGTTACCTCTTCGCCTCTAACGACTCTGCCAcc                                                 <  6:229877/149‑1 (MQ=255)

                                                                                    gTAAAGGCTTTTTAGAAAGCTACTCCACTCCGCTTCGATGTATGCACGCTCCCCAGCATCGGTCACGTTTTCTGATACGACCGGTGCGCCAGGTCTACGAAAGTACGCGGCTGTTACCTCTTCGCCTCTAACGACTCTGCCACCCagag                                            <  5:240284/149‑1 (MQ=255)

                                                                                         ggCTTTTTAGAAAGCTACTCCACTCCGCTTCGATGTATGCACGCTCCCCAGCATCGGTCACGTTTTCTGATACGACCGGTGCGCCAGGTCTACGAAAGTACGCGGCTGTTACCTCTTCGCCTCTAACGACTCTGCCACCCAGAGaaa                                         >  2:198955/1‑147 (MQ=255)

                                                                                            tttttAGAAAGCTACTCCACTCCGCTTCGATGTATGCACGCTCCCCAGCATCGGTCACGTTTTCTGATACGACCGGTGCGCCAGGTCTACGAAAGTACGCGGCTGTTACCTCTTCGCCTCTAACGACTCTGCCACCCAGAGAAATAGAc                                    >  3:73157/1‑149 (MQ=255)

                                                                                                  gAAAGCTACTCCACTCCGCTTCGATGTATGCACGCTCCCCAGCATCGGTCACGTTTTCTGATACGACCGGTGCGCCAGGTCTACGAAAGTACGCGGCTGt                                                                               >  1:236700/1‑100 (MQ=255)

                                                                                                  gAAAGCTACTCCACTCCGCTTCGATGTATGCACGCTCCCCAGCATCGGTCACGTTTTCTGATACGACCGGTGCGCCAGGTCTACGAAAGTACGCGGCTGt                                                                               <  2:236700/100‑1 (MQ=255)

                                                                                                    aaGCTACTCCACTCCGCTTCGATGTATGCACGCTCCCCAGCATCGGTCACGTTTTCTGATACGACCGGTGCGCCAGGTCTACGAAAGTACGCGGCTGTTACCTCTTCGCCTCTAACGACTCTGCCACCCAGAGAAATAGACCACGCATc                            >  3:112809/1‑149 (MQ=255)

                                                                                                              actccGCTTCGATGTATGTACGCTCCCCAGCATCGGTCACGTTTTCTGATACGACCGGTGCGCCAGGTCTACGAAAGTACGCGGCTGTTACCTCTTCGCCTCTAACGActc                                                        <  5:162858/111‑1 (MQ=255)

                                                                                                              actccGCTTCGATGTATGTACGCTCCCCAGCATCGGTCACGTTTTCTGATACGACCGGTGCGCCAGGTCTACGAAAGTACGCGGCTGTTACCTCTTCGCCTCTAACGActc                                                        >  6:162858/1‑111 (MQ=255)

                                                                                                              actccGCTTCGATGTATGCACGCTCCCCAGCATCGGTCACGTTTTCTGATACGACCGGTGCGCCAGGTCTACGAAAGTACGCGGCTGTTACCTCTTCGCCTCTAACGACTCTGCCACCCAGAGAAATAGACCACGCATCGCGAGGATg                   >  7:8812/1‑148 (MQ=255)

                                                                                                               ctccGCTTCGATGTATGCACGCTCCCCAGCATCGGTCACGTTTTCTGATACGACCGGTGCGCCAGGTCTACGAAAGTACGCGGCTGTTACCTCTTCGCCTCTaa                                                              >  8:99377/1‑104 (MQ=255)

                                                                                                               ctccGCTTCGATGTATGCACGCTCCCCAGCATCGGTCACGTTTTCTGATACGACCGGTGCGCCAGGTCTACGAAAGTACGCGGCTGTTACCTCTTCGCCTCTaa                                                              <  7:99377/104‑1 (MQ=255)

                                                                                                                     ttCGATGTATGCACGCTCCCCAGCATCGGTCACGTTTTCTGATACGACCGGTGCGCCAGGTCTACGAAAGTACGCGGCTGTTACCTCTTCGCCTCTAACGACTCTGCCACCCAGAGAAATAGACCACGCATCGCGAGGATGACTAGCCa           >  4:124877/1‑149 (MQ=255)

                                                                                                                        gATGTATGCACGCTCCCCAGCATCGGTCACGTTTTCTGATACGACCGGTGCGCCAGGTCTACGAAAGTACGCGGCTGTTACCTCTTCGCCTCTAACGACTCTGCCACCCAGAGAAATAGACCACGCATCGCGAGGATGACTAGCCAtt         <  4:73157/148‑1 (MQ=255)

                                                                                                                        gATGTATGCACGCTCCCCAGCATCGGTCACGTTTTCTGATACGACCGGTGCGCCAGGTCTACGAAAGTACGCGGCTGTTACCTCTTCGCCTCTAACGACTCTGCCACCCAGAGAAATAGACCACGCATCGCGAGGATGACTAGCCAtt         <  2:89637/148‑1 (MQ=255)

                                                                                                                         aTGTATGCACGCTCCCCAGCATCGGTCACGTTTTCTGATACGACCGGTGCGCCAGGTCTACGAAAGTACGCGGCTGTTACCTCTTCGCCTCTAACGACTCTGCCACCCAGAGAAATAGACCACGCATCGCGAGGATGACTAGCCATtgt       >  4:144966/1‑149 (MQ=255)

                                                                                                                             aTGCACGCTCCCCAGCATCGGTCACGTTTTCTGATACGACCGGTGCGCCAGGTCTACGAAAGTACGCGGCTGTTACCTCTTCGCCTCTAACGACTCTGCCACCCAGAGAAATAGACCACGCATCGCGAGGATGACTAGCCATTGTGca    >  7:240116/1‑148 (MQ=255)

                                                                                                                              tGCACGCTCCCCAGCATCGGTCACGTTTTCTGATACGACCGGTGCGCCAGGTCTACGAAAGTACGCGGCTGTTACCTCTTCGCCTCTAACGACTCTGCCACCCAGAGAAATAGACCACGCATCGCGAGGATGACTAGCCATTGTGcaca  >  8:86759/1‑149 (MQ=255)

GTGAGCCATCATGAGCTGCATTGGCTTATCTTCAGCCGTGAAAATCTTGGATGGTGAGTTAAGCCAACGACCGTCCAAGCGTGAGTAAAGGCTTTTTAGAAAGCTACTCCACTCCGCTTCGATGTATGCACGCTCCCCAGCATCGGTCACGTTTTCTGATACGACCGGTGCGCCAGGTCTACGAAAGTACGCGGCTGTTACCTCTTCGCCTCTAACGACTCTGCCACCCAGAGAAATAGACCACGCATCGCGAGGATGACTAGCCATTGTGCACA  >  NC_002947/4227454‑4227728

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