BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,809,697	T→C	10.5%	E15E (GAA→GAG)	PP_3365 ←	acetolactate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,809,697	0	T	C	10.5%	91.8 / 3.3	38	E15E (GAA→GAG)	PP_3365	acetolactate synthase
Reads supporting (aligned to +/- strand): ref base T (19/15); new base C (2/2); total (21/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 = 4.73e-01

ACCTGCGGACTGCTCGTGACGGGTCAGTACCAGCTTGATCTTCGACTTGCGCAGGGATTCGAGCAGGTCGAGGTTTTCCTCACCGGGAATGCCGAATACATACTCGACACCTTCGTTTTCCAGGCATTGCACAACGACATCGGCGGCCTTGGCCATTTGGGGTACTACCTCAAAATCTTTGTGGGATTGCAGATTTATTAAGCTGCGTGACTGAGATGTTGCAGGATGTCTCAGTGGCGCGCATCGTAGGCCCTCTGGTT  >  NC_002947/3809586‑3809845

acatGCGGACTGCTCGTGACGGGTCAGTACCAGCTTGATCTTCGACTTGCGCAGGGATTCGAGCAGGTCGAGGTTTTCCTCACCGGGAATGCCGAATTCATACTCGACACCTTCGTTTTCCAGGCATTGCACAACGACATCGGCGGCCtt                                                                                                                >  1:1641896/4‑150 (MQ=1)

           gCTCGTGACGGGTCAGTAACAGCTTGATCTTCGACTTGCGCAGGGATTCGAGCAGGTCGAGGTTTTCCTCACCGGGAATGCCGAATACATACTCGACACCTTCGTTTTCCAGGCATTGCACAACGACATCGGCGGCCTTGGCCATTTggg                                                                                                     <  1:1318216/150‑1 (MQ=18)

            ctctTGACGGGTCAGTACCAGCTTGATCTTCGACTTGCGCAGGGATTCGAGCAGGTCGAGGTTTTCCTCACCGGTAATGCCGAATACATACTCGACACCTTCGTTTTCCAGGCATTGCACAACGACATCGGCGGCCTTGGCCATTTgggg                                                                                                    <  1:219154/150‑1 (MQ=21)

             tCGTGACGGGTCAGTACCAGCTTGATCTTCGACTTGCGCAGGGATTCGAGCAGGTCGAGGTTTTCCTCACCGGGAATGCCGAATACATACTCGACACCCTCGTTTTCCAGGCATTGCACAACGACATCGGCGGCCTTGGCCATTTGGGGt                                                                                                   >  2:226901/1‑150 (MQ=25)

             tCGTGACGGGTCAGTACCAGCTTGATCTTCGACTTGCGCAGGGATTCGAGCAGGTCGAGGTTTTCCTCACCGGGAATGCCGAATACATACTCGACACCCTCGTTTTCCAGGCATTGCACAACGACATCGGCGGCCTTGGCCATTTGGGGt                                                                                                   >  2:226917/1‑150 (MQ=25)

                        cAGTACCAGCTTGATCTTCGACTTGCCCCCCGATTCGAGCAGGTCGAGGTTTTCCTCACCGGGAATGCCGAATACATACTCGACACCTTCGTTTTCCAGGCATTGCACAACGACATCGGCGGCCTTGGCCATTTGGGGTACTACCTCaaa                                                                                        >  1:1454578/1‑150 (MQ=255)

                           tACCAGCTTGATCTTCGACTTGCGCAGGGATTCGAGCAGGTCGAGGTTTTCCTCACCGGGAATGCCGAATACATACTCGACACCTTCGTTTTCCAGGCATTGCACAACGACATCGGCGGCCTTGGCCATTTGGGGTACTACCTCAAAATc                                                                                     >  2:171660/1‑150 (MQ=255)

                             ccAGCTTGATCTTCGACTTGCGCAGGGATTCGAGCAGGTCGAGGTTTTCCTCACCGGGAATGCCGAATACATACTCGACACCTTCGTTTTCCAGGCATTGCACAACGACATGGGCGGCCTTGGCCATTTGGGGTACTACCTCAAAATCtt                                                                                   >  1:1163068/1‑150 (MQ=255)

                             ccAGCTTGATCTTCGACTTGCGCAGGGATTCGAGCAGGTCGAGGTTTTCCTCACCGGGAATGCCGAATACATACTCGACACCTTCGTTTTCCAGGCATTGCACAACGACATCGGCGGCCTTGGCCATTTGGGGTACTACCTCAAAATCtt                                                                                   <  2:1163068/150‑1 (MQ=255)

                              cAGCTTGATCTTCGACTTGCGCAGGGATTCGAGCAGGTCGAGGTTTTCCTCACCGGGAATGCCGAATACATACTCGACACCTTCGTTTTCCAGGCATTGCACAACGACATCGGCGGCCTTGGCCATTTGGGGTACTACCTCAAAATCttt                                                                                  >  1:822498/1‑150 (MQ=255)

                              cAGCTTGATCTTCGACTTGCGCAGGGATTCGAGCAGGTCGAGGTTTTCCTCACCGGGAATGCCGAATACATACTCGACACCTTCGTTTTCCAGGCATTGCACAACGACATCGGCGGCCTTGGCCATTTGGGGTACTACCTCAAAATCttt                                                                                  >  2:492449/1‑150 (MQ=255)

                                gCTTGATCTTCGACTTGCGCAGGGATTCGAGCAGGTCGAGGTTTTCCTCACCGGGAATGCCGAATACATACTCGACACCTTCGTTTTCCAGGCATTGCACAACGACATCGGCGGCCTTGGCCATTTGGGGTACTACCTCAAAATCTTtgt                                                                                <  2:822498/150‑1 (MQ=255)

                                 cTTGATCTTCGACTTGCGCAGGGATTCGAGCAGGTCGAGGTTTTCCTCACCGGGAATGCCGAATACATACTCGACACCTTCGTTTTCCAGGCATTGCACAACGACATCGGCGGCCTTGGCCATTTGGGGTACTACCTCAAAATCTTtgtg                                                                               >  1:958647/1‑150 (MQ=255)

                                 cTTGATCTTCGACTTGCGCAGGGATTAAGGCAGGTCGAGGTTTTCCTCACCGGGAATGCCGAATACATACTCGACACCTTCGTTTTCCAGGCATTGCACAACGACATCGGCGGCCTTGGCCATTTGGGGTACTACCTCAAAATCTTtgtg                                                                               >  1:1278398/1‑150 (MQ=255)

                                  ttGATCTTCGACTTGCGCAGGGATTCGAGCAGGTCGAGGTTTTCCTCACCGGGAATGCCGAATACATACTCGACACCTTCGTTTTCCAGGCATTGCACAATGACATCGGCGGCCTTGGCCATTTGGGGtactac                                                                                              <  2:811292/134‑1 (MQ=255)

                                  ttGATCTTCGACTTGCGCAGGGATTCGAGCAGGTCGAGGTTTTCCTCACCGGGAATGCCGAATACATACTCGACACCTTCGTTTTCCAGGCATTGCACAATGACATCGGCGGCCTTGGCCATTTGGGGtactac                                                                                              >  1:811292/1‑134 (MQ=255)

                                     aTCTTCGACTTGCGCAGGGATTCGAGCAGGTCGAGGTTTTCCTCACCGGGAATGCCGAATACATACTCGACACCTTCGTTTTCCAGGCATTGCACAACGACATCGGCGGCCTTGGCCATTTGGGGTACTACCTCAAAATCTTTGTGGGAt                                                                           <  1:492449/150‑1 (MQ=255)

                                      tCTTCGACTTGCGCAGGGATTCGAGCAGGTCGAGGTTTTCCTCACCGGGAATGCCGAATACATACTCGACACCTTCGTTTTCCAGGCATTGCACAACGACATCGGCGGCCTTGGCCATTTGGGGTACTACCTCAAAATCTTTGTGGGAtt                                                                          <  1:364096/150‑1 (MQ=255)

                                       cTTCGACTTGCGCAGGGATTCGAGCAGGTCGAGGTTTTCCTCACCGGGAATGCCGAATACATACTCGACACCTTCGTTTTCCAGGCATTGCACAACGACATCGGCGGCCTTGGCCATTTggg                                                                                                     >  2:413581/1‑122 (MQ=34)

                                       cTTCGACTTGCGCAGGGATTCGAGCAGGTCGAGGGTTTCCTCACCGGGAATGCCGAATACATACTCGACACCTTCGTTTTCCAGGCATTGCACAACGACATCGGCGGCCTTGGCCATTTggg                                                                                                     <  1:413581/122‑1 (MQ=21)

                                        ttCGACTTGCGCAGGGATTCGAGCAGGTCGAGGTTTTCCTCACCGGGAATGCCGAATACATACTCGACACCTTCGTTTTCCAGGCATTGCACAACGACATCGGCGGCCTTGGCCATTTGGGGTACTACCTCAAAATCTTTGTGGGATTGc                                                                        >  1:481715/1‑150 (MQ=255)

                                            aCTTGCGCAGGGATTCGAGCAGGTCGAGGTTTTCCTCACCGGGAATGCCGAATACATACTCGACACCTTCGTTTTCCAGGCATTGCACAACGACATCGGCGGCCTTGGCCATTTGGGGTACTACCTCAAAATCTTTGTGGGATTGCAGAt                                                                    >  1:1500443/1‑150 (MQ=255)

                                                   cAGGGATTCGAGCAGGTCGAGGTTTTCCTCACCGGGAATGCCGAATACATACTCGACACCCTCGTTTTCCAGGCATTGCACAACGACATCGGCGGCCTTGGCCATTTGGGGTACTACCTCAAAATCTTTGTGGGATTGCAGATTTATTaa                                                             <  1:226917/150‑1 (MQ=255)

                                                   cAGGGATTCGAGCAGGTCGAGGTTTTCCTCACCGGGAATGCCGAATACATACTCGACACCCTCGTTTTCCAGGCATTGCACAACGACATCGGCGGCCTTGGCCATTTGGGGTACTACCTCAAAATCTTTGTGGGATTGCAGATTTATTaa                                                             <  1:226901/150‑1 (MQ=255)

                                                              gCAGGTCGAGGTTTTCCTCACCGGGAATGCCGAATACATACTCGACACCTTCGTTTTCCAGGCATTGCACAACGACATCGGCGGCCTTGGCCATTTGGGGTACTACCTCAAAATCTTTGTGGGATTGCAGATTTATTAAGCTGCGTGACt                                                  <  1:1632359/150‑1 (MQ=255)

                                                                aGGTCGAGGTTTTCCTCACCGGGAATGCCGAATACATACTCGACACCTTCGTTTTCCAGGCATTGCACAACGACATCGGCGGCCTTGGCCATTTGGGGTACTACCTCAAAATCTTTGTGGGATTGCAGATTTATTAAGCTGCGTGACTga                                                <  1:171660/150‑1 (MQ=255)

                                                                            tCCTCACCGGGAATGCCGAATACATACTCGACACCTTCGTTTTCCAGGCATTGCACAACGACATCGGCGGCCTTGGCCATTTGGGGTACTACCTCAAAATCTTTGTGGGATTGCa                                                                       <  1:724146/115‑1 (MQ=255)

                                                                            tCCTCACCGGGAATGCCGAATACATACTCGACACCTTCGTTTTCCAGGCATTGCACAACGACATCGGCGGCCTTGGCCATTTGGGGTACTACCTCAAAATCTTTGTGGGATTGCa                                                                       >  2:724146/1‑115 (MQ=255)

                                                                                cACCGGGAATGCCGAATACATACTCGACACCTTCGTTTTCCAGGCATTGCACAACGACATCGGCGGCCTTGGCCATTTGGGGTACTACCTCAAAATCTTTGTGGGATTGCAGATTTATTAAGCTGCGTGACTGAGATGTTGCAGGATGtc                                <  2:1278398/150‑1 (MQ=255)

                                                                                  ccGGGAATGCCGAATACATACTCGACACCTTCGTTTTCCAGGCATTGCACAACGACATCGGCGGCCTTGGCCATTTGGGGTACTACCTCAAAATCTTTGTGGGATTGCAGATTTATTAAGCTGCGTGACTGAGATGTTGCAGGATGtctc                              >  2:891228/1‑150 (MQ=255)

                                                                                             gAATACATACTCGACACCTTCGTTTTCCAGGCATTGCACAACGACATCGGCGGCCTTGGCCATTTGGGGTACTACCTCAAAATCTTTGTGGGATTGCAGATTTATTAAGCTGCGTGACTGAGATGTTGCAGGATGTCTCAGTGgcgc                      <  1:1062882/147‑1 (MQ=255)

                                                                                             gAATACATACTCGACACCTTCGTTTTCCAGGCATTGCACAACGACATCGGCGGCCTAGGCCATTTGGGGTACTACCTCAAAATCTTTGTGGGATTGCAGATTTATTAAGCTGCGTGACTGAGATGTTGCAGGATGTCTCAGTGgcgc                      >  2:1062882/1‑147 (MQ=255)

                                                                                                   atacTCGACACCTTCGTTTTCCAGGCATTGCACAACGACATCGGCGGCCTTGGCCATTTGGGGTACTACCTCAAAATCTTTGTGGGATTGCAGATTTATTAAGCTGCGTGACTGAGATGTTGCAGGATGTCTCAGTGGCGCGCATCGTAg             >  2:886928/1‑150 (MQ=255)

                                                                                                     acTCGACACCTTCGTTTTCCAGGCATTGCACAACGACATCGGCGGCCTTGGCCATTTGGGGTACTACCTCAAAATCTTTGTGGGATTGCAGATTTATTAAGCTg                                                         >  2:127973/1‑104 (MQ=255)

                                                                                                     acTCGACACCTTCGTTTTCCAGGCATTGCACAACGACATCGGCGGCCTTGGCCATTTGGGGTACTACCTCAAAATCTTTGTGGGATTGCAGATTTATTAAGCTg                                                         <  1:127973/104‑1 (MQ=255)

                                                                                                       tCGACACCTTCGTTTTCCAGGCATTGCACAACGACATCGGCGGCCTTGGCCATTTGGGGTACTACCTCAAAATCTTTGTGGGATTGCAGATTTATTAAGCTGCGTGACTGAGATGTTGCAGGATGTCTCAGTg                          <  2:233119/133‑1 (MQ=255)

                                                                                                       tCGACACCTTCGTTTTCCAGGCATTGCACAACGACATCGGCGGCCTTGGCCATTTGGGGTACTACCTCAAAATCTTTGTGGGATTGCAGATTTATTAAGCTGCGTGACTGAGATGTTGCAGGATGTCTCAGTg                          >  1:233119/1‑133 (MQ=255)

                                                                                                              cTTCGTTTTCCAGGCATTGCACAACGACATCGGCGGCCTTGGCCATTTGGGGTACTACCTCAAAATCTTTGTGGGATTGCAGATTTATTAAGCTGCGTGACTGAGATGTTGCAGGATGTCTCAGTGGCGCGCATCGTAGGCCCTCTGGtt  >  1:1347938/1‑150 (MQ=255)

                                                                                                               ttCGTTTTCCAGGCATTGCACAACGACATCGGCGGCCTTGGCCATTTGGGGtactac                                                                                              <  1:737965/57‑1 (MQ=255)

                                                                                                               ttCGTTTTCCAGGCATTGCACAACGACATCGGCGGCCTTGGCCATTTGGGGtactac                                                                                              >  2:737965/1‑57 (MQ=255)

ACCTGCGGACTGCTCGTGACGGGTCAGTACCAGCTTGATCTTCGACTTGCGCAGGGATTCGAGCAGGTCGAGGTTTTCCTCACCGGGAATGCCGAATACATACTCGACACCTTCGTTTTCCAGGCATTGCACAACGACATCGGCGGCCTTGGCCATTTGGGGTACTACCTCAAAATCTTTGTGGGATTGCAGATTTATTAAGCTGCGTGACTGAGATGTTGCAGGATGTCTCAGTGGCGCGCATCGTAGGCCCTCTGGTT  >  NC_002947/3809586‑3809845

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