BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,550,146	G→A	55.8%	T38T (ACC→ACT)	vacJ ←	predicted lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,550,146	0	G	A	55.8%	23.7 / 99.1	120	T38T (ACC→ACT)	vacJ	predicted lipoprotein
Reads supporting (aligned to +/- strand): ref base G (14/39); new base A (46/21); total (60/60)
Fisher's exact test for biased strand distribution p-value = 7.62e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

Only 100 of 120 total aligned reads displayed.

ACGCCAGGCGACAGCGACCGGTCGAACAATATACGGGTCTAATACATTGAAGTTGAAGTTGTACATGGTGCGGTTGAACCCTTCTAACGGGTCAGAACGCCCTTGCTGATCTGTACCGGAACTCGCACACCCC  >  minE/1550080‑1550212

aCGCCAGGCGACAGCGACCGGTCGAACAATATACGGGTCTAATACATTGAAGTTGAAGTTGTACATGGt                                                                  >  1:3291312/1‑69 (MQ=255)

 cGCCAGGCGCCAGCGACCGGTCGAACAATATACGGGTCTAATACATTGAAGTTGAAGTTGTACATAGTg                                                                 >  1:1161913/1‑69 (MQ=255)

 cGCCAGGCGACAGCGACCGGTCGAACAATATACGGGTCTAATACATTGAAGTTGAAGTTGTACATAGTg                                                                 >  1:849923/1‑69 (MQ=255)

 cGCCAGGCGACAGCGACCGGTCGAACAATATACGGGTCTAATACATTGAAGTTGAAGTTGTACATAGTg                                                                 >  1:1574983/1‑69 (MQ=255)

 cGCCAGGCGACAGCGACCGGTCGAACAATATACGGGTCTAATACATTGAAGTTGAAGTTGTACATAGTg                                                                 >  1:1828586/1‑69 (MQ=255)

 cGCCAGGCGACAGCGACCGGTCGAACAATATACGGGTCTAATACATTGAAGTTGAAGTTGTACATAGTg                                                                 >  1:922477/1‑69 (MQ=255)

 cGCCAGGCGACAGCGACCGGTCGAACAATATACGGGTCTAATACATTGAAGTTGAAGTTGTACATAGTg                                                                 >  1:1902593/1‑69 (MQ=255)

 cGCCAGGCGACAGCGACCGGTCGAACAATATACGGGTCTAATACATTGAAGTTGAAGTTGTACATAGTg                                                                 >  1:596903/1‑69 (MQ=255)

 cGCCAGGCGACAGCGACCGGTCGAACAATATACGGGTCTAATACATTGAAGTTGAAGTTGTACATAGTg                                                                 >  1:2462272/1‑69 (MQ=255)

 cGCCAGGCGACAGCGACCGGTCGAACAATATACGGGTCTAATACATTGAAGTTGAAGTTGTACATAGTg                                                                 >  1:1283861/1‑69 (MQ=255)

 cGCCAGGCGACAGCGACCGGTCGAACAATATACGGGTCTAATACATTGAAGTTGAAGTTGTACATAGTg                                                                 >  1:1264376/1‑69 (MQ=255)

 cGCCAGGCGACAGCGACCGGTCGAACAATATACGGGTCTAATACATTGAAGTTGAAGTTGTACATAGTg                                                                 >  1:944683/1‑69 (MQ=255)

 cGCCAGGCGACAGCGACCGGTCGAACAATATACGGGTCTAATACATTGAAGTTGAAGTTGTACATAGTg                                                                 >  1:1194427/1‑69 (MQ=255)

 cGCCAGGCGACAGCGACCGGTCGAACAATATACGGGTCTAATACATTGAAGTTGAAGTTGTACATAGTg                                                                 >  1:2832423/1‑69 (MQ=255)

 cGCCAGGCGACAGCGACCGGTCGAACAATATACGGGTCTAATACATTGAAGTTGAAGTTGTACATAGTg                                                                 >  1:3365805/1‑69 (MQ=255)

 cGCCAGGCGACAGCGACCGGTCGAACAATATACGGGTCTAATACATTGAAGTTGAAGTTGTACATAGTg                                                                 >  1:112004/1‑69 (MQ=255)

 cGCCAGGCGACAGCGACCGGTCGAACAATATACGGGTCTAATACATTGAAGTTGAAGTTGTACATAGTg                                                                 >  1:3214897/1‑69 (MQ=255)

 cGCCAGCCGACAGCGACCGGTCGAACAATATACGGGTCTAATACATTGAAGTTGAAGTTGTACATAGTg                                                                 >  1:439349/1‑69 (MQ=255)

   ccAGGCGACAGCGACCGGTCGAACAATATACGGGTCTAATACATTGAAGTTGAAGTTGTACATGGTGCg                                                               <  1:695383/69‑1 (MQ=255)

    cAGGCGACAGCGACCGGTCGAACAATATACGGGGCTAATACATTGAAGTTGAAGTTGTACATGGTGCgg                                                              <  1:948836/69‑1 (MQ=255)

       gCGACAGCGACCGGTCGAACAATATACGGGTCTAATACATTGAAGTTGAAGTTGTACATGGTGCGGTTg                                                           <  1:1390217/69‑1 (MQ=255)

           cAGCGACCGGTCGAACAATATACGGGTCTAATACATTGAAGTTGAAGTTGTACATGGTGCGGTTGAAcc                                                       <  1:694796/69‑1 (MQ=255)

                aCCGGTCGAACAATATACGGGTCCAATACATTGAAGTTGAAGTTGTACATGGTGCGGTTGAACCCTTCt                                                  <  1:610355/69‑1 (MQ=255)

                     tCGAACAATATACGGGTCTAATACATTGAAGTTGAAGTTGTACATAGTGCGGTTAAACCCTTCTAACgg                                             >  1:3104572/1‑69 (MQ=255)

                     tCGAACAATATACGGGTCTAATACATTGAAGTTGAAGTTGTACATAGTGCGGTTAAACCCTTCTAACgg                                             >  1:620306/1‑69 (MQ=255)

                     tCGAACAATATACGGGTCTAATACATTGAAGTTGAAGTTGTACATAGTGCGGTTAAACCCTTCTAACgg                                             >  1:2664424/1‑69 (MQ=255)

                     tCGAACAATATACGGGTCTAATACATTGAAGTTGAAGTTGTACATAGTGCGGTTAAACCCTTCTAACgg                                             >  1:390147/1‑69 (MQ=255)

                     tCGAACAATATACGGGTCTAATACATTGAAGTTGAAGTTGTACATAGTGCGGTTAAACCCTTCTAACgg                                             >  1:678669/1‑69 (MQ=255)

                     tCGAACAATATACGGGTCTAATACATTGAAGTTGAAGTTGTACATAGTGCGGTTAAACCCTTCTAACgg                                             >  1:2053742/1‑69 (MQ=255)

                     tCGAACAATATACGGGTCTAATACATTGAAGTTGAAGTTGTACATAGTGCGGTTAAACCCTTCTAACgg                                             >  1:1342045/1‑69 (MQ=255)

                     tCGAACAATATACGGGTCTAATACATTGAAGTTGAAGTTGTACATAGTGCGGTTAAACCCTTCTAACgg                                             >  1:1313508/1‑69 (MQ=255)

                     tCGAACAATATACGGGTCTAATACATTGAAGTTGAAGTTGTACATAGTGCGGTTAAACCCTTCTAACgg                                             >  1:276267/1‑69 (MQ=255)

                     tCGAACAATATACGGGTCTAATACATTGAAGTTGAAGTTGTACATAGTGCGGTTAAACCCTTCTAACgg                                             >  1:119204/1‑69 (MQ=255)

                     tCGAACAATATACGGGTCTAATACATTGAAGTTGAAGTTGTACATAGTGCGGTTAAACCCTTCTAACgg                                             >  1:1538414/1‑69 (MQ=255)

                     tCGAACAATATACGGGTCTAATACATTGAAGTTGAAGTTGTACATAGTGCGGTTAAACCCTTCTAACgg                                             >  1:142559/1‑69 (MQ=255)

                     tCGAACAATATACGGGTCTAATACATTGAAGTTGAAGTTGTACATAGTGCGGTTAAACCCTTCTAACgg                                             >  1:1754118/1‑69 (MQ=255)

                     tCGAACAATATACGGGTCTAATACATTGAAGTTGAAGTTGTACATAGTGCGGTTAAACCCTTCTAACgg                                             >  1:2970977/1‑69 (MQ=255)

                     tCGAACAATATACGGGTCTAATACATTGAAGTTGAAGTTGTACATAGTGCGGTTAAACCCTTCTAACgg                                             >  1:170727/1‑69 (MQ=255)

                     tCGAACAATATACGGGTCTAATACATTGAAGTTGAAGTTGTACATAGTGCGGTTAAACCCTTCTAACgg                                             >  1:1033221/1‑69 (MQ=255)

                     tCGAACAATATACGGGTCTAATACATTGAAGTTGAAGTTGTACATAGTGCGGTTAAACCCTTCTAACg                                              >  1:3429790/1‑68 (MQ=255)

                     tCGAACAATATACGGGTCTAATACATTGAAGTTGAAGTTGTACATAGTGCGGTTAAACCCTTCTAACg                                              >  1:3060189/1‑68 (MQ=255)

                         aCAATATACGGGTCTAATACATTGAAGTTGAAGTTGTACATGGTGCGGTTGAACCCTTCTAACGGGt                                           >  1:857442/1‑67 (MQ=255)

                         aCAATATACGGGTCTAATACATTGAAGTTGAAGTTGTACATGGTGCGGTTGAACCCTTCTAACGGGTc                                          >  1:1922784/1‑68 (MQ=255)

                         aCAATATACGGGTCTAATACATTGAAGTTGAAGTTGTACATGGTGCGGTTGAACCCTTCTAACGGGTc                                          >  1:575626/1‑68 (MQ=255)

                            atatACGGGTCTAATACATTGAAGTTGAAGTTGTACATGGTGCGGTTGAACCCTTCTAACGGGTCAGaa                                      >  1:2108636/1‑69 (MQ=255)

                               taCGGGTCTAATACATTGAAGTTGGAGTTGTACATGGTGCGGTTGAACCCTTCTAACGGGTCAGAACg                                    >  1:888379/1‑68 (MQ=255)

                               taCGGGTCTAATACATTGAAGTTGAAGTTGTACATGGTGCGGTTGAACCCTTCTAACGGGTCAGAACGc                                   >  1:788693/1‑69 (MQ=255)

                                 cGGGTCTAATACATTGAAGTTGAAGTTGTACATGGTGCGGTTGAACCCTTCTAACGGGTCAGAACGcc                                  <  1:3416508/68‑1 (MQ=255)

                                 cGGGGCTAATACATTGAAGTTGAAGTTGTACATGGTGCGGTTGAACCCTTCTAACGGGTCAGAACGcc                                  <  1:3059717/68‑1 (MQ=255)

                                           aCATTGAAGTTGAAGTTGTACATGGTGCGGTTGAACCCTTCTAACGGGTCAGAACGCCCTTGCTGATc                        <  1:525986/68‑1 (MQ=255)

                                           aCATTGAAGTTGAAGTTGTACATGGTGCGGTTGAACCCTTCTAACGGGTCAGAACGCCCTTGCTGATc                        <  1:789092/68‑1 (MQ=255)

                                           aCATTGAAGTTGAAGTTGTACATGGTGCGGTTGAACCCTTCTAACGGGTCAGAACGCCCTTGCTGATc                        <  1:2551098/68‑1 (MQ=255)

                                           aCATTGAAGTTGAAGTTGTACATGGTGCGGTTGAACCCTTCTAACGGGTCAGAACGCCCTTGCTGATc                        <  1:1841297/68‑1 (MQ=255)

                                           aCATTGAAGTTGAAGTTGTACATGGTGCGGTTGAACCCTTCTAACGGGTCAGAACGCCCTTGCTGATc                        <  1:361447/68‑1 (MQ=255)

                                            cATTGAAGTTGAAGTTGTACATGGTGCGGTTGAACCCTTCTAACGGGTCAGAACGCCCTTGCTGATCTg                      <  1:2196635/69‑1 (MQ=255)

                                            cATTGAAGTTGAAGTTGTACATGGTGCGGTTGAACCCTTCTAACGGGTCAGAACGCCCTTGCTGATCTg                      <  1:1442311/69‑1 (MQ=255)

                                              ttgaagttgaagTTGTACATGGTGCGGTTGAACCCTTCTAACGGGTCAGAACGCCCTTGCTGATCTGt                     <  1:1714160/68‑1 (MQ=255)

                                              ttgaagttgaagTTGTACATGGTGCGGTTGAACCCTTCTAACGGGTCAGAACGCCCTTGCTGATCTGt                     <  1:3345997/68‑1 (MQ=255)

                                              ttgaagttgaagTTGTACATGGTGCGGTTGAACCCTTCTAACGGGTCAGAACGCCCTTGCTGATCTGt                     <  1:3121613/68‑1 (MQ=255)

                                              ttgaagttgaagTTGTACATGGTGCGGTTGAACCCTTCTAACGGGTCAGAACGCCCTTGCTGATCTGt                     <  1:1171539/68‑1 (MQ=255)

                                              ttgaagttgaagTTGTACATGGTGCGGTTGAACCCTTCTAACGGGTCAGAACGCCCTTGCTGATCTGt                     <  1:3406100/68‑1 (MQ=255)

                                              ttgaagttgaagTTGTACATGGTGCGGTTGAACCCTTCTAACGGGTCAGAACGCCCTTGCTGATCTGt                     <  1:3021414/68‑1 (MQ=255)

                                              ttgaagttgaagTTGTACATGGTGCGGTTGAACCCTTCTAACGGGTCAGAACGCCCTTGCTGATCTGt                     <  1:2495640/68‑1 (MQ=255)

                                              ttgaagttgaagTTGTACATGGTGCGGTTGAACCCTTCTAACGGGTCAGAACGCCCTTGCTGATCTGt                     <  1:740639/68‑1 (MQ=255)

                                              ttgaagttgaagTTGTACATGGTGCGGTTGAACCCTTCTAACGGGTCAGAACGCCCTTGCTGATCTGt                     <  1:1376976/68‑1 (MQ=255)

                                              ttgaagttgaagTTGTACATGGTGCGGTTGAACCCTTCTAACGGGTCAGAACGCCCTTGCTGATCTGt                     <  1:2655166/68‑1 (MQ=255)

                                              ttgaagttgaagTTGTACATGGTGCGGTTGAACCCTTCTAACGGGTCAGAACGCCCTTGCTGATCTGt                     <  1:1363094/68‑1 (MQ=255)

                                              ttgaagttgaagTTGTACATGGTGCGGTTGAACCCTTCTAACGGGTCAGAACGCCCTTGCTGATCTGt                     <  1:2399721/68‑1 (MQ=255)

                                              ttgaagttgaagTTGTACATGGTGCGGTTGAACCCTTCTAACGGGTCAGAACGCCCTTGCTGATCTGTa                    >  1:2613439/1‑69 (MQ=255)

                                               tgaagttgaagtTGTACATGGTGCGGTTGAACCCTTCTAACGGGTCAGAACGCCCTTGCTGATCTGt                     <  1:661535/67‑1 (MQ=255)

                                               tgaagttgaagtTGTACATGGTGCGGTTGAACCCTTCTAACGGGTCAGAACGCCCTTGCTGATCTGt                     <  1:1454873/67‑1 (MQ=255)

                                               tgaagttgaagtTGTACATAGTGCGGTTAAACCCTTCTAACGGGTCAGAACGCCCTTGCTGATCTGt                     <  1:2544125/67‑1 (MQ=255)

                                               tgaagttgaagtTGTACATAGTGCGGTTAAACCCTTCTAACGGGTCAGAACGCCCTTGCTGATCTGt                     <  1:1630989/67‑1 (MQ=255)

                                               tgaagttgaagtTGTACATAGTGCGGTTAAACCCTTCTAACGGGTCAGAACGCCCTTGCTGATCTGt                     <  1:1758684/67‑1 (MQ=255)

                                               tgaagttgaagtTGTACATAGTGCGGTTAAACCCTTCTAACGGGTCAGAACGCCCTTGCTGATCTGt                     <  1:1833118/67‑1 (MQ=255)

                                               tgaagttgaagtTGTACATAGTGCGGTTAAACCCTTCTAACGGGTCAGAACGCCCTTGCTGATCTGt                     <  1:2007886/67‑1 (MQ=255)

                                               tgaagttgaagtTGTACATAGTGCGGTTAAACCCTTCTAACGGGTCAGAACGCCCTTGCTGATCTGt                     <  1:2084711/67‑1 (MQ=255)

                                               tgaagttgaagtTGTACATAGTGCGGTTAAACCCTTCTAACGGGTCAGAACGCCCTTGCTGATCTGt                     <  1:1557076/67‑1 (MQ=255)

                                               tgaagttgaagtTGTACATAGTGCGGTTAAACCCTTCTAACGGGTCAGAACGCCCTTGCTGATCTGt                     <  1:1520661/67‑1 (MQ=255)

                                               tgaagttgaagtTGTACATAGTGCGGTTAAACCCTTCTAACGGGTCAGAACGCCCTTGCTGATCTGt                     <  1:2396033/67‑1 (MQ=255)

                                               tgaagttgaagtTGTACATAGTGCGGTTAAACCCTTCTAACGGGTCAGAACGCCCTTGCTGATCTGt                     <  1:2496759/67‑1 (MQ=255)

                                               tgaagttgaagtTGTACATAGTGCGGTTAAACCCTTCTAACGGGTCAGAACGCCCTTGCTGATCTGt                     <  1:514919/67‑1 (MQ=255)

                                               tgaagttgaagtTGTACATAGTGCGGTTAAACCCTTCTAACGGGTCAGAACGCCCTTGCTGATCTGt                     <  1:2544511/67‑1 (MQ=255)

                                               tgaagttgaagtTGTACATAGTGCGGTTAAACCCTTCTAACGGGTCAGAACGCCCTTGCTGATCTGt                     <  1:861498/67‑1 (MQ=255)

                                               tgaagttgaagtTGTACATAGTGCGGTTAAACCCTTCTAACGGGTCAGAACGCCCTTGCTGATCTGt                     <  1:353339/67‑1 (MQ=255)

                                               tgaagttgaagtTGTACATAGTGCGGTTAAACCCTTCTAACGGGTCAGAACGCCCTTGCTGATCTGt                     <  1:2889724/67‑1 (MQ=255)

                                               tgaagttgaagtTGTACATAGTGCGGTTAAACCCTTCTAACGGGTCAGAACGCCCTTGCTGATCTGt                     <  1:1444291/67‑1 (MQ=255)

                                               tgaagttgaagtTGTACATAGTGCGGTTAAACCCTTCTAACGGGTCAGAACGCCCTTGCTGATCTGt                     <  1:1409170/67‑1 (MQ=255)

                                               tgaagttgaagtTGTACATAGTGCGGTTAAACCCTTCTAACGGGTCAGAACGCCCTTGCTGATCTGt                     <  1:1218516/67‑1 (MQ=255)

                                               tgaagttgaagtTGTACATAGTGCGGTTAAACCCTTCTAACGGGTCAGAACGCCCTTGCTGATCTGt                     <  1:3140897/67‑1 (MQ=255)

                                                      gaagttgTACATGGTGCGGTTGAACCCTTCTAACGGGTCAGAACGCCCTTGCTGATCTGTACCGGAACt            >  1:2324517/1‑69 (MQ=255)

                                                         gttgTACATGGTGCGGTTGAACCCTTCTAACGGGTCAGAACGCCCTTGCTGATCTGTACCGGAACTCg          <  1:3301444/68‑1 (MQ=255)

                                                         gttgTACATGGTGCGGTTGAACCCTTCTAACGGGTCAGAACGCCCTTGCTGATCTGTACCGGAACTCg          <  1:1022534/68‑1 (MQ=255)

                                                         gttgTACATGGTGCGGTTGAACCCTTCTAACGGGTCAGAACGCCCTTGCTGATCTGTACCGGAACTCg          <  1:2993312/68‑1 (MQ=255)

                                                         gttgTACATGGTGCGGTTGAACCCTTCTAACGGGTCAGAACGCCCTTGCTGATCTGTACCGGAACTCg          <  1:3065807/68‑1 (MQ=255)

                                                         gttgTACATGGTGCGGTTGAACCCTTCTAACGGGTCAGAACGCCCTTGCTGATCTGTACCGGAACTCg          <  1:3112253/68‑1 (MQ=255)

                                                           tgTACATGGTGCGGTTGAACCCTTCTAACGGGTCAGAACGCCCTTGCTGATCTGTACCGGAACTCg          <  1:1023839/66‑1 (MQ=255)

                                                                 tGGTGCGGTTGAACCCTTCTAACGGGTCAGAACGCCCTTGCTGATCTGTACCGGAACTCGCCCAcccc  >  1:2345693/1‑68 (MQ=255)

                                                                 tGGTGCGGTTGAACCCTTCTAACGGGTCAGAACGCCCTTGCTGATCTGTACCGGAACTCGCACAcccc  >  1:2853151/1‑68 (MQ=255)

                                                                 tGGTGCGGTTGAACCCTTCTAACGGGTCAGAACGCCCTTGCTGATCTGTACCGGAACTCGCACAcccc  >  1:909599/1‑68 (MQ=255)

ACGCCAGGCGACAGCGACCGGTCGAACAATATACGGGTCTAATACATTGAAGTTGAAGTTGTACATGGTGCGGTTGAACCCTTCTAACGGGTCAGAACGCCCTTGCTGATCTGTACCGGAACTCGCACACCCC  >  minE/1550080‑1550212

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