BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	618,082	A→G	11.5%	V61V (GTA→GTG)	cmk →	cytidylate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	618,082	0	A	G	11.5%	291.6 / 18.0	154	V61V (GTA→GTG)	cmk	cytidylate kinase
Reads supporting (aligned to +/- strand): ref base A (33/103); new base G (4/14); total (37/117)
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.98e-01

Only 100 of 155 total aligned reads displayed.

ATTTATCGCGTACTGGCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGAT  >  minE/618014‑618141

aTTTATCGCGTACTGGCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTa                                                             >  1:2686249/1‑69 (MQ=255)

            cTGGCATTGGCGGCATTACATCACCATGTTGATGTTGCGT‑GGAAGATGCGCTGGTACCGCTGGCATccc                                                <  1:737805/69‑1 (MQ=255)

              ggCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATccc                                                <  1:1692382/68‑1 (MQ=255)

              ggCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATccc                                                <  1:927298/68‑1 (MQ=255)

              ggCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATccc                                                <  1:1200953/68‑1 (MQ=255)

              ggCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATccc                                                <  1:1236277/68‑1 (MQ=255)

              ggCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATccc                                                <  1:431954/68‑1 (MQ=255)

              ggCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATccc                                                <  1:1559930/68‑1 (MQ=255)

              ggCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATccc                                                <  1:1767719/68‑1 (MQ=255)

              ggCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATccc                                                <  1:1505571/68‑1 (MQ=255)

              ggCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCa                                               >  1:2964208/1‑69 (MQ=255)

              ggCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCa                                               >  1:2796475/1‑69 (MQ=255)

              ggCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCa                                               >  1:2207205/1‑69 (MQ=255)

              ggCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTAACGCTGGCATCCCa                                               >  1:1706546/1‑69 (MQ=255)

                   tGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTgg                                          <  1:2681631/69‑1 (MQ=255)

                   tGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTgg                                          <  1:583781/69‑1 (MQ=255)

                   tGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTgg                                          <  1:1402568/69‑1 (MQ=255)

                   tGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTgg                                          <  1:2242161/69‑1 (MQ=255)

                   tGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTgg                                          <  1:1864268/69‑1 (MQ=255)

                   tGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTgg                                          <  1:2536616/69‑1 (MQ=255)

                   tGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTgg                                          <  1:2895788/69‑1 (MQ=255)

                   tGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTgg                                          <  1:3381610/69‑1 (MQ=255)

                   tGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCAGGCATCCCATCTgg                                          <  1:898135/69‑1 (MQ=255)

                            tACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATc                                                  >  1:581458/1‑52 (MQ=255)

                            tACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATc                                                  >  1:2765774/1‑52 (MQ=255)

                            tACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATc                                                  >  1:1969412/1‑52 (MQ=255)

                            tACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATc                                                  >  1:2642756/1‑52 (MQ=255)

                            tACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATc                                                  >  1:2558880/1‑52 (MQ=255)

                            tACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGAATc                                                  >  1:457021/1‑52 (MQ=255)

                            tACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGCCATc                                                  >  1:3106297/1‑52 (MQ=255)

                             aCATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATcc                                                 >  1:1155502/1‑52 (MQ=255)

                             aCATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGttt                                >  1:2146157/1‑69 (MQ=255)

                             aCATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGttt                                >  1:2484340/1‑69 (MQ=255)

                             aCATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGttt                                >  1:1640470/1‑69 (MQ=255)

                             aCATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGttt                                >  1:676037/1‑69 (MQ=255)

                             aCATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGttt                                >  1:1241323/1‑69 (MQ=255)

                             aCATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGttt                                >  1:2905224/1‑69 (MQ=255)

                                    cATGTTGATGTTGCATCGGAAGATGCGCTGGTGCCGCTGGCATCCCATCTGGATGTACGTTTTGTGTca                         <  1:1629153/69‑2 (MQ=255)

                                    cATGTTGATGTTGCATCGGAAGATGCGCTGGTGCCGCTGGCATCCCATCTGGATGTACGTTTTGTGTca                         <  1:220508/69‑2 (MQ=255)

                                    cATGTTGATGTTGCATCGGAAGATGCGCTGGTGCCGCTGGCATCCCATCTGGATGTACGTTTTGTGTca                         <  1:2859246/69‑2 (MQ=255)

                                    cATGTTGATGTTGCATCGGAAGATGCGCTGGTGCCGCTGGCATCCCATCTGGATGTACGTTTTGTGTca                         <  1:3068798/69‑2 (MQ=255)

                                    cATGTTGATGTTGCATCGGAAGATGCGCTGGTGCCGCTGGCATCCCATCTGGATGTACGTTTTGTGTca                         <  1:1156042/69‑2 (MQ=255)

                                    cATGTTGATGTTGCATCGGAAGATGCGCTGGTGCCGCTGGCATCCCATCTGGATGTACGTTTTGTGTca                         <  1:3299411/69‑2 (MQ=255)

                                     atgttgatgttgCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATg                                       <  1:3323854/54‑1 (MQ=255)

                                     atgttgatgttgCATCGGAAGATGCGCTGGTGCCGCTGGCATCCCATCTGGATGTACGTTTTGTGTca                         <  1:447183/68‑2 (MQ=255)

                                     atgttgatgttgCATCGGAAGATGCGCTGGTGCCGCTGGCATCCCATCTGGATGTACGTTTTGTGTca                         <  1:1570318/68‑2 (MQ=255)

                                       gttgatgttgCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGAcc                      <  1:2697321/69‑1 (MQ=255)

                                       gttgatgttgCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGAcc                      <  1:2747128/69‑1 (MQ=255)

                                       gttgatgttgCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGAcc                      <  1:1955333/69‑1 (MQ=255)

                                       gttgatgttgCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGAcc                      <  1:627702/69‑1 (MQ=255)

                                       gttgatgttgCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGAcc                      <  1:3280399/69‑1 (MQ=255)

                                       gttgatgttgCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGAcc                      <  1:1363901/69‑1 (MQ=255)

                                        ttgatgttgCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCa                     >  1:3410157/1‑69 (MQ=255)

                                        ttgatgttgCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCa                     >  1:1127936/1‑69 (MQ=255)

                                          gatgttgCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCaa                    <  1:1892785/68‑1 (MQ=255)

                                              ttgCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCa               <  1:2039385/69‑1 (MQ=255)

                                              ttgCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCa               <  1:1706511/69‑1 (MQ=255)

                                              ttgCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCa               <  1:650636/69‑1 (MQ=255)

                                              ttgCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCa               <  1:934779/69‑1 (MQ=255)

                                               tgCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCa               <  1:2248351/68‑1 (MQ=255)

                                                gCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCt                                            >  1:957782/1‑38 (MQ=255)

                                                gCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTg                                           >  1:1671550/1‑39 (MQ=255)

                                                gCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTg                                           >  1:2942513/1‑39 (MQ=255)

                                                gCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAAt             >  1:1556018/1‑69 (MQ=255)

                                                     ggAAGATGCGCTGGTACCGCTGGCATCCCATCTgg                                          <  1:3140385/35‑1 (MQ=255)

                                                         gATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTg    <  1:85133/69‑1 (MQ=255)

                                                         gATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTg    <  1:3065608/69‑1 (MQ=255)

                                                         gATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTg    <  1:1053898/69‑1 (MQ=255)

                                                         gATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTg    <  1:1020196/69‑1 (MQ=255)

                                                         gATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTg    <  1:1678887/69‑1 (MQ=255)

                                                         gATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTg    <  1:191657/69‑1 (MQ=255)

                                                         gATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTg    <  1:2397567/69‑1 (MQ=255)

                                                         gATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTg    <  1:2295513/69‑1 (MQ=255)

                                                          aTGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGa   <  1:486120/69‑1 (MQ=255)

                                                          aTGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGa   <  1:1340092/69‑1 (MQ=255)

                                                          aTGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGa   <  1:1516963/69‑1 (MQ=255)

                                                          aTGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGa   <  1:726209/69‑1 (MQ=255)

                                                          aTGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGa   <  1:704344/69‑1 (MQ=255)

                                                          aTGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGa   <  1:1822816/69‑1 (MQ=255)

                                                          aTGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGa   <  1:200107/69‑1 (MQ=255)

                                                          aTGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGa   <  1:3124684/69‑1 (MQ=255)

                                                          aTGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGa   <  1:403907/69‑1 (MQ=255)

                                                          aTGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGa   <  1:3202610/69‑1 (MQ=255)

                                                          aTGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGGCGACCAATGGCAATCTGGAAGTGa   <  1:1109085/69‑1 (MQ=255)

                                                           tGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTc                          <  1:2880443/45‑1 (MQ=255)

                                                           tGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTc                          <  1:2429016/45‑1 (MQ=255)

                                                           tGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTc                          <  1:2915416/45‑1 (MQ=255)

                                                           tGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTc                          <  1:2843919/45‑1 (MQ=255)

                                                           tGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTc                          <  1:1439913/45‑1 (MQ=255)

                                                           tGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTc                          <  1:2528082/45‑1 (MQ=255)

                                                           tGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGAt  <  1:1669812/69‑1 (MQ=255)

                                                           tGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGAt  <  1:1179069/69‑1 (MQ=255)

                                                           tGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGAt  <  1:1517077/69‑1 (MQ=255)

                                                           tGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGAt  <  1:633972/69‑1 (MQ=255)

                                                           tGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGAt  <  1:3208911/69‑1 (MQ=255)

                                                           tGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGAt  <  1:247906/69‑1 (MQ=255)

                                                           tGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGAt  <  1:1896901/69‑1 (MQ=255)

                                                           tGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGAt  <  1:487706/69‑1 (MQ=255)

                                                           tGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGAt  <  1:2330221/69‑1 (MQ=255)

                                                           tGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGAt  <  1:2070323/69‑1 (MQ=255)

ATTTATCGCGTACTGGCATTGGCGGCATTACATCACCATGTTGATGTTGCGTCGGAAGATGCGCTGGTACCGCTGGCATCCCATCTGGATGTACGTTTTGTGTCGACCAATGGCAATCTGGAAGTGAT  >  minE/618014‑618141

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